login
Numbers k that are not squarefree such that the difference between sigma(k) and usigma(k) is a square > 0.
1

%I #20 Jun 24 2019 04:38:45

%S 18,28,40,54,68,84,99,120,124,184,204,208,220,284,297,315,372,388,423,

%T 424,475,508,552,616,624,660,765,796,852,900,928,940,945,963,964,1012,

%U 1152,1164,1192,1269,1272,1348,1395,1425,1449,1458,1496,1524,1664,1719,1796,1848,1975

%N Numbers k that are not squarefree such that the difference between sigma(k) and usigma(k) is a square > 0.

%C This sequence is a subsequence of A013929 (nonsquarefree numbers).

%C If n were a squarefree number (A005117), then the difference would be 0.

%C It appears that the terms are of the form q*p^q, p prime.

%C The corresponding squares are: 9, 16, 36, 36, 36, 64, 36, 144, 64, 144, 144, 196, 144,....

%H Amiram Eldar, <a href="/A216259/b216259.txt">Table of n, a(n) for n = 1..10000</a>

%p with(numtheory): for n from 1 to 2000 do :it:=1:s:=0:x:=divisors(n): n1:=nops(x): for k from 1 to n1 do:d:=x[k]:if gcd(d,n/d)=1 then s:=s+d:else fi:od: s1:=sigma(n): if sqrt(s1-s)=floor(sqrt(s1-s)) and s1>s then printf(`%d, `,n):else fi:od:

%t lst={}; usigma[n_] := Block[{d=Divisors[n]}, DivisorSigma[1,n] - Plus@@Select[d,GCD[#,n/#] == 1&]]; Do[If[IntegerQ[Sqrt[usigma[n] && usigma[n] > 0]], AppendTo[lst,n]], {n,2000}]; lst

%Y Cf. A000203 (sigma:sum of divisors), A034448 (usigma:sum of unitary divisors).

%Y Cf. A064212 (sigma+usigma), A013929 (nonsquarefree numbers).

%K nonn

%O 1,1

%A _Michel Lagneau_, Mar 15 2013