OEIS哀悼西蒙斯感谢西蒙斯基金会支持包括OEIS在内的许多科学分支的研究。
登录
OEIS由OEIS基金会的许多慷慨捐赠者.

 

标志
提示
(来自的问候整数序列在线百科全书!)
搜索: a059788-编号:a059798
显示找到的12个结果中的1-10个。 第页12
    排序:相关性|参考文献||被改进的|创建     格式:长的|短的|数据
A075579号 的副本A059788号. +20
0
3, 5, 7, 13, 19, 23, 31, 37, 43, 53, 61, 73, 79, 83, 89, 103, 113, 113, 131, 139, 139, 157 (列表;图表;参考;;历史;文本;内部格式)
抵消
0,1
链接
关键词
死去的
状态
经核准的
A164368号 素数p具有如下性质:如果q是最小素数>p/2,则在p和2q之间存在一个素数。 +10
35
2、11、17、29、41、47、59、67、71、97、101、107、109、127、137、149、151、167、179、181、191、197、227、229、233、239、241、263、269、281、283、307、311、347、349、367、373、401、409、419、431、433、439、461、487、491、503、521、569、571、587、593、599、601、607 (列表;图表;参考;;历史;文本;内部格式)
抵消
1,1
评论
Ramanujan素数具有以下性质:
如果p=素数(n)>2,则所有数字(p+1)/2,(p+3)/2。。。,(素数(n+1)-1)/2是复合的。
序列等于具有此属性的所有素数,无论是否为Ramanujan。
所有Ramanujan素数A104272号都在序列中。
双素数的每一个较小值(A001359号)从11开始,在序列中-弗拉基米尔·舍维列夫2009年8月31日
109是这个序列中的第一个非Ramanujan素数。
生成项的一个非常简单的筛选如下:设p_0=1,对于n>=1,p_n是第n个素数。考虑形式(2p_n,2p_{n+1})的连续区间,n=0,1,2,。。。从每个包含至少一个素数的区间中,我们删除最后一个素值。然后所有剩余的素数构成序列。让我们演示一下这个筛子:对于p_n=1,2,3,5,7,11,。。。考虑区间(2,4)、(4,6)、(6,10)、(10,14)、(14,22)、(22,26)、(26,34)、。从所有素数集中去掉每个区间的最后一个素数,即3,5,7,13,19,23,31,。。。我们得到了2,11,17,29等-弗拉基米尔·舍维列夫,2011年8月30日
这个序列和A194598号是相互包装的序列:
A194598号(1) <=a(1)<=A194598号(2) <=a(2)<=。。。
发件人彼得·穆恩2017年10月30日:(开始)
序列是素数p=prime(k)的列表,因此在素数(k)/2和素数(k+1)/2之间没有素数。将“k”更改为“k-1”,因此将“k+1”更改为A164333号s:它的不同之处在于前缀是一个初始术语3。由此我们得到了一个(n+1)=预素数(1964年(n) )=A151799号(A164333号(n) )。
序列是最大素数不小于2*素数(k)的素数列表,因此,作为集合,它是相对于A000040美元中的一组数字A059788号.
{{2},A166252号A166307号}是一个分区。
(结束)
链接
阿洛伊斯·海因茨,n=1..20000时的n,a(n)表
V.Shevelev,关于含素数的临界小区间,arXiv:0908.2319[math.NT],2009年。
V.Shevelev,Ramanujan和Labos素数及其推广和素数分类,arXiv:0909.0715[math.NT],2009-2011。
V.Shevelev,Ramanujan和Labos素数及其推广和素数分类,J.整数序列。15(2012)第12.5.4条
J.Sondow、J.W.Nicholson和T.D.Noe,Ramanujan Primes:边界、跑动、双杀和空档,arXiv:1105.2249[math.NT],2011(见第5节主要差距)。
J.Sondow、J.W.Nicholson和T.D.Noe,Ramanujan Primes:边界、跑动、双杀和空档,J.整数序列。14(2011)第11.6.2条(见第5节主要差距)。
公式
作为一个集合,这个序列=A000040美元\A059788号=A000040美元\前置词(2*A000040美元) =A000040美元\A151799号(A005843号(A000040美元)). -彼得·穆恩2017年10月30日
例子
2在序列中,因为q=2,所以在2和4之间有一个素数3-N.J.A.斯隆2009年10月15日
MAPLE公司
a: =proc(n)选项记忆;局部q,k,p;
k: =nextprime(‘if’(n=1,1,a(n-1));
do q:=下一个质数(楼层(k/2));
p: =nextprime(k);
如果p<2*q,则打破fi;
k: =p
od;k个
结束时间:
seq(a(n),n=1..55)#阿洛伊斯·海因茨,2011年8月30日
数学
收割[Do[q=NextPrime[p/2];如果[PrimePi[2*q]!=PrimePi[p],Sow[p]],{p,Prime[Range[100]}]][[2,1]]
(*第二个节目:*)
fQ[n_]:=PrimePi[2NextPrime[2]]!=PrimePi[n];
选择[Prime@Range@105,fQ]
黄体脂酮素
(PARI)是(n)=下一素数(n+1)<2*nextprime(n/2)&&isp素数(n)\\查尔斯·格里特豪斯四世2015年4月24日
交叉参考
参考Ramanujan素数,A104272号和相关序列:A164288号A080359号A164294号A193507型A194184号A194186号.
囊性纤维变性。A059788号A164333号A194598号.
囊性纤维变性。A001262号A001567号A062568号A141232号(都与基2的伪素数有关)。
关键词
非n
作者
扩展
定义由澄清和简化乔纳森·桑多2011年10月25日
状态
经核准的
A020900型 最大k,使得k次素数小于2次n次素数。 +10
16
2, 3, 4, 6, 8, 9, 11, 12, 14, 16, 18, 21, 22, 23, 24, 27, 30, 30, 32, 34, 34, 37, 38, 40, 44, 46, 46, 47, 47, 48, 54, 55, 58, 59, 62, 62, 65, 66, 67, 68, 71, 72, 75, 76, 77, 78, 82, 86, 87, 88, 90, 91, 92, 95, 97, 99, 99, 100, 101, 102, 103, 106, 112 (列表;图表;参考;;历史;文本;内部格式)
抵消
1,1
链接
公式
a(n)=A000720号(A100484号(n) )-米歇尔·马库斯2020年2月2日
例子
四素数是7,四素数的两倍是14,小于14的最大素数是13,即第六素数,因此a(4)=6-伯纳德·肖特2020年2月2日
数学
PrimePi[NextPrime[#,-1]]&/@(2素数[Range[70]])(*哈维·P·戴尔2012年7月5日*)
黄体脂酮素
(PARI)a(n)=素数(2*素数(n))\\米歇尔·马库斯2017年10月25日;2020年2月2日
交叉参考
囊性纤维变性。A000720号(π(n)),A100484号(2*prime(n))。
囊性纤维变性。A020901号A059788号A060265型.
囊性纤维变性。A102820号(第一个区别)。
关键词
非n
作者
状态
经核准的
A102820号 2*素数(n)和2*素数(n+1)之间的素数,其中素数(n)是第n个素数。 +10
9
1, 1, 2, 2, 1, 2, 1, 2, 2, 2, 3, 1, 1, 1, 3, 3, 0, 2, 2, 0, 3, 1, 2, 4, 2, 0, 1, 0, 1, 6, 1, 3, 1, 3, 0, 3, 1, 1, 1, 3, 1, 3, 1, 1, 1, 4, 4, 1, 1, 2, 1, 1, 3, 2, 2, 0, 1, 1, 1, 1, 3, 6, 2, 0, 1, 6, 1, 3, 0, 1, 1, 3, 2, 2, 1, 2, 1, 1, 2, 4, 1, 3, 1, 1, 2, 1, 2, 1, 0, 1, 4, 2, 1, 3, 0, 2, 5, 0, 5, 3, 3, 2, 1, 0, 2 (列表;图表;参考;;历史;文本;内部格式)
抵消
1、3
评论
连续偶数半素数之间的素数。[尤里·斯捷潘·格拉西莫夫2010年5月1日]
发件人彼得·穆恩,2023年6月1日:(开始)
的第一个差异A020900型.
A080192号列出与零项相对应的素数(n)。
A104380号(k) 是与k作为项的第一次出现相对应的素数(n)。
如果a(n)非零,A059786号(n) 是最小的A059788号(n+1)a(n)枚举素数中最大的。在中描述的素数树中A290183号,这些素数标记素数(n)的子节点。
猜想:0,1,…的渐近比例,k的。。。是1/3,2/9,2^k/3^(k+1)。
(结束)
链接
莱因哈德·祖姆凯勒(Reinhard Zumkeller),n=1..10000时的n,a(n)表
V.Shevelev,关于含素数的临界小区间,arXiv:0908.2319[math.NT],2009年。
公式
a(n)=A020900型(n+1)-A020900型(n) -彼得·穆恩,2023年6月1日
例子
a(15)=3,因为在第15个和第16个素数的两倍之间有3个素数,即介于2*47和2*53之间。
数学
表[PrimePi[2素数[n+1]]-PrimePi[2素数[n]],{n,150}](*扎克·塞多夫*)
差异[PrimePi[2 Prime[Range[110]]](*哈维·P·戴尔2022年10月29日*)
黄体脂酮素
(哈斯克尔)
a102820 n=a102820_list!!(n-1)
a102820_list=地图(总和。(地图a010051))$
zip带enumFromTo a100484_list(尾部a100484_list)
--莱因哈德·祖姆凯勒2012年4月29日
(PARI)a(n)=素数(2*素数(n+1))-素数\\米歇尔·马库斯2017年9月22日
交叉参考
囊性纤维变性。A104272号A080359号. [弗拉基米尔·舍维列夫,2009年8月24日]
囊性纤维变性。A100484号A010051型.
具有类似定义的序列:A104289号A217564型.
关键词
容易的非n
作者
Ali A.Tanara(Tanara(AT)khayam.ut.ac.ir),2005年2月27日
扩展
来自的更多条款扎克·塞多夫2005年2月28日
状态
经核准的
A174635号 不是Ramanujan素数的素数。 +10
8
3, 5, 7, 13, 19, 23, 31, 37, 43, 53, 61, 73, 79, 83, 89, 103, 109, 113, 131, 137, 139, 157, 163, 173, 191, 193, 197, 199, 211, 223, 251, 257, 271, 277, 283, 293, 313, 317, 331, 337, 353, 359, 379, 383, 389, 397, 421, 443, 449, 457, 463, 467, 479, 499, 509, 521 (列表;图表;参考;;历史;文本;内部格式)
抵消
1,1
评论
的补语A104272号在素数中。与不一样A059788号.
也称为非Ramanujan Primes-约翰·尼克尔森2012年1月29日
链接
多诺万·约翰逊,n=1..10000时的n,a(n)表
J.Sondow,Ramanujan素数与Bertrand公设,arXiv:0907.5232[math.NT],2009-201;阿默尔。数学。《月刊》第116期(2009年),第630-635页-约翰·尼克尔森2012年1月29日
J.Sondow、J.W.Nicholson和T.D.Noe,Ramanujan Primes:边界、跑动、双杀和空档,arXiv:1105.2249[math.NT],2011;J.整数序列。14(2011)第11.6.2条-约翰·尼克尔森2012年1月29日。
数学
nn=100;R=表[0,{nn}];s=0;
Do[If[PrimeQ[k],s++];如果[PrimeQ[k/2],s--];如果[s<nn,R[[s+1]]=k],{k,素数[3 nn]}
];
R=R+1;
补码[Prime[Range[PrimePi[Last[R]]],R](*Jean-François Alcover公司2018年11月5日,之后T.D.诺伊在里面A104272号*)
黄体脂酮素
(Perl)使用理论“:all”;my@n=grep{!is_ramanujan_prime($_)}@{primes(1e3)};说“[@n]”#达娜·雅各布森2016年7月15日
(Perl)使用理论“:all”;我的%r$r{$_}=1表示@{ramanujan_primes(1e7)};比如grep{!exists$r{$}}@{primes(1e7)}#达娜·雅各布森2016年7月15日
交叉参考
囊性纤维变性。A104272号.
关键词
非n
作者
T.D.诺伊2010年11月29日
状态
经核准的
A164917号 通过应用映射x->达到素数(n)的最小步数A060308型(x) 从的任何成员开始A164368号. +10
5
0,1,2,3,0,4,0,1,5,0,1,2,0,6,0,1,0,2,0,3,1,7,1,0,0,2,0,0,1,1,1,0,0,2,8,0,2,0,0,2,1,2,0,0,2,0,0,0,3,1,1,0,1,1,0,0,2,0,3,0,1,2,3,1,1,0,0,2 (列表;图表;参考;;历史;文本;内部格式)
抵消
1,3
评论
从一些素数开始,迭代应用A060308型(或同等产品A059788号)生成递增质数链。
这些链的性质是在素数列表中达到更高的位置,有时会跳过一些更近的素数,几乎是可能的两倍,从而“过度满足”伯特兰的假设。另一方面,A164368号包含素数,这些素数将被包含略高于其值一半的素数的链跳过。该序列显示了从某个成员开始的链向上的距离A164368号我们可以找到素数(n),或者等价地,我们需要多少映射的逆应用才能到达其中的一个成员A164368号如果从质数(n)开始。
注意,通过构造A164368号(k) 从最小素数开始,该素数不是从任何先前开始的任何链的成员A164368号因此,每个素数都存在于这些链中的某个位置,并且从链的起点到达它的步骤a(n)的数量是明确定义的。
链接
V.Shevelev,关于含素数的临界小区间,arXiv:0908.2319[math.NT],2009年。
例子
映射的第一个素链A060308型使用的成员初始化A164368号
2->3->5->7->13->23->43->83->163->317->631->1259->2503->..
11->19->37->73->139->277->547->1093->2179->4357->8713->17419->。。
17->31->61->113->223->443->883->1759->3517->7027->14051->28099->..
29->53->103->199->397->787->1571->3137->6271->12541->25073->。。
41->79->157->313->619->1237->2473->4943->9883->19763->39521->..
47->89->173->337->673->1327->2647->5281->10559->21107->..
表示前4个素数的a(1)到a(4)都在第一条链上,从2开始需要0到3个步骤=A164368号(1). a(5)询问以下步骤的数量A000040美元(5) =11,位于第二条链上,需要0步。
MAPLE公司
A060308型:=proc(n)前缀(2*n+1);结束时间:
isA164368:=进程(p)局部q;q:=下一个质数(楼层(p/2));返回(数值[pi](2*q)-数值[pi](p)>=1);结束时间:
A164368号:=proc(n)选项记住;局部a;如果n=1,则为2;否则a:=下一素数(进程名(n-1));当不是A164368(a)时,执行a:=下一个质数(a);od:返回(a);fi;结束时间:
A164917号:=proc(n)局部p,a,j,q,itr;p:=i素数(n);a:=100000000000000;对于1中的j,do q:=A164368号(j) ;如果q>p,则断裂;fi;itr:=0;而q<p do q:=A060308型(q) ;itr:=itr+1;od;如果q=p,那么如果itr<a,那么a:=itr;fi;fi;od:a;结束时间:
序列(A164917号(n) ,n=1..120)#R.J.马塔尔2009年9月24日
数学
A060308型[n]:=下一素数[2*n+1,-1];
isA164368[p_]:=模块[{q},q=NextPrime[Floor[p/2];返回[PrimePi[2*q]-PrimePi[p]>=1]];
A164368号[编号]:=A164368号[n] =模[{a},如果[n==1,2,a=NextPrime[A164368号[n-1]];而[Not@isA164368[a],a=NextPrime[a]];返回[a]]];
A164917号[n_]:=模[{p,a,j,q,itr},p=素数[n];a=10^15;对于[j=1,True,j++,q=A164368号[j] ;如果[q>p,则中断[]];itr=0;而[q<p,q=A060308型[q] ;itr++];如果[q==p,如果[itr<a,a=itr]];a] ;
表[A164917号[n] ,{n,1120}](*Jean-François Alcover公司2017年12月14日,之后R.J.马塔尔*)
交叉参考
囊性纤维变性。A006992号A104272号A164368号A164288号.
关键词
非n
作者
扩展
编辑,添加和扩展示例R.J.马塔尔2009年9月24日
状态
经核准的
A080191号 素数p,使得p是2*p之前的素数和下一个素数之间的所有素数因子中最大的。 +10
4
2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 61, 67, 73, 79, 83, 89, 97, 103, 109, 113, 127, 131, 137, 139, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 269, 271, 277, 281, 283, 293, 307, 313 (列表;图表;参考;;历史;文本;内部格式)
抵消
1,1
评论
的补语A080192号相对于A000040美元.
发件人弗拉维奥·弗尔南德斯2021年5月26日:(开始)
等价地,素数p使得p是区间[2*p,nextprime(2*p)-1]中数字的所有素数因子中最大的。
对于任何素数p,如果p不是该区间内数字的所有素数因子中最大的(即,如果p并非该序列的项),则该区间内的数字的所有素因子中最大者将是出现在数字2*q中的素数q。
对于所有n,最大素数<2*a(n)是A059788号.(结束)
链接
阿米拉姆·埃尔达尔,n=1..10000时的n,a(n)表
公式
f(precprime(2*p))=p,其中f是由A052248号.
例子
5是一个项,因为7是2*5之前的素数,11是下一个素数,5是8、9和10的所有素数因子中最大的。
数学
选择[Range[300],PrimeQ[#]&&NextPrime[2*#]<2*NextPrime[#]&](*阿米拉姆·埃尔达尔,2020年2月7日*)
黄体脂酮素
(PARI){对于素数(k=2317,p=前素数(2*k);q=下素数(p+1);m=0;对于(j=p+1,q-1,f=因子(j);a=f[矩阵大小(f)[1],1];如果(m<a,m=a);如果(m==k,print1(k,“,”))}
交叉参考
囊性纤维变性。A000040美元A080192号A052248号.
关键词
非n
作者
状态
经核准的
A290183号 a(n)是标记具有素数的根树中高度为n的节点的最小p(A000040美元)作为其非根节点的标签,标记为p的节点的父节点是标记为最大素数小于p/2的节点(如果没有这样的节点,则为根节点)。 +10
4
59, 331, 163, 79, 37, 17, 724777, 941593, 3578683, 73397147, 52457707, 26228849, 2600553509, 1300276753 (列表;图表;参考;;历史;文本;内部格式)
抵消
0,1
评论
标记为p的节点的父节点是与根节点相连的相邻节点。
上面定义的“素数树”与许多较古老的序列有关。标记的节点A055377号(n) 是标记为n的节点的父节点。标记为prime(k)的节点具有A102820号(k) 子节点,除非它没有子节点,否则这些节点都用来自A059786号(k) 至A059788美元(k+1)。叶节点标签为A080192号.深度m的节点是间隔中带有标签的节点[A055496号(m) ,A055496号(m+1))。完整树可以使用A000040美元以表格形式读取,行长度由A102820号前缀为2。
如果节点的高度集有一个最大的有限值k,那么这个序列是有限的,有k+1个连续的定义项。
也就是说,作者的初步评估是,高度n节点的出现与素数k元组的最小素数的出现具有相似性,即:(1)高度n节点标签的出现几乎与素数之间的随机间隔相似;(2) 对于任意n,质数p是这样一个标号的明显概率不会大于log(p)中的多项式;因此,(3)a(n)对所有n来说似乎都存在。
一些初步的经验观察表明,高度n+1的节点出现的频率可能比高度n的节点低5-10倍。
术语a(1)到a(5)来自由标记为17的节点及其后代组成的子树,如下面的示例部分所示。这意味着4个连续的负第一差异,这在序列的后面可能是罕见的。
链接
米歇尔·马库斯和彼得·穆恩,PARI脚本
公式
a(n)=最小值(A055377号^(n+1)\A055377号^(n+2)),其中A055377号^(i) 表示应用后的图像集NA055377号.-澄清了以下建议米歇尔·马库斯2018年11月13日
例子
2*331=662和2*nextprime(331)=2*337=674之间的唯一素数是673,因此标记为331的节点仅是标记为673的节点的父节点。在2*673=1346和2*nextprime(673)=2*677=1354之间没有素数,因此标记为673的节点是叶子。因此,标记为331的节点具有高度1。它是最小的素数,所以a(1)=331。
上面定义的“素数树”开始于:
+----------------------------+
2 3
| |
5 7
+--------------+ +-------+
11 13 17 19
| +-----------+ | +----+
23 29 31 37 41 43
+----+ +----+ +----+----+ | | |
47 53 59 61 67 71 73 79 83 89
由标记为17的节点及其所有5个子节点组成的子树为:
17
|
37
|
79
|
163
|
331
|
673
发件人彼得·穆恩2018年11月22日:(开始)
由标记为a(7)=941593的节点及其后代组成的子树为:
941593
+------+--------+
1883191 1883197
| |
3766387 3766397
|
7532813
+-----------------+-----------------+
15065627 15065639 15065641
| +-----------+ |
30131267 30131279 30131281 30131291
|
60262597
|
120525217
(结束)
黄体脂酮素
(PARI)列表a(nbp)={my(v=素数(nbp,nv,x);对于(n=0,oo,nv=集合(应用(x->precprime(x>>1),v));x=vecmin(Set减号(v,nv)));如果(x>vecmax(nv),中断);打印1(x,“,”);v=nv;);}/*使用nbp=6*10^6表示7个术语,nbp=2*10^7表示8个术语*/\\米歇尔·马库斯2018年11月19日
(PARI)时间(p,n)=(k=1,n,p=前一素数(p>>1));p;
a(n)={p=2;x=时间(p,n);p=下一素数(p+1);q=5;y=时间(q,n+1);q=下一质数(q+1);对于(k=1,oo,if(y!=x,return(x));直到(ny=ntimes(q,n+1))!=x,q=下一素(q+1\\米歇尔·马库斯,2018年12月8日
交叉参考
囊性纤维变性。A080192号是叶节点的列表。
囊性纤维变性。A055377号A059786号A059788美元、和A102820美元与父/子节点关系关联。
囊性纤维变性。A000040美元A055496号A151800型.
关键词
非n更多
作者
彼得·穆恩2017年7月23日
扩展
(8)候选人确认米歇尔·马库斯2018年11月14日
a(9)-a(13)来自米歇尔·马库斯2019年1月2日
状态
经核准的
A370008飞机 a(n)是小于3*prime(n)的最大素数。 +10
5, 7, 13, 19, 31, 37, 47, 53, 67, 83, 89, 109, 113, 127, 139, 157, 173, 181, 199, 211, 211, 233, 241, 263, 283, 293, 307, 317, 317, 337, 379, 389, 409, 409, 443, 449, 467, 487, 499, 509, 523, 541, 571, 577, 587, 593, 631, 661, 677, 683, 691, 709, 719, 751 (列表;图表;参考;;历史;文本;内部格式)
抵消
1,1
链接
例子
5<3*2<7<3*3<11<13<3*5,所以(a(1),a(2),a,a(3))=(5,7,13)。
数学
表[素数[PrimePi[3*Prime[n]]],{n,1200}]
黄体脂酮素
(PARI)a(n)=预素数(3*素数(n))\\米歇尔·马库斯2024年2月10日
交叉参考
关键词
非n
作者
克拉克·金伯利2024年2月9日
状态
经核准的
A370009飞机 a(n)是大于3*prime(n)的最小素数。 +10
1
7, 11, 17, 23, 37, 41, 53, 59, 71, 89, 97, 113, 127, 131, 149, 163, 179, 191, 211, 223, 223, 239, 251, 269, 293, 307, 311, 331, 331, 347, 383, 397, 419, 419, 449, 457, 479, 491, 503, 521, 541, 547, 577, 587, 593, 599, 641, 673, 683, 691, 701, 719, 727, 757 (列表;图表;参考;;历史;文本;内部格式)
抵消
1,1
链接
公式
a(n)=A151800型(A001748号(n) )。
例子
5<3*2<7<3*3<11<13<3*5<17,所以(a(1),a(2),a。
数学
表[NextPrime[Prime[Pi[3*Prime[n]]],{n,1200}]
黄体脂酮素
(PARI)a(n)=下一素数(3*素数(n))\\米歇尔·马库斯2024年2月10日
交叉参考
关键词
非n
作者
克拉克·金伯利2024年2月9日
状态
经核准的
第页12

搜索在0.009秒内完成

查找|欢迎光临|维基|注册|音乐|地块2|演示|索引|浏览|更多|网络摄像头
贡献新序列。或评论|格式|样式表|变换|超级搜索|最近
OEIS社区|维护人OEIS基金会。

许可协议、使用条款、隐私政策。.

上次修改时间:美国东部夏令时2024年5月28日12:54。包含372913个序列。(在oeis4上运行。)