login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A082725 a(n) = n/A100762(n). 4

%I

%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

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified July 8 01:48 EDT 2020. Contains 335502 sequences. (Running on oeis4.)