OFFSET
1,1
COMMENTS
a(n) = p says "In position 'floor(p/10)' is a digit 'p mod 10'."
Each term must be the smallest possible prime not used earlier.
a(1447) = 19173719153, a(3868) = 371379371929.
No further record value up to n = 10^4. Some earlier record values: a(19) = 541, a(20) = 1009, a(35) = 9001, a(110) = 10007, ..., a(142) = 30011, ..., a(278) = 70001, ..., a(474) = 90001, a(523) = 101009, a(657) = 191339, a(902) = 300007, ..., a(1386) = 300221. - M. F. Hasler, Mar 19 2020
LINKS
M. F. Hasler, Table of n, a(n) for n = 1..10000, Mar 19 2020
Eric Angelini, Primes describing digit positions, SeqFan mailing list, March 2, 2020.
M. F. Hasler, PARI/GP program to compute the first N terms of the sequence, Mar 19 2020
M. F. Hasler, Table of n, a(n) for n = 1..30000, Mar 20 2020
EXAMPLE
a(1) = 11 says "In position 1 is a 1" - which is compatible with the term itself. Since any term must have at least two digits, this is certainly the smallest possibility.
a(2) = 41 says "In position 4 is a 1" - which is indeed the last digit of a(2). There is no smaller solution: the term cannot refer to the 2nd nor the 3rd digit of the sequence, since neither 21 nor 33 is prime.
a(3) = 61 says "In position 6 is a 1"; again, there's no smaller solution.
a(4) = 83 says "In position 8 is a 3", and this is again the smallest solution.
a(5) = 113 says "In position 11 is a 3": again the last digit of a(5) itself, and there is no smaller solution.
a(6) = 101 says: "In position 10 is a 1." (This term wasn't possible earlier, but at this position it is.)
a(19) = 541 says "In position 54 there is a 1", which is not yet there: position 54 is the first digit of a(20). So a(20) must start with a digit 1, and the smallest solution is a(20) = 1009, predicting a digit 9 in position 100.
PROG
(PARI) A333085_vec(n, d=[], U=[], F=[], k)={vector(n, i, forprime(p=11, , setsearch(U, p)&& next; k=divrem(p, 10); k[1] > #d + logint(k[1], 10)+1 || k[2] == if( k[1]<=#d, d[ k[1]], digits( k[1] )[ k[1]-#d ]) || next; for(i=1, #F, F[i][1] > #d + logint(p, 10)+1 && break; F[i][2] == digits(p)[ F[i][1]-#d ] || next(2)); d=concat(d, digits(p)); break); while(#F && F[1][1]<=#d, F=F[^1]); k[1]>#d && F=setunion(F, [k]); U=setunion(U, primes([k[1], k[1]+1]*10)); [10, 1]*k)} \\ For n > 500, use the much faster code given in LINKS. - M. F. Hasler, Mar 18 2020
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Eric Angelini and Hans Havermann, Mar 07 2020
EXTENSIONS
Edited by M. F. Hasler, Mar 18 2020
STATUS
approved