|
| |
|
| A214551型 |
| 里德·凯利序列：a（n）=（a（n-1）+a（n-3））/gcd。 |
|
32
|
|
|
|
1, 1, 1, 2, 3, 4, 3, 2, 3, 2, 2, 5, 7, 9, 14, 3, 4, 9, 4, 2, 11, 15, 17, 28, 43, 60, 22, 65, 25, 47, 112, 137, 184, 37, 174, 179, 216, 65, 244, 115, 36, 70, 37, 73, 143, 180, 253, 36, 6, 259, 295, 301, 80, 75, 376, 57, 44, 105, 54, 49, 22, 38, 87, 109, 147
(列表;图表;参考;听;历史;文本;内部格式)
|
|
|
|
抵消
|
0,4
|
|
|
评论
|
就像Narayana的Cows序列A000930号，除了总和被前项的最大公约数（gcd）除。
这是一个强烈的猜测，8和10在这个序列中缺失，但最好有一个证明！请参见A214321型对于推测值。[我在会谈中经常将此称为“里德·凯利的序列”。]-N.J.A.斯隆2017年2月18日
|
|
|
链接
|
N.J.A.斯隆，激励数序列（谈话视频），2021年3月5日。
|
|
|
配方奶粉
|
看起来，非常粗略地说，a（n）~常数*exp（0.123…*n）-N.J.A.斯隆2012年9月7日。有关更精确的估计，请参阅下一条评论。
如果a（n）^（1/n）收敛，则极限应接近1.126（参见链接）-贝诺伊特·克洛伊特2015年11月8日
|
|
|
例子
|
a（14）=9，a（16）=3，因此a（17）=（9+3）/gcd（9,3）=12/3=4。
a（24）=28，a（26）=60，因此a（27）=（28+60）/gcd（28,60）=88/4=22。
|
|
|
MAPLE公司
|
a： =proc（n）a（n）：=`if`（n<3，1，（a（n-1）+a（n-3））/igcd（a（n-1），a（n-3））结束：
|
|
|
数学
|
t={1，1，1}；做[AppendTo[t，（t[-1]]+t[[-3]]）/GCD[t[[-1]]，t[[-3]]]，{100}]
f[l_List]：=附加[l，（l[[-1]]+l[-3]]）/GCD[l[[-1-]]，l[[-3]]]；嵌套[f，{1，1，1}，62](*罗伯特·威尔逊v，2012年7月23日*）
递归表[｛a[0]＝＝a[1]＝＝a[2]＝＝1，a[n]＝＝（a[n-1]+a[n-3]）/GCD[a[n-1]，a[n-3]]｝，a，｛n，70｝](*哈维·P·戴尔2014年5月6日*）
|
|
|
黄体脂酮素
|
（Perl）
使用bignum；
my@seq=（1,1,1）；
打印“1 1\n2 1\n3 1\n”；
对于（我的$i=3；$i<400；$i++）
{
我的$next=（$seq[$i-1]+$seq[$i-3]）/
gcd（$seq[$i-1]，$seq[$i-3]）；
我的$ind=$i+1；
打印“$ind$next\n”；
推送（@seq，$next）；
}
亚gcd{
my（$x，$y）=@_；
（$x，$y）=（$y，$x%$y）而$y；
返回$x；
}
（哈斯克尔）
a214551 n=a214551_列表！！n个
a214551_list=1:1:1:zipWith f a214551 _ list（删除2 a2145512_list）
其中f u v=（u+v）`div`gcd u v
（鼠尾草）
定义A214551Rec（）：
x、 y，z=1，1，1
收益率x
为True时：
x、 y，z=y，z，（z+x）//gcd（z，x）
收益率x
（PARI）第一（n）=我的（v=向量（n+1））；对于（i=1，min（n，3），v[i]=1）；对于（i=4，v，v[i]=（v[i-1]+v[i-3]）/gcd（v[n-1]，v[i-3））；v（v）\\查尔斯·R·Greathouse IV，2017年6月21日
（Python）
从数学导入gcd
定义缺陷（nn）：
alst=[1，1，1]
对于范围（3，nn+1）中的n：
附加alst
返回alst
|
|
|
交叉参考
|
|
|
|
关键词
|
非n,美好的
|
|
|
作者
|
|
|
|
状态
|
经核准的
|
| |
|
|
