OFFSET
1,7
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Michael De Vlieger, Plot f(a(n)) at (x,y) = (n mod m, floor(n/m)) for m = 857 and n = 734449, where f is a color function such that 1 = gray, red indicates primes, gold composite prime powers, green squarefree composites, blue and purple numbers neither squarefree nor prime powers, but purple additionally represents squareful numbers that are not prime powers.
FORMULA
MATHEMATICA
s[n_] := Sqrt[n / Times @@ FactorInteger[n][[;; , 1]]]; s /@ Select[Range[200], AllTrue[FactorInteger[#][[;; , 2]], OddQ] &]
PROG
(PARI) b(n) = {my(f = factor(n)); prod(i = 1, #f~, if(f[i, 2]%2, f[i, 1]^(f[i, 2]-1), 0)); }
lista(kmax) = {my(b1); for(k = 1, kmax, b1 = b(k); if(b1 > 0, print1(sqrtint(b1), ", "))); }
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Nov 17 2023
STATUS
approved