【洛谷p1313】计算系数

（%%%hmr）

计算系数【传送门】

算法呀那个标签：

（越来越懒得写辽）（所以今天打算好好写一写）

---

首先（ax+by）^k的计算需要用到二项式定理：

对于（x+y）^k，有第r+1项的系数为：Tr+1=C_n^ra^n-rb^r

这样对于（ax+by）^k而言，第r+1项的系数就为：a^kb^kC_n^ra^n-rb^r

然而这样算，到就爆掉了呢！

显然不能暴算，然鹅实际上，二项式定理中的系数T,我们可以看成神奇的杨辉三角形：

这样复杂度就降下来了呀，所以又半途而废了

直接带代码：

#include<iostream>
#include<cstdio>
#include<algorithm>
#define mo 10007
using namespace std;
int a,b,k,n,m;
int f[][];
int pow(int a,int k){//快速幂
    int ans=;
    while(k){
        if(k&)ans=ans*a%mo;
        a=a*a%mo;k>>=;
    }return ans;
}
int main()
{
    scanf("%d%d%d%d%d",&a,&b,&k,&n,&m);
    a%=mo;b%=mo;f[][]=;
    for(int i=;i<=k;++i){f[i][]=;
        for(int j=;j<=i;++j)
            f[i][j]=(f[i-][j-]+f[i-][j])%mo;
    }
    printf("%d\n",f[k][n]*pow(a,n)%mo*pow(b,m)%mo);
}

beauty？？

end-