swust oj 1188 计算系数

 给定一个多项式 (ax + by)k ，请求出多项式展开后 x n y m 项的系数。

输出共 1 行，包含一个整数，表示所求的系数，这个系数可能很大， 输出对 10007 取模后的结果。

展开式的公式应该还是记得住，C（n,k）*a^(n-k)*b^k

要对10007取模，在求a的次方的时候不能直接使用pow函数因为会爆掉，所以老老实实乘法取模，至于求组合数，数据的范围是1000，打表既可；

#include<stdio.h>
#include<string.h>
#include<math.h>
int c[2000][2000];
void fun()
{
    int i,j;
    memset(c,0,sizeof(c));
    for(i=0;i<=1005;i++)
     c[0][i]=c[1][i]=1;
    for(i=0;i<=1005;i++)
     c[i][i]=1;
    for(i=0;i<=1005;i++)
     c[i][0]=1;
    for(i=1;i<=1005;i++)
    {
        for(j=1;j<=i;j++)
        {
            if(i!=j)
            c[i][j]=(c[i-1][j]+c[i-1][j-1])%10007;
        }
    }
}
int pow1(int a,int n)
{
    int sum=1;
    for(int i=0;i<n;i++)
    {
        sum=sum*a;
        sum=sum%10007;
    }
    return sum;
}
int main()
{
   int a,b,m,k,n;
   fun();
   while(scanf("%d%d%d%d%d",&a,&b,&k,&n,&m)!=EOF)
   {
       a=a%10007;
       b=b%10007;
      int sum=c[k][n];
      sum=sum%10007;
      sum=sum*pow1(a,n);
      sum=sum%10007;
      sum=sum*pow1(b,m);
      sum=sum%10007;
      printf("%d\n",sum);

   }
return 0;
}