OFFSET
1,2
LINKS
Hugo Pfoertner, Table of n, a(n) for n = 1..10000
MATHEMATICA
nn=100; rads=Select[Range[2, nn], GCD@@FactorInteger[#][[All, 2]]===1&];
Table[Position[rads, Prime[n]][[1, 1]], {n, PrimePi[nn]}]
PROG
(PARI) lista(nn) = Vec(select(x->isprime(x), Vec(select(x->(!ispower(x)&&x>1), [1..nn])), 1)); \\ Michel Marcus, May 04 2018
(Python)
from sympy import prime, mobius, integer_nthroot
def A279984(n): return int((p:=prime(n))-1+sum(mobius(k)*(integer_nthroot(p, k)[0]-1) for k in range(2, p.bit_length()))) # Chai Wah Wu, Oct 12 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Dec 24 2016
STATUS
approved