USACO 玛丽卡(最短路+枚举)

USACO 玛丽卡

麦克找了个新女朋友，玛丽卡对他非常恼火并伺机报复。

    因为她和他们不住在同一个城市，因此她开始准备她的长途旅行。

    在这个国家中每两个城市之间最多只有一条路相通，并且我们知道从一个城市到另一个城市路上所需花费的时间。

    麦克在车中无意中听到有一条路正在维修，并且那儿正堵车，但没听清楚到底是哪一条路。无论哪一条路正在维修，从玛丽卡所在的城市都能到达麦克所在的城市。

    玛丽卡将只从不堵车的路上通过，并且她将按最短路线行车。麦克希望知道在最糟糕的情况下玛丽卡到达他所在的城市需要多长时间，这样他就能保证他的女朋友离开该城市足够远。

编写程序，帮助麦克找出玛丽卡按最短路线通过不堵车道路到达他所在城市所需的最长时间(用分钟表示)。

输入

第一行有两个用空格隔开的数N和M，分别表示城市的数量以及城市间道路的数量。1≤N≤1000，1≤M≤N*(N-1)/2。城市用数字1至N标识，麦克在城市1中，玛丽卡在城市N中。

接下来的M行中每行包含三个用空格隔开的数A，B和V。其中1≤A，B≤N，1≤V≤1000。这些数字表示在A和城市B中间有一条双行道，并且在V分钟内是就能通过。

输出

   输出文件的第一行中写出用分钟表示的最长时间，在这段时间中，无论哪条路在堵车，玛丽卡应该能够到达麦克处，如果少于这个时间的话，则必定存在一条路，该条路一旦堵车，玛丽卡就不能够赶到麦克处。

样例

MARICA.IN

5 7

1 2 8

1 4 10

2 3 9

2 4 10

2 5 1

3 4 7

3 5 10

MARICA.OUT

27

解题报告

处于人道，我们要帮助玛丽卡找出最短路（手动滑稽）。

这一题就是简单的最短路模型。枚举删去的边，用堆优化的dijkstra来计算最短路。

算法复杂度为O(nmlogn)

#include<bits/stdc++.h>
#define Pair pair<int,int>
#define MAXN 1000+10
#define MAXM 600000+1
using namespace std;
int n,m,num,head[MAXN],s,t,pre[MAXN],dis[MAXN],v[MAXM];
int bian[MAXM],ans;
struct Edge{
    int dis,next,to,exi,from;
}edge[MAXM];
void add(int from,int to,int dis)
{
    edge[++num].next=head[from];
    edge[num].to=to;
    edge[num].dis=dis;
    edge[num].from=from;
    head[from]=num;
    edge[num].exi=1;
}
void dij()
{
    memset(dis,0,sizeof(dis));
    memset(pre,0,sizeof(pre));
    memset(v,0,sizeof(v));
    priority_queue<Pair,vector<Pair>,greater<Pair> > h;
    for(int i=1;i<=n;i++) dis[i]=2147483647;
    dis[s]=0;
    h.push(Pair(dis[s],s));
    while(h.size()>0)
    {
        int k=h.top().second;h.pop();
        if(v[k]) continue;
        v[k]=1;
        for(int i=head[k];i;i=edge[i].next)
        if(dis[k]+edge[i].dis<dis[edge[i].to]&&edge[i].exi==1)
        {
            dis[edge[i].to]=dis[k]+edge[i].dis;
            pre[edge[i].to]=i;
            h.push(Pair(dis[edge[i].to],edge[i].to));
        }
    }
}
int main()
{

    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        add(x,y,z);
        add(y,x,z);
    }s=1;
    dij();
    for(int i=pre[n];i;i=pre[edge[i].from])
        bian[++bian[0]]=(i+1)/2;
    for(int i=1;i<=bian[0];i++)
    {
        edge[bian[i]*2].exi=0;
        edge[bian[i]*2-1].exi=0;
        dij();
        edge[bian[i]*2].exi=1;
        edge[bian[i]*2-1].exi=1;
        ans=max(dis[n],ans);
    }
    printf("%d
",ans);
    
    return 0;
}