计算系数

二项式定理

(ax+by)^k%p=sigma(i=0,k) c[i]*x^i*y^(k-i)    c[i]=a^i*b^(k-i)*C(i,k)

要开long long虽然有模数，但在计算加法时可能会溢出使答案不正确

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const long long maxn=1007;
const long long mod=10007;
long long f[maxn+7][maxn+7],ans;
long long a,b,k,n,m;
void maketable(){
    for(long long i=0;i<=maxn;i++)
        for(long long j=1;j<=maxn;j++)
            f[i][j]=1;
    for(long long i=2;i<=k;i++)
        for(long long j=2;j<=i;j++)
            f[i][j]=(f[i-1][j-1]+f[i-1][j])%mod;     
}
long long quick(long long a,long long b,long long p){
    long long ans=1;
    for( ;b;b>>=1,a=(a*a)%p){
        if(b&1==1){
            ans=(ans*a)%p;
        }
    }
    return ans;
}
int main(){
    cin>>a>>b>>k>>n>>m;
    maketable();
    ans=(quick(a,n,mod)*quick(b,m,mod)*f[k][n+1])%mod;
    cout<<ans<<endl;
    return 0; 
}