#include<bits/stdc++.h>
#define ll long long
#define int long long
using namespace std;

const int maxn = 1048576<<1;
const int mod = 998244353;
const int g   = 3;
inline void add(int &a,const int &b){ a+=b;if(a>=mod)a-=mod; }
inline void sub(int &a,const int &b){ a-=b;if(a<0)a+=mod; }

int pw(int x,int k)
{
	int re=1;
	for(;k;k>>=1,x=(ll)x*x%mod) if(k&1)
		re=(ll)re*x%mod;
	return re;
}
int inv(int x){ return pw(x,mod-2); }

int id[maxn],wi[maxn];
void DFT(int f[],const int N,const int sig)
{
	for(int i=1;i<N;i++) if(i<id[i]) swap(f[i],f[id[i]]);
	
	for(int m=2;m<=N;m<<=1)
	{
		int t=m>>1,tt=N/m;
		for(int i=0;i<t;i++)
		{
			int wn= sig==1?wi[i*tt]:wi[N-i*tt];
			for(int j=i;j<N;j+=m)
			{
				int tx= f[j],ty=(ll)f[j+t]*wn%mod;
				f[j]=(tx+ty)%mod;
				f[j+t]=(tx-ty+mod)%mod;
			}
		}
	}
	if(sig==-1)
	{
		int iN=inv(N);
		for(int i=0;i<N;i++) f[i]=(ll)f[i]*iN%mod;
	}
}

int fac[maxn],ifac[maxn];
void pre()
{
	fac[0]=1; for(int i=1;i<maxn;i++) fac[i]=(ll)fac[i-1]*i%mod;
	ifac[maxn-1]=inv(fac[maxn-1]);
	for(int i=maxn-2;i>=0;i--) ifac[i]=(ll)ifac[i+1]*(i+1)%mod;
}
int C(int i,int j)
{
	if(i<j||i<0||j<0) return 0;
	return (ll)fac[i]*ifac[j]%mod*ifac[i-j]%mod;
}

int f[maxn];
int n,m,nm;

signed main()
{
	ios_base::sync_with_stdio(false);
	//freopen("tmp.in","r",stdin);
	
	pre();
	
	cin>>n>>m; nm=n*m;
	for(int i=0;i<=m;i++)
	{
		int cc=(ll)fac[m]*ifac[m-i]%mod*(ll)C(m,i)%mod;
		if(i%2==0) add(f[i],cc);
		else sub(f[i],cc);
	}
	
	int N=1,ln=0;
	while(N<=nm) N<<=1,ln++;
	
	for(int i=1;i<N;i++) id[i]=id[i>>1]>>1|(i&1)<<(ln-1);
	wi[0]=1; wi[1]=pw(g,(mod-1)/N);
	for(int i=2;i<=N;i++) wi[i]=(ll)wi[i-1]*wi[1]%mod;
	DFT(f,N,1);
	
	for(int i=0;i<N;i++) f[i]=pw(f[i],n);
	DFT(f,N,-1);
	
	int ans=0;
	for(int i=0;i<=nm;i++)
	{
	//	if(i%2==0) cerr<<i<<" : "<<f[i]<<endl;
	//	else cerr<<i<<" : "<<(mod-f[i])%mod<<endl;
		add(ans,(ll)f[i]*fac[nm-i]%mod);
	}
	cout<<ans<<endl;
	
	return 0;
}