#include <cstdio>
#include <vector>
#include <algorithm>
#define ALL(vc) (vc).begin(),(vc).end()
using namespace std;
const int N=1000003,P=998244353;
typedef long long ll;
typedef vector<int> vi;
int n,m,k;
int rev[1<<20],cw[1<<20|1];
int ilen,len,bit;
inline void inc(int &x,int v){if((x+=v)>=P) x-=P;}
inline void dec(int &x,int v){if((x-=v)<0) x+=P;}
int qp(int a,int b=P-2){
	int res=1;
	while(b){
		if(b&1) res=(ll)res*a%P;
		a=(ll)a*a%P;b>>=1;
	}
	return res;
}
void init(int _len){
	len=1;bit=-1;
	while(len<_len) len<<=1,++bit;
	int w=qp(3,(P-1)>>(bit+1));
	cw[len]=cw[0]=1;
	for(int i=1;i<len;++i){
		rev[i]=(rev[i>>1]>>1)|((i&1)<<bit);
		cw[i]=(ll)cw[i-1]*w%P;
	}
	ilen=qp(len);
}
struct poly{
	vi f;
	poly(int Len):f(Len){}
	poly(initializer_list<int> Init):f(Init){}
	void NTT(){
		f.resize(len,0);
		for(int i=0;i<len;++i) if(i<rev[i]) swap(f[i],f[rev[i]]);
		for(int i=1,tt=len>>1;i<len;i<<=1,tt>>=1)
			for(int j=0;j<len;j+=(i<<1))
				for(int k=j,t=0;k<(j|i);++k,t+=tt){
					int x=f[k],y=(ll)f[k|i]*cw[t]%P;
					if((f[k]=x+y)>=P) f[k]-=P;
					if((f[k|i]=x-y)<0) f[k|i]+=P;
				}
	}
	void INTT(){
		for(int i=0;i<len;++i) if(i<rev[i]) swap(f[i],f[rev[i]]);
		for(int i=1,tt=len>>1;i<len;i<<=1,tt>>=1)
			for(int j=0;j<len;j+=(i<<1))
				for(int k=j,t=len;k<(j|i);++k,t-=tt){
					int x=f[k],y=(ll)f[k|i]*cw[t]%P;
					if((f[k]=x+y)>=P) f[k]-=P;
					if((f[k|i]=x-y)<0) f[k|i]+=P;
				}
		for(int i=0;i<len;++i) f[i]=(ll)f[i]*ilen%P;
		while(!f.empty()&&!f.back()) f.pop_back();
	}
};
int fac[N],fiv[N];
int C(int n,int k){return (ll)fac[n]*fiv[k]%P*fiv[n-k]%P;}
int a[N],pw[N],res[N];
int main(){
	// freopen("paint.in","r",stdin);
	// freopen("paint.out","w",stdout);
	scanf("%d%d%d",&n,&m,&k);
	int lim=n*m;
	fac[0]=1;
	for(int i=1;i<=lim;++i) fac[i]=(ll)fac[i-1]*i%P;
	fiv[lim]=qp(fac[lim]);
	for(int i=lim;i;--i) fiv[i-1]=(ll)fiv[i]*i%P;
	pw[0]=1;
	for(int i=1;i<=lim;++i) inc(pw[i]=pw[i-1],pw[i-1]);
	for(int i=0;i<=n;++i)
		for(int j=0;j<=m;++j){
			int pos=i*j+(n-i)*(m-j);
			if((i^j)&1){
				a[pos]=(a[pos]-(ll)C(n,i)*C(m,j))%P;
				if(a[pos]<0) a[pos]+=P;
			}
			else a[pos]=(a[pos]+(ll)C(n,i)*C(m,j))%P;
		}
	poly F(lim+1),G(lim+1);
	init(lim<<1|1);
	for(int i=0;i<=lim;++i) F.f[i]=(ll)fac[i]*a[i]%P;
	for(int i=0;i<=lim;++i)
		if(i&1) G.f[i]=P-fiv[i];
		else G.f[i]=fiv[i];
	reverse(ALL(G.f));
	F.NTT();G.NTT();
	poly H(len);
	for(int i=0;i<len;++i) H.f[i]=(ll)F.f[i]*G.f[i]%P;
	H.INTT();
	int Hlen=H.f.size();
	for(int i=lim;i<Hlen;++i) res[lim+lim-i]=(ll)H.f[i]*fiv[i-lim]%P*pw[i-lim]%P;
	int ans=0;
	for(int i=k;i<=lim;++i) ans=(ans+(ll)C(i,k)*res[i])%P;
	ans=(ll)ans*qp(2,P-1-n-m)%P;
	if(ans&&(n&1)) ans=P-ans;
	printf("%d\n",ans);
	return 0;
}