%I#60 2024年5月10日11:08:46
%S 1,0,9,8189,7216012838411523456691210460805529632768,
%电话9830429491288473644236803538944131072002726297673400320,
%电话:419430402548039682264924161324620820384317441819393289069664021474836486424509441932735283257982058496289910292480
%N a(N)是A351048的第一项,有N个素数因子,以重数计算,如果没有该项,则为0。
%C a(n)是第一个数字k(如果有),即A001222(k)=n,A000005(k)除以A003415(k)。
%H Robert Israel,n的表,n=0..1000的a(n)</a>
%对于k>=2,F a(2^k-1)=2^(2^k-1)。
%e a(4)=189,因为189=3^3*7有4个素因子,A000005(189)=8除以A003415(189。
%p儿童:=proc(r)局部L,x,p,q,t,r;
%px:=r[1];
%p L:=r[2];
%p t:=L[-1];
%p p:=t[1];q: =下一素数(p);
%p如果t[2]=1,则t:=[q,1];
%p其他t:=[p,t[2]-1],[q,1]
%p fi;
%p R:=[x*q/p,[op(L[1..-2]),t]];
%p如果nops(L)>=2,则
%p p:=L[-2][1];
%pq:=L[-1][1];
%p如果L[-2][2]=1,则t:=[q,L[-1][2]+1]
%p其他t:=[p,L[-2][2]-1],[q,L[-1][2]+1]
%p fi;
%p R:=R,[x*q/p,[op(L[1..-3]),t]]
%p fi;
%p【R】
%p端程序:
%p f:=进程(n)
%p使用priqueue;
%p局部pq,t,x,V,F;
%p初始化(pq);
%p插入([-2^n,[[2,n]]],pq);
%p do公司
%p V:=提取物(pq);
%px:=-V[1];F: =伏[2];
%p如果(x*add(t[2]/t[1],t=F))mod-mul(t[2]+1,t=F)=0
%p然后返回(x)fi;
%p代表儿童的t(V),插入(t,pq)od;
%p od;
%p端程序:
%p1,0,seq(f(n),n=2..40);
%Y参考A000005、A000122、A003415、A351048。
%K nonn公司
%0、3
%以色列罗贝尔特,2024年5月6日
|