RQNOJ659 计算系数

http://www.rqnoj.cn/problem/659

描述

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

格式

输入格式

共一行，包含5个整数，分别为a，b，k，n，m，每两个整数之间用一个空格隔开。

输出格式

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

样例1

样例输入1[复制]

1 1 3 1 2

样例输出1[复制]

3

限制

1s

提示

对于30%的数据，有0 ≤ k ≤ 10；
对于50%的数据，有a = 1, b = 1；
对于100%的数据，有0 ≤ k ≤ 1000，0 ≤ n, m ≤ k，且n+m = k，0 ≤ a，b ≤ 1,000,000.

来源

NOIp2011提高组Day2第一题

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

题解：求杨辉三角第k行（有k个数的那行）第m个，乘上a^n和b^m的快速幂取模。

特殊处理：if(n+m!=k) puts("0"); ///就是k次幂的各项x^n*y^m的n+m都是等于k的，没有不等于k的项，相当于这项系数是0。（但没有这样的数据，我想多了）

数据较大，最好全用long long。

 //#pragma comment(linker, "/STACK:102400000,102400000")
 #include<cstdio>
 #include<cmath>
 #include<iostream>
 #include<cstring>
 #include<algorithm>
 #include<cmath>
 #include<map>
 #include<set>
 #include<stack>
 #include<queue>
 using namespace std;
 #define ll long long
 #define usint unsigned int
 #define mz(array) memset(array, 0, sizeof(array))
 #define minf(array) memset(array, 0x3f, sizeof(array))
 #define REP(i,n) for(int i=0;i<(n);i++)
 #define FOR(i,x,n) for(int i=(x);i<=(n);i++)
 #define RD(x) scanf("%d",&x)
 #define RD2(x,y) scanf("%d%d",&x,&y)
 #define RD3(x,y,z) scanf("%d%d%d",&x,&y,&z)
 #define WN(x) printf("%d\n",x);
 #define RE  freopen("D.in","r",stdin)
 #define WE  freopen("1biao.out","w",stdout)

 const int MOD=;

 int f[][];
 int a,b,k,x,y;
 ll pow_mod(int x,int k){
     ll re=,t=x;
     while(k>){
         if(k&==){re*=t;re%=MOD;}
         t*=t;
         t%=MOD;
         k>>=;
     }
     return re;
 }

 int main() {
     ll i,j;
     mz(f);
     f[][]=;
     for(i=; i<=; i++) {
         f[i][]=;
         for(j=; j<=i; j++)
             f[i][j]=(f[i-][j-]+f[i-][j])%MOD;
     }
     while(scanf("%d%d%d%d%d",&a,&b,&k,&x,&y)!=EOF) {
         if(x+y!=k) puts("");
         else {

             printf("%lld\n",(((((ll)f[k][y]) * pow_mod(a,x)) % MOD) * pow_mod(b,y)) % MOD);
         }
     }
     return ;
 }