#include<bits/stdc++.h>
using namespace std;
using ll=long long;
using ull=unsigned long long;
const int mod=998244353;
const int N=1<<22;
using arr=int[N+5];
int add(int x,int y){ return (x+=y)>=mod?x-mod:x; }
int sub(int x,int y){ return (x-=y)<0?x+mod:x; }
int red(int &x){ return x+=(x>>31)&mod; }
int ksm(ll x,int tp,int s=1){
	for(;tp;x=x*x%mod,tp>>=1) if(tp&1) s=x*s%mod;
	return s;
}
arr fac,ifac,inv;
void prep(){
	fac[0]=ifac[0]=1;
	for(int i=1;i<=N;i++) fac[i]=1ll*fac[i-1]*i%mod;
	ifac[N]=ksm(fac[N],mod-2);
	for(int i=N;i;i--) ifac[i-1]=1ll*ifac[i]*i%mod;
	for(int i=1;i<=N;i++) inv[i]=1ll*ifac[i]*fac[i-1]%mod;
}
int C(int n,int m){
    return 1ll*fac[n]*ifac[m]%mod*ifac[n-m]%mod;
}
namespace poly{
	int L,iv;
	arr w,rev;
	arr pa,pb;
	#define szf sizeof(int)
	void cl(arr a){ memset(a,0,L*szf); }
	void r_prep(){
		int L=N>>1;
		int val=ksm(3,(mod-1)/N);
		w[L]=1;
		for(int i=1;i<L;i++) w[i+L]=1ll*w[i+L-1]*val%mod;
		for(int i=L-1;i;i--) w[i]=w[i<<1];
	}
	void pre_n(int n){
		L=1; while(L<n) L<<=1;
		iv=mod-(mod-1)/L;
		for(int i=1;i<L;i++) rev[i]=(rev[i>>1]>>1)|(i&1?L>>1:0);
	}
	void FFT(arr t,bool ok=1){
		int x,y;
		for(int i=1;i<L;i++)
			if(i<rev[i]) swap(t[i],t[rev[i]]);
		for(int i=1;i<L;i<<=1)
			for(int j=0;j<L;j+=i<<1)
				for(int l=0;l<i;l++){
					x=t[j+l],y=1ll*t[i+j+l]*w[i+l]%mod;
					t[i+j+l]=sub(x,y),t[j+l]=add(x,y);
				}
		if(!ok){
			reverse(t+1,t+L);
			for(int i=0;i<L;i++) t[i]=1ll*t[i]*iv%mod;
		}
	}
	void NTT(arr a){ FFT(a); }
	void INTT(arr a){ FFT(a,0); }
	void NTT(arr a,arr b){ memcpy(b,a,L*szf); FFT(b); }
	void INTT(arr a,arr b){ memcpy(b,a,L*szf); FFT(b,0); }

