%I#17 2023年4月11日07:17:33

%序号1,3,4,15,76,9123113317141942249749635379123972322426395

%N数N，其中R（N）+4*10^层（N/2）是素数，其中R[N]=（10^N-1）/9（单位：A002275）。

%C对于n>1，对应的素数是A105992的子集：近重单位素数。

%C在基数10中，R（n）+4*10^ floor（n/2）有上限（n/2。对于奇数n，这是一个回文（也称为翼素数，参见A077783：奇数项的子序列），对于偶数n，数字5正好位于数字的中间。

%C另请参阅变体A331868，其中数字5正好在中间的右边。

%H Brady Haran和Simon Pampena，<a href=“https://www.youtube.com/watch？v=HPfAnX5blO0网址“>《Glitch Primes and Cyclops Numbers》，数字爱好者视频（2015）。

%e对于n=1，R（1）+4*10^楼层（1/2）=5是素数。

%e对于n=3，R（3）+4*10^楼层（3/2）=151是素数。

%e对于n=4，R（4）+4*10^楼层（4/2）=1511为素数。

%e对于n=15，R（15）+4*10^楼层（15/2）=1111111 51111111是素数。

%t选择[Range[0，2500]，PrimeQ[（10^#-1）/9+4*10^Floor[#/2]]&]

%o（PARI）用于（n=09999，ispseudoprime（p=10^n\9+4*10^（n\2））&&print1（n“，”））

%Y参见A105992（近重单位素数）、A002275（重单位）、A004023（重单位指数）、A011557（10的幂）。

%Y参见A331862、A331861、A33186、A331666（数字0、2、3或4的变体，而不是5），A331868（楼层（n/2-1）而不是楼层（n/2）的变体）。

%Y参考A077783（奇数项）。

%K nonn，基础，硬，更多

%O 1,2号机组

%A _M.F.Hasler，2020年2月9日

%E a（12）-a（14），来自Michael S.Branicky，2023年2月3日

%E a（15）-a（16），来自_Michael S.Branicky_，2023年4月11日