%I #37 Jan 14 2024 12:11:13
%S 1,14,30,105,248,264,418,714,1485,3080,3135,3596,3828,3956,4064,5396,
%T 6678,8636,10098,12648,20026,20790,21318,22152,23374,24882,25714,
%U 26040,35074,35343,39105,41656,43890,44660,49938,55154,56134,56536,61344,71145,74613,86304,87087,94944
%N Integers n such that sigma(n)/phi(n) is a perfect square.
%C From _Robert Israel_, Dec 12 2017: (Start)
%C Intersection of A011257 and A020492.
%C If x and y are coprime members of the sequence, then x*y is in the sequence.
%C Contains all members of A133028 except 3. (End)
%H Giovanni Resta, <a href="/A293391/b293391.txt">Table of n, a(n) for n = 1..10000</a> (first 349 terms from Robert Israel)
%H J. A. B. Dris and C. Leibovici, <a href="https://math.stackexchange.com/q/2462552">Why is this sequence not in the OEIS?</a>, October 8 2017.
%H ProofWiki, <a href="https://proofwiki.org/wiki/Integers_whose_Ratio_between_Sigma_and_Phi_is_Square">Integers whose ratio between sigma and phi is square</a>, misses the second term, 14 as of Dec 2017.
%F a(n) = sigma(n)/phi(n) = m^2, for some integer m.
%e sigma(14)=3*8=24, phi(14)=14*(1/2)*(6/7)=6, sigma(14)/phi(14)=2^2, so 14 is in the list.
%p for n from 1 to 100000 do
%p r := numtheory[sigma](n)/numtheory[phi](n) ;
%p if issqr(r) then
%p printf("%d,",n) ;
%p end if;
%p end do: # _R. J. Mathar_, Dec 07 2017
%t Select[Range[10^5], IntegerQ@ Sqrt[DivisorSigma[1, #]/EulerPhi[#]] &] (* _Michael De Vlieger_, Dec 08 2017 *)
%o (PARI) isok(n) = my(q=sigma(n)/eulerphi(n)); issquare(q) && (denominator(q) == 1); \\ _Michel Marcus_, Dec 07 2017; corrected Sep 21 2019
%Y Cf. A000010, A000203, A011257, A020492, A133028, A289336, A289412, A292422.
%K nonn
%O 1,2
%A _Jose Arnaldo Bebita Dris_, Oct 08 2017