bzoj1179: [Apio2009]Atm tarjan+spfa

最新推荐文章于 2019-03-12 00:16:02 发布

EMber _

最新推荐文章于 2019-03-12 00:16:02 发布

阅读量512

点赞数

分类专栏： bzoj tarjan 最短路

本文链接：https://blog.csdn.net/qq_35866453/article/details/78343813

版权

bzoj 同时被 3 个专栏收录

452 篇文章 0 订阅

订阅专栏

最短路

29 篇文章 0 订阅

订阅专栏

tarjan

7 篇文章 0 订阅

订阅专栏

博客介绍了如何解决Apio2009比赛中关于有向图的Atm问题，通过tarjan和spfa算法求解从源点s出发到特定点的最长路径。博主分享了初始尝试中遇到的错误，即误将最短路问题当作最长路问题，导致样例未通过。在简化题意后，问题变得清晰，博主指出该问题其实是一个基础的最长路径问题。

摘要由CSDN通过智能技术生成

题意：一个有向图求从s出发，到某一些规定点的最长路径。
一开始不化模型然后直接无脑dfs，这就是不把最短路当最短路，煞笔一个我，样例都过不去才发现问题。
然后简化了一下题意。。这不就裸的么。没什么好说的。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=5e5+5;
int n,m,a[N];
int head[N],next[N<<1],sta[N],top,go[N<<1],tot,s;
queue<int> q;
bool bz[N],vis[N];

int ans,dis[N],val[N];
struct node
{
    int x,y;
}e[N];
int sz=0,scc=0,bl[N],dfn[N],low[N];
inline void add(int x,int y)
{
    go[++tot]=y;
    next[tot]=head[x];
    head[x]=tot;
}

inline void tarjan(int x)
{
    dfn[x]=low[x]=++sz;
    bz[x]=1,sta[++top]=x;
    for(int i=head[x];i;i=next[i])
    {
        int v=go[i];
        if (!dfn[v])
        {
            tarjan(v);
            low[x]=min(low[x],low[v]);
        }
        else if (bz[v])low[x]=min(low[x],dfn[v]);
    }
    if (dfn[x]==low[x])
    {
        scc++;
        int y=0;
        while (y!=x)
        {
            y=sta[top--];
            bl[y]=scc;
            val[scc]+=a[y];
            bz[y]=0;
        }
    }
}
inline void spfa(int s)
{
    s=bl[s];
    dis[s]=val[s];
    q.push(s),vis[s]=1;
    while (!q.empty())
    {
        int x=q.front();
        q.pop();
        for(int i=head[x];i;i=next[i])
        {
            int v=go[i];
            if (dis[x]+val[v]>dis[v])
            {
                dis[v]=dis[x]+val[v];
                if (!vis[v])
                {
                    vis[v]=1;
                    q.push(v);
                }
            }

        }vis[x]=0;
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    fo(i,1,m)
    {
        scanf("%d%d",&e[i].x,&e[i].y);
        add(e[i].x,e[i].y);
    }
    fo(i,1,n)scanf("%d",&a[i]);
    int k;
    scanf("%d%d",&s,&k);
    fo(i,1,n)
    if (!dfn[i])tarjan(i);
    tot=0;memset(head,0,sizeof(head));
    fo(i,1,m)
    {
        int x=e[i].x,y=e[i].y;
        if (bl[x]!=bl[y])add(bl[x],bl[y]);
    }
    spfa(s);
    fo(i,1,k)
    {
        int x;
        scanf("%d",&x);
        ans=max(ans,dis[bl[x]]);
    }
    printf("%d\n",ans);
}