组合数学打表程序(矩阵快速幂中也经常用到):
#include<stdio.h>
#define MAXN 15
#define MOD 10000
int c[15][15];
void calc_comb()
{
c[0][0]=1;
for(int i=1;i<MAXN;i++)
{
c[i][0] = c[i][i] = 1;
for(int j=1;j<i;j++)
c[i][j] = (c[i-1][j]+c[i-1][j-1])%MOD;
}
}
int main()
{
int a,b;
calc_comb();
while(~scanf("%d%d",&a,&b))
{
printf("%d\n",c[a][b]);
}
return 0;
}