ACdream 1083 有向无环图dp

题目链接：点击打开链接

人民城管爱人民

Time Limit: 12000/6000 MS (Java/Others) Memory Limit: 128000/64000 KB (Java/Others)

SubmitStatistic Next Problem

Problem Description

一天GG正在和他的后宫之一的MM在外面溜达，MM突然说了一句，“我想吃鸡蛋灌饼”……当他们吃的正high的时候，城管出现了！作为传说中的最强军事力量，卖鸡蛋灌饼的小贩在他们面前也只算是战力为的5的渣滓，一秒钟就被秒杀了……

在这场屠杀中，GG和他的后宫本来只是围观群众，但是不幸的是，城管看到了GG胃里的鸡蛋灌饼，他们要逮捕GG！但是GG显然不能让他们如愿，于是GG带着后宫开始了往大运村的逃亡之旅。

整个地图有n个路口，灌饼摊在0号路口，大运村在n-1号路口。有m条只能单向通过的道路连接这n个路口，每条道路用一个正整数表示走过需要的时间。整个地图没有环路，但两个路口之间可能有多条通路。现在GG希望以最短的时间到大运村，但不幸的是，城管为了抓住他动用了卫星对他进行空中跟踪，并且会在某一时刻空降到某一条道路上进行封锁（封锁会在瞬间完成，可惜动静太大了GG也能在第一时间知道哪条道路被封锁了），之后这条路就无法通过了。在整个行动中只会出现一次空降，而且不会在GG经过这条道路的时候进行封锁，也就是说，不会在GG在某条路上走了一半的时候封锁这条路。而且，城管们希望尽可能的延缓GG到达大运村的时间。

现在GG希望知道，自己多久能到达大运村，方便安排之后和其他后宫的约会。

注意双方是以博弈的思想来进行选择，即GG希望时间最短，城管希望时间最长，而且他们都非常聪明会做出最佳的选择。

Input

输入第一行为数据组数T（T<=30）。

每组数据第一行包含两个整数n,m（2<=n<=10000，1<=m<=100000），表示路口数和道路数。之后m行描述了所有的道路，每行有三个整数u,v,w（0<=u,v<n，0<w<=1000），表示路口u到路口v有一条需要w时间走过的道路。

Output

对于每组数据输出一个整数，表示GG最后到达大运村需要的时间。如果GG无法到达大运村，输出-1。

Sample Input

25 60 1 11 2 12 4 11 4 30 3 23 4 13 40 1 10 1 21 2 31 2 4

Sample Output

45

Source

Dshawn

思路：

1、首先这个dp一定是逆向拓扑序(把所有边反向进行拓扑排序），这样我们在计算u点时，通过有向边(u->v)把u点的状态从v点转移过来，能保证v点一定是已计算过的

2、用dis[i]表示i点到终点的最短路，dp[i]表示i点到终点的删边最短路(这里所谓的删边，删除的边是在i-终点的路径上）

3、那么我们设GG站在u点，通过u->v的边把状态从v转移过来。

4、计算dis[u]比较容易，就是min(dis[v]+edge[i].dis)

5、计算dp[u]：根据定义dp[u] 是删边最短路，则删除的边是在u-终点的路径上，那么显然有2种情况

  1）删除u-v这种与u相连的边

  2）删除v-终点上的边

对于2）：那么GG有主动权，显然是选择edge[i].dis + dp[v] 中最小的（我们设这个最小值为smin）

对于1）：军队拥有主动权，军队一定选择删掉一条边，那么删完以后GG自然还是选择走最短路，也就是GG只能选择edge[i].dis+dis[v] 中次小值（设这个值为dmin

这里注意的是，军队拥有的主动权是可以任意选择一条边，所以dp[u]=max(dmin, smin)

#include<stdio.h>#include<iostream>#include<algorithm>#include<string.h>#include<set>#include<queue>#include<math.h>#include<map>using namespace std;#define N 10050#define M 100500#define inf 10000000struct Edge{int from, to, dis, nex;}edge[M<<1];int head[N], edgenum;void add(int u, int v, int d){Edge E={u,v,d,head[u]};edge[edgenum] = E;head[u] = edgenum++;}int n,m;int sp[N], in[N];void topsort(){//把图进行拓扑序queue<int>q;for(int i = 0; i < n; i++)if(in[i]==0)q.push(i);int top = 0;while(!q.empty()){int u = q.front(); q.pop();sp[top++] = u;for(int i = head[u]; ~i; i = edge[i].nex)if(i&1){int v = edge[i].to;in[v]--;if(!in[v])q.push(v);}}}int dp[N][2], a[M];//dp[u][0]表示1到u的最短路 dp[u][1]表示1到u的次小最短路void work(int u){if(u == n-1){ dp[u][0] = dp[u][1] = 0; return ; }int top = 0, smin = inf;for(int i = head[u]; ~i; i = edge[i].nex)if(!(i&1)){int v = edge[i].to;dp[u][0] = min(dp[u][0], dp[v][0]+edge[i].dis);smin = min(smin, dp[v][1]+edge[i].dis);a[top++] = dp[v][0] + edge[i].dis;}if(top<2){dp[u][1] = inf; return ;}int fir = a[0], sec = a[1];if(fir>sec)swap(fir,sec);for(int i = 2; i < top; i++){if(a[i]<=fir)sec = fir, fir = a[i];else sec = min(sec, a[i]);}dp[u][1] = max(smin, sec);}void init(){memset(head, -1, sizeof head);edgenum = 0;}int main(){int T, i, j, u, v, d;scanf("%d",&T);while(T--){init();memset(in, 0, sizeof in);scanf("%d %d",&n,&m);while(m--){scanf("%d %d %d",&u,&v,&d);add(u,v,d);add(v,u,d);in[u]++;}topsort();for(i = 0; i < n; i++)dp[i][0] = dp[i][1] = inf;for(i = 0; i < n; i++)work(sp[i]);if(dp[0][1]==inf)dp[0][1] = -1;printf("%d\n",dp[0][1]);}return 0;}