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

 

标志
提示
(来自的问候整数序列在线百科全书!)
A252868型 的无平方版本A252867型. 4
1, 2, 3, 10, 21, 5, 6, 35, 22, 15, 14, 33, 7, 30, 77, 26, 105, 13, 42, 65, 66, 91, 11, 70, 143, 210, 187, 39, 55, 78, 385, 34, 165, 182, 51, 110, 273, 85, 154, 195, 119, 330, 17, 231, 170, 429, 238, 715, 102, 455, 374, 1365, 38, 1155, 442, 57, 770, 663, 95, 462, 1105, 114, 1001, 255 (列表图表参考历史文本内部格式)
抵消
1,2
评论
如果n<=3,则a(n)=n,否则前面没有出现的第一个无平方数与a(n-2)至少有一个公因数,但与a(n-1)没有公因数。无平方数按其在中的出现次数排序A019565号.
它们表示的整数集与A252867型是的,但在表示中使用了无平方数的因式分解。
这是无平方数的排列。[我认为这目前只是一种猜测-N.J.A.斯隆2015年1月10日]
链接
David L.Applegate、Hans Havermann、Bob Selcoe、Vladimir Shevelev、N.J.A.Sloane和Reinhard Zumkeller,黄石公园排列,arXiv预印arXiv:1501.016692015和J.国际顺序。18 (2015) 15.6.7.
配方奶粉
a(n)=A019565号(A252867型(n-1)
数学
(*b)=A019565号,c=A252867型*)
b[n_]:=倍@@Prime[Flatten[Position[#,1]]&[Reverse[IntegerDigits[n,2]];
c[n_]:=c[n]=如果[n<3,n,对于[k=3,True,k++,如果[FreeQ[Array[c,n-1],k],如果[BitAnd[k,c[n-2]]>=1&&BitAnd[k,c[n-1]]==0,返回[k]]];
a[n]:=b[c[n-1]];
数组[a,64](*Jean-François Alcover公司2018年10月3日*)
黄体脂酮素
(PARI)invecn(v,k,x)=对于(i=1,k,如果(v[i]==x,返回(i)));0
无平方(n)=局部(r=1,i=1);while(n>0,if(n%2,r*=prime(i));i++;n=2);第页
alist(n)=局部(v=向量(n,i,i-1),x);对于(k=4,n,x=3;而(invecn(v,k-1,x)||!位和(v[k-2],x)||位和(v[k-1],x),x++);v[k]=x);向量(n,i,无平方(v[i]))
(Python)
从运算符导入mul
从functools导入reduce
从sympy导入质数
定义A019565号(n) :
….如果n>0,则返回reduce(mul,(枚举(bin(n)[:1:-1])中i,v的素数(i+1),如果v==“1”)),否则返回1
A252868型_列表,l1,l2,s,b=[1,2,3],2,1,3,set()
对于范围内的_(10**4):
….i=秒
….为真时:
如果不是(b或i&l1中的i)和i&l2:
............A252868型_列表.附录(A019565号(i) )
…………l2,l1=l1,i
…………b.添加(i)
当s在b:
…………..删除
………….s+=1
…………中断
……..i+=1#柴华武2014年12月25日
交叉参考
关键词
非n
作者
状态
经核准的

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

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

上次修改时间:2024年4月19日11:31 EDT。包含371792个序列。(在oeis4上运行。)