OFFSET
1,1
EXAMPLE
sigma(8) mod (8+7) = 15 mod 15 = 0.
MAPLE
with(numtheory); P:=proc(q, h) local n; for n from 1 to q do
if n+h>0 then if type(sigma(n)/(n+h), integer) then print(n); fi; fi; od; end: P(10^9, 7);
MATHEMATICA
Select[Range[10^6], Mod[DivisorSigma[1, #], # + 7] == 0 &] (* Michael De Vlieger, Jul 05 2016 *)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Paolo P. Lava, Jul 05 2016
EXTENSIONS
a(6)-a(7) from Giovanni Resta, Jul 05 2016
STATUS
approved