%I #20 Aug 23 2022 10:17:54
%S 11,47,251,9431,191,19457,280627,2213,1006253,9129563,66945301,
%T 184171621,726512053,2732087209,10206934519,59883612989,25650350371
%N a(n) is the first prime that is the start of a sequence of exactly n consecutive primes that are in A354589.
%e a(3) = 251 because the 4 primes starting with 251 are 251, 257, 263, 269, and 251, 257 and 263 are in A354589 but 269 is not; and 251 is the least prime for which this works.
%p P:= select(isprime, [seq(i,i=1..10^7,2)]):
%p nP:= nops(P);
%p P10:= P mod 10:
%p s:= 0: V:= Vector(10):
%p for i from 1 to nP-3 do
%p if convert(P10[i..i+3],set) = {1,3,7,9} then s:= s+1
%p else
%p if s > 0 and s <= 10 and V[s] = 0 then V[s]:= P[i-s] fi;
%p s:= 0
%p fi
%p od:
%p convert(V, list);
%Y Cf. A007652, A354589.
%K nonn,more,base
%O 1,1
%A _J. M. Bergot_ and _Robert Israel_, Aug 18 2022
%E a(11)-a(14) from _Michael S. Branicky_, Aug 18 2022
%E a(15)-a(17) from _Michael S. Branicky_, Aug 20 2022