|
| |
|
|
A107649
|
|
Numbers n such that (10^(2n+1)+72*10^n-1)/9 is prime.
|
|
44
| | |
|
|
|
OFFSET
| 1,2
|
|
|
COMMENTS
| n is in the sequence iff the palindromic number 1(n).9.1(n) is prime (dot between numbers means concatenation). If n is in the sequence then n is not of the forms 3m, 6m+5, 22m+3, 22m+7, etc. (the proof is easy).
|
|
|
REFERENCES
| C. Caldwell and H. Dubner, "Journal of Recreational Mathematics", Volume 28, No. 1, 1996-97, pp. 1-9.
|
|
|
LINKS
| Patrick De Geest, World!Of Numbers, Palindromic Wing Primes (PWP's)
Makoto Kamada, Prime numbers of the form 11...11911...11
Index entries for primes involving repunits.
|
|
|
FORMULA
| a(n) = (A077795(n)-1)/2.
|
|
|
EXAMPLE
| 26 is in the sequence because (10^(2*26+1)+72*10^26-1)/9=1(26).9.1(26)
= 11111111111111111111111111911111111111111111111111111 is prime.
|
|
|
MATHEMATICA
| Do[If[PrimeQ[(10^(2n + 1) + 72*10^n - 1)/9], Print[n]], {n, 3000}]
|
|
|
PROG
| (PARI) for(n=0, 1e4, if(ispseudoprime(t=(10^(2n+1)+72*10^n)\9), print1(t", "))) \\ Charles R Greathouse IV, Jul 15 2011
|
|
|
CROSSREFS
| Cf. A004023, A077775-A077798, A107123-A107127, A107648, A107649, A115073, A183174-A183187.
Sequence in context: A098443 A052775 A137964 * A052763 A084211 A114496
Adjacent sequences: A107646 A107647 A107648 * A107650 A107651 A107652
|
|
|
KEYWORD
| nonn
|
|
|
AUTHOR
| Farideh Firoozbakht (mymontain(AT)yahoo.com), May 19 2005
|
|
|
EXTENSIONS
| Edited by Ray Chandler (rayjchandler(AT)sbcglobal.net), Dec 28 2010
|
| |
|
|
