A082725 a(n) = n/A100762(n).

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

Author

N. J. A. Sloane, Nov 17 2008

Keywords

• nonn

Links

• Antti Karttunen, Table of n, a(n) for n = 1..16384

Crossrefs

Cf. A100762, A100549, A141586.

Programs

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 *)

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