#include<bits/stdc++.h>
using namespace std;
#define int long long
const int M=998244353;
const int q=3;
const int qq=332748118;
int n,m,f[4000005],tr[4000005];
int qpow(int x,int y)
{
	int ans=1;
	while(y)
	{
		if(y&1)ans=ans*x%M;
		x=x*x%M;
		y>>=1;
	}
	return ans;
}
void ntt(int lim,int *f,int typ)
{
	for(int i=0;i<lim;i++)
		tr[i]=(tr[i>>1]>>1)|((i&1)?(lim>>1):0);
	for(int i=0;i<lim;i++)
		if(i<tr[i])swap(f[i],f[tr[i]]);
	for(int o=2,k=1;o<=lim;k<<=1,o<<=1)
	{
		int wn=qpow(typ?q:qq,(M-1)/o);
		for(int i=0;i<lim;i+=o)
		{
			int w=1;
			for(int j=0;j<k;j++)
			{
				int a=f[i+j],b=f[i+j+k];
				f[i+j]=(a+b*w)%M;
				f[i+j+k]=(a-b*w%M+M)%M;
				w=w*wn%M;
			}
		}
	}
	if(typ)return;
	int inv=qpow(lim,M-2);
	for(int i=0;i<lim;i++)
		f[i]=f[i]*inv%M;
}
int jc[4000005],inv[4000005],jcinv[4000005];
signed main()
{
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	cin>>n>>m;
	jc[0]=jcinv[0]=1;
	inv[1]=1;
	for(int i=2;i<=n*m;i++)
		inv[i]=(M-M/i)*inv[M%i]%M;
	for(int i=1;i<=m*n;i++)
		jc[i]=i*jc[i-1]%M,jcinv[i]=inv[i]*jcinv[i-1]%M;
	for(int i=0;i<=m;i++)
		f[i]=jc[m]*jc[m]%M*jcinv[m-i]%M*jcinv[m-i]%M*jcinv[i]%M;
	int lim=1;
	while(lim<n*m+1)lim<<=1;
	ntt(lim,f,1);
	for(int i=0;i<lim;i++)
		f[i]=qpow(f[i],n);
	ntt(lim,f,0);
	int ans=0;
	for(int i=0,t=1;i<=n*m;i++,t=t*(M-1)%M)
	{
		ans=(ans+t*f[i]%M*jc[n*m-i])%M;
	}
	cout<<ans;
}