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

最新推荐文章于 2022-07-15 08:56:49 发布

ikrvxt

最新推荐文章于 2022-07-15 08:56:49 发布

阅读量637

点赞数 1

分类专栏： # 状压DP # 容斥文章标签：深度优先动态规划算法容斥搜索

本文链接：https://blog.csdn.net/Emm_Titan/article/details/123483512

版权

状压DP 同时被 2 个专栏收录

18 篇文章 0 订阅

订阅专栏

容斥

9 篇文章 1 订阅

订阅专栏

本文详细解析了洛谷P3160题目的解决方案，主要涉及动态规划和容斥原理的应用。在不超过28!种可能的情况下，通过非暴力计算避免了暴力填充的复杂度。文章首先介绍了问题背景，然后阐述了动态规划的状态转移方程，接着讨论了如何处理局部极小值点的填充约束，并指出了一种错误情况——非法的局部极小值点。最后，通过搜索修正错误并给出完整代码，实现了有效的解题策略。

摘要由CSDN通过智能技术生成

problem

luogu-P3160

solution

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

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

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

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

$dp(s,i+1)\leftarrow dp(s,i)*(cnt-i)$ 。

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

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

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

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

具体可见下方代码。

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

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

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

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

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

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

所以 $k$ 并不会多大。

我们完全可以搜索做。

在这里插入图片描述

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;
}