登录
a(0)=0,a(1)=1;如果a(n-1)+a(n-2)不是半素数,则a(n)=a(n-1)+a(n-3);否则,a(n)是a(n-1)+a(n-2)的最大素因子。
5

%I#19 2017年11月13日10:13:58

%S 0,1,1,2,3,5,8,13,7,20,27,47,37,84,11,19,30,7,37,44,81125103228,

%电话331,4337413951365223388511820033121512482454618706103,

%电话:383486,79113192,61,23,841071911493401635036661678331000183323351667499

%N a(0)=0,a(1)=1;如果a(n-1)+a(n-2)不是半素数,则a(n)=a(n-1)+a(n-3);否则a(n)是a(n-1)+a(n-2)的最大素数。

%C不带半素数的斐波那契数A000045的模拟。

%C序列是无界的吗?(参见有界的对偶序列A214156。)

%H Peter J.C.Moses,n表,a(n)表示n=0..499</a>

%p A214094:=程序(n)

%p选项记忆;

%p如果n<=1,则

%p n;

%p其他

%p a:=进程名(n-1)+进程名(n-2);

%p如果numtheory[bigomega](a)=2,则

%p max(op(数值[因子集](a));

%p其他

%p返回a;

%p end if;

%p end if;

%p end程序:#_R.J.Mathar_,2013年2月18日

%t A214094[0]:=0;

%t A214094[1]:=1;

%t A214094[n_]:=A214094[n]=如果[PrimeOmega[#]==2,最后[Most[Divisors[#]]],#]&[A214094[0-1]+A214094[2]];

%t表[A214094[n],{n,0,99}](*_Peter J.C.Moses_,2013年2月18日*)

%tnxt[{a_,b_}]:={b,如果[PrimeOmega[a+b]==2,因子整数[a+b][-1,1]],a+b]};NestList[nxt,{0,1},70][[All,1]](*哈维·P·戴尔,2017年11月13日*)

%Y参考A001358,A214156。

%K nonn公司

%0、4

%2013年2月16日,A_Vladimir Shevelev