登录
A251756型
a(0)=0;对于n>0,a(n)是列表中尚未包含的最小整数,其复合公共因子为a(n-1)。
2
0, 4, 8, 12, 6, 18, 9, 27, 36, 16, 20, 10, 30, 15, 45, 54, 24, 28, 14, 42, 21, 63, 72, 32, 40, 44, 22, 66, 33, 99, 81, 90, 48, 52, 26, 78, 39, 117, 108, 56, 60, 50, 25, 75, 100, 64, 68, 34, 102, 51, 153, 126, 70, 35, 105, 84, 76, 38, 114, 57, 171, 135, 120, 80
抵消
0,2
评论
这个序列似乎包含了每个复合数。
这些值被分组到从原点延伸的五条线附近,各自的斜率约为{0.608、0.912、1.22、1.82、2.74}={1、1.5、2、3、4.5}*0.608。(如中所示A098550号这些线不是很直。) -M.F.哈斯勒2014年12月14日
链接
莱因哈德·祖姆凯勒(Reinhard Zumkeller),n=0..10000时的n,a(n)表
David L.Applegate、Hans Havermann、Bob Selcoe、Vladimir Shevelev、N.J.A.Sloane和Reinhard Zumkeller,《黄石置换》,arXiv预印本arXiv:1501.01669,2015和J.国际顺序。 18 (2015) 15.6.7.
数学
g[a_List]:=块[{k=4},而不是[CompositeQ[GCD[a[-1]],k]]]||MemberQ[a,k],k++];追加[a,k]];嵌套[g,{0},63](*L.埃德森·杰弗里2014年12月8日(之后罗伯特·威尔逊v) *)
黄体脂酮素
(PARI)invecn(v,k,x)=对于(i=1,k,如果(v[i]==x,返回(i)));0
alist(n)=局部(v=向量(n),x,g);v[1]=4;对于(k=2,n,x=4;while(invecn(v,k-1,x)(g=gcd(v[k-1],x))==1|isprime(g),x++);v[k]=x);v(v)
(哈斯克尔)
导入数据。列表(删除)
a251756 n=a251756_列表!!(n-1)
a251756_list=0:f 0 a002808_list,其中
f x zs=g zs,其中
g(y:ys)|d==1||a010051'd==1=gys
|否则=y:f y(删除y zs)
其中d=gcd x y
--莱因哈德·祖姆凯勒2014年12月8日
(Python)
从gmpy2导入gcd,is_prime
51756英镑_列表,l,s,b=[0],0,1,{}
对于范围(10**3)内的_:
i=秒
为True时:
如果b中没有i:
m=gcd(i,l)
如果不是(m==1或is_prime(m)):
A251756型_列表.附加(i)
l、 b[i]=i,正确
而s在b中:
b.流行音乐
s+=1
打破
i+=1#柴华武2014年12月8日
关键词
非n
作者
状态
经核准的