uoj#422. 【集训队作业2018】小Z的礼物(MIn-Max容斥+插头dp)

时间:2021-03-05 21:52:55

题面

传送门

题解

好迷……

很明显它让我们求的是\(Max(S)\),我们用\(Min-Max\)容斥,因为\(Min(S)\)是很好求的,只要用方案数除以总方案数算出概率,再求出倒数就是期望了

然而如果爆搜枚举子集的话复杂度是\(O(2^{cnt})\)的

发现总共的方案数只有\(2*n*m-n-m\)种,而且\(n\)非常小,我们可以考虑插头\(dp\)

设\(f_{i,S,k}\)表示做到了第\(i\)列,插头的状态为\(S\),覆盖方案数为\(k\)时的方案总数,并且这个里面已经考虑了容斥系数

然后直接转移就是了

//minamoto
#include<bits/stdc++.h>
#define R register
#define fp(i,a,b) for(R int i=a,I=b+1;i<I;++i)
#define fd(i,a,b) for(R int i=a,I=b-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
using namespace std;
const int N=1505,P=998244353;
inline int add(R int x,R int y){return x+y>=P?x+y-P:x+y;}
inline int dec(R int x,R int y){return x-y<0?x-y+P:x-y;}
inline int mul(R int x,R int y){return 1ll*x*y-1ll*x*y/P*P;}
char s[N][N];int f[2][(1<<6)+5][N],inv[N],g[N];
int n,m,t,G,sum,tmp,res;
int main(){
// freopen("testdata.in","r",stdin);
scanf("%d%d",&n,&m),G=(1<<n)-1;
fp(i,0,n-1)scanf("%s",s[i]);
sum=n*m*2-n-m;
inv[0]=inv[1]=1;fp(i,2,sum)inv[i]=mul(P-P/i,inv[P%i]);
f[0][0][0]=P-1,t=0;
fp(j,0,m-1)fp(i,0,n-1){
memset(f[t^1],0,sizeof(f[t^1]));
fp(S,0,G)fp(k,0,sum)if(f[t][S][k]){
int T=S&(G^(1<<i));
f[t^1][T][k]=add(f[t^1][T][k],f[t][S][k]);
if(s[i][j]=='*'){
T|=1<<i;
tmp=(i&&!(S&(1<<(i-1))))+(j&&!(S&(1<<i)))+(i<n-1)+(j<m-1);
f[t^1][T][k+tmp]=add(f[t^1][T][k+tmp],P-f[t][S][k]);
}
}
t^=1;
}
fp(S,0,G)fp(i,1,sum)res=add(res,mul(f[t][S][i],inv[i]));
res=mul(res,sum);printf("%d\n",res);
return 0;
}