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

思路

我们考虑设fS表示S里面的这些题目都有区分度的方案数。
那么我们可以知道fS =
Σ︀
S′⊂S (−1)|S′|gS′ ,其中S′是S的一个
子集,而gS的含义是有多少种方案使得S中的所有题都没有区
分度.
最后我们考虑如何计算gS。我们设计一个hS ,S是一个长度
为m的三进制数,如果第i位为2表示这一位无所谓，然后如果
这一位为0/1则表示这一位选择0/1的数的个数.
设S′是S中不为2的位置集合,则gS′+ = 2hS − 1,最后再给每
个gS加上空集就行。
如何计算hS?如果S中没有2,则可以直接算出来，否则我们随
便找一个2,将这个2分别变成0/1得
到S′, S′′,hS = hS′ + hS′′ ,按照0 + 1的个数从小往大枚举就
行，两步的复杂度都是O(3M),总复杂度为O(3M + NM).
模拟

代码

#include<bits/stdc++.h>
#define N 300005
#define Mo 998244353
using namespace std;
 
int f[15000005],n,pw3[N],pw2[N],g[1<<15],m,size[N],Cas;
 
char c[N];
 
void dfs(int u,int _2,int _3,int _4)
{
	if (u==m)
	{
		if (_4!=-1)
			f[_3]=f[_3-pw3[m-1-_4]]+f[_3-2*pw3[m-1-_4]];
		(g[_2]+=pw2[f[_3]]-1)%=Mo;
	}
	else
	{
		dfs(u+1,(_2<<1)+1,_3*3,_4);
		dfs(u+1,(_2<<1)+1,_3*3+1,_4);
		dfs(u+1,(_2<<1),_3*3+2,u);
	}
}
 
int main()
{
	pw2[0]=1; 
	for(int i=1; i<N;++i) pw2[i]=pw2[i-1]*2%Mo;
	for(int i=1; i<(1<<15);++i) size[i]=size[i/2]+(i&1); 
	pw3[0]=1; 
	for(int i=1; i<20;++i) pw3[i]=pw3[i-1]*3;
	scanf("%d%d",&n,&m);
	for(int i=1; i<=n;++i)
	{
		scanf("%s",c); int tmp=0;
		for(int j=0;j<m;++j)  tmp=tmp*3+(c[j]=='1');
		f[tmp]++;
	}
	dfs(0,0,0,-1); int yjy=0;
	for(int i=0; i<(1<<m);++i) (g[i]+=1)%=Mo;
	for(int i=1; i<(1<<m);++i)
	{
		int tmp=0;
		for(int j=i;;j=(j-1)&i)
		{
			if (size[j]&1) tmp=(tmp-g[j]+Mo)%Mo;
			else tmp=(tmp+g[j])%Mo;
			if (!j) break;
		}
		yjy^=tmp;
	}
	printf("%d\n",yjy);
}