HDU 3001 状压(三进制）

最新推荐文章于 2024-12-10 21:50:35 发布

beihai2013

最新推荐文章于 2024-12-10 21:50:35 发布

阅读量398

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/beihai2013/article/details/51136426

----动态规划---- 专栏收录该内容

46 篇文章

订阅专栏

本文深入探讨了状压DP（状态压缩动态规划）的应用，特别关注于如何利用三进制表示来优化访问次数的问题。通过实例分析，详细解释了算法的实现过程和背后的逻辑思维，旨在提升读者在解决复杂动态规划问题时的技巧。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

/*
状压DP。
三进制表示访问几次
*/

#include <bits/stdc++.h>
using namespace std;
#define inf ((long long)1000000007 * 1000000007)
#define LL long long
const int MAXN = 10 + 2;
const int MAXM = 120000;
LL dp[MAXM][MAXN];
LL g[MAXN][MAXN];
int thr[MAXN];
int n, m;
void init()
{
    thr[0] = 1;
    for(int i = 1 ; i < MAXN ; i++) thr[i] = thr[i - 1] * 3;
}
int main()
{
    init();
    while(scanf("%d%d", &n, &m) != EOF) {
        memset(g, -1, sizeof g);
        for(int i = 0 ; i < n ; i++) g[i][i] = 0;
        for(int i = 0 ; i < m ; i++) {
            int u, v;
            LL w;
            scanf("%d%d%I64d", &u, &v, &w); u--, v--;
            if(g[u][v] == -1 || g[u][v] > w) {
                g[u][v] = g[v][u] = w;
            }
        }

        memset(dp, -1, sizeof dp);
        for(int i = 0 ; i < n ; i++) dp[thr[i]][i] = 0;
        for(int i = 0 ; i < thr[n] ; i++) {
            for(int j = 0 ; j < n ; j++) {
                if((i / thr[j]) % 3 == 0 || dp[i][j] == -1) continue;
                for(int k = 0 ; k < n ; k++) {
//                    if(i == 1) {
//                        printf("ts = %d, i / thr[k] % 3 = %d, j = %d\n", i + thr[k], (i / thr[k]) % 3, j);
//                        printf("thr[k] = %d, k = %d, g[j][k] = %I64d\n", thr[k], k, g[j][k]);
//                    }
                    if((i / thr[k]) % 3 == 2 || g[j][k] == -1) continue;
                    int ts = i + thr[k];
//                    if(i == 1)
//                    printf("ts = %d\n", ts);
                    if(dp[ts][k] == -1 || dp[ts][k] > dp[i][j] + g[j][k]) {
//                        printf("i = %d, j = %d, ts = %d, k = %d\n", i, j, ts, k);
                        dp[ts][k] = dp[i][j] + g[j][k];
                    }
//                    system("pause");
                }
            }
        }
        LL ans = -1;
        for(int i = 0 ; i < n ; i++) {
            LL temp = dp[thr[n] - 1][i];
            if(temp != -1)
                if(temp < ans || ans == -1) ans = temp;
        }
        printf("%I64d\n", ans);
    }
    return 0;
}