[bzoj]1179 ATM

1179: [Apio2009]Atm

Time Limit: 15 Sec  Memory Limit: 162 MB
Submit: 4645  Solved: 2054
[Submit][Status][Discuss]

Description

Siruseri 城中的道路都是单向的。不同的道路由路口连接。按照法律的规定， 在每个路口都设立了一个 Siruser

i 银行的 ATM 取款机。令人奇怪的是，Siruseri 的酒吧也都设在路口，虽然并不是每个路口都设有酒吧。Bandit

ji 计划实施 Siruseri 有史以来最惊天动地的 ATM 抢劫。他将从市中心 出发，沿着单向道路行驶，抢劫所有他

途径的 ATM 机，最终他将在一个酒吧庆 祝他的胜利。使用高超的黑客技术，他获知了每个 ATM 机中可以掠取的

现金数额。他希 望你帮助他计算从市中心出发最后到达某个酒吧时最多能抢劫的现金总数。他可 以经过同一路口

或道路任意多次。但只要他抢劫过某个 ATM 机后，该 ATM 机 里面就不会再有钱了。 例如，假设该城中有 6 个

路口，道路的连接情况如下图所示：

市中心在路口 1，由一个入口符号→来标识，那些有酒吧的路口用双圈来表示。每个 ATM 机中可取的钱数标在了

路口的上方。在这个例子中，Banditji 能抢 劫的现金总数为 47，实施的抢劫路线是：1-2-4-1-2-3-5。

Input

第一行包含两个整数N、M。N表示路口的个数，M表示道路条数。

接下来M行，每行两个整数，这两个整数都在1到N之间，

第i+1行的两个整数表示第i条道路的起点和终点的路口编号。

接下来N行，每行一个整数，按顺序表示每个路口处的ATM机中的钱数。

接下来一行包含两个整数S、P，S表示市中心的编号，也就是出发的路口。P表示酒吧数目。

接下来的一行中有P个整数，表示P个有酒吧的路口的编号

N, M<=500000。每个ATM机中可取的钱数为一个非负整数且不超过4000。

输入数据保证你可以从市中心沿着Siruseri的单向的道路到达其中的至少一个酒吧。

Output

输出一个整数，表示Banditji从市中心开始到某个酒吧结束所能抢劫的最多的现金总数。

Sample Input

6 7
1 2
2 3
3 5
2 4
4 1
2 6
6 5
10
12
8
16
1 5
1 4
4
3
5
6

Sample Output

47

Tarjan缩点之后，跑一遍spfa求最长路就行了。

中间因为变量，初始化问题死了好多次，我会在代码里写出我给自己挖的坑。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<vector>
#include<stack>
#include<map>
using namespace std;
#define mem(a,b) memset(a,b,sizeof(a))
#define ll long long
#define inf 1000000000
#define maxn 500000+100
#define maxm 500000+100
struct node
{
    int to,next;
} edge[maxm];
int n,m,head[maxn],vis[maxn],dfn[maxn],low[maxn],cdu[maxn],num[maxn],cnt,timi,stack1[maxn],top,cut,rdu[maxn];
int mon[maxn],jiu[maxn],ss,p,from[maxn],go[maxn],s1[maxn],d[maxn],q[maxn],c[maxn];
stack <int> s;
inline int read()
  {
      int x=0,f=1;char ch=getchar();
      while(ch<'0'||ch>'9') {if(ch=='-') f=-1;ch=getchar();}
      while(ch>='0'&&ch<='9') {x=10*x+ch-'0';ch=getchar();}
      return x*f;
  }
void init()
{
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(head,-1,sizeof(head));
    memset(vis,0,sizeof(vis));
    memset(num,0,sizeof(num));
    memset(cdu,0,sizeof(cdu));
    memset(rdu,0,sizeof(rdu));
    mem(jiu,0);
    mem(c,0);
    //while(!s.empty()) s.pop();
    cnt=0;
    timi=1;
    top=0;
    cut=0;
}
void addedge(int u,int v)
{
    edge[cnt].to=v;
    edge[cnt].next=head[u];
    head[u]=cnt;
    cnt++;
}
void tarjan(int u)
{
    dfn[u]=timi;
    low[u]=timi;
    timi++;
    s.push(u);
    vis[u]=1;
    for(int i=head[u]; i!=-1; i=edge[i].next)
    {
        int v=edge[i].to;
        if(!dfn[v])
        {
            tarjan(v);
            low[u]=min(low[u],low[v]);
        }
        else
        {
            if(vis[v])
                low[u]=min(low[u],dfn[v]);
        }
    }
    if(low[u]==dfn[u])
    {
        cut++;
        int x=s.top();
        while(x!=u)
        {
            vis[x]=0;
            num[x]=cut;
            s1[cut]+=mon[x];
            //printf("cut=%d s1=%d \n",cut,s1[cut]);
            if(jiu[x]) c[cut]=1;
            s.pop();
            x=s.top();
        }
        num[x]=cut;
        s1[cut]+=mon[x];//我今天又在这里死了，你不能因为人家在外面你就不管人家啊。。瞎浪
        if(jiu[x]) c[cut]=1;
        //printf("cut=%d s1=%d \n",cut,s1[cut]);
        vis[x]=0;
        s.pop();
    }
}
void spfa()
{
    for(int i=1;i<=cut;++i) d[i]=0;
    mem(vis,0);
    int l=0,r=1,x,y;
    
    q[1]=num[ss];d[num[ss]]=s1[num[ss]];
    //cout<<q[1]<<endl;
    //cout<<l<<" "<<r<<endl;
    while(l!=r)
    {
        x=q[++l];
        //cout<<head[x]<<endl;
        if(l==maxn) l=0;vis[x]=0;
        for(int i=head[x];i!=-1;i=edge[i].next)
        {
         // cout<<"i="<<i<<endl;
            if(d[x]+s1[y=edge[i].to]>d[y])
            {
                d[y]=d[x]+s1[y];
               // cout<<"y="<<y<<endl;
                if(!vis[y])
                {
                    vis[y]=1;q[++r]=y;if(r==maxn) r=0;
                }
            }
        }
    }
}
int main()
{
   init();//第一次跑不出来死在这里，没有初始化，今天在这个地方死第二次了（微笑
   int x,y;
   n=read();m=read();
   for(int i=1;i<=m;++i)
   {
       int a,b;
       a=read();b=read();
       from[i]=a;go[i]=b;
       addedge(a,b);
   }
   for(int i=1;i<=n;++i) mon[i]=read();
   ss=read();p=read();
   for(int i=1;i<=p;++i)
   {
       int h=read();jiu[h]=1;
   }
   for(int i=1;i<=n;++i)
    if(!dfn[i]) tarjan(i);
   mem(head,-1);
   cnt=0;
   for(int i=1;i<=m;++i)
   {
       if(num[x=from[i]]!=num[y=go[i]]) addedge(num[x],num[y]);//最后发现死在了这里，我真的是zz，缩点后不用缩点。。
   }
   spfa();
   int ans=0;
   //for(int i=1;i<=n;++i) cout<<num[i]<<" "<<mon[i]<<endl;
   //cout<<cut<<endl;
   //for(int i=1;i<=cut;++i) cout<<s1[i]<<endl;
   for(int i=1;i<=cut;++i) if(c[i]) ans=max(ans,d[i]);
   cout<<ans<<endl;
}