A028988 Repdigit + 1 is prime.
1, 2, 4, 6, 22, 66, 88, 222, 666666, 22222222, 66666666, 666666666, 22222222222, 66666666666, 88888888888888, 88888888888888888, 66666666666666666666, 66666666666666666666666, 88888888888888888888888888888888888
Offset: 1
Links
- Michel Marcus, Table of n, a(n) for n = 1..39
- Patrick De Geest, Repdigits: WONplate 83
Programs
-
Mathematica
a[n_]:=NestList[FromDigits[Append[{#},n]]&,n,34]; Union[Join[{1},Flatten[Table[Select[a[n],PrimeQ[#+1]&],{n,2,8,2}]]]] (* Jayanta Basu, May 30 2013 *) Select[Flatten[Table[FromDigits[PadRight[{},n,k]],{n,40},{k,9}]],PrimeQ[ #+1]&] (* Harvey P. Dale, Apr 22 2018 *)
-
PARI
u=30; for(n=1, u, r=(10^n-1)/9; for(a=1, 9, m=r*a; if(ispseudoprime(m+1), print1(m, ", ")))) \\ Michel Marcus, Jan 26 2023; after A028987
Extensions
Offset corrected by Arkadiusz Wesolowski, Aug 14 2011
Comments