luogu P1073 最优贸易 NOIP2009 题解

题目描述

C 国有 n 个大城市和 m 条道路，每条道路连接这 n 个城市中的某两个城市。任意两个城市之间最多只有一条道路直接相连。这 m 条道路中有一部分为单向通行的道路，一部分为双向通行的道路，双向通行的道路在统计条数时也计为 1 条。

C 国幅员辽阔，各地的资源分布情况各不相同，这就导致了同一种商品在不同城市的价格不一定相同。但是，同一种商品在同一个城市的买入价和卖出价始终是相同的。

商人阿龙来到 C 国旅游。当他得知同一种商品在不同城市的价格可能会不同这一信息之后，便决定在旅游的同时，利用商品在不同城市中的差价赚回一点旅费。设 C 国 n 个城市的标号从 1~ n，阿龙决定从 1 号城市出发，并最终在 n 号城市结束自己的旅行。在旅游的过程中，任何城市可以重复经过多次，但不要求经过所有 n 个城市。阿龙通过这样的贸易方式赚取旅费：他会选择一个经过的城市买入他最喜欢的商品――水晶球，并在之后经过的另一个城市卖出这个水晶球，用赚取的差价当做旅费。由于阿龙主要是来 C 国旅游，他决定这个贸易只进行最多一次，当然，在赚不到差价的情况下他就无需进行贸易。

假设 C 国有 5 个大城市，城市的编号和道路连接情况如下图，单向箭头表示这条道路为单向通行，双向箭头表示这条道路为双向通行。

假设 1~n 号城市的水晶球价格分别为 4，3，5，6，1。

阿龙可以选择如下一条线路：1->2->3->5，并在 2 号城市以 3 的价格买入水晶球，在 3号城市以 5 的价格卖出水晶球，赚取的旅费数为 2。

阿龙也可以选择如下一条线路 1->4->5->4->5，并在第 1 次到达 5 号城市时以 1 的价格买入水晶球，在第 2 次到达 4 号城市时以 6 的价格卖出水晶球，赚取的旅费数为 5。

现在给出 n 个城市的水晶球价格，m 条道路的信息（每条道路所连接的两个城市的编号以及该条道路的通行情况）。请你告诉阿龙，他最多能赚取多少旅费。

输入输出格式

输入格式：

第一行包含 2 个正整数 n 和 m，中间用一个空格隔开，分别表示城市的数目和道路的数目。

第二行 n 个正整数，每两个整数之间用一个空格隔开，按标号顺序分别表示这 n 个城市的商品价格。

接下来 m 行，每行有 3 个正整数，x，y，z，每两个整数之间用一个空格隔开。如果 z=1，表示这条道路是城市 x 到城市 y 之间的单向道路；如果 z=2，表示这条道路为城市 x 和城市y 之间的双向道路。

输出格式：

输出文件 trade.out 共 1 行，包含 1 个整数，表示最多能赚取的旅费。如果没有进行贸易，则输出 0。

解法：

其实如果没有环的话，就是个DAG啊。直接跑一遍bfs记录走过的路上最小的买入价格和对于每到一个点时候的利润的最大值一直更新到最后即可。（显然最小买入价格有传递性，如果在到这个节点之前可以以这个价格买入，那么到这个节点之后仍然可以。但是最大卖出价格是没有的，所以保存的第二个值是最大的利润，对于每一个点只计算以当前点上面的价格卖出所产生的利润。）
那么考虑如何消掉环即可。有个经典的东西叫做Tarjan（强连通分量）
可以把环缩成点（经典题是受欢迎的牛）
然后就做完啦！
不过代码实现十分繁琐。
最后改成了用STL栈实现Tarjan以及vector存边才通过，我也不是很懂为什么。
（哦顺便这个方法比那个两次SPFA跑的快一些 但是我可能是因为用vector存边和用STL栈的缘故，所以跑的并不快。）

代码：

#include<iostream>
#include<cstdio>
#include<queue>
#include<vector>
#include<stack>
using namespace std;
const int N=100005;
const int INF=2147483647;
int n,m;
vector<int> G[N],E[N];
int scnt=0,qltn=0;
int low[N],dfn[N];
int c[N],mx[N],mm[N],co[N],f[N];
bool instack[N];
stack<int> sta;
void tarjan(int now)
{
  low[now]=dfn[now]=++scnt;
  sta.push(now);
  instack[now]=true;
  int sizz=G[now].size();
  for(int i=0;i<sizz;i++)
  {
    int v=G[now][i];
    if(!dfn[v])
    {
      tarjan(v);
      low[now]=min(low[now],low[v]);
    }
    else if(instack[v])
      low[now]=min(low[now],dfn[v]);
  }
  if(dfn[now]==low[now])
  {
    qltn++;mm[qltn]=INF;mx[qltn]=0;
    for(int i=sta.top();;i=sta.top())
    {
      co[i]=qltn,mm[qltn]=min(mm[qltn],c[i]),mx[qltn]=max(mx[qltn],c[i]);
      sta.pop(),instack[i]=false;if(i==now) break;
    }
  }
}
queue<int> q;bool inq[N];
void bfs()
{
  for(int i=1;i<=qltn;i++) f[i]=mx[i]-mm[i];
  q.push(co[1]),inq[co[1]]=true;
  while(!q.empty())
  {
    int tt=q.front();q.pop();inq[tt]=false;
    int sizz=E[tt].size();
    for(int i=0;i<sizz;i++)
    {
      int xx=E[tt][i];
      mm[xx]=min(mm[xx],mm[tt]);
      f[xx]=max(max(f[tt],f[xx]),mx[xx]-mm[xx]);
      if(!inq[xx]) q.push(xx),inq[xx]=true;
    }
  }
}
int main()
{
  scanf("%d%d",&n,&m);
  for(int i=1;i<=n;i++)
  {
     scanf("%d",&c[i]);
  }
  for(int i=1;i<=m;i++)
  {
    int x,y,z;
    scanf("%d%d%d",&x,&y,&z);
    if(z==1) G[x].push_back(y);
    else
    {
      G[x].push_back(y);
      G[y].push_back(x);
    }
  }
  for(int i=1;i<=n;i++)
  {
    if(!dfn[i]) tarjan(i);
  }
  for(int i=1;i<=n;i++)
  {
      int sizz=G[i].size();
      for(int j=0;j<sizz;j++)
      {
          int v=G[i][j];
          if(co[i]==co[v]) continue;
              E[co[i]].push_back(co[v]);
      }
  }
  //cout<<qltn<<endl;
  bfs();
  printf("%d",f[co[n]]);
  return 0;
}

这种把一个带环的图转换成DAG并在上面跑DP的方法真是妙啊！