OFFSET
1,12
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..16384
FORMULA
PROG
(PARI)
A209229(n) = (n && !bitand(n, n-1));
A334870(n) = if(issquare(n), sqrtint(n), my(c=core(n), m=n); forprime(p=2, , if(!(c % p), m/=p; break, m*=p)); (m));
(PARI)
\\ Much faster, A302777 like in above:
A048675(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; };
A052126(n) = if(1==n, n, (n/vecmax(factor(n)[, 1])));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jun 08 2020
STATUS
approved