【BZOJ 2669】局部极小值

局部极小值题目链接
蓄水池题目链接
(双倍经验)

题目描述

有一个n行m列的整数矩阵，其中1到nm之间的每个整数恰好出现一次。如果一个格子比所有相邻格子（相邻是指有公共边或公共顶点）都小，我们说这个格子是局部极小值。给出所有局部极小值的位置，你的任务是判断有多少个可能的矩阵。

Sol

诈一看不可做…
回归问题本质
假设啥要求都没有,我们知道答案就是$(n\ast m)!$ ,这是因为一开始的数有$n\ast m$种选择,然后之后的方案数依次递减1

换言之,假设我们一个个放数,我们只要知道了每一次放数的可能的方案数,那么我们就能求出总方案了

发现一个性质:
由于局部极小值的要求是周围的数都比它大,所以假设我们从小往大放数(倒着也类似),那么如果局部极小值处还没有填上数 , 那么它周围的一切格子都不能填上数 , 也就是说我们一次次放数的时候 , 假设我们知道了局部极小值处的填数状况 , 那么我们很容易就能算出这一步的填数方案数

考虑到$4\ast 7$的大小最多只会有$8$个可能的局部极小值 , 一个显然的状压dp就出来了

我们设$dp[i][j]$表示填到了数 $i$ , 局部极小值处填数状态为 $j$ 的方案数
答案就是 $dp[n\ast m][11...11]$

对于每一个状态,我们先把它还能填数的不是局部极小值的位置数算出来(下面记为$ret[j]$),之后就可以转移了
转移分两种情况:

1. 不放在局部极小值处
   这时状态$j$不变,假设有$k$个蓄水池已经填了数,方案数为 $ret[j]-(i-1)+k$
   这表示的就是有多少个没有放数的并且可以放数的非局部极小值格子
2. 放在局部极小值处
   那就枚举放在哪一个上,改变状态直接加上方案数

然而这样会算出不合法的方案 , 因为题目要求不能有其他非局部极小值的地方成为局部极小值
也就是说我们的最后的方案中有的方案局部极小值多了

那么就容斥一下 , 暴搜哪些地方多出来了 , 套上公式算就行了

代码:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#define Set(a,b) memset(a,b,sizeof(a))
using namespace std;
int n,m;
bool c[10][10];
const int N=30;
const int MAXN=1<<10;
const int mod=12345678;
typedef long long ll;
struct chi{
	int x,y;
}a[N];int num=0;
int id[10][10];
int dp[30][MAXN];
int ans=0;
int tot;
inline void upd(int &x,int y){x+=y;if(x>=mod) x-=mod;}
const int dx[9]={0,1,-1,-1,1,0,0,-1,1};
const int dy[9]={0,1,-1,1,-1,-1,1,0,0};
int ret[MAXN];
bool vis[10][10];
inline int cal(int X){
	Set(vis,0);
	for(register int i=1;i<=num;++i){
		register int x=a[i].x,y=a[i].y;
		if(X&1) vis[x][y]=1;
		else{
			for(register int p=1;p<=8;++p){
				register int j=x+dx[p],k=y+dy[p];
				vis[j][k]=1;
			}
		}
		X>>=1;
	}
	register int cnt=0;
	for(register int i=1;i<=n;++i)for(register int j=1;j<=m;++j) if(!vis[i][j]&&!id[i][j]) ++cnt;
	return cnt;
}
#define lowbit(a) ((a)&(-a))
inline int Count(int x){int r=0;while(x) ++r,x-=lowbit(x);return r;}
inline void solve(int flag)
{
	Set(dp,0);
	register int UP=1<<num;
	for(register int i=0;i<UP;++i) ret[i]=cal(i);
	dp[0][0]=1;
	for(register int i=1;i<=tot;++i){
		
		for(register int j=0;j<UP;++j){
			register int ct=Count(j);
			upd(dp[i][j],1ll*dp[i-1][j]*(ret[j]-i+1+ct)%mod);//下到一个非局部极小值位置
			for(register int k=1,z=1;k<=num;++k,z<<=1){//下到一个局部极小值位置
				if(z&j) continue;
				upd(dp[i][j|z],dp[i-1][j]);
			}
		}
		
	}
	if(flag) dp[tot][UP-1]=mod-dp[tot][UP-1];
	upd(ans,dp[tot][UP-1]);
}
void dfs(int x,int y,int tot)
{
	if(y>m) ++x,y=1;
	if(x>n) return solve(tot&1);
	if(id[x][y]) return dfs(x,y+1,tot);
	bool flag=1;dfs(x,y+1,tot);
	for(register int i=1;i<=8;++i){
		register int p=x+dx[i],q=y+dy[i];
		if(id[p][q]) {flag=0;break;}
	}
	if(flag) {
		id[x][y]=++num,a[num]=(chi){x,y},dfs(x,y+1,tot+1);
		id[x][y]=0,--num;
	}
}
int main()
{
	scanf("%d %d",&n,&m);tot=n*m;
	for(register int i=1;i<=n;++i){
		for(register int j=1;j<=m;++j){
			char ch=getchar();
			while(ch!='.'&&ch!='X') ch=getchar();
			if(ch=='X') {
				c[i][j]=1,a[++num]=(chi){i,j};id[i][j]=num;
				if(c[i-1][j]||c[i-1][j-1]||c[i][j-1]) return puts("0"),0;
			}
		}
	}
	if(!num) return puts("0"),0;
	dfs(1,1,0);
	printf("%d\n",ans);
}