|
#Iglesias Z Krist中的eq（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）；
结束条件：；
结束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)
FracR:=过程（Phalf）
如果类型（Phalf，“even”），则
（bb（法尔夫，法尔夫）-A045683号（Phalf））/2；
其他的
0;
结束条件：；
结束进程：
#伊格莱西亚斯·克里斯特（Iglesias Z Krist）中的等式（24）。221 (2006)
A011951型：=进程（n）
局部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；
结束进程：
序列(A011951型（2*n），n=5..40）#R.J.马塔尔2024年4月15日
| 