|
%I #26 Feb 10 2019 05:19:19
%S 286330897,10858338851,12869802851,15845166851,29837412851,
%T 45480846799,56676324799,56676324863,68105187851,73915118861,
%U 114737845853,129282912851,154648223809,155738371853,207036953861,271077075851,358515148853,373169411809,373169411861,395705343799
%N Primes such that iteration of A062028 (n + its digit sum) yields 7 primes in a row.
%C The first 15 terms are immediately calculated from A320878(1..200) using the formula.
%H Lars Blomberg, <a href="/A320879/b320879.txt">Table of n, a(n) for n = 1..445</a> (Terms < 10^14)
%H Carlos Rivera, <a href="http://www.primepuzzles.net/puzzles/puzz_163.htm">Puzzle 163. P+SOD(P)</a>
%F A320879 = { n in A320878 | A062028(n) in A320878 } = { n = A320878(k) | A062028(n) = A320878(k+1) }.
%o (PARI) is_A320879(n)=isprime(n=A062028(n))&& is_A320878(n) \\ If possible, use this to select terms from a sufficiently large precomputed array A320878:
%o A320879 = select( is_A320879, A320878) \\ or, in that case, the more efficient:
%o A320879 = select( p->setsearch(A320878, A062028(p)), A320878) \\ or: vecextract(A320878, select(i->A062028(A320878[i])==A320878[i+1], [1..#A320878-1]))
%Y Cf. A047791, A048519, A062028 (n + digit sum of n).
%Y Cf. A048523 .. A048527.
%Y a(1) = A090009(8) = start of first chain of 8 primes under iteration of A062028.
%Y Subsequence of A320878; A320880 is a subsequence.
%K nonn,base
%O 1,1
%A _Zak Seidov_ and _M. F. Hasler_, Nov 08 2018
%E a(16)-a(20) from _Lars Blomberg_, Feb 10 2019
| 