c++ 分情况排列组合计数

#include <cstring>
#include <iostream>

using namespace std;
typedef long long LL;

const int N=2005;
const int mod=100003;
int fact[N], infact[N];

int qmi(int a, int k){
    int ans=1;
    while(k){
        if(k&0x01) ans=(LL)ans*a%mod;
        a=(LL)a*a%mod;
        k>>=1;
    }
    return ans;
}

int C(int a, int b){
    if(a<b) return 0;
    return (LL)fact[a]*infact[a-b]*infact[b]%mod;
}

int P(int a, int b){
    if(a<b) return 0;
    return (LL)fact[a]*infact[a-b]%mod;
}

int main(){
    fact[0]=infact[0]=1;
    for(int i=1; i<N; ++i){
        fact[i]=(LL)fact[i-1]*i%mod;
        infact[i]=(LL)infact[i-1]*qmi(i, mod-2)%mod;
    }
    int a, b, c, d, k;
    cin>>a>>b>>c>>d>>k;
    LL ans=0;
    for(int i=0; i<=k; ++i){
        ans=(ans%mod+(LL)C(b, i)*P(a, i)%mod*(LL)C(d, k-i)*P(a+c-i, k-i)%mod)%mod;
    }
    cout<<ans<<endl;
    return 0;
}