[BZOJ]1097: [POI2007]旅游景点atr

　FGD想从成都去上海旅游。在旅途中他希望经过一些城市并在那里欣赏风景，品尝风味小吃或者做其他的有趣
的事情。经过这些城市的顺序不是完全随意的，比如说FGD不希望在刚吃过一顿大餐之后立刻去下一个城市登山，
而是希望去另外什么地方喝下午茶。幸运的是,FGD的旅程不是既定的，他可以在某些旅行方案之间进行选择。由于
FGD非常讨厌乘车的颠簸，他希望在满足他的要求的情况下，旅行的距离尽量短，这样他就有足够的精力来欣赏风
景或者是泡MM了^_^.整个城市交通网络包含N个城市以及城市与城市之间的双向道路M条。城市自1至N依次编号，道
路亦然。没有从某个城市直接到它自己的道路，两个城市之间最多只有一条道路直接相连，但可以有多条连接两个
城市的路径。任意两条道路如果相遇，则相遇点也必然是这N个城市之一，在中途，由于修建了立交桥和下穿隧道
，道路是不会相交的。每条道路都有一个固定长度。在中途，FGD想要经过K(K<=N-2)个城市。成都编号为1，上海
编号为N,而FGD想要经过的N个城市编号依次为2,3,…,K+1.举例来说，假设交通网络如下图。FGD想要经过城市2,3,
4,5，并且在2停留的时候在3之前，而在4,5停留的时候在3之后。那么最短的旅行方案是1-2-4-3-4-5-8,总长度为1
9。注意FGD为了从城市2到城市4可以路过城市3,但不在城市3停留。这样就不违反FGD的要求了。并且由于FGD想要
走最短的路径，因此这个方案正是FGD需要的。

思路：状压DP应该很容易想到吧，我的代码常数太大了，还是不要看吧

#include <cstdio>
#include <queue>
using namespace std;
inline int min(int a,int b)
{
    return a<b?a:b;
}
inline char tc(void)
{
    static char fl[1000000],*A=fl,*B=fl;
    return A==B&&(B=(A=fl)+fread(fl,1,1000000,stdin),A==B)?EOF:*A++;
}
inline int read(void)
{
    int a=0;static char c;
    while((c=tc())<'0'||c>'9');
    while(c>='0'&&c<='9')
        a=a*10+c-'0',c=tc();
    return a; 
}
struct Y
{
    int to,cost,next;
}w[400001];
int now,n,m,k,ans=2e9,head[200001],tot,dis[22][200001],cnt,a[22],dp[21][1048756];
char b[200001];
queue<int>que;
void spfa(int x)
{
    while(!que.empty())
        que.pop();
    for (int i=1;i<=n;++i)
        dis[x][i]=1073741823;
    b[x]=1,que.push(x),dis[x][x]=0;
    while(!que.empty())
    {
        now=que.front(),que.pop(),b[now]=0;
        for (int i=head[now];i;i=w[i].next)
            if(dis[x][w[i].to]>dis[x][now]+w[i].cost)
            {
                dis[x][w[i].to]=dis[x][now]+w[i].cost;
                if(!b[w[i].to])
                    b[w[i].to]=1,que.push(w[i].to);
            }
    }
    return ;
}
int main(void)
{
    register int i,j,l,x,y,z;
    n=read(),m=read(),k=read();
    for (i=1;i<=m;++i)
        x=read(),y=read(),z=read(),w[++tot].to=y,w[tot].cost=z,w[tot].next=head[x],head[x]=tot,w[++tot].to=x,w[tot].cost=z,w[tot].next=head[y],head[y]=tot;
    for (i=1;i<=k+1;++i)
        spfa(i);
    for (cnt=read(),i=1;i<=cnt;++i)
        x=read(),y=read(),a[y]+=1<<x-2;
    for (i=1;i<(1<<k);++i)
        for (j=1;j<=k;++j)
            dp[j][i]=1073741823;
    for (i=1;i<(1<<k);++i)
        for (j=1;j<=k;++j)
            if(i&(1<<j-1))
            {
                if(a[j+1]==0&&(1<<j-1)==i)
                    dp[j][i]=min(dp[j][i],dis[1][j+1]);
                for (l=1;l<=k;++l)
                    if((i&(1<<l-1))&&(a[j+1]&i)==a[j+1]&&l!=j)
                        dp[j][i]=min(dp[j][i],dp[l][i-(1<<j-1)]+dis[l+1][j+1]);
            }
    for (i=1;i<=k;++i)
        ans=min(ans,dp[i][(1<<k)-1]+dis[i+1][n]);
    if(k!=0)
        printf("%d",ans);
    else
        printf("%d",dis[1][n]);
    return 0;
}