ZYB的染色计划

题目

题目描述
ZYBZYB是一名资深特级OIOI教师。在他的题库里，一共有MM道判断题。现在，他正在教授NN名OIOI选手，作为一名特级教师，他一眼就看出了每个学生对每道问题的答案是YESYES还是NONO.

现在他希望出一场堂测，他会选择这些判断题当中的一个非空子集(2M−12M−1种可能)。他希望这场堂测能有一些区分度。定义一道题目有区分度，当且仅当至少有一个学生回答YESYES,并且至少有一个学生回答NONO，而一场堂测有区分度，当且仅当所有的题目都是有区分度的。

然而，这些学生有可能会咕咕咕。每个学生只有50%50%的可能出现在这堂课上进行堂测。现在ZYBZYB想知道，对于每个不同的子集，有多大的概率这场堂测是有区分度的呢?

为了方便起见，假设这2M−12M−1个概率分别是Ans1,Ans2,…,Ans2M−1Ans1,Ans2,…,Ans2M−1,你只需要输出(Ans1×2Nmod998244353)xor(Ans2×2Nmod998244353)xor…xor(Ans2M−1×2Nmod998244353)(Ans1×2Nmod998244353)xor(Ans2×2Nmod998244353)xor…xor(Ans2M−1×2Nmod998244353)的值.

输入格式
第一行两个整数n,mn,m.

接下来nn行,每行一个长度为mm的0101字符串,其中第ii位为11表示该学生对这个问题回答YESYES,00表示回答NONO.

输出格式
一行输出所求值。

样例1
输入

2 2
01
10
输出

1
解释 对于每道题目，当且仅当所有人都到时才有区分度.故Ans1=Ans2=Ans3=1/4Ans1=Ans2=Ans3=1/4,带入可知是11.

样例2
输入

4 2
00
01
10
11
输出

7
解释 对于第一题,当只有1,21,2或只有3,43,4到时是没区分度的,所以可能性是(16-7)/16,第二题也是如此.

若需要两题都有区分度,则需要1,41,4同时到或2,32,3同时到,则(1,4),(2,3)以及任意超过33个人到都是合法解。可能性为7/167/16.

于是，答案为9 xor 9 xor 7=7

数据范围
对于30%30%的数据,M≤4M≤4.

对于60%60%的数据,M≤10M≤10.

对于100%100%的数据,N≤100000,M≤15

思路

我们考虑如何对一组染色方案判断是否可行。可以记录一个
插头dpi ,j ,S表示到(i , j)为止,上面n + 1个插头的状态为S是否
可行。只需要判断dpn,m,0是否可行就行.
这个形式让我们想起了dp of dp,我们只需要把这个长度
为2n+1的dp数组都放到状态里就行了,让fi ,j ,c,S来表示
到(i , j)这个格子,前面n个格子的颜色
是c(0 − 26 − 1),128个dp状态为S的方案数。
我们对这个dp做一个小优化，将状态变为26个.也就是说我
们只记录每个位置当前存不存在插头。对于当前第j − 1个位
置的插头，我们暂时不分它是右插头还是下插头，而当处
理j时,如果j − 1有插头,可以让它成为右插头，如果不用就当
下插头用了。这样优化之后实际这个S就比较少了，不会超
过2000个。

代码

#include <bits/stdc++.h>
#define int long long
#define For(it,c) for(__typeof(c.begin())it=c.begin();it!=c.end();it++)
using namespace std;
const int Mod=998244353;
typedef unsigned long long ULL;
map<pair<int,ULL>,int> f[2];
int n,m;
signed main()
{
	cin>>m>>n;
	f[0][make_pair(0,1)]=1;
	bool yjy=0;
	for(int i=1;i<=n;i++)
	{
		for(int j=0;j<m;j++)
		{
			yjy=!yjy;
			f[yjy].clear();
			For(it,f[!yjy])
			{
				int sum=it->second,col=it->first.first;
				ULL aii=it->first.second;
				for(int k=0;k<2;k++)
				{
					ULL naii=0;
					for(int bit=0;bit<(1<<m);bit++)
					{
						if(aii>>bit&1)
						{
							int upcol=col>>j&1,upbit=bit>>j&1;
							int lecol,lebit;
							if(j)lecol=col>>(j-1)&1,lebit=bit>>(j-1)&1;
							if(upbit)
							{
								if(k!=upcol)naii|=1ULL<<(bit xor (1<<j));
							}
							else
							{
								if(j && k!=lecol && lebit)naii|=1ULL<<(bit xor (1<<(j-1)));
								naii|=1ULL<<(bit xor (1<<j));
							}
						}
					}
					if(naii)
					{
						f[yjy][make_pair(col xor ((col>>j&1)<<j) xor (k<<j),naii)]+=sum;
						f[yjy][make_pair(col xor ((col>>j&1)<<j) xor (k<<j),naii)]%=Mod;
					}
				}
			}
		}
	}
	int ans=0;
	For(it,f[yjy])
	{
		if(it->first.second&1)
		{
			ans+=it->second;
			ans%=Mod;
		}
	}
	cout<<ans<<endl;
	return 0;
}