|
%I #12 Nov 15 2018 14:08:43
%S 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,
%T 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,
%U 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
%N a(n) = n/A100762(n).
%H Antti Karttunen, <a href="/A082725/b082725.txt">Table of n, a(n) for n = 1..16384</a>
%t {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 *)
%o (PARI)
%o A100549(n) = if(1==n,1,prime(primepi(1+vecmax(factor(n)[,2]))));
%o 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)));
%o A082725(n) = (n/A100762(n)); \\ _Antti Karttunen_, Nov 11 2018
%Y Cf. A100762, A100549, A141586.
%K nonn
%O 1,3
%A _N. J. A. Sloane_, Nov 17 2008
|
