%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)