#include <bits/stdc++.h>
using namespace std;
const int N=20100,Mod=10003;
long long fact[N]={1},in[N]={0,1},res;
long long C(int n,int m){return fact[n]*in[fact[m]]*in[fact[n-m]]%Mod;}
long long P(int n,int m){return fact[n]*in[fact[n-m]]%Mod;}
int main(){
	int a,b,c,d,k;cin>>a>>b>>c>>d>>k;
	for(int i=1;i<=1000;i++) fact[i]=fact[i-1]*i%Mod;
	for(int i=2;i<=1000;i++) in[i]=(Mod-Mod/i)*in[Mod%i]%Mod;
	for(int i=0;i<=min(a,b);i++) res+=C(b,i)*P(a,i)%Mod*C(d,k-i)*P(a+c-i,k-i)%Mod;
	return cout<<res<<endl,0;
}