#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
constexpr int dx[]={1,0,0,-1},dy[]={0,1,-1,0};
constexpr int Dx[]={1,1,0,-1,-1,-1,0,1},Dy[]={0,1,1,1,0,-1,-1,-1};
constexpr int mod=998244353,Mod=1e9+7,inf=Mod;
constexpr ll linf=(ll)inf*inf;
typedef pair<ll,int>P;
#define m_p make_pair
struct fastio{
	fastio(){
		cin.tie(0);
		ios::sync_with_stdio(0);
		cout<<fixed<<setprecision(20);
	};
}fio;
template<class T,class U>
bool chmax(T&a,const U&b){
	if(a<b){
		a=b;
		return 1;
	}
	return 0;
}
template<class T,class U>
bool chmin(T&a,const U&b){
	if(a>b){
		a=b;
		return 1;
	}
	return 0;
}
//*
#include<atcoder/all>
using namespace atcoder;
//*/
ll mypow(int x,ll n){
	if(n==0)
		return 1;
	if(n%2==1)
		return mypow(x,n-1)*x%mod;
	ll res=mypow(x,n/2);
	return res*res%mod;
}
ll inv(int x){
	return mypow(x,mod-2);
}
ll dp1[5001][5001],dp2[5001][5001];
int main(){
	int n,k;
	cin>>n>>k;
	int in=inv(n);
	dp1[0][0]=1;
	dp2[0][0]=n;
	for(int i=0;i<k;i++){
		for(int j=0;j<=min(n,i);j++){
			if(j>0)
				(dp1[i+1][j-1]+=dp1[i][j]*in%mod*j)%=mod;
			if(j<n)
				(dp1[i+1][j+1]+=dp1[i][j]*in%mod*(n-j))%=mod;
			if(j>0)
				(dp2[i+1][j-1]+=(dp2[i][j]+dp1[i][j])*in%mod*j)%=mod;
			if(j<n)
				(dp2[i+1][j+1]+=dp2[i][j]*in%mod*(n-j))%=mod;
		}
	}
	int ans=0;
	for(int i=0;i<=n;i++)
		(ans+=dp2[k][i])%=mod;
	cout<<ans<<endl;
}