2018.11.06 bzoj1097: [POI2007]旅游景点atr（最短路+状压dp）

传送门

预处理出不能在每个点停留之后才停留的点的状态。

对kkk个点都跑一次最短路存下来之后只需要简单状压一下就能过了吐槽原题空间64MB蒟蒻无能为力

然后用fillfillfill赋极大值的时候当memsetmemsetmemset用赋成了很小的值导致挑了一个晚上都不知道哪里错了。

于是f[i][j]f[i][j]f[i][j]表示现在停留在iii这个点并且已经停留过的状态为jjj，随便转移一波就做完了。

代码：

#include<bits/stdc++.h>
using namespace std;
inline int read(){
    int ans=0;
    char ch=getchar();
    while(!isdigit(ch))ch=getchar();
    while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
    return ans;
}
const int N=2e4+5,M=2e5+5,K=22;
int n,m,k,up,dis[K][N],first[N],cnt=0,f[K][1<<(K-1)],ban[K];
bool vis[N];
struct edge{int v,w,next;}e[M<<1];
inline void add(int u,int v,int w){e[++cnt].v=v,e[cnt].w=w,e[cnt].next=first[u],first[u]=cnt;}
inline void dijkstra(int s){
    fill(dis[s],dis[s]+n+1,0x3f3f3f3f),fill(vis,vis+n+1,0);
    priority_queue<pair<int,int> >q;
    q.push(make_pair(dis[s][s]=0,s));
    while(!q.empty()){
        int x=q.top().second;
        q.pop();
        if(vis[x])continue;
        vis[x]=1;
        for(int i=first[x];~i;i=e[i].next){
            int v=e[i].v;
            if(dis[s][v]>dis[s][x]+e[i].w)dis[s][v]=dis[s][x]+e[i].w,q.push(make_pair(-dis[s][v],v));
        }
    }
}
inline int dfs(int pos,int sta){
    if(f[pos][sta])return f[pos][sta];
    if(sta==up)return f[pos][sta]=dis[pos][n-1];
    f[pos][sta]=0x3f3f3f3f;
    for(int i=1;i<=k;++i)if((sta&ban[i])==ban[i])f[pos][sta]=min(f[pos][sta],dfs(i,sta|(1<<(i-1)))+dis[pos][i]);
    return f[pos][sta];
}
int main(){
    memset(first,-1,sizeof(first)),n=read(),m=read(),k=read(),up=(1<<k)-1;
    for(int i=1,u,v,w;i<=m;++i)u=read()-1,v=read()-1,w=read(),add(u,v,w),add(v,u,w);
    for(int i=read(),u,v;i;--i)u=read()-1,v=read()-1,ban[v]|=1<<(u-1);
    for(int i=0;i<=k;++i)dijkstra(i);
    cout<<dfs(0,0);
    return 0;
}