|
%I #22 Jan 29 2022 13:04:59
%S 1,1,2,1,2,2,2,1,1,2,2,2,2,2,2,1,2,1,2,2,2,2,2,2,1,2,2,2,2,2,2,1,2,2,
%T 2,3,2,2,2,2,2,2,2,2,2,2,2,2,3,5,2,2,2,2,2,2,2,2,2,2,2,2,2,1,2,2,2,2,
%U 2,2,2,1,2,2,2,2,2,2,2,2,1,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,2,1,2,2,2,2,2
%N Smallest common prime-divisor of phi(n) = A000010(n) and sigma_2(n) = A001157(n); a(n)=1 if no common prime-divisor exists.
%H Antti Karttunen, <a href="/A082071/b082071.txt">Table of n, a(n) for n = 1..16384</a>
%F a(n) = A020639(gcd(A000010(n), A001157(n))). - _Antti Karttunen_, Nov 03 2017
%t Array[If[CoprimeQ[#1, #2], 1, Min@ Apply[Intersection, Map[FactorInteger[#][[All, 1]] &, {#1, #2}]]] & @@ {EulerPhi@ #,
%t DivisorSigma[2, #]} &, 105] (* _Michael De Vlieger_, Nov 03 2017 *)
%o (PARI)
%o A020639(n) = if(1==n,n,vecmin(factor(n)[, 1]));
%o A082071(n) = A020639(gcd(eulerphi(n),sigma(n,2))); \\ _Antti Karttunen_, Nov 03 2017
%Y Cf. A000010, A001157, A020639.
%Y Cf. also A082065, A082067, A082068, A082069, A082070, A082072.
%K nonn
%O 1,3
%A _Labos Elemer_, Apr 07 2003
%E Values corrected by _R. J. Mathar_, Jul 09 2011
%E More terms from _Antti Karttunen_, Nov 03 2017
%E Changed "was found" to "exists" in definition. - _N. J. A. Sloane_, Jan 29 2022
|
