象棋可能出现的摆法总数

C语言不支持大数，所以我用数组来模拟，导致整个程序很混乱，可读性很差。但是C语言中似乎没有什么更好的解决方式，我又不会别的语言，只能勉为其难了。

#include<stdio.h>
#include<math.h>
#define KK for(k=0;k<9;++k)
#define CARRY(a) for(k=0;k<8;++k){if(a[k]>99999){a[k+1]+=a[k]/100000;a[k]%=100000;}}
#define FOR(i,top) for(i=0;i<top;++i)
//摆放顺序：将士象和未过河三七路兵、未过河一五九路兵、过河兵、车马炮 
int c(unsigned int y,unsigned int x){//组合数 
	if(x){
		int t=1,q=1,k; 
		for(k=1;k<=x;++k)
			t*=(y-k+1),q*=k;
		t/=q;
		return t;
	}
	else return 1;
}
void prt(int x){//输出时带0 
	if(x==0)printf(" 00000");
	else if(x<10)printf(" 0000%d",x);
	else if(x<100)printf(" 000%d",x);
	else if(x<1000)printf(" 00%d",x);
	else if(x<10000)printf(" 0%d",x);
	else printf(" %d",x);
}
int main(){
	int k;//k表示位数 
	int i,j;
	int tem[9]={0};//累加用的临时变量
	int s[9]={0};//s[]表示象棋可能出现的摆法总数
	int a1[7][31][9]={0};
	//a1[i][j][]表示在双方已摆放士、相、兵、车、马、炮的数量总和为j，还剩i对车、马、炮情况未知时未知棋子的摆法总数 
	for(j=0;j<31;++j)
		a1[0][j][0]=1;	
	for(i=1;i<7;++i){
		for(j=0;j<31-2*i;++j){
			KK a1[i][j][k]=c(88-j,2)*a1[i-1][j+2][k]+(88-j)*a1[i-1][j+1][k]+a1[i-1][j][k];
			CARRY(a1[i][j])//进位 
		}
	}
	printf("步骤1完成\n");
	
	
	int C[10][6][2]={0};
	/*这个数列存在的意义是拆分组合数，防止100000进制计数法出现紊乱。C[i][j][0]*C[i][j][1]=c(44-i,j) 
	 这里是这个程序最失败的地方，但是我找不到更好的解决办法。*/ 
	for (i=0;i<10;++i){
	for(j=0;j<6;++j){
		if(j<4){			
			C[i][j][0]=c(44-i,j);
			C[i][j][1]=1;
		}
		if(j==4){
			C[i][j][0]=(41-i)*(42-i);
			C[i][j][1]=(43-i)*(44-i);
			for(k=2;k<=4;++k){//这里的k不是位数，而是一个暂时的循环变量
				if (C[i][j][0]%k==0)
					C[i][j][0]/=k;
				else  C[i][j][1]/=k;
			}
		}
		if(j==5){
			C[i][j][0]=(40-i)*(41-i)*(42-i);
			C[i][j][1]=(43-i)*(44-i);
			for(k=2;k<=5;++k){
				if (C[i][j][0]%k==0)
					C[i][j][0]/=k;
				else  C[i][j][1]/=k;
			}
		}	
	}
	}
	printf("步骤2完成\n");
	
	
	int r1,r2,r3,b1,b2,b3;
	int a2[5][5][6][6][9]={0};
	/*a2[r1][b1][r2][b2][]表示表示已摆放红、黑双方的士、相数量之和分别为r1、b1， 
	红、黑双方分别有r2、b2个未过河兵时未知棋子的摆法总数,
	r3、b3分别表示红、黑双方过河兵的数量。*/ 
	FOR(r1,5){FOR(b1,5){FOR(r2,6){FOR(b2,6){
		FOR(r3,6-r2){FOR(b3,6-b2){
			KK tem[k]=a1[r1+r2+r3+b1+b2+b3][6][k];
			
			KK tem[k]*=C[r1+r2][b3][0];CARRY(tem);
			KK tem[k]*=C[r1+r2][b3][1];CARRY(tem);
			KK tem[k]*=C[b1+b2][r3][0];CARRY(tem);
			KK tem[k]*=C[b1+b2][r3][1];CARRY(tem);
			
			KK a2[r1][b1][r2][b2][k]+=tem[k];
			CARRY(a2[r1][b1][r2][b2]);
		}}
	}}}}
	printf("步骤3完成\n");
	
	
	int r4,r5,r6,b4,b5,b6;
	int a3[5][5][3][3][9]={0};
	/*f[r4][b4][r5][b5][]表示已摆放红、黑双方各自的士相数量之和分别为r4、b4，
	红、黑双方未过河三七路兵的数量分别为r5、b5时未知棋子的摆法总数,
	r6、b6分别表示红、黑双方未过河一五九路兵的数量。 */ 
	FOR(r4,5){FOR(b4,5){FOR(r5,6){FOR(b5,6){
		FOR(r6,4){FOR(b6,4){
			KK a3[r4][b4][r5][b5][k]+=c(3,r6)*c(3,b6)*pow(2,r6)*pow(2,b6)*a2[r4][b4][r5+r6][b5+b6][k];
			CARRY(a3[r4][b4][r5][b5]);
		}}
		printf("*");
	}}
	printf("%d		%d\n",r4,b4);
	}}
	printf("步骤4完成\n");
	
	
	int x7;
	int a4[3][3][3]={0};//a[x7][x8][x9][]表示帅、士、相、未过河三七路兵的组合,x7、x8、x9分别表示士、相、未过河三七路兵的数量 
	for(x7=0;x7<3;++x7){
		a4[x7][0][0]=c(5,x7)*(9-x7);
		a4[x7][0][1]=c(5,x7)*(36-4*x7);
		a4[x7][0][2]=c(5,x7)*(36-4*x7);
		a4[x7][1][0]=c(5,x7)*(62-7*x7);
		a4[x7][1][1]=c(5,x7)*(1230-126*x7);
		a4[x7][1][2]=c(5,x7)*(212-24*x7);
		a4[x7][2][0]=c(5,x7)*(183-21*x7);
		a4[x7][2][1]=c(5,x7)*(626-72*x7);
		a4[x7][2][2]=c(5,x7)*(529-61*x7);
	} 
	printf("步骤5完成\n");
	
	
	int r7,r8,r9,b7,b8,b9;//r7、r8、r9、b7、b8、b9依次表示红士、红相、红方未过河三七路兵、黑士、黑相、黑方未过河三七路兵的数量  
	FOR(r7,3){FOR(r8,3){FOR(r9,3){FOR(b7,3){FOR(b8,3){FOR(b9,3){
		KK tem[k]=a3[r7+r8][b7+b8][r9][b9][k];
		KK tem[k]*=a4[r7][r8][r9];CARRY(tem);
		KK tem[k]*=a4[b7][b8][b9];CARRY(tem);
		KK s[k]+=tem[k];CARRY(s);
	}}}}}}
	printf("步骤6完成\n");
	
	
	printf("\n象棋可能出现的摆法总数：\n");
	for(k=8;k>=0;--k)prt(s[k]);
	return 0;