OFFSET
1,1
COMMENTS
Or: Primes p such that there is exactly one prime in each decade [10d-9, 10d] for 10p-9 <= d < 10p. - M. F. Hasler, May 02 2019
EXAMPLE
4073 is not in the sequence because 407203 and 407207 are both prime and in the same decade. 211217 is not in the sequence because 21121691 and 21121697 are both prime and in the same decade. 5046053 is not in the sequence because 504605291 and 504605293 are both prime and in the same decade. - R. J. Mathar, May 02 2019
MAPLE
isA103608 := proc(n)
local p, dec ;
if not isprime(n) then
false;
else
p := 100*(n-1) ;
for dec from 0 to 9 do
p := nextprime(p) ;
if modp(floor(p/10), 10) <> dec then
return false;
end if;
end do:
p := nextprime(p) ;
if p > 100*n then
true ;
else
false;
end if;
end if;
end proc:
for i from 1 do
p := ithprime(i) ;
if isA103608(p) then
printf("%d, \n", p) ;
end if;
end do: # R. J. Mathar, May 02 2019
PROG
(PARI) select( is_A103608(p)={for(k=10*p-9, 10*p, #primes([10*k-9, 10*k])==1||return); isprime(p)}, primes(10^5)) \\ M. F. Hasler, May 02 2019
CROSSREFS
KEYWORD
nonn,less
AUTHOR
R. J. Mathar, _M. Hasler_, May 02 2019
STATUS
approved