OFFSET
1,1
COMMENTS
a(5) > 10^8 if it exists. - Felix Fröhlich, Jul 01 2016
No more terms < 6.5*10^14. - Jud McCranie, Dec 02 2019
EXAMPLE
sigma(4) mod (4+3) = 7 mod 7 = 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, 3);
MATHEMATICA
Select[Range[10^6], Mod[DivisorSigma[1, #], # + 3] == 0 &] (* Michael De Vlieger, Jul 01 2016 *)
PROG
(PARI) is(n) = Mod(sigma(n), n+3)==0 \\ Felix Fröhlich, Jul 01 2016
(Magma) [n: n in [1..2*10^6] | SumOfDivisors(n) mod (n+3) eq 0 ]; // Vincenzo Librandi, Jul 02 2016
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Paolo P. Lava, Jun 28 2016
STATUS
approved