OFFSET
1,1
COMMENTS
Tested up to n = 10^6.
EXAMPLE
Anti-divisors of 58 are 3, 4, 5, 9, 13, 23, 39. Their sum is 96.
The only anti-divisor of 96 is 64.
Again, anti-divisors of 64 are 3, 43. Their sum is 46. Finally, anti-divisors of 46 are 3, 4, 7, 13, 31. Their sum is 58 and 58 / 58 = 1.
MAPLE
with(numtheory); P:=proc(q, h) local a, i, j, k, n;
for n from 5 to q do a:=n; for i from 1 to h do
k:=0; j:=a; while j mod 2 <> 1 do k:=k+1; j:=j/2; od;
a:=sigma(2*a+1)+sigma(2*a-1)+sigma(a/2^k)*2^(k+1)-6*a-2; od;
if type(a/n, integer) then print(n); fi; od; end: P(10^6, 4);
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Oct 01 2013
EXTENSIONS
Offset corrected and a(33)-a(42) from Donovan Johnson, Jan 09 2014
STATUS
approved