	void Mult(arr a,arr b,arr c,int n){
		pre_n(n);
		NTT(a,pa); NTT(b,pb);
		for(int i=0;i<L;i++)
			c[i]=1ll*pa[i]*pb[i]%mod;
		INTT(c);
		fill(c+n,c+L,0);
		cl(pa); cl(pb);
	}
	void Inv(arr a,arr b,int n){
		int i,j,l;
		pre_n(n); l=L;
		cl(b); *b=ksm(*a,mod-2);
		for(i=2;i<=l;i<<=1){
			pre_n(i);
			NTT(a,pa); NTT(b,pb);
			for(j=0;j<i;j++) pa[j]=1ll*pa[j]*pb[j]%mod;
			INTT(pa); fill(pa,pa+(i>>1),0); NTT(pa);
			for(j=0;j<i;j++) pa[j]=1ll*(mod-pa[j])*pb[j]%mod;
			INTT(pa); j=i>>1;
			memcpy(b+j,pa+j,j*szf);
		}
		fill(b+n,b+l,0);
		cl(pa); cl(pb);
	}
	void Inv2(arr a,arr b,arr c,int n){
		static arr d;
		Inv(a,d,n);
		Mult(b,d,c,n*2);
		fill(c+n,c+2*n,0);
		cl(d);
	}
	void diff(arr t,int n){
		for(int i=1;i<n;i++)
			t[i-1]=1ll*i*t[i]%mod;
		t[n-1]=0;
	}
	void inte(arr t,int n){
		for(int i=n-1;i;i--)
			t[i]=1ll*t[i-1]*inv[i]%mod;
		t[0]=0;
	}
	void Ln(arr a,arr b,int n){
		pre_n(n); cl(b);
		memcpy(b,a,n*szf);
		diff(b,n);
		Inv2(a,b,b,n);
		inte(b,n);
	}
	void Exp(arr a,arr b,int n){
		static arr e;
		int i,j,l;
		pre_n(n);l=L;
		cl(b); *b=1;
		for(i=2;i<=l;i<<=1){
			Ln(b,e,i);
			for(j=0;j<i;j++) e[j]=sub(a[j],e[j]);
			Mult(e,b,e,i); j=i>>1;
			memcpy(b+j,e+j,j*szf);
		}
		fill(b+n,b+l,0);
		cl(e);
	}
}
namespace Get_val{
	using namespace poly;
	arr G[20],F[20];
	void pre_g(int k,int l,int p,arr vx){
		int *g=G[k]+p;
		if(!k){
			g[0]=1; g[1]=sub(0,vx[l]);
			return;
		}
		int L=1<<k;
		pre_g(k-1,l,p,vx); pre_g(k-1,l+L/2,p+L,vx);
		Mult(G[k-1]+p,G[k-1]+p+L,g,L);
		g[L]=sub(g[0],1); g[0]=1;
	}
	void solve(int k,int l,int p,arr ans){
		int *f=F[k]+p;
		if(!k){
			ans[l]=f[0];
			return;
		}
		int L=1<<k;
		int *fl=F[k-1]+p,*fr=F[k-1]+p+L;
		int *gl=G[k-1]+p,*gr=G[k-1]+p+L;
		static arr vf;
		pre_n(L);
		NTT(f,vf);
		reverse(gl+1,gl+L); NTT(gl);
		reverse(gr+1,gr+L); NTT(gr);
		for(int i=0;i<L;i++){
			fl[i]=1ll*vf[i]*gr[i]%mod;
			fr[i]=1ll*vf[i]*gl[i]%mod;
		}
		INTT(fl); fill(fl+L/2,fl+L,0);
		INTT(fr); fill(fr+L/2,fr+L,0);
		solve(k-1,l,p,ans);
		solve(k-1,l+L/2,p+L,ans);
	}
	void work(int n,int m,arr f,arr vx,arr vy){
		int up=(int)log2(max(n,m)-1)+1;
		pre_g(up,1,0,vx);
		int *st=F[up];
		static arr ig;
		int L=1<<up+1;
		Inv(G[up],ig,L/2);
		reverse(ig+1,ig+L);
		Mult(f,ig,st,L);
		fill(st+L/2,st+L,0);
		solve(up,1,0,vy);
	}
}
int n,m;
arr f,g;
int main(){
    prep();
    poly::r_prep();
    scanf("%d %d",&n,&m);
    for(int i=0;i<=m;i++){
        f[i]=1ll*C(m,i)*fac[m]%mod*ifac[m-i]%mod;
        if(i&1) f[i]=mod-f[i];
    }
    poly::Ln(f,g,n*m+1);
    for(int i=0;i<=n*m;i++) g[i]=1ll*g[i]*n%mod;
    memset(f,0,sizeof f);
    poly::Exp(g,f,n*m+1);
    int ans=0;
    for(int i=0;i<=n*m;i++)
        ans=(ans+1ll*fac[n*m-i]*f[i])%mod;
    printf("%d\n",ans);
    return 0;
}