标签：oez int CF 奇偶性 read ksm Problem Unnamed oee

Description

计算满足下列条件的 \([m\times n]0/1\) 矩阵的数量：

• 每行后 \(n-k\) 列至少有一个 \(1\)

• 每行互不相同

• \(1\sim p\) 列共有奇数个 \(1\)，\(k+1\sim k+q\) 列共有奇数个 \(1\)

• 行之间无序

\(1\le k<n\le 10^9,m\le 10^6\)

Solution

先不考虑最后一个限制，最后除以 \(m!\) 来无序化

奇偶性的限制可以使用烂大街的套路：最后一行来抵消前面行的奇偶性需求，前面的行任意选

那么需要处理此时不满足第一和第二条限制的问题，将不合法方案数减去

设 \(f_m\) 表示行数为 \(m\) 时的答案，\(g_m\) 表示行数为 \(m\) 且不一定满足限制 \(1\) 时的答案

前者转移可以用 \(\displaystyle\binom{2^n-2^k}{m-1}(m-1)!\) 减掉出现相同的情况数 \((m-1)(2^n-2^k-(n-2))f_{m-2}\) 再减掉不满足最后一行后 \(n-k\) 列至少有 \(1\) 个 \(1\) 的方案数

那么非常精彩地，现在最后一行的后 \(n-k\) 列没有任何 \(1\) 了，前 \(m-1\) 行的后 \(n-k\) 列奇偶性已经满足，所以这是 \(g_{m-1}\)

后者的转移由于不考虑前 \(k\) 列的合法性，所以总方案数要附加最后一行前 \(k\) 列乱选的 \(2^k\) 的系数，出现相同还是 \((m-1)(2^n-2^k-(n-2))g_{m-2}\)

不满足有 \(1\) 的部分也是符合 \(g_{m-1}\) 的定义，但是由于最后一行前 \(k\) 列可以乱选，那么也要附上 \(2^k\) 作为系数

时间复杂度 \(\Theta(m)\)

Code

const int N=1e6+10;
int f[N],g[N];
int n,m,k,oee,oez;
signed main(){
	n=read(); k=read(); m=read(); 
	oee=read(),oez=read();
	int s=ksm(2,k),tot=del(ksm(2,n),ksm(2,k));
	int fac=1,P=1;
	f[0]=!oez&&!oee;
	g[0]=(bool)!oez;
	f[1]=(bool)oez;
	g[1]=f[1]*s;
	for(int i=2;i<=m;++i){
		int cho=del(tot,i-2);
		ckmul(P,cho);
		ckmul(fac,i);
		f[i]=del(P,g[i-1]);
		ckdel(f[i],mul(i-1,mul(f[i-2],cho)));
		g[i]=mul(s,del(P,g[i-1]));
		ckdel(g[i],mul(i-1,mul(g[i-2],cho)));
	}
	print(mul(f[m],ksm(fac,mod-2)));
	return 0;
}

标签：oez,int,CF,奇偶性,read,ksm,Problem,Unnamed,oee
来源： https://www.cnblogs.com/yspm/p/16339502.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。