login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 


Indices of primes in A057137.
11

%I #21 Apr 16 2024 14:02:57

%S 171,277,367,561,567,18881

%N Indices of primes in A057137.

%C Digits are in ascending order beginning with 1 and after 9 comes 0.

%C Indices of primes in A057137.

%C All terms must end in 1 or 7: A057137(n) is even when n is even, and divisible by 3 iff n == 0, 2, 3, 5, 6, 8 or 9 (mod 10). - _M. F. Hasler_, Apr 14 2024

%D Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 61, 298.

%e a(1) = 12345678901234567890...01234567890...012345678901 = A057137(171) is the first prime term in A057137.

%t fQ[n_] := PrimeQ@ FromDigits@ Mod[Range@n, 10]; lst = {}; Do[ If[fQ@n, AppendTo[lst, n]; Print@n], {n, 10000}]; lst

%t IntegerLength[Select[Table[FromDigits[PadRight[{},n,{1,2,3,4,5,6,7,8,9,0}]],{n,1,1001,2}],PrimeQ]] (* _Harvey P. Dale_, Feb 07 2024 *)

%o (PARI) N=0;for(n=1,600,if(ispseudoprime(N=10*N+n%10),print1(n", "))) \\ _Charles R Greathouse IV_, May 10 2014 (Comment: Surprisingly, this is faster than calling ispseudoprime() only when n ends in 1 or 7, even when much larger N's are considered, e.g., up to 3000. - _M. F. Hasler_, Apr 14 2024)

%Y Cf. A006055, A057137, A120828, A120820, A120821, A120822, A120823, A120824, A120825, A120826, A120827.

%K base,hard,nonn

%O 1,1

%A _Robert G. Wilson v_, Jul 05 2006

%E 18881 from Arjen K. Lenstra, Feb 20 2012

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified September 22 06:15 EDT 2024. Contains 376097 sequences. (Running on oeis4.)