1. 程式人生 > >【BZOJ 1922】【SDOI 2010】大陸爭霸(有限制的最短路)

【BZOJ 1922】【SDOI 2010】大陸爭霸(有限制的最短路)

(誰告訴我這是個分層最短路來著的???)

進入一個城市必須要這個城市所有結界都打破才能進

那我們可以邊炸邊走 也就是說 我們可以維護d1陣列:走到的時間(結界可能沒炸完) d2陣列:可進入的時間(結界都被炸完了) d陣列:真實的到達時間

容易發現d[i]=max(d1[i],d2[2])

那我們就用dijkstra維護d陣列 每次都要更新當前點所保護的城市的可進入時間 以及當前點連的城市的到達時間 當一個城市的結界全部被炸了 就入堆

注意 這種有限制的最短路最好用dijkstra 因為每個點只會走一次

#include<bits/stdc++.h>
#define N 3005
#define M 70005
using namespace std;
template<class T>
inline void read(T &x)
{
    x=0;
    static char ch=getchar();
    while(ch<'0'||ch>'9')   ch=getchar();
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
}
struct Edge
{
    int to,next,val;
}edge[M];
int n,m,tot,first[N],a[N][N],pro[N],p[N];//pro:每個城市保護多少城市   p:每個城市受多少城市保護 
inline void addedge(int x,int y,int z)
{
    tot++; edge[tot].to=y; edge[tot].next=first[x]; edge[tot].val=z; first[x]=tot;
}
typedef pair<int,int> Pair;
priority_queue<Pair,vector<Pair>,greater<Pair> > heap;
bool vis[N];
int d1[N],d2[N],d[N];//d1:到達時間 d2:可進入時間 
void Dijkstra()
{
    memset(d,63,sizeof(d)); memset(d1,63,sizeof(d1));
    heap.push(make_pair(0,1)); 
    d1[1]=0; d[1]=0;
    while(!heap.empty())
    {
        int now=heap.top().second; heap.pop();
        if(vis[now]) continue; 
        vis[now]=true;
        for(int u=first[now];u;u=edge[u].next)
        {
            int v=edge[u].to;
            if(d[now]+edge[u].val<d1[v])    //到達時間 
            {
                d1[v]=d[now]+edge[u].val;
                d[v]=max(d1[v],d2[v]);
                if(!p[v])   heap.push(make_pair(d[v],v));
            }
        }
        for(int i=1;i<=pro[now];i++)    //更新保護城市 
        {
            int v=a[now][i];
            p[v]--,d2[v]=max(d2[v],d[now]),d[v]=max(d1[v],d2[v]);
            if(!p[v])   heap.push(make_pair(d[v],v));
        }
    }
    cout<<d[n]<<endl;   
}
int main()
{
    read(n); read(m);
    for(int i=1,u,v,w;i<=m;i++) 
    {
        read(u); read(v); read(w);
        addedge(u,v,w);
    }
    for(int i=1;i<=n;i++)
    {
        read(p[i]);
        for(int j=1,x;j<=p[i];j++)
        {
             read(x);
             a[x][++pro[x]]=i;
        }
    }
    Dijkstra();
    return 0;
}