登录
a(1)=a(2)=2;此后,a(n)=gpf(1+a(n-1)a(n-2)),其中gpf是最大素因子。
1

%I#24 2018年8月3日15:44:53

%S 2,2,5,11,7,13,23,5,29,73353859842392761856717126948763,

%电话:122613509210218446725300792659346387041447334098941853251,

%电话:148080802321968921649521449033686346008003007716781536311893316247302477109839787560145369313122241776027171580811013839862377046165714073

%N a(1)=a(2)=2;此后,a(n)=gpf(1+a(n-1)a(n-2)),其中gpf是最大素因子。

%tα={2,2};执行[AppendTo[a,FactorInteger[1+a[[n-1]]a[[n-2]][[-1,1]],{n,3,26}];a(*迈克尔·德弗里格,2015年8月5日*)

%tnxt[{a_,b_}]:={b,因子整数[a*b+1][[-1,1]]};NestList[nxt,{2,2},30][[All,1]](*哈维·P·戴尔,2018年8月3日*)

%o(鼠尾草)

%o定义gpf(n):

%o返回(系数(n)[-1])[0]

%o定义A259828vec(m):#m>2=f

%o f=2

%o v=[2,2]

%o对于范围(f,m)内的i:

%o v附加(gpf(1+v[i-1]*v[i-2])

%o返回v

%o(PARI)gpf(n)=我的(f=系数(n)[,1]);f[#f];

%o第一个(m)=我的(v=向量(m));v[1]=2;v[2]=2;对于(i=3,m,v[i]=gpf(1+v[i-1]*v[i-2]));v;

%Y参考A260725、A000945。

%K nonn公司

%O 1,1号机组

%A _Anders Hellström_,2015年8月5日

%E 2015年8月5日_Michael De Vlieger_的条款a(20)-a(26)