1. 程式人生 > >bozj2669(容斥+狀壓dp)

bozj2669(容斥+狀壓dp)

一張圖最多8個區域性最小值。

dp[i][j]表示正在填從小到大第i個數,區域性最小值所在的位置已被填的情況為j時的方案數

p[i] 為區域性最小值所在的位置已被填的情況為j時,所有可以填數的位置數。(x為沒被填的區域性最小值,x的周圍都不能填,如果填了x就不是區域性最小值了!所以才會存在可以填數的位置數)

當i 不為區域性最小值時 dp[i][j] += dp[i-1][j]*(p[i] - i + 1); 

當i 為區域性最小值時,列舉i填的是狀態j中的哪一個區域性最小值 dp[i][j] += dp[i-1][j^(1<<(k-1))];

在第二種情況時,是區域性最小值的點一定是區域性最小值,不是區域性最小值的點可能會變成區域性最小值,需要剪掉多算的部分。

在dfs時容斥一下=、=

http://www.cto800.com/view/50553177130254010924.html(這個寫的特別好)

#include<iostream>
#include<cstdio>
#include<cstring> 
using namespace std;
const int mod = 12345678;
int n, m, zl[15][2] = {{1,1},{1,0},{-1,0},{-1,-1},{0,-1},{0,1},{-1,1},{1,-1},{0,0}};
char chu[5][10];
int ans, p[260], f[50][260], zhan[10][8], vis[10][10];
int dp()
{
	memset(p,0,sizeof(p));
	memset(f,0,sizeof(f));
   int top = 0;
   for(int i = 1; i <= n; i++)
    for(int j = 1; j <= m; j++)
    {
    	if(chu[i][j] == 'X')
    	{
    		top++;
    		zhan[top][0] = i;
    		zhan[top][1] = j;
    	}
    }
   for(int k = 0; k < (1 << top); k++)
   {
   	memset(vis,0,sizeof(vis));
   	for(int j = 1; j <= top; j++)
   		if((k & (1 << (j - 1))) == 0) vis[zhan[j][0]][zhan[j][1]] = 1;
   	for(int i = 1; i <= n; i++)
   	 for(int j = 1; j <= m; j++)
   	 	for(int z = 0; z <= 9; z++)
   	 	{
   	 		if(z == 9){ p[k] ++; break;} 
   	 	   if(vis[i + zl[z][0]][j + zl[z][1]] == 1) break;
		}
   }
   f[0][0] = 1;
   for(int i = 1; i <= n * m; i++)
      for(int j = 0; j < (1 << top); j++)
      {
        f[i][j] = 1ll * f[i - 1][j] * max(0, p[j] - i + 1) % mod;
		for(int k = 1; k <= top; k++)
		{
		  if(j & 1 << (k-1)) f[i][j] = (f[i][j] + f[i-1][j^(1<<(k-1))]) % mod;
		} 
      }
   return f[n*m][(1<<top) - 1];
}
void dfs(int x, int y, int t)
{
	if(y == m+1) {
		dfs(x+1,1,t);	return ;
	}
	if(x == n+1) 
	{
		int ha = -1;
		if(t % 2 == 0) ha = 1;
		ans = (ans + 1ll * ha * dp() + mod) % mod;
		return ;
	}
   dfs(x, y+1, t);
   for(int i = 0; i <= 9; i++)
   {
   	if(i == 9)
   	{
   	chu[x][y] = 'X';
   	dfs(x,y+1,t+1);
   	chu[x][y] = '.';
   	break;
   	}
	if(chu[x + zl[i][0]][y + zl[i][1]] == 'X')
   	 break;
   }
}
int main()
{
	cin >> n >> m;
	for(int i = 1; i <= n; i++)
	 scanf("%s",chu[i]+1);
	for(int i = 1; i <= n; i++)
	 for(int j = 1; j <= m; j++)
	 {
	 	if(chu[i][j] == 'X') 
	 	for(int k = 0; k < 8; k++)
	     if(chu[i+zl[k][0]][j+zl[k][1]] == 'X')
	     {
	        printf("0");
			return 0; 
	     }
	 }
	dfs(1,1,0);
	cout << (ans + mod) %mod;
	return 0;
}