bzoj1097[POI2007]旅游景点atr spfa+状压DP

题意：给出一个图，要求一定经过一些点，而且经过顺序有要求，点的范围是2到k+1，问最短路径。

这题卡常卡的我蛋都碎了都没卡过，最后懒得卡了。
其实很简单的一道题目如果不卡常的话，先对于1..k+1跑spfa以后，把题目要求的遍历顺序加入状态中，直接转移就可以了，挺显然的。
代码：(TLE)

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=4e5+5;
const int inf=0x7fffffff;
int n,m,k;
typedef long long ll;
int val[N];
int head[20005],next[N],go[N];
int dis[30][20005],d[20005];
int bin[21],state[20005];
int f[30][1<<21];
bool vis[20005];
int read()
{
    int x=0,f=1; char ch=getchar();
    while (ch<'0' || ch>'9') {if (ch=='-') f=-1; ch=getchar();}
    while (ch>='0' && ch<='9') {x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
inline void spfa(int s)
{
    queue<int>q;
    for(int i=1;i<=n;i++)d[i]=inf;
    d[s]=0;
    q.push(s),vis[s]=1;
    while (!q.empty())
    {
        int x=q.front();
        q.pop();
        for (int i=head[x];i;i=next[i])
        {
            int v=go[i];
            if (d[v]>d[x]+val[i])
            {
                d[v]=d[x]+val[i];
                if (!vis[v])
                {
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
        vis[x]=0;
    }
    for(int i=1;i<=n;i++)dis[s][i]=d[i];
}
int tot=0;
inline void add(int x,int y,int z)
{
    go[++tot]=y;
    val[tot]=z;
    next[tot]=head[x];
    head[x]=tot;
}
inline int mn(int a,int b)
{
    if (a==-1)return b;
    if (a>b)return b;
    else return a;
}
int main()
{
    n=read(),m=read(),k=read();
    fo(i,1,m)
    {
        int x=read(),y=read(),z=read();
        add(x,y,z);
        add(y,x,z);
    }
    fo(i,1,k+1)spfa(i);
    int q=read();
    fo(i,1,q)
    {
        int x=read(),y=read();
        state[y]|=(1<<(x-2));
    }
    memset(f,-1,sizeof(f));
    f[1][0]=0;
    int tot=1<<k;
    fo(s,0,tot-1)
    {
        fo(i,1,k+1)
        if(f[i][s]!=-1)
        {
            fo(j,2,k+1)
            if ((s&state[j])==state[j])
            f[j][s|(1<<(j-2))]=mn(f[j][s|(1<<(j-2))],f[i][s]+dis[i][j]);        
        }
    }
    int ans=inf;
    fo(i,1,k+1)
    if (f[i][tot-1]!=-1)ans=mn(ans,f[i][tot-1]+dis[i][n]);
    printf("%d\n",ans);
    return 0;
}