problem

luogu-P3160

solution

这么小的数据范围，非暴力不状压。暴力 O(28!)O(28!)O(28!) 呵呵呵可以拉走了。

我们不妨从小到大填数字，这样如果局部极小值点还没有填的话，周围的九宫格就一定不能被填。

设 dp(s,i):dp(s,i):dp(s,i): 局部极小值点的是否填了数字的情况 sss，已经填完了 [1,i][1,i][1,i] 以内的数字。

将 i+1i+1i+1 填入非极小值点影响区域。

dp(s,i+1)←dp(s,i)∗(cnt−i)dp(s,i+1)\leftarrow dp(s,i)*(cnt-i)dp(s,i+1)←dp(s,i)∗(cnt−i)。

cnt:cnt:cnt: 除去现在还没有被填的极小值以及其九宫格内的格子，剩下的个子数量，显然要空出 n−cntn-cntn−cnt 个格子去满足极小值点周围的要求。

不能简单的用 9×9\times9× 极小值点个数，显然两个极小值的九宫格有可能重叠部分。
将 i+1i+1i+1 填入极小值点。

那么该极小值一旦被填，就可以释放周围九宫格的限制。

dp(s∣(1<<j),i+1)←dp(s,i)dp(s|(1<<j),i+1)\leftarrow dp(s,i)dp(s∣(1<<j),i+1)←dp(s,i)。

具体可见下方代码。

但是这样 dpdpdp 我们发现一个致命的问题。

我们会计算到『某些九宫格内中心点不是 X\text{X}X，但中心点又恰好满足局部极小值的要求』的情况。

而包含这种情况的方案在题目中是不合法的。题目已经给出了所有局部极小值点，其余点就必须不是。

这里我们就想到了容斥。不好意思我没有想到

我们可以额外多钦定了 kkk 个极小值点，那么算出来的 dpdpdp 应该前面配一个 (−1)k(-1)^k(−1)k 容斥出最后的答案。

而可以同时存在的局部极小值点个数并没有多少，顶多 888 个。

所以 kkk 并不会多大。

我们完全可以搜索做。

code

#include <bits/stdc++.h>
using namespace std;
#define int long long
#define mod 12345678
int n, m, ans;
int px[10], py[10];
char ch[10][10];
int vis[10][10];
int f[1 << 10][30];bool inside( int x, int y ) {if( x < 0 or x >= n or y < 0 or y >= m ) return 0;else return 1;
}int calc() {memset( f, 0, sizeof( f ) );int tot = 0;for( int i = 0;i < n;i ++ )for( int j = 0;j < m;j ++ )if( ch[i][j] == 'X')px[tot] = i, py[tot] = j, ++ tot;f[0][0] = 1;for( int s = 0;s < (1 << tot);s ++ ) {memset( vis, 0, sizeof( vis ) );for( int i = 0;i < tot;i ++ )if( ! (s >> i & 1) )for( int x = -1;x <= 1;x ++ )for( int y = -1;y <= 1;y ++ )if( inside( px[i] + x, py[i] +  y) )vis[px[i] + x][py[i] + y] = 1;int cnt = n * m;for( int i = 0;i < n;i ++ )for( int j = 0;j < m;j ++ )cnt -= vis[i][j];for( int i = 0;i <= cnt;i ++ )if( f[s][i] ) {( f[s][i + 1] += f[s][i] * ( cnt - i ) ) %= mod;for( int j = 0;j < tot;j ++ )if( ! (s >> j & 1) )( f[s | (1 << j)][i + 1] += f[s][i] ) %= mod;}}return f[(1 << tot) - 1][n * m];
}void dfs( int x, int y, int k ) {if( x >= n ) return ( ans += k * calc() ) %= mod, void();if( y >= m ) dfs( x + 1, 0, k );else {dfs( x, y + 1, k );bool flag = 1;for( int i = -1;i <= 1;i ++ )for( int j = -1;j <= 1;j ++ )if( inside( x + i, y + j ) and ch[x + i][y + j] == 'X' )flag = 0;if( flag ) {ch[x][y] = 'X';dfs( x, y + 1, -k );ch[x][y] = '.';}}
}signed main() {scanf( "%lld %lld", &n, &m );for( int i = 0;i < n;i ++ )scanf( "%s", ch[i] );for( int i = 0;i < n;i ++ )for( int j = 0;j < m;j ++ )if( ch[i][j] == 'X' ) for( int x = -1;x <= 1;x ++ )for( int y = -1;y <= 1;y ++ )if( ( x or y ) and inside( i + x, j + y ) and ch[i + x][j + y] == 'X' )return puts("0"), 0;dfs( 0, 0, 1 );printf( "%lld\n", ( ans + mod ) % mod );return 0;
}

