%I#10 2017年12月3日07:30:14
%S 1,1,2,3,4,7,10,11,20,31,42,73110183292473762123519923209,
%电话:51988407136042201156145762593238150863244100394963954,
%电话:103401716730702707089438158708724111467398185546393002203648576675
%N基于斐波那契的递归。
%H E.S.Rowland,<a href=“http://arxiv.org/abs/0710.3217“>自然素数生成递归</A>,arXiv:071.3217[math.NT]
%F a(n)=a(n-1)+a(n-2)+gcd(n,a(n-1))-gcd(n,a(n-2))。
%F a(n)~c*phi^n,其中phi=A001622=(1+sqrt(5))/2是黄金比率,c=0.3434866160389779937344617212678945874922532000472607933856634329169…-_Vaclav Kotesovec_,2017年12月3日
%p A139759:=proc(n)选项记忆;如果n<=1,则为1;否则an_1:=A139759(n-1);an_2:=A139759(n-2);an_1+an_2+gcd(n,an_1)-gcd(n,an_2);fi;结束:序列(A139759(n),n=0..60);#_R.J.Mathar,2008年5月20日
%ta[0]=a[1]=1;a[n]:=a[n]=a[n-1]+a[n-2]+GCD[n,a[n-1]]-GCD[n,a[n-2]];
%t表[a[n],{n,0,60}](*Jean-François Alcover_,2017年12月3日*)
%Y参考A000045。
%K容易,不是
%0、3
%A_Ctibor O.Zizka,2008年5月20日
%E更多来自R.J.Mathar_的条款,2008年5月20日
%E将引用转换为链接-R.J.Mathar_,2009年10月30日
|