题链：

http://www.lydsy.com/JudgeOnline/problem.php?id=2669

题解：

容斥，DP，DFS

先看看 dp 部分：
首先呢，X的个数不会超过 8个。
个数很少，所以考虑状压，把需要填 X的那几个位置状压为二进制10表示对应的那个X位置是否已经填数。
同时填的数互不重复，考虑从小填到大。
令 cnt[S] 表示除了不在集合 S 里的 X 位置及其周围的位置，剩下的位置个数。
定义 dp[i][S]表示从小到大填数填完了i这个数，且已经填了的 S 这个集合里的 X 位置的方案数。
转移：依次去填数 1~N*M，每次有两种选择：
1).把这个数填在 某个 X 位置(枚举一个 k表示第 k个 X 位置填当前数)
dp[i][S]+=dp[i-1][S^(1<<(k-1))]
2).把这个数填在非 X 位置，那么填的位置有 cnt[S]-(i-1) 种。
dp[i][s]+=dp[i-1][s]*(cnt[s]-(i-1)) (好好理解一下这个转移)
这样 dp 可以保证那些给出的 X 位置一定是局部最小值，
因为第二种转移的填数位置都不能填在还没有填数的 X 位置的周围。
所以就完了么？
当然还没有，尽管我们保证了给出的 X 位置一定是局部最小值，
但是没有保证非 X位置一定不是非局部最小值。即,求出来的 dp[N*M][all_S(全集)]的意思是至少all_S集合里的 X位置为局部最小值的方案数。
所以容斥如下：
ANS = 至少多填了0个局部最小值的方案数(dp[N*M][all_S])
          -至少多填了1个局部最小值的方案数
         +至少多填了2个局部最小值的方案数
          -....+ ....
这些用于容斥的方案数的求法：
DFS 搜索出哪些非 X 位置还可以改为 X ,
然后对于每一种新的填法，去跑一遍上述的dp即可求得对应的方案数。

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#define _ % mod
#define filein(x) freopen(#x".in","r",stdin);
#define fileout(x) freopen(#x".out","w",stdout);
using namespace std;
const int mv[9][2]={{0,0},{-1,0},{-1,1},{0,1},{1,1},{1,0},{1,-1},{0,-1},{-1,-1}};
const int mod=12345678;
char mp[10][10];
int N,M,ANS;
int solve(){static bool vis[10][10];static int dp[30][1<<8],cnt[1<<8],x[10],y[10],tot,tmp;tot=0; memset(dp,0,sizeof(dp));for(int i=1;i<=N;i++)for(int j=1;j<=M;j++) if(mp[i][j]=='X') tot++,x[tot]=i,y[tot]=j;for(int s=0;s<1<<tot;s++){tmp=0; memset(vis,0,sizeof(vis));for(int i=1;i<=tot;i++) if(!(s&(1<<(i-1))))for(int k=0;k<9;k++)vis[x[i]+mv[k][0]][y[i]+mv[k][1]]=1;for(int i=1;i<=N;i++)for(int j=1;j<=M;j++) if(!vis[i][j]) tmp++;cnt[s]=tmp;}dp[0][0]=1;for(int i=1;i<=N*M;i++)for(int s=0;s<1<<tot;s++){dp[i][s]=(1ll*dp[i][s]+1ll*dp[i-1][s]*max(cnt[s]-(i-1),0)_)_;for(int k=1;k<=tot;k++) if(s&(1<<(k-1)))dp[i][s]=(1ll*dp[i][s]+dp[i-1][s^(1<<(k-1))])_;}return dp[N*M][(1<<tot)-1];
}
void dfs(int x,int y,int t){if(y==M+1){dfs(x+1,1,t);return;}if(x==N+1){int tmp=solve();if(t&1) tmp=(-1ll*tmp+mod)_;ANS=((1ll*ANS+tmp)_+mod)_;return;}dfs(x,y+1,t);bool fg=1;for(int k=0;k<9;k++) if(mp[x+mv[k][0]][y+mv[k][1]]=='X') fg=0;if(fg){mp[x][y]='X';dfs(x,y+1,t+1);mp[x][y]='.';}
}
int main()
{scanf("%d%d",&N,&M);for(int i=1;i<=N;i++) scanf("%s",mp[i]+1);dfs(1,1,0);printf("%d",ANS);return 0;
}