OFFSET
1,1
PROG
(PARI) a(n) = gcd(apply(x->x-1, divisors(2*n+1)));
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Nov 10 2015
STATUS
approved