[NOIP2009 提高组] 最优贸易（分层图）

最新推荐文章于 2025-08-12 07:00:00 发布

原创最新推荐文章于 2025-08-12 07:00:00 发布 · 324 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

本文介绍了一种使用分层图解决洛谷P1073水晶球问题的方法。通过构建三张平行图来模拟交易状态的变化，并运用SPFA算法求解最大收益。

题目链接：https://www.luogu.com.cn/problem/P1073

题目描述

简要题解

分层图应用题
这道题目实际上要求我们求出从1->n的旅途中(每个点可以多次询问)至多进行一次水晶球的买卖赚取差价。以此求最多能赚的钱，那么我们可以将整个旅途划分为以下三个状态。

当前没有进行过任何买入和卖出操作
当前已经进行过买入操作（即手里有水晶球，但未卖出）
当前已经将手里的水晶球卖出。

那么我们可以按照这三个状态，构造出三张平行图（分层图），首先由题目中的边建出原图，并将原图“复制”
出三份,接下去我们分析一下，三张图(也就是三个状态)的转化关系。

对于第一张图：

一直不进行操作,继续走下去
在当前点买入一颗水晶球,此时转化到第二个状态,跳转到第二张图

对于第二张图：

一直不进行操作,继续走下去
在当前点卖出手中的水晶球，此时转化到第三个状态，跳转到第三张图

对于第三张图:

一直不进行操作,一直到结束旅行

分层图设计：

首先明确分层图的设计，我们假设第一张图的点为x,那么x在第二张图中的点为x+n,在第三张图为x+2*n…
我们可以受到启发，当有一条a->b的边时,我们进行如下的加边操作。

add(a,b,0)//表示从a->b,即从第一张图走到第一张图，不进行操作，代价为0
add(a+n,b+n,0)//表示从a+n->b+n,即从第二张图走到第二张图，不进行操作，代价为0
add(a+2n,b+2n,0)//表示从a+2n->b+2n,即从第三张图走到第三张图，不进行操作，代价为0
add(a,b+n,-w[a])//表示从a->b+n,即在a点进行了买入水晶球操作，进入第二张图,代价为-w[a]
add(a+n,b+2n,w[a])//表示从a+n->b+2n,即在a点进行了卖出水晶球操作,进入了第三张图，代价为w[a]

注意的是，我们可以设置一个0号点表示旅途结束，即我们需要add(n,0,0),add(2n,0,w[n]),add(3n,0,0);
我们在三张图的n号点,我们都可以选择结束旅途。第一张图结束旅途代价为0,第二张图结束旅途代价为wn，第三张图结束旅途代价也为0。

建立好分层图后，只需要从1跑一边最短路算法找出最后贡献最大的一条路(有负权,可以跑spfa)
最后答案即为dist[0]

题解代码

#include <bits/stdc++.h>
using namespace std;
typedef pair<int,int>PII;
const int N=500050*6;
int h[300010],e[N],ne[N],sz[100010],w[N];
int n,m,idx,dist[300010],st[300030];
void add(int a,int b,int c)
{
    e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}
void fun(int a,int b)
{
    add(a,b,0);
    add(a,b+n,-sz[a]);
    add(a+n,b+n,0);
    add(a+2*n,b+2*n,0);
    add(a+n,b+2*n,sz[a]);
}
void spfa()
{
    memset(dist,-0x3f,sizeof dist);
    dist[1]=0;
    st[1]=1;
    queue<int>q;
    q.push(1);
    while(q.size())
    {
        int u=q.front();
        q.pop();
        st[u]=0;
        for(int i=h[u];i!=-1;i=ne[i])
        {
            int j=e[i];
            if(dist[j]<dist[u]+w[i])
            {
                dist[j]=dist[u]+w[i];
                if(!st[j])
                {
                    st[j]=1;
                    q.push(j);
                }
            }
        }
    }
}
int main(void)
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    memset(h,-1,sizeof h);
    cin>>n>>m;
    for(int i=1;i<=n;i++)cin>>sz[i];
    for(int i=1;i<=m;i++)
    {
        int a,b,c;
        cin>>a>>b>>c;
        if(c==1)fun(a,b);
        else
        {
            fun(a,b);
            fun(b,a);
        }
    }
    add(n,0,0);
    add(n*2,0,sz[n]);
    add(n*3,0,0);
    spfa();
    cout<<dist[0];
    //system("pause");
    return 0;
}