P3160 [CQOI2012]局部极小值（dfs+状压dp）

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
char a[10];
int m,n,min[6][10];//描述整个矩阵
int num,x[30],y[30];//描述坑的个数、位置
int w=12345678;//膜
int dx[10]={-1,-1,-1,0,0,1,1,1,0};//移动位置
int dy[10]={-1,0,1,-1,1,-1,0,1,0};
int vis[6][10],f[29][(1<<8)+10],hi[1<<9];//dp用到的东西
int dp(){
    int i,j,k;
    memset(f,0,sizeof(f));
    f[0][0]=1;
    for(i=0;i<(1<<num);i++){//预处理出每个状态i对应的可填点数量
        hi[i]=n*m;
        memset(vis,0,sizeof(vis));
        for(j=0;j<num;j++)if(!(i&(1<<j)))for(k=0;k<9;k++)vis[x[j]+dx[k]][y[j]+dy[k]]=1;
        for(j=1;j<=n;j++)for(k=1;k<=m;k++)if(vis[j][k])hi[i]--;
    }
    for(i=1;i<=n*m;i++){//枚举填哪个数
        for(j=0;j<(1<<num);j++){//枚举状态
            if(hi[j]-i+1>0)f[i][j]+=f[i-1][j]*(hi[j]-i+1),f[i][j]%=w;
            for(k=0;k<num;k++)if((1<<k)&j)f[i][j]+=f[i-1][j^(1<<k)],f[i][j]%=w;
        }
    }
    return f[n*m][(1<<num)-1];
}
int cnt=0;
int P=0;
int dfs(int X,int Y){
    if(Y==m+1)X++,Y=1;
    if(X==n+1)
    {
        if(cnt%2==1)P-=dp();else P+=dp();
        P%=w;P+=w;P%=w;
        return 0;
    }
     dfs(X,Y+1);
 int i;
    for(i=0;i<9;i++)if(min[X+dx[i]][Y+dy[i]])return 1;

    //如果没return说明这个地方是个可能变成坑的地方，那就把它变成坑dfs一下
    x[num]=X;y[num++]=Y;
    min[X][Y]=1;
    cnt++;
    dfs(X,Y+1);
    min[X][Y]=0;num--;//别忘了变回去
    cnt--;
    return  1;
}
int main(){
    int i,j;
    scanf("%d%d",&n,&m);
    for(i=1;i<=n;i++){
        scanf("%s",a+1);
        for(j=1;j<=m;j++)if(a[j]=='X'){
            min[i][j]=1;
            x[num]=i;y[num++]=j;
        }
    }
    for(i=0;i<num;i++)for(j=0;j<i;j++)if(abs(x[i]-x[j])<2&&abs(y[i]-y[j])<2)return printf("0"),0;
    if(!num)return printf("0"),0;
    dfs(1,1);
    printf("%d",P);
    return 0;
}