最小奖励
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;
}