【UOJ#422】【集訓隊作業2018】小Z的禮物(min-max容斥,輪廓線dp)


【UOJ#422】【集訓隊作業2018】小Z的禮物(min-max容斥,輪廓線dp)

題面

UOJ

題解

毒瘤xzy,怎么能搬這種題當做WC模擬題QwQ
一開始開錯題了,根本就不會做。
后來發現是每次任意覆蓋相鄰的兩個,那么很明顯就可以套\(min-max\)容斥。
要求的就是\(max(All)\),而每個集合的\(min\)是很好求的。
如果直接暴力枚舉集合復雜度就是\(2^{cnt}cnt\)
仔細想想每個子集我們要知道的是什么,只需要知道子集大小來確定前面的容斥系數,還需要知道覆蓋子集的方案數,這樣就可以知道\(min\)的概率,倒數就是期望了。
而覆蓋子集的方案數不會超過\(2*n*m-n-m\),顯然要比\(2^{cnt}\)優秀。
所以枚舉覆蓋子集的方案數來\(dp\),至於子集大小之和容斥系數相關,而容斥系數只有正負\(1\),所以直接乘進去一起轉移就好了,不需要單獨存一維狀態。
考慮每次新加入一個點之后的覆蓋方案,只需要知道當前位置四周是否已經存在於子集當中,那么直接狀壓當前的輪廓線就好了。
設狀態\(f[S][k]\)表示覆蓋方案數為\(k\),輪廓線為\(S\)時的方案數,容斥系數已經考慮進去。
顯然當前位置可以不選,那么直接轉移。
如果當前位置可以選入子集,那么乘上系數\(-1\),同時修改覆蓋方案數以及輪廓線的狀態轉移。
最后按照\(min-max\)容斥統計答案即可。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
#define MOD 998244353
void add(int &x,int y){x+=y;if(x>=MOD)x-=MOD;}
char g[200][200];
int n,m,inv[1500],sum,pw,nw,ans,S;
int f[2][1<<6][1200];
int main()
{
	scanf("%d%d",&n,&m);S=1<<n;sum=2*n*m-n-m;
	for(int i=1;i<=n;++i)scanf("%s",g[i]+1);
	inv[0]=inv[1]=1;for(int i=2;i<1500;++i)inv[i]=1ll*inv[MOD%i]*(MOD-MOD/i)%MOD;
	f[0][0][0]=MOD-1;
	for(int i=1;i<=m;++i)
		for(int j=1;j<=n;++j)
		{
			pw=nw;nw^=1;memset(f[nw],0,sizeof(f[nw]));
			for(int T=0;T<S;++T)
				for(int k=0;k<=sum;++k)
					if(f[pw][T][k])
					{
						int nT=T&((S-1)^(1<<(j-1)));
						add(f[nw][nT][k],f[pw][T][k]);
						if(g[j][i]=='*')
						{
							nT|=1<<(j-1);int pls=0;
							if(j>1&&!(T&(1<<(j-2))))++pls;
							if(i>1&&!(T&(1<<(j-1))))++pls;
							if(i<m)++pls;if(j<n)++pls;
							add(f[nw][nT][k+pls],MOD-f[pw][T][k]);
						}
					}
		}
	for(int T=0;T<S;++T)
		for(int i=1;i<=sum;++i)
			add(ans,1ll*f[nw][T][i]*inv[i]%MOD);
	ans=1ll*ans*sum%MOD;printf("%d\n",ans);
	return 0;
}


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM