bozj2669(容斥+狀壓dp)
阿新 • • 發佈:2018-11-08
一張圖最多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; }