洛谷P3387

3 篇文章 0 订阅
3 篇文章 0 订阅

思路就是跑tarjan缩点,把同一个联通分量里的点全部缩成一个点,所有的权值都记在这个点上,然后根据tarjan的结果,在不同联通分量之间重新建图,联通分量里连出去的边现在全部归结到剩余的这个点上,连进来的边也同理。然后统计那些入度为0的点,它们肯定能跑出最好的结果,因为入度不为0的点肯定会被那些入度为0的点跑到。对于这些点每个点跑一次spfa,统计答案即可。实际操作中一定要小心…很容易出错…

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<vector>
#include<stack>
#include<queue>
#include<cstring>
using namespace std;
//tarjan缩点然后重新建图,将不再一个强联通分量中的点重新连边,每一个联通分量变为一个点,点权为其中点之和,然后直接spfa遍历那些入度为0的点(因为显然入度不为0的点出发不如入度为0的点出发)
int sccindex=0,belong[10005],scccnt=0,low[10005],dfn[10005],indegree[10005];//belong记录各个点属于哪个联通分量,indegree记录缩点后各个联通块的入度
bool instack[10005];
stack<int>sta;
vector<int>G1[10005];//原始图
vector<int>G2[10005];//缩点后的图
int dq[10005],n ,m;//原始点权
int newdq[10005];//缩点之后的点权
void tarjan(int v)
{
    low[v]=dfn[v]=++sccindex;
    sta.push(v);instack[v]=true;
    for(auto j=G1[v].begin();j!=G1[v].end();j++){
        if(!dfn[*j]){
            tarjan(*j);
            low[v]=min(low[v],low[*j]);
        }
        else if(instack[*j]&&dfn[*j]<low[v])
            low[v]=dfn[*j];
    }
    if(low[v]==dfn[v]){
        scccnt++;int j;
        do{
            j=sta.top();sta.pop();instack[j]=false;belong[j]=scccnt;newdq[scccnt]+=dq[j];//将联通分量中的点权加和
        }while(j!=v);
    }
}
int spfa(int x)
{
    int dis[10005];memset(dis,0,sizeof(dis));
    bool inque[10005];memset(inque,0,sizeof(inque));
    dis[x]=newdq[x];//将每条边的长度认为是终点的点权,所以起点距离加上起点的点权
    queue<int>que;que.push(x);int res=dis[x];
    while(!que.empty()){
        int t=que.front();que.pop();inque[t]=false;
        for(auto j=G2[t].begin();j!=G2[t].end();j++){
            if(dis[*j]<dis[t]+newdq[*j]){
                dis[*j]=dis[t]+newdq[*j];
                res=max(res,dis[*j]);
                if(!inque[*j]){
                    que.push(*j);inque[*j]=true;
                }
            }
        }
    }
    //for(int i=1;i<=scccnt;i++)res=max(res,dis[i]);
    return res;
}
int main()
{
    int i,j,k;
    cin>>n>>m;
    for(i=1;i<=n;i++)cin>>dq[i];
    for(i=1;i<=m;i++){
        scanf("%d%d",&j,&k);
        G1[j].push_back(k);
    }
    for(i=1;i<=n;i++)
        if(!dfn[i])
            tarjan(i);
    /*for(i=1;i<=scccnt;i++){
        for(j=1;j<=n;j++)
            if(belong[j]==i)cout<<j<<' ';
        cout<<endl;
    }*/
    for(i=1;i<=n;i++){
        for(auto j=G1[i].begin();j!=G1[i].end();j++){
            if(belong[i]!=belong[*j]){
                G2[belong[i]].push_back(belong[*j]);//如果两点不在一个联通分量中,那么在缩点后的图中添加一条边
                indegree[belong[*j]]++;//统计入度
            }
        }
    }
    //for(i=1;i<=scccnt;i++)cout<<indegree[i]<<' ';
    int ans=0;
    for(i=1;i<=scccnt;i++)
        if(!indegree[i])
            ans=max(ans,spfa(i));
    cout<<ans<<endl;

    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值