OFFSET
1,1
COMMENTS
Least k such that k*A020522(n)-1 is prime. - Michel Marcus, Apr 13 2021
LINKS
Robert Israel, Table of n, a(n) for n = 1..2090
EXAMPLE
1*2^6*(2^6-1) - 1 = 4031 = 29*139
2*2^6*(2^6-1) - 1 = 8063 = 11*733
3*2^6*(2^6-1) - 1 = 12095 = 5*2419
4*2^6*(2^6-1) - 1 = 16127, which is prime, so a(6)=4.
MAPLE
f:= proc(n) local c, k;
c:= 2^n*(2^n-1);
for k from 1 do if isprime(c*k-1) then return k fi od
end proc:
map(f, [$1..100]); # Robert Israel, Apr 12 2021
PROG
(PARI) a(n) = my(k=1); while(!isprime(k*2^n*(2^n-1) - 1), k++); k; \\ Michel Marcus, Apr 13 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Nov 18 2004
STATUS
approved