[CQOI2012] 局部极小值（状压DP + 容斥 + 搜索）相关推荐

bzoj2560串珠子状压dp+容斥(?)
2560: 串珠子 Time Limit: 10 Sec Memory Limit: 128 MB Submit: 515 Solved: 348 [Submit][Status][Discuss ...
[蓝桥杯][算法提高VIP]Sharing Chocolate（状压dp记忆化搜索）
题目描述每天,巧克力在它的许多形式上被全世界数百万人分享.它是一个真正普遍的糖果,实际上在世界上每个国家都能得到. 你发现唯一比吃巧克力更好的事情是把它分享给朋友.不幸的是,你的朋友非常挑剔,有着不 ...
[九省联考2018]一双木棋chess（状压DP+博弈论+Minmax搜索）
洛谷题目传送门解题思路 nnn和mmm都很小,很容易想到状压DP 题目描述的规则大致是长这个样子也就是说填的位置单调递减,我们大胆一点,设一个10进制的状压,表示每一列填到了那个位置,那么这个数要 ...
[蓝桥杯][算法提高VIP]金陵十三钗（状压dp记忆化搜索）
题目描述在电影<金陵十三钗>中有十二个秦淮河的女人要自我牺牲代替十二个女学生去赴日本人的死亡宴会.为了不让日本人发现,自然需要一番乔装打扮.但由于天生材质的原因,每个人和每个人之间的相似 ...
【BZOJ-2669】局部极小值状压DP + 容斥原理
2669: [cqoi2012]局部极小值 Time Limit: 3 Sec Memory Limit: 128 MB Submit: 561 Solved: 293 [Submit][Stat ...
牛客国庆集训派对Day2: E. 数据排序（状压DP+记忆化搜索）
E. 数据排序题目描述机器学习通常需要用到大量的人工标注好的数据进行训练.现在有这么一个数据集,有 N 个张照片,每张照片中都有一个模特.某个研究员想要训练一个机器学习算法,能够根据照片对模特的魅 ...
BZOJ 2560: 串珠子 (状压DP+枚举子集补集+容斥)
(Noip提高组及以下),有意者请联系Lydsy2012@163.com,仅限教师及家长用户. 2560: 串珠子 Time Limit: 10 Sec Memory Limit: 128 MB Su ...
BZOJ 2669 Luogu P3160 [CQOI2012]局部极小值 (容斥原理、DP)
题目链接 (bzoj) https://www.lydsy.com/JudgeOnline/problem.php?id=2669 (luogu) https://www.luogu.org/prob ...
【uoj#37/bzoj3812】[清华集训2014]主旋律状压dp+容斥原理
题目描述求一张有向图的强连通生成子图的数目对 $10^9+7$ 取模的结果. 题解状压dp+容斥原理设 $f[i]$ 表示点集 $i$ 强连通生成子图的数目,容易想到使用总方案数 $2^{sum ...

[CQOI2012] 局部极小值（状压DP + 容斥 + 搜索）

problem

solution

code

[CQOI2012] 局部极小值（状压DP + 容斥 + 搜索）相关推荐

最新文章

热门文章