【bzoj2669】【cqoi2012】【局部极小值】【状压dp】

最新推荐文章于 2020-03-15 13:19:10 发布

sunshinezff

最新推荐文章于 2020-03-15 13:19:10 发布

阅读量1.6k

点赞数

分类专栏： dp

本文链接：https://blog.csdn.net/sunshinezff/article/details/51029181

版权

dp 专栏收录该内容

63 篇文章 3 订阅

订阅专栏

Description

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

给出所有局部极小值的位置，你的任务是判断有多少个可能的矩阵。

Input

输入第一行包含两个整数n和m（1<=n<=4, 1<=m<=7），即行数和列数。以下n行每行m个字符，其中“X”表示局部极小值，“.”表示非局部极小值。

Output

输出仅一行，为可能的矩阵总数除以12345678的余数。

Sample Input

3 2
X.
..
.X

Sample Output

题解：

考虑将1-n*m依次往矩阵里填。

因为局部极小值的数量不会超过8,所以我们可以将每个局部极小值是否填写压成一维状态。

设f[i][j]表示填到第i个数,局部极小值的填写状态为j的方案数。

转移时预处理数组p[j]表示当局部极小值的填写状态为j时有哪些位置可以填数。

枚举当前数是否填在局部极小值的位置上。

f[i][j]=f[i-1][j]*(p[j]-i+1)+f[i-1][k]

但是这样会把一些不是局部极小值的位置当成局部极小值。

所以我们容斥一下即可。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#define P 12345678
using namespace std;
int dx[9]={1,0,-1,0,1,-1,1,-1,0},dy[9]={0,1,0,-1,1,-1,-1,1,0};
int n,m,ans,i,l,f[30][1<<8],p[1<<8],vis[8][8],a[10][2];
char ch[8][8];
int dp(){
  memset(f,0,sizeof(f));
  memset(p,0,sizeof(p));
  int top(0);
  for(int i=1;i<=n;i++)
   for (int j=1;j<=m;j++)
    if (ch[i][j]=='X')
      a[++top][0]=i,a[top][1]=j; 
  for (int k=0;k<(1<<top);k++){
    memset(vis,0,sizeof(vis));
    for (int i=1;i<=top;i++)
     if (~k&(1<<i-1)) vis[a[i][0]][a[i][1]]=1; 
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++){
      for (l=0;l<9;l++)
       if (vis[i+dx[l]][j+dy[l]]) break; 
      if (l==9) p[k]++;
     }
  }
  f[0][0]=1;
  for (int i=1;i<=n*m;i++)
   for (int j=0;j<(1<<top);j++){
    (f[i][j]+=(long long)f[i-1][j]*max(p[j]-i+1,0))%=P;
    for (int k=1;k<=top;k++)
      if(j&(1<<k-1)) (f[i][j]+=f[i-1][j^(1<<k-1)])%=P; 
   }
  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){(ans+=dp()*(t&1?-1:1))%=P;return;}
  dfs(x,y+1,t);
  for (i=0;i<9;i++)
   if (ch[x+dx[i]][y+dy[i]]=='X') break;
  if (i==9){ch[x][y]='X';dfs(x,y+1,t+1);ch[x][y]='.';} 
}
int main(){
 scanf("%d%d",&n,&m);
 for (int i=1;i<=n;i++) scanf("%s",ch[i]+1); 
 for (int i=1;i<=n;i++)
  for (int j=1;j<=m;j++)
   if (ch[i][j]=='X')
    for (int k=0;k<8;k++)
     if (ch[i+dx[k]][j+dy[k]]=='X')
      return printf("0"),0;
 dfs(1,1,0);
 cout<<(ans+P)%P<<endl;
}