教皇与僭称教皇

教皇与僭称教皇

---

题目描述

稳固了后方之后，光下巴开始寻求加冕为神圣罗马帝国皇帝，$1153$ 年发动了第一次远
征意大利。
在进军途中，他接到了教皇尤金三世的求救信，原来是教士阿诺德领导人民夺取了城市政权，选举了罗马元老院成员和执政官，主张教会放弃领地过使徒的简朴生活。
这真是天赐的机会，金天成加速进军镇压了这次起义，以拯救者的身份进入罗马城。
但新任教皇哈德良四世是个不懂变通的倔老头，他在为金天成加冕时，坚持要金天成按照惯例为教皇牵马、扶镫，激怒了金天成，结果使庆典会场变成了战场，近千人被杀。
但不管过程如何狼狈，金天成已通过加冕获得了“神圣罗马帝国皇帝”的称号（在他之前的德皇，一直只称为“罗马帝国皇帝”，从他开始又冠上“神圣”二字） 。
好吧，现在出题人要管管这不知道如何狼狈过程了，教堂的区域很少，只有 $N$ 块区域。
然后某些区域之间有通路。这 $N$ 块区域中每块区域都有敌军，$N$ 个区域之间有 $M$ 条通道。
金天成大帝现在需要以闪电般的速度将这 $N$ 个区域中敌军迅速歼灭。
金天成大帝的军队比教皇军队要强很多，他们到达某个区域就能瞬间歼灭该区域的军队。
但是走第 $i$ 条通道花的时间为 $T_i$ 。
由于某种特殊的作者编造不出来的原因导致金天成大帝不能回到以前曾经过的区域，所以金天成大帝希望有一种方案，使得他能最快歼灭教皇的军队，如果不存在方案，则输出 -1。
一开始金天成大帝所在区域为 $1$ 号区域。

---

输入格式

第一行两个数字 $N$，$M$。
接下来 $M$ 行，每行三个数字 $A$,$B$,$C$，表示 $A$ 到 $B$ 有一条通路，耗时为 $C$。

---

输出格式

输出一行，一个数字，表示金天成大帝干掉教皇军队的最少时间。

---

样例输入

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

---

样例输出

18

---

样例解释

从 $1$ 开始一条路走到黑就 OK。

---

数据范围

$30\%$ 的数据：$1≤N≤8$，$1≤M≤28$。
$100\%$ 的数据：$1≤N≤16$，$1≤M≤120$。$1≤C≤10^8$

---

Solution

TSP问题，用状态压缩可以解决。

---

Code

#include <iostream>
#include <cstdio>
#include <cstring>

#define LL long long
#define Min(x,y) ((x)<(y)?(x):(y))
#define INT64_MAX 9223372036854775807LL

using namespace std;

int n,m,cnt,head[20],nxt[400],data[400],flow[400];
LL ans=INT64_MAX,f[20][70000];

void add(int x,int y,int z){
    nxt[cnt]=head[x];data[cnt]=y;flow[cnt]=z;head[x]=cnt++;
    nxt[cnt]=head[y];data[cnt]=x;flow[cnt]=z;head[y]=cnt++;
}

LL dp(int x,int about){
    if(x==1)return f[x][about]=0;
    if(f[x][about]!=-1)return f[x][about];
    LL minx=INT64_MAX;
    for(int i=head[x];i!=-1;i=nxt[i]){
        int d=data[i];
        if((about&(1<<(d-1)))&&(d!=1||((about^(1<<(x-1)))==1))){
            LL tmp=dp(d,about^(1<<(x-1)));
            if(tmp==INT64_MAX)continue;
            if(minx==-1||minx>tmp+flow[i])
                minx=tmp+flow[i];
        }
    }
    return f[x][about]=minx;
}

int main(){

    freopen("pope.in","r",stdin);   
    freopen("pope.out","w",stdout);

    memset(head,-1,sizeof head);
    memset(f,-1,sizeof f);

    scanf("%d%d",&n,&m);
    const int end=(1<<n)-1;

    for(int i=1;i<=m;i++){
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        add(x,y,z);
    }

    for(int i=2;i<=n;i++){
        LL tmp=dp(i,end);
        ans=Min(ans,tmp);
    }

    if(end==1)ans=Min(ans,dp(1,end));

    if(ans==INT64_MAX)ans=-1;

    printf("%lld\n",ans);

    return 0;
}