OFFSET
1,3
COMMENTS
2*a(n) + 1: 1, 3, 627, 61707, 85985, 495507, 1169979, 8261497, ...
MATHEMATICA
Do[ k=2*n+1; f=PowerMod[ 2, n, k ] + PowerMod[ n, 2, k ]; If[ IntegerQ[ f/k ], Print[ n ] ], {n, 0, 10^7} ] (* Metin Sariyar, Sep 21 2019 *)
PROG
(Magma) [n: n in [0..100000] | Denominator((2^n+n^2)/(2*n+1)) eq 1];
(PARI) is(n) = Mod(2, 2*n+1)^n==-n^2 \\ Felix Fröhlich, Sep 20 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Sep 20 2019
EXTENSIONS
More terms from Felix Fröhlich, Sep 20 2019
Offset 1 from Michel Marcus, Jul 02 2021
STATUS
approved