1. 程式人生 > >[luogu3959 noip2017] 寶藏 (狀壓dp)

[luogu3959 noip2017] 寶藏 (狀壓dp)

return mem inline line 最小 noip void code clas

Description

技術分享圖片

Input

第一行兩個用空格分離的正整數 n,m ,代表寶藏屋的個數和道路數。

接下來 m 行,每行三個用空格分離的正整數,分別是由一條道路連接的兩個寶藏 屋的編號(編號為 1?n ),和這條道路的長度 v 。

Output

一個正整數,表示最小的總代價。

Sample Input 1

4 5
1 2 1
1 3 3
1 4 1
2 3 4
3 4 1

Sample Output 1

4

Sample Input 2

4 5
1 2 1
1 3 3
1 4 1
2 3 4
3 4 2

Sample Output 2

5

HINT

技術分享圖片

Solution

n最大12顯然狀壓至於思路 直接暴力枚舉ヽ( ̄▽ ̄)?居然就過了。。

註意:
兩點之間有重邊,要取最小的 \(40*1\) QAQ

有零邊,要建鄰接矩陣 \(60*1\) QAQ

Code

//By Menteur_Hxy
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define M(a,b) memset(a,(b),sizeof(a))
#define F(i,a,b) for(register int i=(a);i<=(b);i++)
using namespace std;
typedef long long LL;

int read() {
    int x=0,f=1; char c=getchar();
    while(!isdigit(c)) {if(c=='-')f=-f;c=getchar();}
    while(isdigit(c)) x=(x<<1)+(x<<3)+c-48,c=getchar();
    return x*f;
}

const int INF=0x3f3f3f3f;
int n,m;
int dis[15],edg[15][15],in[15][15];
LL dp[1<<15],ans=INF; 

void dfs(int x) {
    F(i,0,n-1) if(x&(1<<i)) 
        F(j,0,n-1) if(in[i][j] && i!=j && !(x&(1<<j))) 
            if(dp[x|(1<<j)]>dp[x]+dis[i]*edg[i][j]) {
                int reg=dis[j];
                dp[x|(1<<j)]=dp[x]+dis[i]*edg[i][j];
                dis[j]=dis[i]+1;
                dfs(x|(1<<j));
                dis[j]=reg;
            }
}

int main() {
    n=read(),m=read();
    M(edg,0x3f);
    F(i,1,m) {
        int u=read()-1,v=read()-1,w=read();
        if(w<edg[u][v]) in[u][v]=in[v][u]=1,edg[u][v]=edg[v][u]=w;
    }
    F(i,0,n-1) {
        M(dis,0x3f); M(dp,0x3f);
        dis[i]=1; dp[1<<i]=0;
        dfs(1<<i);
        ans=min(ans,dp[(1<<n)-1]);      
    }
    printf("%lld",ans);
    return 0;
}

[luogu3959 noip2017] 寶藏 (狀壓dp)