这个网站是通过捐款来支持的。OEIS基金会.

M. F. Hasler/A0399

从奥伊斯维基
跳转到:导航搜索

A0399=素数,使得每一个不同的数字排列是复合的(包括排列与前导零点):

2, 3, 5,7, 11, 19,23, 29, 41,43, 47, 53,59, 61, 67,83, 89, 151,211, 223, 227,229, 233, 257,263, 269, 353,383, 409, 431,433, 443, 449,487, 499, 523,487, 499, 523,γ,γ,γ,γ,γ

此序列的代码:(添加到OEIS.Org/A0399《01》第2018卷

A039866-行(n)={*/*,计算长度为n*/i的条件(n>1),最后一个数字仅适用于长度>1 */*,局部(d=EVE(Vec(“0245681379”))/*所有可能的数字,最后可能的结尾数字:相当于vCeSo[(0,9),(a,b)-IS-素数(10 +a)-IS-素数(10 +b))*/*,u=vEctov(n,i,10 ^(n-1))/*幂10:点积dU比OffDiges(d)*/y,Next PiMm()= /*的稍微快和短,将“全局”变量p改为“1…n**/*”的“下一个”置换(i=2,n,/*i=我们将减少的索引位置)。(123)……213,132……(*=)(t=p[i] - 1)其次,/*我们能减少这个索引而不得到0吗?否则转到下一个位置*/*(set)(set(p[i++n.]),t)这个左边的索引可用吗?*[t]=d[p[i] ] /*这会改变数字吗?*(或),T—γ(下)(2)/*…否则尝试下一个更小;如果我们运行到0 */i,转到下一个位置;i i n/*最右边的位置?*/π}位测试(650,d[t])/*然后,数字必须在1379×/π返回/*For(一个数字以外的1379个到达右边)我们完成:所有后续的Prm将是复合*//;现在,所有的标准都满足了,如果我们用较小的t替换p[i]并重新排序其余的索引到p[i] */* p=CONAT的左边([SETIMUUS(set(p [1…i]),[t])/*重新排序索引(t(t)到左边*/y,t,p[i+.n.n]);/*然后附加t(=新P[i]),和其他索引向右(不变)*/yr返回(1)^ / *结束为(i);也结束NeXPrimf():如果所有排列都完成*/*,L=列表()/*列表结果*/,f/*“找到”*/,p/*排列*/d,d/*以增加顺序*/y)/*结束“局部”变量和程序;现在跟随主循环*//;FoVEC(i=向量(n,j,[7 ^(j=n),10)] /*选择数字:索引1对应于数字“0”、“2”、…、“7”、“9”;只生成递增向量i,而第n=n个索引必须是>7,即数字1,3或7或9 */*,vEcSUM(D=VECUSED(D),i)/*对应于索引列表i的*/*)%(3××)的下一个/*和数字必须不可被3 */整除;f f=0 / *迄今为止没有发现任何*;;p=[1…n] /*初始排列=ID*/;直到()NEXPURME()/*在完成循环体之后,转到下一个相关排列;如果所有都完成,则结束。*(/ p)/*(VP(D,P)/*以排列顺序P**/*U/*转换为数字*//**素数?如果我们已经找到另一个素数:Fo & &nf(2)/*,如果我们已经找到了另一个素数,则下一个数的选择*/yf*f=p/*另一个集合“找到“等于当前置换*/*))/*结束()(*)/*;f f/*,如果我们已经找到了一个素数*/&d&[f](1)] *第一个数字非零*/&&&ListPoT(L,VEC摘取物(D,F)*U)/*对应于列表*/*的数,1/*只考虑增加的指数*/*);集(L)/*以递增顺序*/*,/*-OR(n=1)返回:返回1位素数(2,3,5,7)*/*素数(4)/*相当于[2,3,5,7],但OEIS在每个“之后”插入一个空间,因此显式表取更多的空间*/*)/*结束,如果n>1*/}返回n个数字的所有项。- M. F. Hasler,JUL 01 2018