登录
邮编:304752
从a(1)=1开始,对于n>1,a(n)=前面未包括的a(n-1)的最小除数,否则a。
2
1, 2, 6, 3, 12, 4, 20, 5, 10, 30, 15, 60, 420, 7, 14, 42, 21, 84, 28, 140, 35, 70, 210, 105, 630, 9, 18, 72, 8, 24, 120, 40, 240, 16, 48, 336, 56, 168, 840, 280, 1680, 80, 480, 32, 96, 672, 112, 560, 3360, 160, 960, 64, 192, 1344, 224, 1120, 6720, 320, 1920, 128, 384, 2688, 448, 2240, 13440, 640, 3840, 256, 768, 5376, 896, 4480, 26880, 1280, 7680,512
抵消
1,2
评论
这是的变体A113552号.
发件人迈克尔·德弗利格2018年5月20日:(开始)
在下表中,我们注意到一个循环的子循环为41。
设e=地板(n/8)。如果条件为假,则写倍数k;如果d不在a中,则写除数d的奇偶性。我们可以将a(n)表示为最小四个素数的乘积,如下所示。
n(mod 8)k或d 2 3 5 7
-------------------------------------------
0 5 2 ^(e-2)5 7
1 6 2 ^(e-1)3 5 7
2偶数2^(e-1)5
3 6 2 ^(e-1)3 5
4偶数2^e
5 3 2 ^e 3
6 7 2 ^e 3 7
7偶数2^(e-1)7
推测:
1.所有项只能被最小的4个素数的某些组合整除。
2.幂2^e,正整数e,位于n={1,2,6,29,34,44,52,60,68,…};第一个差异是{1,4,23,5,10,8,8,8,…},然后是8。
3.对于n>41,使得n(mod 8)=4,a(n)=2^((n-4)/8)。
4.对于n>26,所有条件均为偶数。索引{1、4、8、11、14、17、21、24、26}处的奇数项为{1、3、5、15、7、21、35、105、9}。(结束)
链接
Antti Karttunen和Michael De Vlieger,n=1时的n,a(n)表。.16384
例子
在a(27)=18=2*3^2之后,下一项a(28)既不是2*18=2^2*3^ 2,也不是3*18=2*3 ^ 3既是2又是18的除法。但4不除以18,并且4*18=72尚未用于序列中,因此a(28)=72。
MAPLE公司
lim:=60:使用(numtheory):membera:=proc(val)global a,n:local j:对于j从1到n-1 do if(a[j]=val)然后返回true:fi:od:return false:end:a[1]:=1:对于n从2到lim do d:=sort([除数(a[n-1])[]]):s:=true:对于k从1到nops(d)do如果(不是membera(d[k]))),那么a[n]:=d[k]:s:=false:break:fi:od if(s)然后,对于2中的j,如果(非成员(j,d)和非成员(j*a[n-1])),则a[n]:=j*a[n-1]:break:fi:od:fi:od:seq(a[n',n=1..lim)#纳撒尼尔·约翰斯顿,2011年5月10日,最初为A113552号
#第二个Maple项目:
b: =过程(n)是(n=1)结束:
a: =proc(n)选项记忆;局部j,l,i,m;
j: =a(n-1):l:=排序([numtheory[divisors](j)[]]);
对于i到nops(l),如果不是b(l[i]),则执行
则b(l[i]):=真;返回l[i]
光纤;
对于m,而m在l或b(m*j)中做od;
b(m*j):=真;米*j
结束:a(1):=1:
seq(a(n),n=1..100)#阿洛伊斯·海因茨2018年5月22日
数学
f[s_]:=附加[s,d=除数[s[-1]]];如果[Complement[d,s]!={},补码[d,s][[1],k=2;而[Mod[s[-1]],k]==0||MemberQ[s,k*s[-1]],k++];k*s[[-1]]];嵌套[f,{1},60](*罗伯特·威尔逊v,2006年8月20日,最初为A113552号*)
黄体脂酮素
(PARI)
up_to=(2^14)+1;
v304752=矢量(up_to);
m_occurrences=地图();
k=0;prev=1;对于(n=1,up_to,fordiv(prev,d,如果(!mapisdefined(m_occurrences,d),v304752[n]=d;mapput(m_occurrences,d,n);断裂);如果(!v304752[n],m=1;try=prev;while(!(prev%m)|mapisdefined(m_occurrences,try),m++;try=上一个*m);mapput(m_occurrences,v304752[n]=try,n);上一版本=v304752[n]);
邮编:304752(n) =v304752【n】;
(PARI)邮编:304752(n,a=1,list=list(a)/*设置为0以仅获得a(n)*/,U=[])={for(i=2,n,U=setunion(U,[a]);fordiv(a,d,setsearch(U,d)||[a=-d,break]);if(a>0,for(m=2,oo,a%m&&!setsearch(U,m*a)&&(a*=m)&break),a=-a);list&listput(list,a);/*a%2&print f(“a(%d)=%d,”,i,a)*/);if(list,list,a)}\\M.F.哈斯勒2020年12月26日
交叉参考
不同于A113552号首次在n=28时,其中a(28)=72,而A113552号(28) = 90.
关键词
非n
状态
经核准的