bzoj1179（Atm）

---恢复内容开始---

1179: [Apio2009]Atm

Time Limit: 15 Sec  Memory Limit: 162 MB

Description

Input

第一行包含两个整数N、M。N表示路口的个数，M表示道路条数。接下来M行，每行两个整数，这两个整数都在1到N之间，第i+1行的两个整数表示第i条道路的起点和终点的路口编号。接下来N行，每行一个整数，按顺序表示每个路口处的ATM机中的钱数。接下来一行包含两个整数S、P，S表示市中心的编号，也就是出发的路口。P表示酒吧数目。接下来的一行中有P个整数，表示P个有酒吧的路口的编号

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

HINT

50%的输入保证N, M<=3000。所有的输入保证N, M<=500000。每个ATM机中可取的钱数为一个非负整数且不超过4000。输入数据保证你可以从市中心沿着Siruseri的单向的道路到达其中的至少一个酒吧。

这道题写了好久，终于写出来啦（哭。。。。。）

这道题思路大概是这样的，用adjacency list存完边后，先做一个tarjan并染色的预处理，然后把所有相同颜色的城市看成一个强连通分量，把这个强连通分量的颜色设置为新的地图编号为color的一个城市，然后将记录下来的原先的城市的道路，一一判断是否在新的map中是一座城市，若不是，那重新建立道路。（当然也得重新设置酒吧的位置）

最后跑一遍spfa求出每个酒吧的最大值即可（因为把强连通分量看成了一个城市了，所以不会再产生环）。

还是一样，贴代码（其实后来发现有一点可以省略不写的，就是重新设置好城市以后，把酒吧原来所在城市的的编号改成现在的就行，，还多开了一个array）

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 ;
 using namespace std;
 struct node {int next,aim;};
 node f[maxn],e[maxn];
 int point,head[maxn],Low[maxn],stack[maxn],DFN[maxn],dep,top,cl[maxn],color,queue[maxn],n,m,newmoney[maxn];
 bool vis[maxn],bar[maxn],newbar[maxn];
 int x1[maxn],y1[maxn],money[maxn],start,t;
 long long ans,d[maxn];
 void add(node *array,int x,int y)
 {
   point++;
   array[point].next=head[x];
   head[x]=point;
   array[point].aim=y;
 }
 void tarjan(int x)
 {
     DFN[x]=Low[x]=++dep;//为节点u设定次序编号和Low初值
      stack[++top]=x;
      vis[x]=true; //将节点u压入栈中
     ;i=e[i].next)//枚举每一条边
         if (!DFN[e[i].aim])
         {//如果节点v未被访问过
             tarjan(e[i].aim);//继续向下找
             Low[x]=min(Low[x],Low[e[i].aim]);
         }
         else if (vis[e[i].aim])//如果节点v还在栈内
                 Low[x]=min(Low[x],DFN[e[i].aim]);
     if (DFN[x]==Low[x])//如果节点u是强连通分量的根
     {
      vis[x]=)cl[x]=++color;
      while(stack[top]!=x)
      {
       cl[stack[top]]=color;
       vis[stack[top]]=false;
       top--;
      }
      top--;
     }

 }
 void spfa(int x)
 {
  ,tail=;
  queue[Head]=x;
  while(Head<tail)
  {
   ;i=f[i].next)
   {
    int u=f[i].aim;
    if(d[queue[Head]]+newmoney[u]>d[u])
    {
     d[u]=d[queue[Head]]+newmoney[u];
     if(!vis[u])queue[tail++]=u;
     vis[u]=true;
    }
   }
   vis[queue[Head]]=false;
   Head++;
  }
 }
 int main()
 {
 scanf("%d%d",&n,&m);
 ;i<=m;i++)
 {
  scanf("%d%d",&x1[i],&y1[i]);
  add(e,x1[i],y1[i]);
 }
 ;i<=n;i++)
 {
  scanf("%d",&money[i]);
 }
 scanf("%d",&start);
 scanf("%d",&t);
 ;i<=t;i++)
 {
  int l;
  scanf("%d",&l);
  bar[l]=true;
 }
 ;i<=n;i++)
 if(!DFN[i])tarjan(i);
 memset(head,,sizeof(head));
 point=;
 ;i<=m;i++)
 {
  int p=cl[x1[i]],q=cl[y1[i]];
  if(p!=q)add(f,p,q);
 }
 ;i<=n;i++)
 newmoney[cl[i]]+=money[i];
 start=cl[start];
 ;i<=n;i++)
 {
 if(bar[i]){
  newbar[cl[i]]=true;
 }
 }
 memset(vis,,sizeof(vis));
 d[start]=newmoney[start];
 vis[start]=true;
  spfa(start);
  ;i<=color;i++)
  if(newbar[i])ans=max(ans,d[i]);
  printf("%lld",ans);
  ;
 } 

---恢复内容结束---