%I#31 2022年12月29日09:10:52

%S 3,31653415914159141591592653288419717950288415935897，

%电话9323846264314159265358992653535897932323846264338327，

%电话：841971693993751862803482534211789793238462643338425342117067982138327950288419716938972384626383279

%N第一个N位素数出现在Pi的十进制展开式A000796中。

%C对于接下来的n位素数，请参阅b文件（链接）。

%C序列A198344给出了这些素数Pi的位置。

%H Michael S.Branicky，n表，n=1..1000的a（n）（M.F.Hasler的术语1..162）

%o（PARI）默认值（realprecision，2000）；A104841（n）={for（c=0，default（realprecision）-n-2，Pi\.1^c%10&ispseudoprime（p=Pi\.1*（n+c-1）%10^n）&return（p））；error（“请增加默认值（realpecision）以计算A10484（“n”）。”）}\\_M.F.Hasler_，2011年10月23日

%o（Python）

%o从sympy导入S，isprime

%o pi=“3”+str（S.pi.n（10**5））[2:]#或从文件加载数据

%o def A104841_A198344（n）：如果pi[i]，则返回next（（（p，i+1），对于范围（len（pi）-n）中的i=“0”和isprime（p:=int（pi[i:i+n]）），“位数不足”）

%o打印（[A104841_A198344（n）[0]表示范围（1，21）内的n）#_Michael S.Branicky_2022年12月28日

%Y参考A104819-A104842。

%Y参考A000796，A198344。

%K nonn，基础

%O 1,1号机组

%A _Zak Seidov，2005年3月27日