#include<bits/stdc++.h>
#define LL long long
#define LLL __int128
#define uint unsigned
#define ldb long double
#define uLL unsigned long long
using namespace std;
typedef vector<int> poly;
const int Mod=998244353,G=3,img=86583718,winv=932051910;// winv = 2^64 mod Mod
const int Lim=64,inv2=(Mod+1)/2,inv3=(Mod+1)/3,rinv=Mod-2;
const LL Mod2=(LL)Mod*Mod;
poly w0={1},W0;vector<poly>Tr;
inline int fadd(int x,const int&y){return (x+=y-2*Mod)>=0?x:x+2*Mod;}
inline int fsub(int x,const int&y){return (x-=y)>=0?x:x+2*Mod;}
inline int reduce(LL x){return (x+(LL)((uint)x*rinv)*Mod)>>32;}
inline int qpow(int x,int y,int z=1){
	for(;y;(y>>=1)&&(x=(LL)x*x%Mod))if(y&1)z=(LL)z*x%Mod;return z;
}
inline int BSGS(int x){
	unordered_map<int,int>mp;
	int S=sqrtl(Mod)+1,j=1,q=1;
	for(int i=0;i<S;++i,j=(LL)j*G%Mod)mp[(LL)j*x%Mod]=i;
	for(int i=0;i<=S;++i,q=(LL)q*j%Mod)if(mp.count(q)&&i*S>=mp[q])return i*S-mp[q];
	return -1;
}
inline poly Add(poly P,poly Q){
    if(P.size()<Q.size())swap(P,Q);
    for(int i=Q.size();i--;)(P[i]+=Q[i])>=Mod&&(P[i]-=Mod);
    return P;
}
inline poly Sub(poly P,poly Q){
    if(P.size()<Q.size())P.resize(Q.size());
    for(int i=Q.size();i--;)(P[i]-=Q[i])<0&&(P[i]+=Mod);
    return P;
}
template<int T>inline void dft(poly&P){
	const int m=P.size(),n=m/T,n2=n/2;
    int L=w0.size();
	if(L*2<n)for(w0.resize(n2);L*2<n;L<<=1){
		const int cr=qpow(G,Mod/(L<<2));
		for(int i=0;i<L;++i)w0[i+L]=(LL)cr*w0[i]%Mod;
	}
	if(W0.size()<w0.size()){
		const int pr=W0.size(),ed=w0.size();
		W0.resize(ed);
		for(int i=pr;i<ed;++i)W0[i]=reduce((LL)winv*w0[i]);
	}
	for(int k=n,st1,st2;st1=k*T,st2=st1>>1,k>1;k>>=1)
		for(int i=0,ie=n/k;i<ie;++i)
			for(int p=i*st1,pe=p+st2;p<pe;++p){
				const int q=p+st2,z=reduce((LL)P[q]*W0[i]);
				P[q]=fsub(P[p],z),P[p]=fadd(P[p],z);
			}
    for(int&i:P)i>=Mod&&(i-=Mod);
}
template<int T>inline void idft(poly&P){
	const int m=P.size(),n=m/T,ni=qpow(n,Mod-2),n2=n/2;
	poly rev(n);
	for(int i=1;i<n;++i)rev[i]=(rev[i>>1]>>1)+(i&1)*n2;
	for(int&i:P)i=(LL)i*ni%Mod;
	for(int i=1;i<n;++i)if(i<rev[i])swap_ranges(P.begin()+T*i,P.begin()+T*(i+1),P.begin()+T*rev[i]);
	for(int i=1;i<n2;++i)swap_ranges(P.begin()+T*i,P.begin()+T*(i+1),P.begin()+T*(n-i));
	dft<T>(P);
	for(int i=1;i<n;++i)if(i<rev[i])swap_ranges(P.begin()+T*i,P.begin()+T*(i+1),P.begin()+T*rev[i]);
}
template<int T>inline void mul_mod(poly&P,poly&Q){
	const int m=P.size(),n=m/T;
	dft<T>(P),dft<T>(Q);
	for(int i=0;i<n;++i){
		vector<uLL>ans(T<<1);
		for(int j=0;j<T;++j){
			for(int k=0;k<T;++k)
				ans[j+k]+=(LL)P[i*T+j]*Q[i*T+k];
			if(!((~j)&7))for(int k=j;k<j+T;++k)
				if(ans[k]>=(Mod2<<3))ans[k]-=Mod2<<3;
		}
		const int c=((i&1)?Mod-w0[i>>1]:w0[i>>1]);
		for(int j=T;j--;)P[i*T+j]=(ans[j]+ans[j+T]%Mod*c)%Mod;
	}
	idft<T>(P);
}
template<size_t... m>inline void solve(index_sequence<m...>,int n,poly&P,poly&Q){
	static void (*ptrs[])(poly&,poly&)={&mul_mod<m+1>...};
	ptrs[n-1](P,Q);
}
inline poly Mul(poly P,poly Q){
    if(P.empty()||Q.empty())return poly();
	int pn=P.size(),qn=Q.size(),rn=pn+qn-1,b=1;
	while((Lim<<b)<rn)++b;
	const int T=((rn-1)>>b)+1,m=T<<b;
	P.resize(m),Q.resize(m);
	solve(make_index_sequence<Lim>{},T,P,Q);
	return P.resize(rn),P;
}
signed main(){
    cin.tie(0)->sync_with_stdio(0);
    int n,m,k,ans=0;
    cin>>n>>m>>k;
    poly frc(n*m+1),ivf(n*m+1);
    frc[0]=1;
    for(int i=1;i<=n*m;++i)frc[i]=1ll*frc[i-1]*i%Mod;
    ivf[n*m]=qpow(frc[n*m],Mod-2);
    for(int i=n*m;i;--i)ivf[i-1]=1ll*ivf[i]*i%Mod;
    const auto binom=[&](const int&n,const int&m){
        return n<m||m<0?0:1ll*frc[n]*ivf[m]%Mod*ivf[n-m]%Mod;
    };
    poly F(n*m),G(n*m);
    for(int i=0;i<=n*m;++i){
        if(k>=i)F[i]=1ll*(i&1?Mod-ivf[i]:ivf[i])*ivf[k-i]%Mod;
        if(n*m>=i+k)G[i]=1ll*ivf[i]*ivf[n*m-k-i]%Mod;
    }
    F=Mul(F,G);
    for(int x=0;x<=n;++x)
        for(int y=0;y<=m;++y){
            const int k=x*(m-y)+(n-x)*y;
            ans=(ans+(x+y&1?Mod-1ll:1ll)*F[k]%Mod*frc[k]%Mod*frc[n*m-k]%Mod*binom(n,x)%Mod*binom(m,y))%Mod;
        }
    for(int i=1;i<=n+m;++i)((ans&1)?ans+=Mod:ans)>>=1;
    return cout<<ans,0;
}