12pts求调玄关

#include<bits/stdc++.h>
using namespace std;
const int mod=998244351;
int n,m,t[510],dp[50][510][510]={1},fac[510]={1},ans;
string s;
int mp(int x){
	return x&1;
}
signed main(){
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);
	cin>>n>>m>>s;
	s='?'+s;
	for(int i=1;i<=n;i++){
		int x;
		cin>>x;
		t[x]++;
		fac[i]=i*fac[i-1]%mod;
	}
	for(int i=1;i<=n;t[i]+=t[i-1],i++);
	for(int i=0;i<=n;i++){
		for(int j=0;j<=i;j++)
			for(int k=0;k<=i;dp[mp(i+1)][j][k]=0,k++);
		if(s[i+1]=='0')
			for(int j=0;j<=i;j++)
				for(int k=0;k<=i;dp[mp(i+1)][j+1][k]=(dp[mp(i+1)][j+1][k]+dp[mp(i)][j][k])%mod,k++);
		else
			for(int j=0;j<=i;j++)
				for(int k=0;k<=i;k++){
					int x=t[j]-k;
					dp[mp(i+1)][j][k]=(dp[mp(i+1)][j][k]+dp[mp(i)][j][k])%mod;
					dp[mp(i+1)][j][k+1]=(dp[mp(i+1)][j][k+1]-x*dp[mp(i)][j][k])%mod;
					dp[mp(i+1)][j+1][k+1]=(dp[mp(i+1)][j+1][k+1]+x*dp[mp(i)][j][k])%mod;
				}
	}
	for(int j=0;j<=n-m;j++)
		for(int k=0;k<=n;ans=(ans+dp[mp(n)][j][k]*fac[n-k]%mod)%mod,k++);
	cout<<ans;
	return 0;
}