#-----------
# nの根基（radical）を求める
# 引数として素数データ配列、素数データ数が必要
#-----------
rad=function(n,arr_pn,maxi){
  i=1
  flag1=0
  flag2=0
  b=1
  while(i<=maxi && flag1==0){ #素数でまわす
    if(n<arr_pn[i,1]*arr_pn[i,1]){ #素因数がない範囲に入った場合
      if(flag2==0){ #自分自身がこれまで出てきた素因数かどうか
        b=b*n
        return(b)
      }
    }
    r=n%%arr_pn[i,1] #素数で割れるかどうか
    if(r==0){ #素数で割れる場合
      if(flag2==0){
        b=b*arr_pn[i,1] #根基の素因数とする
        flag2=1
      }
      n=n/arr_pn[i,1] #素因数で割る
    }else{
      i=i+1 #次の素数に移る
      flag2=0
    }
  }
}