最小奖励（拓扑、广搜、spfa）

最小奖励

         Jzabc更是一个狂热的旅游爱好者。这次他想去一个诡异的地方。这个地方有n个村庄，编号为1到n。此刻他在1号村庄，他想去n号村庄，这n个村庄之间有m条单向道路。通过某条道路时，可能会花费一些钱作为“买路钱”，可是通过一些道路时，不仅不需要交“买路钱”，而且会得到一些奖励（这就是诡异的地方？）。当然两个村庄之间可能有多条直接相连的道路。但是每一条路只能通过一次。这些村庄有这样的特性，从任何一个村庄出发，沿着任一条路径走都不会回到出发点。找到一条路径从1号村庄到n号村庄后，他需要计算一共得到多少奖励，一共交了多少“买路钱”。如果得到的总奖励钱数大于交的“买路钱”数，那么称走这一条路径可以得到奖励；相反，如果前者小于后者，那么称走这一条路径需要花费。如果两者相等，那么Jzabc不会选择这一条路（我也不知道他为什么不选这一条路径）。不会出现所有的路径两者都相等。他又需要你的帮助，让你找一条路径，使他得的奖励最小（为什么不是最大的奖励？），并输出最小奖励。如果找不到一条路径能使他得到奖励，那么就找一条路径使他得到的花费最大（为什么就不是最小花费？），并输出最大花费。

【输入】

第一行两个正整数n，m，中间用空格隔开。n是村庄的总数，m是道路总数。

以下m行，每行三个正整数x，y，w，描述一条道路的信息，中间用空格隔开。x是此条道路的起点，y是终点。若w为正，则是通过此条道路得到的奖励钱数；若w为负，则是通过此条道路交的“买路钱”数；w不为0，并且|w|≤10。

【输出】

一个整数。若有最小奖励，则输出。否则输出最大花费。

【样例输入1】

3 4

1 2 2

2 3 -1

2 3 3

1 3 -1

【样例输出1】

1

【样例输入2】

3 4

1 2 2

2 3 -3

2 3 -5

1 3 -2

【样例输出2】

3

【数据范围】

20%的数据，n≤20；m≤200；

50%的数据，n≤50；m≤2000；

100%的数据，n≤100；m≤20000

算法分析：本题需要求的实际上是最小奖励（优先）和最大花费（次先），根据题中的性质可以很轻易地知道这是一个有向无环图，因此很容易想到和拓扑排序有关的算法。

最小奖励：使用minaw[i]记录从起点到i点最小奖励（可以小于零及花费大于奖励）。又要保证尽量小又要大于零（不是常规的某个最大或最小），因此引入一个数组maxaw[i]表示从i到终点可以获得的最大奖励（此时就是一个裸的极值算法），对于与从 i 可以走到的点 j ，当minaw[i] + maxaw[j] + i->j 边的权值三者之和大于零时才更新minaw（因为如果到 j 点后走最大的奖励路径都是花费多的话便不满足条件）。利用拓扑排序中的入度进行更新，最后判断终点的minaw是否大于零（大于零则输出，否则进入最大花费）。

最大花费：这个算法就是裸的极值求法，先将maxco[ ]赋予极大值，然后走一遍拓扑排序来求最小值，输出终点maxco的绝对值即可。

代码较长并且采用的算法有bfs、spfa、拓扑（跑了三次图，其实三次用一种都可以），因此仅供参考。

#include<queue>
#include<cmath>
#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define name "minaw"
using namespace std;
const int MAXN=101,MAXM=20001;
int n,m,maxaw[MAXN],minaw[MAXN],maxco[MAXN];
int cnt,save[MAXN],enter1[MAXN],enter2[MAXN];
struct edgeb{
	int to,v,next;
}edge[MAXM];
vector <int> to[MAXN],val[MAXN];
void add(int x,int y,int v)
{
	edge[++cnt].to=y;
	edge[cnt].v=v;
	edge[cnt].next=save[x];
	save[x]=cnt;
	to[y].push_back(x);
	val[y].push_back(v);
	enter1[y]++;
	enter2[x]++;
}
void bfs()//预处理maxaw 
{
	queue <int> d;
	d.push(n);
	memset(maxaw,-0x7f,sizeof(maxaw));
	maxaw[n]=0;
	while(!d.empty())
	{
		int x=d.front();
		d.pop();
		for(int i=0;i<to[x].size();i++)
		{
			int y=to[x][i];
			maxaw[y]=max(maxaw[y],maxaw[x]+val[x][i]);
			enter2[y]--;
			if(enter2[y]==0)
			d.push(y);
		}
	}
}
void topdp()//最小奖励 
{
	queue <int> d;
	d.push(1);
	memset(minaw,0x7f,sizeof(minaw));
	minaw[1]=0;
	while(!d.empty())
	{
		int x=d.front();
		d.pop();
		for(int i=save[x];i;i=edge[i].next)
		{
			int y=edge[i].to;
			enter1[y]--;
			if(minaw[x]+maxaw[y]+edge[i].v>0)//!!!:更新minaw[y]的必要条件 
			{
				minaw[y]=min(minaw[y],minaw[x]+edge[i].v);
				if(enter1[y]==0)
				d.push(y);
			}
		}
	}
}
void spfa()//最大花费 
{
	for(int i=1;i<=cnt;i++)edge[i].v=-edge[i].v;
	queue <int> d;
	d.push(1);
	memset(maxco,-0x6f,sizeof(maxco));
	maxco[1]=0;
	while(!d.empty())
	{
		int x = d.front();
		d.pop();
		for(int i=save[x];i;i=edge[i].next)
		{
			int y=edge[i].to;
			if(maxco[y]<maxco[x]+edge[i].v)
			{
				maxco[y]=maxco[x]+edge[i].v;
				d.push(y);
			}
		}
	}
}
int main()
{
	freopen(name ".in","r",stdin);
	freopen(name ".out","w",stdout);
	cin>>n>>m;
	for(int i=1;i<=m;i++)
	{
		int a,b,c;
		scanf("%d%d%d",&a,&b,&c);
		add(a,b,c);
	}
	bfs();
	topdp();
	if(minaw[n]>0&&minaw[n]<=200000){
		cout<<minaw[n]<<endl;
		return 0;
	}
	spfa();
	cout<<abs(maxco[n])<<endl;
	return 0;
}