登录
OEIS由OEIS基金会的许多慷慨捐赠者.

 

标志
提示
（来自的问候整数序列在线百科全书!)
搜索： a336923-编号：a336922
显示找到的9个结果中的1-9个。 第页1
    排序：关联|参考文献|数|被改进的|创建     格式：长的|短的|数据
A359579型 的Dirichlet逆A336923飞机，其中A336923飞机（n） =1，如果σ（2n）-σ（n）是2的幂，否则为0。 +20
三
1, -1, -1, 0, 0, 1, -1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, -1, 0, 0, 1, 0, 0, 0, 0, 0, -1, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, -1 (列表;图表;参考;听;历史;文本;内部格式)
抵消
1
链接
Antti Karttunen，n=1..107163的n，a（n）表
配方奶粉
a（1）=1，对于n>1，a（n）=-和{d|n，d<n}A336923飞机（n/d）*a（d）。
如果e=1，与a（2^e）=-1相乘，如果e>1，与0相乘，对于素数p>2，a（p^e）=(-A036987号（p） ）^e.-更正人阿米拉姆·埃尔达尔和安蒂·卡图恩2023年1月24日
对于所有n>=1，abs（a(A056652号（n） ）=abs（a（2）*A056652号（n） ））＝1。
对于所有n>=1，abs（a(A219174型（n） ）=1，如果A219174型（n） 不是4的倍数。
数学
f[p_，e_]：=如果[2^整数指数[p+1，2]==p+1，（-1）^e，0]；f[2，e_]：=如果[e==1，-1，0]；a[1]=1；a[n_]：=倍@@f@@FactorInteger[n]；数组[a，100](*阿米拉姆·埃尔达尔2023年1月24日*）
黄体脂酮素
（PARI）
A209229型（n） =（n&&！比特（n，n-1））；
A359579型（n） ={my（f=因子（n））；prod（k=1，#f~，如果（2==f[k，1]，-（1==f[k，2]）(-A209229型（1+f[k，1]））^f[k、2]）；}；
交叉参考
囊性纤维变性。A036987号,A056652号,A219174型,A336923飞机.
另请参阅A359578型.
关键词
签名,复数
作者
安蒂·卡图恩2023年1月8日
状态
经核准的
A331410型 a（n）是从n开始并使用映射k->k+k/p达到2次幂所需的迭代次数，其中p是k的最大素因子。 +10
54
0, 0, 1, 0, 2, 1, 1, 0, 2, 2, 2, 1, 2, 1, 3, 0, 3, 2, 3, 2, 2, 2, 2, 1, 4, 2, 3, 1, 4, 3, 1, 0, 3, 3, 3, 2, 4, 3, 3, 2, 3, 2, 3, 2, 4, 2, 2, 1, 2, 4, 4, 2, 4, 3, 4, 1, 4, 4, 4, 3, 2, 1, 3, 0, 4, 3, 4, 3, 3, 3, 3, 2, 5, 4, 5, 3, 3, 3, 3, 2, 4, 3, 3, 2, 5, 3, 5, 2, 5, 4, 3, 2, 2, 2, 5, 1, 3, 2, 4, 4, 5, 4, 3, 2, 4 (列表;图表;参考;听;历史;文本;内部格式)
抵消
1,5
评论
设f（n）=A000265号（n） 是n的奇数部分。设p是k的最大素因子，设k=p*m。假设k不是2的幂，即p>2，则f（k）=p*f（m）。迭代次数为k->k+k/p=p*m+m=（p+1）*m。因此，p*f（m）->f（p+1。因为对于p>2，f（p+1）<p，每次迭代中的奇数部分都会减少，直到变为1，也就是说，直到我们达到2的幂-阿米拉姆·埃尔达尔2020年2月19日
在迭代的任何一步都可以使用k的任何奇数素数因子，结果都是一样的。因此，就像A329697型，这也是完全加性序列-安蒂·卡图恩2020年4月29日
当且仅当a（n）等于A005087号（n） 则σ（2n）-σ（n）是2的幂。（请参见A336923飞机,A046528号). -安蒂·卡图恩2021年3月16日
链接
Antti Karttunen，n=1..65537时的n，a（n）表
迈克尔·德弗利格，a（n）的带注释扇形二叉树标记索引n并应用颜色代码，其中黑色表示a（n）=0，红色表示a（n）=1，洋红色表示n=1.16383时a（n”）的最大值。
配方奶粉
发件人安蒂·卡图恩2020年4月29日：（开始）
这是一个完全可加序列：对于奇素数p，a（2）=0，a（p）=1+a（p+1），a（m*n）=a（m）+a（n），如果m，n>1。
a（2n）=a(A000265号（n） ）=a（n）。
如果A209229型（n） ==1，a（n）=0，否则a（n+A052126号（n） ），或相等，1+a（n+（n/A078701型（n） ）。
a（n）=A334097型（n）-A334098型（n） ●●●●。
一个(A122111号（n） ）=A334108型（n） ●●●●。
（结束）
a（n）=A334861飞机（n）-A329697型（n） ●●●●-安蒂·卡图恩2020年5月14日
a（n）=a(A336467飞机（n） ）+A087436号（n）=A336921型（n）+A087436号（n） ●●●●-安蒂·卡图恩2021年3月16日
例子
15的轨迹是[15,18,24,32]，需要3次迭代才能达到32。因此，a（15）=3。
数学
a[n_]：=-1+长度@NestWhileList[#+#/FactorInteger[#][[-1，1]]&，n，#/2^IntegerExponent[#，2]！=1 &]; 数组[a，100](*阿米拉姆·埃尔达尔，2020年1月16日*）
黄体脂酮素
（岩浆）f:=func<n|n+n div p，其中p是Max（素数除数（n））>；g： =func<n|n eq 1或Max（素数除数（n））eq 2>；a： =[]；对于[1..1000]中的n，做k:=n；s： =0；而不是g（k）做s：=s+1；k： =f（k）；结束while；附加（~a，s）；结束；a//马吕斯·A·伯蒂2020年1月19日
（PARI）A331410型（n） =如果（！位和（n，n-1），0，1+A331410型（n+（n/vecmax（因子（n）[，1]）））\\安蒂·卡图恩2020年4月29日
（PARI）A331410型（n） ={my（k=0）；while（比特（n，n-1），k++；my（f=因子（n）[，1]）；n+=（n/f[2-（n%2）]）；（k）；}\\安蒂·卡图恩2020年4月29日
（PARI）A331410型（n） ={my（f=因子（n））；和（k=1，#f~，如果（2==f[k，1]，0，f[k+A331410型（1+f[k，1]））；}\\安蒂·卡图恩2020年4月30日
交叉参考
囊性纤维变性。A000265号,A005087号,A006530号（最大素因子），A052126号,A078701型,A087436号,A329662型（记录位置和每个n的首次出现），A334097型,A334098型,A334108型,A334861飞机,A336467飞机,A336921型,A336922飞机,A336923飞机(A046528号).
Cf.阵列A335430型、及其行A335431型,A335882型、以及A335874飞机.
另请参阅A329697型（使用映射k->k-k/p时的类似序列），A335878飞机.
另请参阅A330437型,A335884型,A335885美元,A336362飞机,A336363飞机用于其他类似迭代。
关键词
非n
作者
阿里·萨达2020年1月16日
扩展
数据段扩展至a（105）安蒂·卡图恩2020年4月29日
状态
经核准的
A062731号 2*n的除数之和。 +10
32
3, 7, 12, 15, 18, 28, 24, 31, 39, 42, 36, 60, 42, 56, 72, 63, 54, 91, 60, 90, 96, 84, 72, 124, 93, 98, 120, 120, 90, 168, 96, 127, 144, 126, 144, 195, 114, 140, 168, 186, 126, 224, 132, 180, 234, 168, 144, 252, 171, 217, 216, 210, 162, 280, 216, 248, 240, 210 (列表;图表;参考;听;历史;文本;内部格式)
抵消
1,1
评论
a（n）也是2*n划分为相等部分的所有部分的总数-奥马尔·波尔2021年2月14日
链接
N.J.A.斯隆，n=1..20000时的n，a（n）表[哈里·史密斯的前1000条条款]
配方奶粉
a（n）=A000203号（2*n）-R.J.马塔尔2011年4月6日
a（n）=A000203号（n）+A054785号（n） ●●●●-R.J.马塔尔2020年5月19日
发件人瓦茨拉夫·科特索维奇，2022年8月7日：（开始）
Dirichlet g.f.：zeta（s）*zeta（s-1）*（3-2^（1-s））。
求和{k=1..n}a（k）~5*Pi^2*n^2/24。（结束）
数学
lst={}；执行[AppendTo[lst，DivisorSigma[1，n]]，{n，2，6！，2}]；第一次(*弗拉基米尔·约瑟夫·斯蒂芬·奥尔洛夫斯基2008年9月20日*）
除数Sigma[1,2*范围[60]](*哈维·P·戴尔2022年6月8日*）
黄体脂酮素
（PARI）矢量（66，n，sigma（2*n，1））
（PARI）对于（n=1000，写入（“b062731.txt”，n，“”，西格玛（2*n））\\哈里·史密斯2009年8月9日
（MuPAD）编号：：sigma（2*n）$n=0..81//零入侵拉霍斯2008年5月13日
（岩浆）[SumOfDivisors（2*n）：n in[1..70]]//文森佐·利班迪2014年10月31日
交叉参考
西格玛（k*n）：A000203号（k=1），A144613号（k=3），A193553号（k=4，偶数等分），A283118号（k=5），A224613型（k=6），A283078型（k=7），A283122号（k=8），A283123号（k=9）。
囊性纤维变性。A008438号,A074400型,A182818号,239052英镑（奇数平分），A326124型（部分金额），A054784号,215947英镑,A336923飞机,A346870型,A346878飞机,A346880型,A355750型.
第2行，共行A319526型第2列（&R）A216626型.第1行，共行A355927飞机.
浅对角线（2n，n）A265652型。另请参阅A244658型.
关键词
容易的,非n
作者
杰森·厄尔斯2001年7月11日
扩展
零点移除和偏移校正奥马尔·波尔2009年7月17日
状态
经核准的
A054784号 整数n，使得σ（2n）-σ（n）是2的幂，其中σ是n的除数之和。 +10
9
1, 2, 3, 4, 6, 7, 8, 12, 14, 16, 21, 24, 28, 31, 32, 42, 48, 56, 62, 64, 84, 93, 96, 112, 124, 127, 128, 168, 186, 192, 217, 224, 248, 254, 256, 336, 372, 381, 384, 434, 448, 496, 508, 512, 651, 672, 744, 762, 768, 868, 889, 896, 992, 1016, 1024, 1302, 1344, 1488 (列表;图表;参考;听;历史;文本;内部格式)
抵消
1,2
评论
如果n是梅森素数乘以2的幂的无平方积，那么sigma（2n）-sigma。
反之亦然。这个序列中的所有数字都有这个形式-伊凡·内雷廷2016年8月12日
发件人安蒂·卡图恩，2021年9月1日：（开始）
数字k，使其奇数除数之和[A000593号（k） ]是2的幂。
数字k的奇数部分[A000265号（k） ]在中A046528号.
（结束）
链接
伊万·内雷廷，n，a（n）表，n=1.10000
配方奶粉
数字n是这样的A000203号（2*n）-A000203号（n） 对于一些w，=2^w。
和{n>=1}1/a（n）=2*Product_{p inA000668号}（1+1/p）=2*A306204型= 3.1711177758... . -阿米拉姆·埃尔达尔2023年1月11日
例子
对于n=12，σ（2n）=σ（24）=1+2+3+4+6+8+12+24=60，σ。所以σ（2n）-σ（n）=60-28=32=2^5是2的幂，因此12在序列中-迈克尔·波特2016年8月15日
MAPLE公司
N： =10^6:#获取所有术语<=N
M： =选择（isprime，[seq（2^i-1，i=select（isprim，[$2.ilog2（N+1）]）]）：
R： =映射（t->seq（2^i*t，i=0..floor（log[2]（N/t））），映射（转换，组合：-powerset（M），`*`））：
排序（转换（R，列表））#罗伯特·伊斯雷尔2016年8月12日
数学
排序@选择[外侧压扁[次数，p2=2^范围[0，11]，次数@@#&/@子集@选择[p2-1，PrimeQ]]，#<=最大值@p2 &] (*伊凡·内雷廷2016年8月12日*）
选择[Range[1500]，IntegerQ[Log2[DivisorSigma[1，2#]-DivisorSigma[1，#]]&](*哈维·P·戴尔，2019年4月23日*）
黄体脂酮素
（PARI）
A209229型（n） =（n&&！比特（n，n-1））；
是A054784（n）=A209229型（sigma（n>>估值（n，2））\\安蒂·卡图恩2021年8月28日
交叉参考
囊性纤维变性。A000203号,A000265号,A000396号（甚至术语构成子序列），A000593号,A000668号,A046528号,A063883号,A209229型,A306204型,A331410型,A336923飞机（特征函数）。
零的位置A336922飞机中0和1的位置A336361型.
另请参阅A003401号.
关键词
非n
作者
拉博斯·埃利默2000年5月22日
状态
经核准的
A336477飞机 a（n）=1，如果规则n-gon可以用尺子（或者更准确地说，是一条没有标记的直尺）和罗盘构造，否则为0。 +10
9
1, 1, 1, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0 (列表;图表;参考;听;历史;文本;内部格式)
抵消
1
评论
对于所有n>=1，a（n）=1=>A295297号（n） =0。
链接
Antti Karttunen，n=1..65537时的n，a（n）表
特征函数的索引项
配方奶粉
a（n）=1，如果A005087号（n） 等于A329697型（n） [即，如果A336469飞机（n） =0]，否则为0。
a（n）=A209229型(A000010号（n） ）-维林·亚涅夫和安蒂·卡图恩2021年3月2日
与a（2^e）=1相乘，对于奇素数p，a（p^e）=A209229型（p-1）如果e=1，如果e>1，则为0-安蒂·卡图恩2023年1月6日
黄体脂酮素
（PARI）
A329697型（n） =如果（！位和（n，n-1），0，1+329697英镑（n-（n/vecmax（因子（n）[，1]）））；
A336477飞机（n） =（ω（n>>估值（n，2））==A329697型（n） ）；
（PARI）
A209229型（n） =（n&&！比特（n，n-1））；
A336477飞机（n） ={my（f=因子（n））；prod（k=1，#f~，（2==f[k，1]||A209229型（f[k，1]-1）*（1==f[k（2]））；}\\安蒂·卡图恩2023年1月6日
交叉参考
的特征函数A003401号.
囊性纤维变性。A000010号,A000265号,A001221号,A005087号,A209229型,A295297号,A329697型,A336469,A359578型（Dirichlet逆）。
另请参阅A336471型,A336923飞机（梅森素数的类似序列）。
关键词
非n,复数
作者
安蒂·卡图恩2020年7月25日
扩展
关键字：mult添加者安蒂·卡图恩2023年1月6日
状态
经核准的
A336922飞机 a（n）=A331410型（n）-A005087号（n） ●●●●。 +10
7
0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 0, 2, 1, 2, 1, 0, 1, 1, 0, 3, 1, 2, 0, 3, 1, 0, 0, 1, 2, 1, 1, 3, 2, 1, 1, 2, 0, 2, 1, 2, 1, 1, 0, 1, 3, 2, 1, 3, 2, 2, 0, 2, 3, 3, 1, 1, 0, 1, 0, 2, 1, 3, 2, 1, 1, 2, 1, 4, 3, 3, 2, 1, 1, 2, 1, 3, 2, 2, 0, 3, 2, 3, 1, 4, 2, 1, 1, 0, 1, 3, 0, 2, 1, 2, 3, 4, 2, 2, 1, 1 (列表;图表;参考;听;历史;文本;内部格式)
抵消
1,17
链接
Antti Karttunen，n=1..65537的n，a（n）表
配方奶粉
a（n）=A331410型（n）-A005087号（n） ●●●●。
a（n）=A336921型（n）+A046660号(A000265号（n） ）。
黄体脂酮素
（PARI）
A005087号（n） =ω（n>>估值（n，2））；
A331410型（n） =如果（！位和（n，n-1），0，1+A331410型（n+（n/vecmax（因子（n）[，1]）））；
A336922飞机（n） =(A331410型（n）-A005087号（n） ）；
交叉参考
囊性纤维变性。A000265号,A005087号,A054784号（零位置），A331410型,A336467美元,A336921型,A336923飞机.
另请参阅A336469飞机.
关键词
非n
作者
安蒂·卡图恩2020年8月9日
状态
经核准的
A336391型 词法上最早的无限序列，使得a（i）=a（j）=>A331410型（i）=A331410型（j） 和A336158型（i）=A336158型（j） ，对于所有i，j>=1。 +10
6
1, 1, 2, 1, 3, 2, 2, 1, 4, 3, 3, 2, 3, 2, 5, 1, 6, 4, 6, 3, 7, 3, 3, 2, 8, 3, 9, 2, 10, 5, 2, 1, 5, 6, 5, 4, 10, 6, 5, 3, 6, 7, 6, 3, 11, 3, 3, 2, 4, 8, 12, 3, 10, 9, 12, 2, 12, 10, 10, 5, 3, 2, 13, 1, 12, 5, 10, 6, 5, 5, 6, 4, 14, 10, 15, 6, 5, 5, 6, 3, 16, 6, 6, 7, 17, 6, 17, 3, 14, 11, 5, 3, 7, 3, 17, 2, 6, 4, 11, 8, 14, 12, 6, 3, 18 (列表;图表;参考;听;历史;文本;内部格式)
抵消
1,3
评论
有序对的限制增长序列变换[A331410型（n） ，A336158型（n） ]。
对于所有i，j：
A336390型（i）=A336390型（j） =>a（i）=a（j）
a（i）=a（j）=>A336921型（i）=A336921型（j） ，
a（i）=a（j）=>A336922飞机（i）=A336922飞机（j） =>A336923飞机（i）=A336923飞机（j） ●●●●。
链接
Antti Karttunen，n=1..65537的n，a（n）表
黄体脂酮素
（PARI）
up_to=65537；
rgs_transform（invec）={my（om=Map（），outvec=vector（length（invec）），u=1）；对于（i=1，长度（invesc），如果（mapisdefined（om，invec[i]），my（pp=mapget（om，invec[i]；
A000265号（n） =（n>>估值（n，2））；
A046523号（n） ={my（f=vecsort（factor（n）[，2]，4），p）；prod（i=1，#f，（p=nextprime（p+1））^f[i]）；}；\\发件人A046523号
A336158型（n）=A046523号(A000265号（n） ）；
A331410型（n） =如果（！位和（n，n-1），0，1+A331410型（n+（n/vecmax（因子（n）[，1]）））；
Aux336391（n）=[A331410型（n） ，A336158型（n） ]；
v336391=rgs_transform（向量（up_to，n，Aux336391（n）））；
A336391型（n） =v336391【n】；
交叉参考
囊性纤维变性。A331410型,A336467飞机,A336390型,A336921型,A336922飞机,A336923飞机.
另请参阅A336471型.
关键词
非n
作者
安蒂·卡图恩2020年8月10日
状态
经核准的
A364251型 如果n的形式为q*（2^k），则a（n）=1，其中q是梅森素数之一(A000668号)且k>=0，否则a（n）=0。 +10
4
0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0 (列表;图表;参考;听;历史;文本;内部格式)
抵消
1
链接
Antti Karttunen，n=1..65537时的n，a（n）表
特征函数的索引项
配方奶粉
a（n）=[A331410型（n） ==1]，其中[]是艾弗森括号。
a（n）=[A364260型（n） ==1]。
a（n）<=A336923飞机（n）<=A364252型（n） ●●●●。
黄体脂酮素
（PARI）
isA000668（n）=（i素数（n）&&！比特（n，1+n））；
A364251型（n） =是A000668（n>>估值（n，2））；
交叉参考
的特征函数A335431型.
囊性纤维变性。A000668号,A331410型,A364260型.
另请参阅A336923飞机,A364252型.
关键词
非n
作者
安蒂·卡图恩2023年7月16日
状态
经核准的
A364252型 如果n除了2和/或梅森素数之外没有其他素数因子，则a（n）＝1，否则为0。 +10
三
1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 1 (列表;图表;参考;听;历史;文本;内部格式)
抵消
1
链接
Antti Karttunen，n=1..100000时的n，a（n）表
特征函数的索引项
配方奶粉
a（p^e）=[p==2]的完全乘法+A036987号（p） ，其中[]是艾弗森括号。
a（n）=[A336467飞机（n） ==1]。
a（n）>=A336923飞机（n） >=A364251型（n） ●●●●。
黄体脂酮素
（PARI）
A209229型（n） =（n&&！比特（n，n-1））；
A364252型（n） ={my（f=因子（n））；prod（k=1，#f~，（2==f[k，1]||A209229型（f[k，1]+1）；}；
交叉参考
囊性纤维变性。A036987号,A209229型,A336467飞机,A336923飞机,A364251型.
的特征函数A219174型.
关键词
非n,复数
作者
安蒂·卡图恩2023年7月16日
状态
经核准的
第页1

搜索在0.007秒内完成

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

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

上次修改时间：2024年6月23日15:51 EDT。包含373651个序列。（在oeis4上运行。）