#伊格莱西亚斯·克里斯特(Iglesias Z Krist)中的等式(6)。221 (2006)
b:=过程(p,q)
局部d;
a:=0;
对于d从1到min(p,q)do
如果modp(p,d)=0且modp(q,d)=0,则
ph:=地板(p/2/d);
qh:=楼层(q/2/d);
a:=a+numtheory[mobius](d)*二项式(ph+qh,ph);
结束if;
结束do:
a;
结束进程:
#伊格莱西亚斯·克里斯特(Iglesias Z Krist)中的等式(17)。221 (2006)
bt:=进程(p,q)
如果类型(p+q,“奇数”),则
b(p,q);
其他的
0;
结束条件:;
结束进程:
#更正了Iglesias Z Krist中的等式(15)。221(2006),d(p/2)和d(q/2)
bbtemp:=进程(p,q)
局部d、ph、qh;
a:=0;
对于d从1到min(p,q)do
如果modp(p,2*d)=0且modp(q,2*d)=0,则
ph:=p/2/d;
qh:=q/2/d;
a:=a+numtheory[mobius](d)*二项式(ph+qh,ph);
结束条件:;
结束do:
a;
结束进程:
#伊格莱西亚斯·克里斯特(Iglesias Z Krist)中的等式(16)。221 (2006)
bb:=进程(p,q)
如果type(p,‘even’)和type(q,‘even'),则
(bbtemp(p,q)-bt(p/2,q/2))/2;
其他的
0 ;
结束条件:;
结束进程:
#伊格莱西亚斯·克里斯特(Iglesias Z Krist)中的等式(25)。221 (2006)
压裂反应:=过程(Phalf)
如果类型(Phalf,“even”),则
其他的
0;
结束条件:;
结束进程:
#伊格莱西亚斯·克里斯特(Iglesias Z Krist)中的等式(24)。221 (2006)
局部a、p、q、p;
P:=n;
a:=FracR(P/2);
对于从0到P的q do
p:=p-q;
如果modp(p-q,3)=0且p<q,则
a:=a+bb(p,q);
结束条件:;
结束do:
a;
结束进程: