%I#25 2021年7月3日10:57:27

%编号：4,12,6038125221909415128612377921034144848407746766869330

%N A069090的项数，精确到N位。

%H Barry Carter，<a href=“https://github.com/barrycarter/bcapps/blob/master/QUORA/bc-primes.m“>Mathematica处理第一个n，其中随机滚动的数字将产生素数</a>

%p A276707：=程序（n）

%p局部a，k；

%p a：=0；

%p k：=下一素数（10^（n-1））；

%p，而k<10^n do

%p如果是A069090（k），则

%pa:=a+1；

%p end if；

%p k：=下一素数（k）；

%p端do：

%p a；

%结束程序：#R.J.Mathar_，2016年12月15日

%o（perl/tcsh）“bzcat A069090.b.txt.bz2|perl-anle'打印长度（$F[1]）'|sort|uniq-c”带https://github.com/barrycarter/bcapps/blob/master/QUORA/A069090.b.txt.bz2

%o（Python）

%o来自sympy import primerage，isprime

%o定义正常（p）：

%o s=str（p）

%o如果len（s）==1：返回True

%o返回范围（1，len（s））中i的所有（非isprime（int（s[：i]））

%o定义a（n）：返回和（质数范围内p的ok（p）（10**（n-1），10**n））

%o打印（[a（n）代表范围（1，7）中的n）]#_Michael S.Branicky_，2021年7月3日

%o（Python）#跳过错误前缀的更快版本

%o从sympy导入isprime，nextprime

%o定义a（n）：

%o如果n==1：返回4

%o p，c=下一素数（10**（n-1）），0

%o当p<10**n时：

%o s，失败=str（p），错误

%o对于范围（1，n）中的i：

%o ti=整数（s[：i]）

%o如果isprime（ti）：失败=i；打破

%o如果失败：p=下一个素数（（ti+1）*10**（n-i））

%o其他：p，c=下一素数（p），c+1

%o返回c

%o打印（[a（n）代表范围（1，8）中的n）]#_Michael S.Branicky_，2021年7月3日

%Y参考A069090。

%K nonn，基础，更多

%O 1,1号机组

%A _巴里·卡特，2016年9月15日

%E a（9）-a（11），来自Michael S.Branicky，2021年7月3日