cp's OEIS Frontend

A058250 GCD of n-th primorial number and its totient.

%I A058250 #39 Aug 16 2024 21:12:31
%S A058250 1,1,2,2,6,30,30,30,30,330,2310,2310,2310,2310,2310,53130,690690,
%T A058250 20030010,20030010,20030010,20030010,20030010,20030010,821230410,
%U A058250 821230410,821230410,821230410,13960916970,739928599410,739928599410
%N A058250 GCD of n-th primorial number and its totient.
%H A058250 Michael De Vlieger, <a href="/A058250/b058250.txt">Table of n, a(n) for n = 0..1000</a>
%F A058250 a(n) = gcd(A002110(n), A000010(A002110(n))) = gcd(A002110(n), A005867(n)).
%F A058250 a(n) = A005867(n) / A038110(n+1).  For example: For n = 4: a(4) = 48 / 8 = 6. - _Jamie Morken_, Apr 12 2019
%e A058250 a(6) = gcd(30030,5760) = 30.
%p A058250 [seq(igcd(product(ithprime(k), k=1..m), product(ithprime(k)-1, k=1..m)), m=1..50)];
%t A058250 GCD[#,EulerPhi[#]]&/@Rest[FoldList[Times,1,Prime[Range[30]]]] (* _Harvey P. Dale_, Dec 19 2012 *)
%t A058250 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 *)
%o A058250 (PARI) a(n) = my(pr=prod(k=1, n, prime(k))); gcd(pr, eulerphi(pr)); \\ _Michel Marcus_, Apr 13 2019
%Y A058250 Cf. A002110, A000010, A005867, A038110.
%K A058250 nonn,look
%O A058250 0,3
%A A058250 _Labos Elemer_, Dec 05 2000
%E A058250 a(0) = 1 inserted by _Michael De Vlieger_, Apr 13 2019