|
%I #16 Oct 13 2019 02:32:00
%S 1,2,36,108,128,225,288,441,450,600,864,882,1089,1176,1521,1944,2000,
%T 2178,2601,2904,3042,3240,3249,4056,4761,5202,6144,6498,6561,6936,
%U 7569,8649,8664,9522,10000,11760,12321,12696,13122,15129,15138,16641,17298
%N Numbers k such that k/A000005(k) is a square.
%C Subsequence of A033950.
%H Amiram Eldar, <a href="/A145450/b145450.txt">Table of n, a(n) for n = 1..5479</a> (terms below 10^9)
%e 36/A000005(36) = 36/9 = 4 = 2^2, hence 36 is in the sequence.
%t Select[Range[20000],IntegerQ[Sqrt[#/DivisorSigma[0,#]]]&] _Harvey P. Dale_, May 21 2012
%o (PARI) {m=18000; for(n=1, m, s=sigma(n,0); if(n%s==0&&issquare(n/s), print1(n,",")))}
%Y Cf. A000005 (sigma_0, number of divisors of n), A033950 (refactorable numbers), A140480 (RMS numbers).
%K nonn
%O 1,2
%A _Ctibor O. Zizka_, Oct 10 2008
%E Edited and extended by _Klaus Brockhaus_, Oct 15 2008
|
