%I #40 Sep 08 2022 08:46:04
%S 1,3,14,12,88,14,116,15,190,88,989,35,477,116,209,105,6901,190,7067,
%T 88,196,989,6439,35,15049,477,2754,172,10207,209,4976,336,989,6901,
%U 1189,190,10877,7067,477,248,13529,377,44461,989,418,6439,79523,105,10244,15049
%N Least number k such that n divides gcd(sigma(k), phi(k)) (A009223).
%C For each n there are infinitely many numbers k for which n divides sigma(k) and phi(k). - _Marius A. Burtea_, Mar 28 2019
%H Giovanni Resta, <a href="/A222713/b222713.txt">Table of n, a(n) for n = 1..10000</a> (terms 1-388 from Marius A. Burtea, 389-2808 from David A. Corneth)
%e Given A009223 = 1, 1, 2, 1, 2, 2, 2, 1, 1, 2, 2, 4, 2, 6, 8, 1, 2, 3, ...
%e 1 first divides A009223(1); 2 first divides A009223(3); 3 first divides A009223(14)=6.
%t Array[Block[{i = 1}, While[Mod[GCD[DivisorSigma[1, i], EulerPhi@ i], #] != 0, i++]; i] &, 50] (* _Michael De Vlieger_, Mar 28 2019 *)
%o (PARI) A009223_hunt(x)=local(n=0,g);while(n++,g=A009223(n);if(g%x,,return(n)));
%o for(x=1,50,print1(A009223_hunt(x)", "))
%o (Magma) [Min([n: n in [1..300000] | IsIntegral(SumOfDivisors(n)/m) and IsIntegral(EulerPhi(n)/m) ]): m in [1..70]]; // _Marius A. Burtea_, Mar 28 2019
%o (Magma) v:=[];
%o for n in [1..60] do
%o m:=1;
%o while not EulerPhi(m) mod n eq 0 or not SumOfDivisors(m) mod n eq 0 do
%o v[n]:=0;
%o m:=m+1;
%o end while;
%o v[n]:=m;
%o end for;
%o v // _Marius A. Burtea_, Mar 30 2019
%Y Cf. A009223, A222714.
%K nonn
%O 1,2
%A _Phil Carmody_, Mar 01 2013
|