%I #24 Mar 14 2020 05:09:05
%S 15,26,44,56,95,96,119,122,124,140,143,194,215,216,236,287,304,364,
%T 386,407,495,511,527,551,556,560,575,639,740,752,764,780,791,794,815,
%U 871,900,935,936,992,1004,1036,1116,1159,1196,1199,1232,1295,1328,1346
%N Numbers k such that the sum of unitary proper divisors of k is a square > 1.
%C A unitary divisor of n is a divisor d of n such that gcd(d, n/d) = 1.
%H Amiram Eldar, <a href="/A063936/b063936.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1000 from Harry J. Smith)
%F A034460(a(n)) > 1 and A010052(A034460(a(n))) = 1. - _Reinhard Zumkeller_, Aug 15 2012
%e The unitary divisors of 15 are 1,3,5,15 and then the unitary aliquot part is 9 which is a square.
%t us[1] = 0; us[n_] := Times @@ (1 + Power @@@ FactorInteger[n]) - n; Select[Range[1500], (s = us[#]) > 1 && IntegerQ@Sqrt[s] &] (* _Amiram Eldar_, Mar 14 2020 *)
%o (PARI) us(n) = sumdiv(n,d, if(gcd(d,n/d)==1,d));
%o j=[]; for(n=1,3000, if(us(n)-n > 1 && issquare(us(n)-n),j=concat(j,n))); j
%o (PARI) us(n) = sumdiv(n, d, if(gcd(d, n/d)==1, d))
%o { n=0; for (m=1, 10^9, u=us(m) - m; if (issquare(u) && u > 1, write("b063936.txt", n++, " ", m); if (n==1000, break)) ) } \\ _Harry J. Smith_, Sep 03 2009
%o (Haskell)
%o import Data.List (findIndices)
%o a063936 n = a063936_list !! (n-1)
%o a063936_list = map (+ 1) $
%o findIndices (\x -> x > 1 && a010052 x == 1) a034460_list
%o -- _Reinhard Zumkeller_, Aug 15 2012
%Y Cf. A034448, A063937.
%K easy,nonn
%O 1,1
%A _Felice Russo_, Aug 31 2001
%E More terms from _Jason Earls_, Sep 04 2001
|