蒜头君的多项式

#include <iostream>
#include <bits/stdc++.h>
using namespace std;
int C[1001][1001];
int mod = 10007;
void init(){
	for(int i=0;i<=1000;i++){
		C[i][0]=1;
		C[i][i]=1;
	}
	for(int i=2;i<=1000;i++){
		for(int j=1;j<i;j++){
			C[i][j]=(C[i-1][j-1]+C[i-1][j])%mod; // 帕斯卡恒等式 
		}
	}		
}
int main()
{
 	int p,q,k,a,b;
 	int s;
	init();
 	while(cin>>p>>q>>k>>a>>b){ 
 	s = C[k][b];
	 for(int i=1;i<=a;i++){
 		s=(s*(p%mod))%mod;
	 }
	 for(int i=1;i<=b;i++){
 		s=(s*(q%mod))%mod;
	 }
	cout<<s<<'\n';	 
	 }
    return 0;
}