OFFSET
1,1
COMMENTS
The sequence is not a permutation of prime numbers.
E.g., after calculating 2001 terms of the sequence, the first absent primes are 1973,3719,3917,7193,9137,9173,9371. It's evident that these numbers will never appear in the sequence because any last term of the sequence should use at least one of digits 1,3,7,9.
The first nine terms {2, 3, 5, 7, 11,23, 17, 29,13] coincide with A068863(1..9).
The only fixed points are a(n) for n={1, 2, 3, 4, 5, 7, 12, 13, 17, 19} are {2, 3, 5, 7, 11, 17, 37, 41, 59, 67} that is for these n's a(n)=prime(n)=A000040(n).
a (100*k) for k = 1,20: {443, 1193, 1741, 1621, 4567, 6047, 5851, 6491, 7151, 7559, 9349, 10601, 11119, 11699, 13001, 11839, 14107, 16111, 15073, 16487}.
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..10000
Zak Seidov, Graph of first 2000 terms.
MATHEMATICA
a = {2}; While[ Length[a] < 100, d = IntegerDigits@ Last@ a; p = 2; While[ Intersection[ IntegerDigits@p, d] != {} || MemberQ[a, p], p = NextPrime@ p]; AppendTo[a, p]]; a (* Giovanni Resta, Mar 21 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Mar 20 2017
STATUS
approved