P3160 [CQOI2012]局部极小值

题目

P3160 [CQOI2012]局部极小值

一眼就是状压，接下来就不知道了\(qwq\)

做法

我们能手玩出局部小值最多差不多是\(8,9\)个的样子，\(dp_{i,j}\)为填满\(1~i\)数字，局部小值的状态为\(j\)

第\(k\)个局部极小值填\(i\)：\(dp[i][j]=(dp[i][j]+dp[i-1][j^(1<<k-1)])%p\)

不填在局部极小值，显然有些地方不能填\(i\)的，首先还没填的局部极小值不填，其周围也不能填(填\(i\)后后面再填比不符合局部极小值)

我们预处理出每种状态不能填时的位置：\(dp[i][j]=(dp[i][j]+dp[i-1][j]*max(num[j]-i+1,0))%p\)

一顿操作后发现WA了，我们得到的局部极小值可能不是恰好是给出的\(X\)(更多)

容斥原理：Ans=至少多0个极小值-至少多1个极小值+至少多两个极小值......

理解：至少多0个\((x,x+1,x+2,x+3...)-k(x+1,x+2,x+3...)\)其中\(k\)为某些值的系数

填了至少多0个极小值后，有多填的部分那就得减去至少多1个极小值，至少多1个极小值的位置也有很多，就会有一些位置减掉的系数为(k>1)，就要又加上

\(dfs\)时加上能被多余的极小值填上的地方换成'X'，然后多次做\(dp\)

My complete code

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long LL;
const int maxn=31;
const int p=12345678;
const int dx[8]={-1,-1,-1,0,1,1,1,0},dy[8]={-1,0,1,1,1,0,-1,-1};
int n,m,tot;
LL ans;
int x[maxn],y[maxn],num[1<<9],dp[maxn][1<<9];
bool visit[maxn][maxn];
char s[maxn][maxn];
inline int Solve(){
	tot=0;
	for(int i=1;i<=n;++i)
	    for(int j=1;j<=m;++j)
	        if (s[i][j]=='X')
	            x[++tot]=i,
				y[tot]=j;
	int Up=1<<tot;
	for(int i=0;i<Up;i++){
	 	int cnt(0);
	 	memset(visit,0,sizeof(visit));
		for(int j=1;j<=tot;++j)
	 	    if(!((i>>(j-1))&1)){
	 	  	    visit[x[j]][y[j]]=1;
	 	  	    for(int k=0;k<8;++k){
	 	  	 	    int xx=x[j]+dx[k];
	 	  	 	    int yy=y[j]+dy[k];
	 	  	 	    if(xx>0&&yy>0&&xx<=n&&yy<=m)
	 	  	 	        visit[xx][yy]=1;
			    }
		    }
		for(int j=1;j<=n;++j)
		    for(int k=1;k<=m;++k)
		        if(visit[j][k])
				    ++cnt;
		num[i]=n*m-cnt;
	}
	memset(dp,0,sizeof(dp));
	dp[0][0]=1;
	for(int i=1;i<=n*m;++i)
	    for(int j=0;j<Up;++j){
	  	    dp[i][j]=(dp[i][j]+dp[i-1][j]*max(num[j]-i+1,0))%p;
	  	    for(int k=1;k<=tot;k++)
	  	        if(j&(1<<(k-1)))
	  	            dp[i][j]=(dp[i][j]+dp[i-1][j^(1<<k-1)])%p;
	    }
	return dp[n*m][Up-1];
}
void Dfs(int x,int y,int k){
	if(y==m+1){
	    Dfs(x+1,1,k);
	 	return;
	}
	if(x==n+1){
		ans=(ans+(((k&1)==0)?1:-1)*Solve()+p)%p;
	 	return;
	}
	Dfs(x,y+1,k);
	if(s[x][y]!='X'){
	    bool f=true;
	    for(int i=0;i<8;i++)
	        if(s[x+dx[i]][y+dy[i]]=='X'){
	  	        f=false;
	  	        break;
	        }
	    if(f){
	        s[x][y]='X',
	        Dfs(x,y+1,k+1),
	        s[x][y]='.';
	    }
	}
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
    	scanf(" %s",s[i]+1);
	Dfs(1,1,0);
	printf("%lld\n",ans);
	return 0;
}
/*
5 5
.....
...X.
.X...
.....
....X
*/