|
（PARI）
α（n）={
如果（n<2，返回（最大值（0，n-1））；
本地（nm1=n-1，
mi=m=天花板（nm1/2），
r=地板（对数（m）/对数（2）），
i、 fi，α=0，a）；
对于步骤（i=1,2*r+1，2，
mi/=2；
fi=（1+2^i）\ 3；
α+=fi*地板（0.5+mi）；
);
α*=2；
如果（nm1%2，\\调整为偶数n
a=系数（2*m）[1,2]-1；
α-=（1+2^（1+2*a））\3；
);
回报（α）；
}
（J）
注意：。向量上的J函数
注意：。注意大参数的舍入错误
注意：。确定到~1e8
alphav=：3:0
n=.<：年
如果.+/ntlo=。n>0 do。
n=。ntlo#n
m=.>.-：n个
r=.<。2^.米
f=.<。3%~2+2^2*>：i.>./>：第页
z=。0
mi=。米
for_i.i.#f do。
z=。z+（i{f）*<0.05+mi=.mi%2
结束。
nzer=。(+/ @: (0=>./\)@:|.)“1@：#：m
ntlo#^：_1 z-（2|n）*<.-：nzer{f
否则。
ntlo公司
结束。
)
注意：。例如alphav 1 3 5 100 2 8 33
| 