rad(p)をpの異なる素因数の積としてR(p)=2^rad(p)かな
素因数分解は必須のような気がする