%I#16 2023年4月6日08:30:09

%序号2,14,42144200302724287181056824438

%N数k，使k*（10^k-1）+1是素数。

%C经ECM认证的较大值。不超过6000。

%o（PARI）is（n）=ispseudoprime（n*（10^n-1）+1）\\查尔斯·格里特豪斯IV_，2017年6月13日

%o（Python）

%o来自sympy import isprime

%o定义afind（极限，启动=1）：

%o k，pow10=启动，10**启动

%o对于范围内的k（启动，极限+1）：

%o如果是素数（k*（pow10-1）+1）：打印（k，end=“，”）

%o k+=1

%o功率10*=10

%o afind（500）#_Michael S.Branicky_，2021年8月26日

%Y参考A110807。

%K更多，nonn

%O 1，1

%A _Jason Earls，2005年8月18日

%E a（7）-a（8）摘自_Ryan Propper，2006年9月20日

%E a（9）摘自Michael S.Branicky，2021年8月26日

%E a（10），来自_Michael S.Branicky，2023年4月5日