[go: up one dir, main page]

login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A058250 GCD of n-th primorial number and its totient. 8
1, 1, 2, 2, 6, 30, 30, 30, 30, 330, 2310, 2310, 2310, 2310, 2310, 53130, 690690, 20030010, 20030010, 20030010, 20030010, 20030010, 20030010, 821230410, 821230410, 821230410, 821230410, 13960916970, 739928599410, 739928599410 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
LINKS
FORMULA
a(n) = gcd(A002110(n), A000010(A002110(n))) = gcd(A002110(n), A005867(n)).
a(n) = A005867(n) / A038110(n+1). For example: For n = 4: a(4) = 48 / 8 = 6. - Jamie Morken, Apr 12 2019
EXAMPLE
a(6) = gcd(30030,5760) = 30.
MAPLE
[seq(igcd(product(ithprime(k), k=1..m), product(ithprime(k)-1, k=1..m)), m=1..50)];
MATHEMATICA
GCD[#, EulerPhi[#]]&/@Rest[FoldList[Times, 1, Prime[Range[30]]]] (* Harvey P. Dale, Dec 19 2012 *)
Fold[Append[#1, {#1, #2, GCD[#1, #2]} & @@ {#4 #1, #2 (#4 - 1)} & @@ Append[#1[[-1]], #2]] &, {{1, 1, 1}}, Prime@ Range[29]][[All, -1]] (* Michael De Vlieger, Apr 25 2019 *)
PROG
(PARI) a(n) = my(pr=prod(k=1, n, prime(k))); gcd(pr, eulerphi(pr)); \\ Michel Marcus, Apr 13 2019
CROSSREFS
Sequence in context: A004304 A326907 A270487 * A179929 A278258 A067644
KEYWORD
nonn,look,changed
AUTHOR
Labos Elemer, Dec 05 2000
EXTENSIONS
a(0) = 1 inserted by Michael De Vlieger, Apr 13 2019
STATUS
approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 29 08:01 EDT 2024. Contains 375510 sequences. (Running on oeis4.)