bzoj 2669: [cqoi2012]局部极小值

题意
有一个n行m列的整数矩阵，其中1到nm之间的每个整数恰好出现一次。如果一个格子比所有相邻格子（相邻是指有公共边或公共顶点）都小，我们说这个格子是局部极小值。
给出所有局部极小值的位置，你的任务是判断有多少个可能的矩阵。
###题解：
挺好的一道题。
状压dp+容斥。
既然极小，考虑从小到大填。
容易发现，极小值的位置一定不超过8，所以可以状压各个极小值位置是否填写。
$f[i][t]$表示填入了前i个数，极小值填写状态是t，的方案数。
先与处理当填写状态为t时，有多少个位置可以填入当前数（包括之前已经填入的极小值位置，方便转移，具体见代码），那么就可以转移了。
枚举填入那个极小值位置，或者其他位置。

for(int i=1;i<=n*m;i++)
		for(int t=0;t<(1<<tp);t++)
		{
			(f[i][t]+=f[i-1][t]*max(p[t]-i+1,0LL)%mod)%=mod;
			for(int j=1;j<=tp;j++)
				if(t&(1<<(j-1))) (f[i][t]+=f[i-1][t^(1<<(j-1))])%=mod;
		}

但是还有个问题，就是可能不是极小值的地方变成了极小值，所以容斥一发。
code：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
const LL mod=12345678;
int yh[10]={-1,-1,-1,0,1,1,1,0,0};
int yl[10]={-1,0,1,1,1,0,-1,-1,0};
LL f[35][1<<8],ans=0,p[1<<8];
int n,m;
char s[6][9];
int X[10],Y[10];
int dp()
{
	memset(f,0,sizeof(f));memset(p,0,sizeof(p));
	int tp=0;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			if(s[i][j]=='X') X[++tp]=i,Y[tp]=j;
	for(int t=0;t<(1<<tp);t++)
	{
		bool v[6][9];
		memset(v,false,sizeof(v));
		for(int i=1;i<=tp;i++)
			if(~t&(1<<(i-1))) v[X[i]][Y[i]]=true;
		for(int i=1;i<=n;i++)
			for(int j=1;j<=m;j++)
			{
				int k;
				for(k=0;k<9;k++) if(v[yh[k]+i][yl[k]+j]) break;
				if(k>8) p[t]++;
			}
	}
	//for(int t=0;t<(1<<tp);t++) printf("%d ",p[t]);printf("\n");
	f[0][0]=1;
	for(int i=1;i<=n*m;i++)
		for(int t=0;t<(1<<tp);t++)
		{
			(f[i][t]+=f[i-1][t]*max(p[t]-i+1,0LL)%mod)%=mod;
			for(int j=1;j<=tp;j++)
				if(t&(1<<(j-1))) (f[i][t]+=f[i-1][t^(1<<(j-1))])%=mod;
		}
	//printf("%d %lld\n",tp,f[n*m][(1<<tp)-1]);
	return f[n*m][(1<<tp)-1];
}
void dfs(int x,int y,int cnt)
{
	if(y==m+1) {dfs(x+1,1,cnt);return;}
	if(x==n+1)
	{
		(ans+=dp()*(cnt&1?-1:1))%=mod;
		return;
	}
	dfs(x,y+1,cnt);
	int k;
	for(k=0;k<9;k++) if(s[yh[k]+x][yl[k]+y]=='X') break;
	if(k>8) {s[x][y]='X';dfs(x,y+1,cnt+1);s[x][y]='.';}
}
int main()
{
	scanf("%d %d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%s",s[i]+1);
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++)
		if(s[i][j]=='X')
			for(int k=0;k<8;k++)
				if(s[i+yh[k]][j+yl[k]]=='X') return puts("0"),0;
	dfs(1,1,0);
	printf("%lld",(ans+mod)%mod);
}