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

18, 28, 40, 54, 68, 84, 99, 120, 124, 184, 204, 208, 220, 284, 297, 315, 372, 388, 423, 424, 475, 508, 552, 616, 624, 660, 765, 796, 852, 900, 928, 940, 945, 963, 964, 1012, 1152, 1164, 1192, 1269, 1272, 1348, 1395, 1425, 1449, 1458, 1496, 1524, 1664, 1719, 1796, 1848, 1975

Offset

1,1

Comments

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

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

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

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

Links

Amiram Eldar, Table of n, a(n) for n = 1..10000

Maple

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:

Mathematica

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

Crossrefs

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

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

Sequence in context: A093648 A171221 A363659 * A117101 A063840 A180117

Adjacent sequences: A216256 A216257 A216258 * A216260 A216261 A216262

Keyword

nonn

Author

Michel Lagneau, Mar 15 2013

Status

approved