#include<bits/stdc++.h>
using namespace std;
const int mod=998244353,G=3,GI=(mod+1)/3;
#define vi vector<int>
#define pb push_back
int qpow(int a,int b){
	int c=1;for(;b;b>>=1){
		if(b&1)c=1ll*a*c%mod;
		a=1ll*a*a%mod;
	}
	return c;
}
int L,lb[1<<22];
void init(int n){
	L=1;while(L<=n)L<<=1;
	for(int i=0;i<L;i++)lb[i]=(lb[i>>1]>>1)|((i&1)?(L>>1):0);
}
void NTT(vi &a,int fl){
	for(int i=0;i<L;i++)if(i>lb[i])swap(a[i],a[lb[i]]);
	for(int o=1;o<L;o<<=1){
		int Wn=qpow(fl?G:GI,(mod-1)/(o<<1));
		for(int i=0;i<L;i+=(o<<1))for(int j=0,w=1;j<o;j++,w=1ll*w*Wn%mod){
			int x=a[i+j],y=1ll*w*a[i+j+o]%mod;
			a[i+j]=(x+y)%mod,a[i+j+o]=(x-y)%mod;
		}
	}
	if(!fl){
		int I=qpow(L,mod-2);
		for(int i=0;i<L;i++)a[i]=1ll*a[i]*I%mod;
	}
}
vi mul(vi a,vi b){
	int n=(int)a.size()+(int)b.size()-1;init(n);
	a.resize(L),b.resize(L),NTT(a,1),NTT(b,1);
	for(int i=0;i<L;i++)a[i]=1ll*a[i]*b[i]%mod;
	NTT(a,0);a.resize(n);
	return a;
}
int n,m,jc[1010000],ij[1001000];
int C(int a,int b){return 1ll*jc[a]*ij[b]%mod*ij[a-b]%mod;}
int main(){
    scanf("%d%d",&n,&m);
    jc[0]=1;
    for(int i=1;i<=n*m;i++)jc[i]=1ll*i*jc[i-1]%mod;
    ij[n*m]=qpow(jc[n*m],mod-2);
    for(int i=n*m;i;i--)ij[i-1]=1ll*i*ij[i]%mod;
    vector<int>f,g;
    g={1};
    for(int i=0;i<=m;i++)f.pb((((m-i)&1)?-1ll:1ll)*C(m,i)*ij[i]%mod);
    for(int z=n;z;z>>=1){
        if(z&1)g=mul(g,f);
        f=mul(f,f);
    }
    int ans=0;
    for(int i=0;i<=n*m;i++)(ans+=1ll*jc[i]*g[i]%mod)%=mod;
    ans=1ll*ans*qpow(jc[m],n)%mod;
    return printf("%d",(ans+mod)%mod),0;
}