#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#include<bits/stdc++.h>
#include<atcoder/all>
using namespace std;
using namespace atcoder; 
using ll=long long;
using mint=modint998244353;
const ll inf=(1e9+7)*(1e9+7);
int dx[4]={1,0,0,-1},dy[4]={0,1,-1,0};
const int mod=998244353;
using P=pair<ll,int>;
mint dp[5001][5002],fact[5002];
int main(){
	int n,m;
	cin>>n>>m;
	vector<int>a(n);
	for(int i=0;i<n;i++)
		cin>>a[i];
	int cnt=a[0]+a[n-1];
	for(int i=1;i<n;i++)
		if(a[i]!=a[i-1])
			cnt++;
	fact[0]=1;
	for(int i=1;i<=n+1;i++)
		fact[i]=fact[i-1]*i;
	mint comb=fact[n+1]/(fact[cnt]*fact[n+1-cnt]);
	int memo1[5002],memo2[5002],memo3[5003];
	for(int j=0;j<=n+1;j+=2){
		if(j>=2)
			memo1[j]=j*(j-1)/2;
		memo2[j]=j*(n+1-j);
		if(j+2<=n+1)
			memo3[j]=(n+1-j)*(n+1-j-1)/2;
	}
	dp[0][0]=1;
	for(int i=0;i<m;i++){
		for(int j=0;j<=n+1;j+=2){
			if(j>=2)
				dp[i+1][j-2]+=dp[i][j]*memo1[j];
			dp[i+1][j]+=dp[i][j]*memo2[j];
			if(j+2<=n+1)
				dp[i+1][j+2]+=dp[i][j]*memo3[j];
		}
	}
	cout<<(dp[m][cnt]/comb).val()<<endl;
}