%I#31 2023年5月19日15:21:26

%S 1,2,3,5,6,9,4,7,10,17,12,11,13,8,19,27,21,26,20,16,14,28,22,18,47,30，

%电话：31,23,34,37,41,45,49,33,36,58,24,62,39,56,42,93,54,25,51,53,15,70,72，

%U 73,46,50,64,97,55,57171,96,79,81,66,71132,89121,29,61,60177,32

%N自然数的重排，使得每个部分积+1都是素数。

%C From _Robert G.Wilson v_，2012年12月24日：（开始）

%C记录：1、2、3、5、6、9、10、17、19、27、28、47、49、58、62、93、97、171、177、184、221、243、470、512、573、768、856、999、1028、1226、1659、2522。。。

%C逾期记录：1、2、3、4、7、8、14、15、29、32、35、59、75。。。（完）

%H Michael S.Branicky，n表，n的a（n）表示n=1..1600（术语751..954来自Bill McEachen，1..750来自Robert G.Wilson v）

%F猜想：n/log（n）<<a（n）<<n*log（n_托马斯·奥多夫斯基，2017年8月9日

%t f[s_List]：=块[{k=1，p=Times@@s}，而[MemberQ[s，k]||！素数Q[k*p+1]，k++]；附加[s，k]]；巢穴[f，{1}，69]（*_Robert G.Wilson v_，2012年12月24日*）

%o（PARI）v=[1]；n=1；而（n<100，s=1+n*prod（i=1，#v，v[i]）；if（isprime（s））&&！vecsearch（vecsort（v），n），v=concat（v，n）；n=0）；n++）；v\\_Derek Orr_，2015年6月16日

%o（Python）

%o从gmpy2导入is_prime

%o从itertools导入岛

%o def agen（startp=1，startset=set（））：#术语生成器

%o a set，p，mink=启动集，启动p，1

%o为True时：

%o an=水貂

%o当a在aset中或不在is_prime中时（p*an+1）：an+=1

%o产生an；增加（a）；p*=一个

%o当水貂在aset中时：AST.丢弃（水貂）；水貂+=1

%o打印（列表（islice（agen），70））#_Michael S.Branicky_，2023年5月19日

%Y参考A073674、A083769。

%K nonn公司

%O 1,2号机组

%A _Amarnath Murthy，2002年8月11日

%E 2003年2月1日，来自_Sascha Kurz的更多条款