poj 3160

Father Christmas flymouse
题意：飞鼠想给大家发礼物，收到礼物后大家给飞鼠的满足程度都已数量化（有正有负）。每个人都住一个房间，房间之间的路是单向的，飞鼠可以选任一个点为起点，路可以重复走，房间不能重复访问，但可以经过而不访问。现在问飞鼠得到的满足程度最大会是多少？

分析：负分的房间肯定不访问。强连通分支可以缩点，缩成的点的分数取所有正分房间的总和，缩点后的图就变成若干个单向连通分支，由出度为零的点倒推即可算得最大值。

#include<cstdio>
#include<vector>
#include<stack>
#include<algorithm>
using namespace std;
int low[30000],num[30000],comfort[30000],total[30000],ans[30000],counter,set[30000],sn,n,m,in[30000],max_ans;
vector<int> group[30000],arc[30000],reverse_arc[30000];
bool visited[30000],instack[30000];
stack<int> connect,start;
void tarjan(int u)
{
    int i,v;
    connect.push(u);
    instack[u]=true;
    visited[u]=true;
    low[u]=num[u]=++counter;
    for(i=0;i<arc[u].size();i++)
    {
        v=arc[u][i];
        if(!visited[v])
        {
            tarjan(v);
            if(low[u]>low[v])
                low[u]=low[v];
        }
        else if(instack[v] && low[u]>num[v])
            low[u]=num[v];
    }
    if(num[u]==low[u])
    {
        total[sn]=0;
        do
        {
            v=connect.top();
            connect.pop();
            set[v]=sn;
            group[sn].push_back(v);
            instack[v]=false;
            if(comfort[v]>0)
                total[sn]+=comfort[v];
        }while(u!=v);
        ans[sn]=total[sn];
        sn++;
    }
}
void backtrace(int u)//on condiction of indegree[u]==0
{
    int i,v;
    for(i=0;i<reverse_arc[u].size();i++)
    {
        v=reverse_arc[u][i];
        if(ans[v]<ans[u]+total[v])
            ans[v]=ans[u]+total[v];
        in[v]--;
        if(in[v]==0)
            start.push(v);
    }
    if(max_ans<ans[u])
        max_ans=ans[u];
}
int main()
{
    int i,u,v;
    while(~scanf("%d%d",&n,&m))
    {
        for(i=0;i<n;i++)
        {
            visited[i]=false;
            instack[i]=false;
            scanf("%d",&comfort[i]);
        }
        for(i=0;i<m;i++)
        {
            scanf("%d%d",&u,&v);
            arc[u].push_back(v);
        }
        counter=sn=0;
        for(i=0;i<n;i++)
        {
            if(!visited[i])
                tarjan(i);
        }
        //reverse_arc
        for(i=0;i<sn;i++)
            in[i]=0;
        for(u=0;u<n;u++)
        {
            for(i=0;i<arc[u].size();i++)
            {
                v=arc[u][i];
                if(set[v]!=set[u] && find(reverse_arc[set[v]].begin(),reverse_arc[set[v]].end(),set[u])==reverse_arc[set[v]].end())
                {
                    reverse_arc[set[v]].push_back(set[u]);
                    in[set[u]]++;
                }
            }
        }
        for(i=0;i<sn;i++)
        {
            if(in[i]==0)
                start.push(i);
        }
        max_ans=0;
        while(!start.empty())
        {
            u=start.top();
            start.pop();
            backtrace(u);
        }
        for(i=0;i<sn;i++)
        {
            reverse_arc[i].clear();
            group[i].clear();
        }
        for(i=0;i<n;i++)
            arc[i].clear();
        printf("%d\n",max_ans);
    }
    return 0;
}