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