bzoj 1179 Atm

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1179

题解：

　　一道比较综合的图论题

　　直接讲正解：
　　如果这个图G中存在某个强连通分量，那么这个强连通分量中的所有ATM即可视为都被抢到，所有的酒吧都可视为重点，并且也可以从这个强连通分量的任何结点出发继续向外扩展

　　所以先做一遍Tarjan，找出强连通分量，然后重新构图，把每个强连通分量缩成一个点，此点的权值即为原先强连通分量里所有点权之和，判断此点中有没有酒吧，再将原先所有连接强连通分量的边连接在这个点

　　最后做SPFA，找出权值最大的路径

　　这道题时间限制15s，空间162MB，所以一般不会TLE或MLE

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAXN 500010
#define MAXM 500010
int n,m,time,cnt,head1[MAXN],DFN[MAXN],Low[MAXN],stack[MAXN],top,a1[MAXN],belong[MAXN];
bool vis[MAXN],bar1[MAXN];
struct node1
{
    int v,next;
}edge1[MAXM];
//以上变量基本是在进行tarjan（及之前）使用 
int heads[MAXN],mm,a2[MAXN],s,p,d[MAXN],q[MAXN],head,tail,ans;
bool viss[MAXN];
struct node2
{
    int v,next;
}edge2[MAXM];
//以上变量是在进行SPFA时使用 
//a存放点权，bar记录是否有酒吧 
bool bar2[MAXN];
void add(int x,int y)
{
    edge1[++cnt].next=head1[x];
    head1[x]=cnt;
    edge1[cnt].v=y;
}
void buildmap(int k)//重新建图、连边 
{
    for(int i=1;i<=n;i++)
    {
        if(belong[i]==k)
        {
            for(int j=head1[i];j;j=edge1[j].next)
            {
                if(belong[edge1[j].v]==k)continue;
                edge2[++m].next=heads[k];
                heads[k]=m;
                edge2[m].v=belong[edge1[j].v];
            }
        }
    }
}
void tarjan(int u)
{
    DFN[u]=Low[u]=++time;
    vis[u]=true;
    stack[++top]=u;
    for(int i=head1[u];i;i=edge1[i].next)
    {
        int v=edge1[i].v;    
        if(!DFN[v]) 
        {
            tarjan(v);
            Low[u]=min(Low[u],Low[v]);
        }
        else if(vis[v])Low[u]=min(Low[u],DFN[v]);
    }
    if(DFN[u]==Low[u]) 
    {
        int i;
        cnt++;
        do
        {
            i=stack[top--];
            vis[i]=false;
            belong[i]=cnt;
            a2[cnt]+=a1[i];
            if(bar1[i])bar2[cnt]=true;//先更新点权和是否有酒吧 
        }while(u!=i);
    }
}
void SPFA()
{
    head=1;tail=2;
    q[1]=s;
    viss[s]=1;
    while(head<tail)
    {
        for(int i=heads[q[head]];i!=0;i=edge2[i].next)
        {
            if(d[q[head]]+a2[edge2[i].v]>d[edge2[i].v])//注意，这里不是边权，是点权 
            {
                d[edge2[i].v]=d[q[head]]+a2[edge2[i].v];
                if(!viss[edge2[i].v])
                {
                    q[tail++]=edge2[i].v;
                    viss[edge2[i].v]=true;
                }
            }
        }
        viss[q[head]]=false;
        head++;
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        add(x,y);
    }
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a1[i]);
    }
    scanf("%d%d",&s,&p);
    for(int i=1;i<=p;i++)
    {
        int x;
        scanf("%d",&x);
        bar1[x]=true;
    }
    cnt=0;
    m=0;
    for(int i=1;i<=n;i++)
        if(!DFN[i])tarjan(i);
    for(int i=1;i<=cnt;i++)
        buildmap(i);
    s=belong[s];
    for(int i=1;i<=cnt;i++)
    {
        d[i]=a2[i];
    }
    SPFA();
    for(int i=1;i<=cnt;i++)
    {
        if(bar2[i])ans=max(ans,d[i]);
    }
    printf("%d",ans);
    return 0;
}