#include<bits/stdc++.h>
#define fo(i,l,r) for(int i=(l);i<=(r);++i)
#define fd(i,l,r) for(int i=(l);i>=(r);--i)
#define fu(i,l,r) for(int i=(l);i<(r);++i)
#define ci const int
#define ll long long
using namespace std;
ci mod=998244353,N=4e5+5;
int fac[N],dfac[N],inv[N];
inline void add(int &x,ci v){
	x+=v,x-=x<mod?0:mod;
}
inline void sub(int &x,ci v){
	x-=v,x+=x<0?mod:0;
}
ll qk(ll x,int y=mod-2){
	ll ans=1;
	while(y){
		if(y&1)ans=ans*x%mod;
		x=x*x%mod,y>>=1;
	}
	return ans;
}
ll C(ci n,ci m){
	if(n<m||m<0)return 0;
	return (ll)fac[n]*dfac[m]%mod*dfac[n-m]%mod;
}
int n,k,w[8],co[8],cnt[8];
int trs[4][4][N],ed[4][N],cur[8],tmp[8];
int dp[N][4];
void Pre(){
	for(int S=0;S<1<<k-2;++S){
		memset(cur,0,sizeof(cur));
		cur[S]=1;
		for(int i=1;i<=n;++i){
			memset(tmp,0,sizeof(tmp));
			for(int s=0;s<1<<k-2;++s)
				for(int t=0;t<1<<k-1;++t){
					int nx=s|t;
					if(nx&1)add(tmp[nx>>1],(ll)cur[s]*co[t]%mod);
					else add(trs[S][nx>>1][i],(ll)cur[s]*co[t]%mod*dfac[i]%mod);
					if(nx==0)add(ed[S][i],(ll)cur[s]*co[t]%mod*dfac[i]%mod);
					if(nx==1)add(ed[S][i],(ll)cur[s]*co[t]%mod*dfac[i+1]%mod);
					if(nx==2)add(ed[S][i],(ll)cur[s]*co[t]%mod*dfac[i]%mod*dfac[2]%mod);
					if(nx==3)add(ed[S][i],(ll)cur[s]*co[t]%mod*dfac[i+2]%mod);
					if(nx==4)add(ed[S][i],(ll)cur[s]*co[t]%mod*dfac[i]%mod*dfac[2]%mod);
					if(nx==5)add(ed[S][i],(ll)cur[s]*co[t]%mod*dfac[i+1]%mod*dfac[2]%mod);
					if(nx==6)add(ed[S][i],(ll)cur[s]*co[t]%mod*dfac[i]%mod*dfac[3]%mod);
					if(nx==7)add(ed[S][i],(ll)cur[s]*co[t]%mod*dfac[i+3]%mod);
				}
			memcpy(cur,tmp,sizeof(cur));
		}
	}
}
ci g=3,ig=qk(3);
int R[N];
void NTT(int *a,ci n,ci tp){
	for(int i=0;i<n;++i)R[i]=(R[i>>1]>>1)|((i&1)?(n>>1):0);
	for(int i=0;i<n;++i)if(i<R[i])swap(a[i],a[R[i]]);
	for(int len=2;len<=n;len<<=1){
		ci gn=qk(tp==1?g:ig,(mod-1)/len);
		for(int l=0;l<n;l+=len)
			for(int k=l,al,ar,cur=1;k<l+(len>>1);++k)
				al=a[k],ar=a[k+(len>>1)],
				a[k]=(al+(ll)ar*cur)%mod,
				a[k+(len>>1)]=(al-(ll)ar*cur%mod+mod)%mod,
				cur=(ll)cur*gn%mod;
	}
	if(tp==-1){
		ci inv=qk(n);
		for(int i=0;i<n;++i)a[i]=(ll)a[i]*inv%mod;
	}
}
namespace poly{
	const int N=2097153,mo=998244353;
	int n,m,a[N],b[N],G[N],Ginv[N],rev[N];
	ll ksm(ll x,int y){ll s=1;for(;y;y>>=1,x=x*x%mo)if(y&1)s=s*x%mo;return s;}
	void init(int n)
	{
		for(int i=1;i<n;i<<=1)
		{
			int g0=ksm(3,(mo-1)/(i<<1)),g0i=ksm(g0,mo-2),g=1,gi=1;
			fu(j,i,i+i) G[j]=g,Ginv[j]=gi,g=1ull*g*g0%mo,gi=1ull*gi*g0i%mo;
		}
	}
	void ntt(int n,int *a,int pd)
	{
		fu(i,0,n) if(i<rev[i]) swap(a[i],a[rev[i]]);
		for(int i=1,x;i<n;i<<=1)
			for(int *j=a;j<a+n;j+=(i<<1))
			{
				for(int *k=j,*buf=(pd?G:Ginv)+i;k<j+i;++k,++buf)
				{
					x=1ull*k[i]*(*buf)%mo;
					if((k[i]=*k+mo-x)>=mo) k[i]-=mo;
					if((*k+=x)>=mo) *k-=mo;
				}
			}
		if(!pd)
		{
			int invn=ksm(n,mo-2);
			fu(i,0,n) a[i]=1ull*a[i]*invn%mo;
		}
	}
	int gt(int deg){
		int n=1,s=0;while(n<deg) n<<=1,s++;
		fu(i,0,n) rev[i]=(rev[i>>1]>>1)|((i&1)<<(s-1));
		init(n);
		return n;
	}
	void nt(int *f,ci n){
		ntt(n,f,1);
	}
	void nt2(int *f,int *g,int *h,ci n){
		fu(i,0,n) h[i]=1ull*f[i]*g[i]%mo;
		ntt(n,h,0);
	}
	// void NTT(int n1,int n2,int deg,int *f,int *g)
	// {
	// 	int n=1,s=0;while(n<deg) n<<=1,s++;
	// 	init(n);
	// 	for(int i=n1;i<n;++i)f[i]=0;
	// 	for(int i=n2;i<n;++i)g[i]=0;
	// 	fu(i,0,n) rev[i]=(rev[i>>1]>>1)|((i&1)<<(s-1));
	// 	ntt(n,f,1);
	// 	ntt(n,g,1);
	// 	fu(i,0,n) f[i]=1ull*f[i]*g[i]%mo;
	// 	ntt(n,f,0);
	// }
}
int A[4][N],B[4][N],c[N];
void Solve(ci l,ci r){
	if(l==r){
		if(l==1)++dp[1][0];
		return;
	}
	ci mid=l+r>>1;
	Solve(l,mid);
	if(r-l<=150){
		for(int i=l;i<=mid;++i)
			for(int s=0;s<1<<k-2;++s)
				if(dp[i][s])
					for(int j=mid+1;j<=r;++j)
						for(int t=0;t<1<<k-2;++t)
							add(dp[j][t],(ll)dp[i][s]*trs[s][t][j-i]%mod);
	}
	else{
		int n1=mid-l+1,n2=r-l;
		int G=poly::gt(n1+n2-1);
		for(int s=0;s<1<<k-2;++s){
			int n1=0;
			for(int i=0;i<G;++i)A[s][i]=0;
			for(int i=l;i<=mid;++i)A[s][n1++]=dp[i][s];
			poly::nt(A[s],G);
		}
		for(int s=0;s<1<<k-2;++s)
			for(int t=0;t<1<<k-2;++t){
				int n2=0;
				for(int i=0;i<G;++i)B[t][i]=0;
				for(int i=1;i<=r-l;++i)B[t][n2++]=trs[s][t][i];
				poly::nt(B[t],G);
				poly::nt2(A[s],B[t],c,G);
				for(int i=mid+1;i<=r;++i)add(dp[i][t],c[i-l-1]);
			}
		// for(int s=0;s<1<<k-2;++s)
		// 	for(int t=0;t<1<<k-2;++t){
		// 		int n1=0,n2=0;
		// 		for(int i=l;i<=mid;++i)A[n1++]=dp[i][s];
		// 		for(int i=1;i<=r-l;++i)B[n2++]=trs[s][t][i];
		// 		poly::NTT(n1,n2,n1+n2-1,A,B);
		// 		for(int i=mid+1;i<=r;++i)add(dp[i][t],A[i-l-1]);
		// 	}
	}
	Solve(mid+1,r);
}
int main(){
	fac[0]=fac[1]=dfac[0]=dfac[1]=inv[1]=1;
	for(int i=2;i<N;++i)
		fac[i]=(ll)fac[i-1]*i%mod,
		inv[i]=(ll)(mod-mod/i)*inv[mod%i]%mod,
		dfac[i]=(ll)dfac[i-1]*inv[i]%mod;
	scanf("%d%d",&n,&k);
	for(int s=0;s<1<<k-1;++s)scanf("%d",&w[s]),cnt[s]=cnt[s>>1]^(s&1);
	int st=clock();
	for(int A=0;A<1<<k-1;++A){
		ci B=((1<<k-1)-1)^A;		
		for(int t=B;;t=(t-1)&B){
			(cnt[t]?sub:add)(co[t|A],w[A]);
			if(!t)break;
		}
	}
	Pre();
	int ans=0;
	n=n-k+1;
	Solve(1,n);
	n=n+k-1;
	for(int i=1;i<=n-k+1;++i)for(int s=0;s<1<<k-2;++s)add(ans,(ll)dp[i][s]*ed[s][n-k+1-i+1]%mod);
	ans=(ll)ans*fac[n]%mod;
	printf("%d",ans);
	return 0;
}