%I#10 2020年10月18日02:49:59
%S 2,5,12,29,7819950813553592958922575270579194228539961,
%电话:15076024228745119139403369044395181182272003821776082524,
%电话:22198231756363074656
%A022837第N个低点的N指数。
%C序列中的“低点”是指小于前一个术语的术语(对于初始术语,跳过此条件),然后是两个或更多的增加。
%H Rémy Sigrist,适用于A324787的PARI程序</a>
%p Riecaman:=过程(a,s,M)
%p#从s开始,加上或减去a[n],得到M项。如果a有w项,可以得到M=w+1项。
%p局部b,M2,n,t;
%p如果whattype(a)<>list,则返回ERROR(“第一个参数应该是列表”);fi;
%p如果a[1]=0,则错误(“a[1]不应为零”);fi;
%p M2:=最小值(nops(a),M-1);
%p b:=[s];t: =秒;
%p表示n从1到M2 do
%p如果a[n]>t,则t:=t+a[n],否则t:=t-a[n';fi;b: =[运算(b),t];日期:
%p b;结束;
%p块:=proc(a,S)局部b,c,d,M,L,n;
%p#给定一个列表a,其前导项具有索引S,return[b,c,d],其中b列出a中低点的索引,c列出a在低点的值,d列出低点之间的游程长度。
%p b:=[];c: =[];d: =[];五十: =1;
%p#如果[1]是低点?
%p n:=1;
%p如果((a[n+1]>a[n])和(a[n+2]>a[n+1])),则
%p b:=[op(b),n+S-1];c: =[op(c),a[n]];d: =[op(d),n-L];五十: =n;fi;
%n的p从2到nops(a)-2 do
%p#如果a[n]是低点?
%p如果((a[n-1]>a[n])和(a[n+1]>a[n])以及(a[n+2]>a[n+1]),则
%p b:=[op(b),n+S-1];c: =[op(c),a[n]];d: =[op(d),n-L];五十: =n;fi;od;
%p[b,c,d];结束;
%p0:=[seq(ithprime(n),n=1..100001)]:
%p q1:=里卡曼(p0,1100000):
%p块(q1,0);#产生[当前序列,A324788、A324789]
%o(PARI)见链接部分。
%Y参见A022837、A324788、A324889、A324790。
%Y如果基本序列(A022837)以0开始,而不是1,我们将得到A008348、A309226、A324782、A324883、A309225。
%K nonn,更多
%0、1
%A _N.J.A.Sloane,2019年9月4日
%E修改定义,使偏移量为0.-_N.J.A.Sloane,2019年10月2日
%E a(12)-a(22),来自Rémy Sigrist,2020年10月18日
|