%I #20 Jan 29 2022 13:10:59
%S 1,1,2,7,2,2,2,5,13,2,2,7,2,2,2,31,2,13,2,7,2,2,2,5,31,2,5,7,2,2,2,7,
%T 2,2,2,13,2,5,2,5,2,2,2,7,13,2,2,31,19,31,2,7,2,5,2,5,5,5,2,7,2,2,13,
%U 127,2,2,2,7,2,2,2,13,2,2,31,7,2,2,2,31,11,2,2,7,2,2,5,5,2,13,2,7,2,2,2,7,2
%N Greatest common prime-divisor of sigma_1(n)=A000203(n) and sigma_2(n)=A001157(n); a(n)=1 if no common prime-divisor exists.
%H Antti Karttunen, <a href="/A082066/b082066.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = A006530(A179931(n)). - _Reinhard Zumkeller_, Jul 10 2011
%t ffi[x_] := Flatten[FactorInteger[x]] lf[x_] := Length[FactorInteger[x]] ba[x_] := Table[Part[ffi[x], 2*w-1], {w, 1, lf[x]}] f1[x_] := DivisorSigma[1, n]; f2[x_] := DivisorSigma[2, x] Table[Max[Intersection[ba[f1[w]], ba[f2[w]]]], {w, 1, 128}]
%t (* Second program: *)
%t Table[Last[Apply[Intersection, FactorInteger[Map[DivisorSigma[#, n] &, {1, 2}]][[All, All, 1]]] /. {} -> {1}], {n, 109}] (* _Michael De Vlieger_, May 22 2017 *)
%o (PARI) gpf(n)=if(n>1,my(f=factor(n)[,1]);f[#f],1)
%o a(n)=gpf(gcd(sigma(n),sigma(n,2))) \\ _Charles R Greathouse IV_, Feb 19 2013
%o (Python)
%o from sympy import primefactors, gcd, divisor_sigma
%o def a006530(n): return 1 if n==1 else primefactors(n)[-1]
%o def a(n): return a006530(gcd(divisor_sigma(n), divisor_sigma(n, 2))) # _Indranil Ghosh_, May 22 2017
%Y Cf. A006530, A001157, A000203, A082061-A082065.
%K nonn
%O 1,3
%A _Labos Elemer_, Apr 07 2003
%E Changed "was found" to "exists" in definition. - _N. J. A. Sloane_, Jan 29 2022
|