OFFSET
1,1
LINKS
Rémy Sigrist, Table of n, a(n) for n = 1..8192
EXAMPLE
2*5-1 = 9 is 1001 in binary. Looking at the binary numbers that end with 1001: 1001 = 9 in decimal is composite; 11001 = 25 in decimal is composite. But 101001 = 41 in decimal is prime. So a(5) = 41. - Corrected by Rémy Sigrist, Feb 05 2020
MAPLE
isA105888 := proc(p, n) local pdgs, n21dgs ; pdgs := convert(p, base, 2) ; n21dgs := convert(2*n-1, base, 2) ; if nops(n21dgs) > nops(pdgs) then return false; else verify( [op(1..nops(n21dgs), n21dgs)], [op(1..nops(n21dgs), pdgs)], 'sublist') ; end if; end proc: A105888 := proc(n) p := 2 ; while not isA105888(p, n) do p := nextprime(p) ; end do ; p ; end proc: seq(A105888(n), n=1..80) ; # R. J. Mathar, Dec 06 2009
MATHEMATICA
pr=-16; Select[Prime[Range[200]], MultiplicativeOrder[pr, # ] == #-1 &]
PROG
(PARI) a(n) = my (m=2*n-1); forstep (p=m, oo, 2^#binary(m), if (isprime(p), return (p)))
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Leroy Quet, Aug 08 2009
EXTENSIONS
Extended beyond a(10) by R. J. Mathar, Dec 06 2009
STATUS
approved