poj 3160 Father Christmas flymouse(强连通+dp)

原文：http://blog.acmj1991.com/?p=697

题意：给你一个有向图。有向图中每个点都有权值（可能为负）而且能走回头路但权值只能加一次，求权值最大的一条路。

思路：由于能走回头路所以必须先用强连通去环缩点，重新构图。然后必须到达一点时使其权值之和最大==>d[v]=max(dp[v],dp[u]+w[v])(v为子节点，u为父节点)。这里必须把所有的父节点都走完，所以用topsort

这题其实开始没想做，但是xs让做…..于是…..主要是很久没写强连通的题了，什么tarjan一直认为比较麻烦，但是现在谢谢也很容易，而且完全是自己想自己写的，虽然比以前刷题慢了一些但是自己想自己写的效果就是不一样

#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;
 
#define maxN 30010
int qhead[maxN],head[maxN],deg[maxN],wi[maxN],cost[maxN];
int dfn[maxN],low[maxN],vis[maxN],id[maxN],stack[maxN];
int n,res,index,top,num,maxx=0;
struct node{
    int v,next;
}q[6*maxN],p[6*maxN];
 
void init()
{
    memset(head,-1,sizeof(head));
    memset(qhead,-1,sizeof(qhead));
    memset(deg,0,sizeof(deg));
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(vis,0,sizeof(vis));
    memset(cost,0,sizeof(cost));
    maxx=top=res=0;
    num=index=1;
}
int min(int x,int y){return x<y?x:y;}
int max(int x,int y){return x>y?x:y;}
void add(int u,int v){
    p[num].v=v;
    p[num].next=head[u],head[u]=num++;
}
void dfs(int s)
{
    dfn[s]=low[s]=index++;
    stack[++top]=s;
    vis[s]=1;
    for(int i=qhead[s];i!=-1;i=q[i].next)
    {
        int v=q[i].v;
        if(!dfn[v]){
            dfs(v);
            low[s]=min(low[s],low[v]);
        }else if(vis[v])
            low[s]=min(low[s],low[v]);
    }
    if(dfn[s]==low[s])
    {
        int m;
        res++;
        do{
            m=stack[top--];
            vis[m]=0;
            id[m]=res;
            cost[res]+=wi[m];
        }while(m!=s);
    }
}
void tarjan()
{
    for(int i=0;i<n;i++)
        if(!dfn[i])
            dfs(i);
}
void build()
{
    for(int u=0;u<n;u++)
        for(int i=qhead[u];i!=-1;i=q[i].next){
            int v=q[i].v;
            if(id[u]!=id[v]){
                add(id[u],id[v]);
                deg[id[v]]++;
            }
        }
    top=0;
    memset(wi,0,sizeof(wi));
    for(int i=1;i<=res;i++)
        if(deg[i]==0){
            stack[++top]=i;
            wi[i]=cost[i];
            maxx=max(maxx,cost[i]);
        }
}
void topsort()
{
    while(top>0){
        int u=stack[top--];
        for(int i=head[u];i!=-1;i=p[i].next)
        {
            int v=p[i].v;
            if(--deg[v]==0)stack[++top]=v;
            wi[v]=max(wi[v],wi[u]+cost[v]);
            maxx=max(maxx,wi[v]);
        }
    }
}
 
int main()
{
    int m,u,v;
    while(~scanf("%d%d",&n,&m))
    {
        init();
        for(int i=0;i<n;i++){
            scanf("%d",&wi[i]);
            if(wi[i]<0)wi[i]=0;
        }
        for(int i=0;i<m;i++){
            scanf("%d%d",&u,&v);
            q[i].v=v;
            q[i].next=qhead[u],qhead[u]=i;
        }
        tarjan();
        build();
        topsort();
        printf("%d\n",maxx);
    }
}