A097703 Numbers j such that m = 216*j + 108 satisfies sigma(m) != 2*usigma(m).
1, 4, 7, 10, 12, 13, 16, 19, 22, 24, 25, 28, 31, 34, 37, 40, 43, 46, 49, 52, 55, 58, 60, 61, 62, 64, 67, 70, 73, 76, 79, 82, 84, 85, 87, 88, 91, 94, 97, 100, 103, 106, 109, 112, 115, 118, 121, 122, 124, 127, 130, 133, 136, 137, 139, 142, 144, 145, 148, 151, 154, 157
Offset: 1
Keywords
Links
- Amiram Eldar, Table of n, a(n) for n = 1..10000
Programs
-
Mathematica
usigma[n_] := Block[{d = Divisors[n]}, Plus @@ Select[d, GCD[ #, n/# ] == 1 &]]; Complement[ Range[157], (Select[ Range[37000], DivisorSigma[1, # ] == 2usigma[ # ] &] - 108)/216] (* Robert G. Wilson v, Aug 28 2004 *)
-
PARI
is(k) = {my(f = factor(216*k + 108)); sigma(f) != 2 * prod(i = 1, #f~, 1 + f[i,1]^f[i,2]);} \\ Amiram Eldar, Aug 31 2024
Comments