login
A082725
a(n) = n/A100762(n).
4
1, 1, 3, 1, 5, 3, 7, 1, 1, 5, 11, 1, 13, 7, 15, 1, 17, 1, 19, 5, 21, 11, 23, 1, 25, 13, 1, 7, 29, 15, 31, 1, 33, 17, 35, 1, 37, 19, 39, 5, 41, 21, 43, 11, 5, 23, 47, 1, 49, 25, 51, 13, 53, 1, 55, 7, 57, 29, 59, 5, 61, 31, 7, 1, 65, 33, 67, 17, 69, 35, 71, 1, 73, 37, 25, 19, 77, 39, 79, 1, 1, 41
OFFSET
1,3
LINKS
MATHEMATICA
{1}~Join~Table[Function[{q, P}, n/Times @@ Power @@@ Select[q, First@ # <= P &]] @@ {#, Prime@ PrimePi[1 + Max@ #[[All, -1]] ]} &@ FactorInteger[n], {n, 2, 82}] (* Michael De Vlieger, Nov 13 2018 *)
PROG
(PARI)
A100549(n) = if(1==n, 1, prime(primepi(1+vecmax(factor(n)[, 2]))));
A100762(n) = if(1==n, 1, my(u = A100549(n), f=factor(n)); prod(i=1, #f~, if(f[i, 1]<=u, f[i, 1]^f[i, 2], 1)));
A082725(n) = (n/A100762(n)); \\ Antti Karttunen, Nov 11 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Nov 17 2008
STATUS
approved