洛谷 题解 P3627 【[APIO2009]抢掠计划】

图论 tarjan缩点+最短路 的一道题

• tarjan求强连通份量（为之后缩点打下良好的基础）

(若是不会tarjan的请点击这儿)

你须要的东西：
(1)、dfn[]，表示这个点在dfs时是第几个被搜到的。
(2)、low[]，表示这个点以及其子孙节点连的全部点中dfn最小的值
(3)、stack[]，表示当前全部可能能构成是强连通份量的点。
(4)、vis[]，表示一个点是否在stack中。
(5)、color[],记录每个点强连通份量的编号。
(6)、deep,记录dfs树的深度
inline void tarjan(int now)
{
    dfn[now]=++deep;
    low[now]=deep;
    vis[now]=1;
    st.push(now);
    for(int i=0;i<ver[now].size();i++)
    {
        int x=ver[now][i];
        if(!dfn[x])
        {
            tarjan(x);
            low[now]=min(low[now],low[x]);
        }
        else
        {
            if(vis[x])
                low[now]=min(low[now],low[x]);
        }
    }
    if(dfn[now]==low[now])
    {
        color[now]=++sum;
        vis[now]=0;
        while(st.top()!=now)
        {
            color[st.top()]=sum;
            vis[st.top()]=0;
            st.pop();
        }
        st.pop();
    }
}

• 缩点（去除图中的环）

//重点：这里建新图是依托强连通份量的编号来建的
for(int i=1;i<=n;i++)
    {
        for(int j=0;j<ver[i].size();j++)
        {
            int x=ver[i][j];
            if(color[i]!=color[x])//若是不是属于同一个强连通份量中，那么就合并
            {
                g[color[i]].push_back(color[x]);
           //千万不能写成g[i].push_back(x);坑死我了
            }
        }
    }
//这一部分代码还能够适当优化...（想想）

• 对于点权与酒馆的一些处理（方便求最短路）

for(int i=1;i<=n;i++)
    {
        ww[color[i]]+=w[i];//将这个强连通份量中全部的点权所有加起来
        if(tf[i])tf[color[i]]=1;//只要这个强连通份量中有一个结点有酒馆，那么就设定为有酒馆
    }

• 求最短路模板（然而实际是最长路）

//模板不作解释
inline void spfa()
{
    d[color[s]]=ww[color[s]];
    queue<int>q;
    q.push(color[s]);
    while(q.size())
    {
        int now=q.front();
        q.pop();
        for(int i=0;i<g[now].size();i++)
        {
            int x=g[now][i];
            if(d[now]+ww[x]>d[x])
            {
                d[x]=d[now]+ww[x];
                q.push(x);
            }
        }
    }
}

因此...

在全部有酒馆的节点中选一个最大值输出就行了

for(int i=1;i<=sum;i++)
    {
        //cout<<d[i]<<" ";
        if(tf[i])ans=max(ans,d[i]);
    }

完整代码

#include<bits/stdc++.h>
using namespace std;
const int MAXN=500000+10;
int n,m;
vector<int>ver[MAXN];
vector<int>g[MAXN];
int w[MAXN],ww[MAXN];
bool tf[MAXN];
int s,p,ans=0;
int dfn[MAXN],color[MAXN],low[MAXN];
int deep,sum;
bool vis[MAXN];
int d[MAXN];
stack<int>st;
inline int read()
{
    int tot=0;
    char c=getchar();
    while(c<'0'||c>'9')
        c=getchar();
    while(c>='0'&&c<='9')
    {
        tot=tot*10+c-'0';
        c=getchar();
    }
    return tot;
}
inline void tarjan(int now)
{
    dfn[now]=++deep;
    low[now]=deep;
    vis[now]=1;
    st.push(now);
    for(int i=0;i<ver[now].size();i++)
    {
        int x=ver[now][i];
        if(!dfn[x])
        {
            tarjan(x);
            low[now]=min(low[now],low[x]);
        }
        else
        {
            if(vis[x])
                low[now]=min(low[now],low[x]);
        }
    }
    if(dfn[now]==low[now])
    {
        color[now]=++sum;
        vis[now]=0;
        while(st.top()!=now)
        {
            color[st.top()]=sum;
            vis[st.top()]=0;
            st.pop();
        }
        st.pop();
    }
}
inline void spfa()
{
    d[color[s]]=ww[color[s]];
    queue<int>q;
    q.push(color[s]);
    while(q.size())
    {
        int now=q.front();
        q.pop();
        for(int i=0;i<g[now].size();i++)
        {
            int x=g[now][i];
            if(d[now]+ww[x]>d[x])
            {
                d[x]=d[now]+ww[x];
                q.push(x);
            }
        }
    }
}
int main()
{
    //freopen("testdata.in","r",stdin);
    n=read();m=read();
    for(int i=1;i<=m;i++)
    {
        int x=read(),y=read();
        ver[x].push_back(y);
    }
    for(int i=1;i<=n;i++)
        w[i]=read();
    s=read();p=read();
    for(int i=1;i<=p;i++)
    {
        int x=read();
        tf[x]=1;
    }
    for(int i=1;i<=n;i++)
    {
        if(!dfn[i])
            tarjan(i);
    }
    /*cout<<endl;
    for(int i=1;i<=n;i++)
        cout<<color[i]<<" ";
    cout<<endl;*/
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<ver[i].size();j++)
        {
            int x=ver[i][j];
            if(color[i]!=color[x])
            {
                //cout<<i<<" "<<x<<" "<<color[i]<<" "<<color[x]<<endl;
                g[color[i]].push_back(color[x]);
            }
        }
    }
    for(int i=1;i<=n;i++)
    {
        ww[color[i]]+=w[i];
        if(tf[i])tf[color[i]]=1;
    }
    /*cout<<color[s]<<endl;
    cout<<endl;
    for(int i=1;i<=sum;i++)
    {
        for(int j=0;j<g[i].size();j++)cout<<i<<" "<<g[i][j]<<"\n";
    }
    cout<<endl;
    for(int i=1;i<=sum;i++)
        cout<<ww[i]<<" ";cout<<endl;
    for(int i=1;i<=sum;i++)
        cout<<tf[i]<<" ";cout<<endl;
    cout<<ww[color[s]]<<endl;*/
    spfa();
    for(int i=1;i<=sum;i++)
    {
        //cout<<d[i]<<" ";
        if(tf[i])ans=max(ans,d[i]);
    }
    //cout<<endl;
    cout<<ans<<endl;
    return 0;
}