OFFSET
1,1
LINKS
Pierre CAMI, Table of n, a(n) for n = 1..10000
EXAMPLE
2^1+5*2^1-1 = 11 prime so a(1) = 5.
2^2+5*2^1-1 = 13 prime so a(2) = 5.
PROG
PFGW64 from Primeform group and SRYPTIFY
command : pfgw64 -f in.txt
in.txt file :
SCRIPT
DIM kk
DIM nn, 0
DIM mm
DIMS tt
OPENFILEOUT myfil, prem.txt
LABEL loopn
SET nn, nn+1
IF nn>10000 THEN END
IF nn%2==0 THEN SET mm, nn/2
IF nn%2==1 THEN SET mm, nn/2+1
SET kk, 2
LABEL loopk
SET kk, kk+1
SETS tt, %d, %d, %d, %d\ ; nn; kk; p(kk); mm
PRP (2^(nn-mm)+p(kk))*2^mm-1, tt
IF ISPRP THEN GOTO a
IF ISPRIME THEN GOTO a
GOTO loopk
LABEL a
WRITE myfil, tt
GOTO loopn
(Haskell)
a206553 n = head [p | p <- drop 2 a000040_list,
a010051 (2^n + p*2^(div (n+1) 2) - 1) == 1]
-- Reinhard Zumkeller, Feb 10 2012
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Feb 09 2012
STATUS
approved