#附件1
Hra:=proc（r:：integer，a:：intiger，q:：integer）
二项（r-1，a-1）*（q-3）+二项（r-1，a）；
%*（q-3）^（r-a-1）；
结束进程：
Jra:=proc（r:：integer，a:：intiger，q:：integer）
二项式（r-2，a-2）*（q-3）^2+2*二项式；
%*（q-3）^（r-a-2）；
结束进程：
#附件2，I_m
局部q，a，f；
q:=9；
a:=0；
f:=1+（-1）^（r+a）+（1+（-1-）^a）*（1-（-1））^r）*楼层（（q-3）/2）/2；
Jra（r，a，q）+二项式（2，r-a）+f*Hra（floor（r/2），floor（a/2），q）；
%/4 ;
结束进程：
