bzoj 1179(强连通+spfa)
1179: [Apio2009]Atm
Time Limit: 15 Sec Memory Limit: 162 MBSubmit: 2358 Solved: 972
[ Submit][ Status][ Discuss]
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
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的单向的道路到达其中的至少一个酒吧。
解题思路：先强连通，再用spfa（可把边权改为负数再求最短路）
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cmath>
using namespace std;
const int maxn=500001;
int len,len1,S,tcc,head,size,n,m,tail;
int to[500001],to1[500001],next[500001],next1[500001];
int h[500001],h1[500001],dfn[500001],twf[500001],dui[500001],zhi[500001];
int q[500001],pan[500001],v[500001],dis[500001];
bool b[500001];
inline int read()
{
char y; int x=0,f=1; y=getchar();
while (y<'0' || y>'9') {if (y=='-') f=-1; y=getchar();}
while (y>='0' && y<='9') {x=x*10+int(y)-48; y=getchar();}
return x*f;
}
void insert(int x,int y)
{
++len; to[len]=y; next[len]=h[x]; h[x]=len;
}
void insert1(int x,int y)
{
++len1; to1[len1]=y; next1[len1]=h1[x]; h1[x]=len1;
}
void dfs(int now)
{
++tcc; ++head; q[head]=now; dfn[now]=twf[now]=tcc; pan[now]=1;
int u=h[now];
while (u!=0)
{
if (pan[to[u]]==0)
{
dfs(to[u]);
}
if (pan[to[u]]==1)
{
twf[now]=min(twf[now],twf[to[u]]);
}
u=next[u];
}
if (dfn[now]==twf[now])
{
++size;
while (q[head]!=now)
{
--head;
dui[q[head+1]]=size; zhi[size]+=v[q[head+1]];
pan[q[head+1]]=2;
}
--head; dui[now]=size; zhi[size]+=v[now]; pan[now]=2;
}
}
int main()
{
freopen("large-tree-0.in","r",stdin);
freopen("main.out","w",stdout);
n=read(); m=read();
for (int i=1;i<=m;++i)
{
int x=read(); int y=read();
insert(x,y);
}
for (int i=1;i<=n;++i)
{
v[i]=read();
}
S=read();
size=0; head=0; tcc=0;
memset(dfn,0,sizeof(dfn));
memset(twf,0,sizeof(twf));
memset(pan,0,sizeof(pan));
for (int i=1;i<=n;++i)
if (pan[i]==0)
{
dfs(i);
}
for (int i=1;i<=n;++i)
{
int u=h[i];
while (u!=0)
{
int u1=dui[i]; int u2=dui[to[u]];
if (u1!=u2)
{
insert1(u1,u2);
}
u=next[u];
}
}
memset(dis,0x7f,sizeof(dis));
memset(b,true,sizeof(b));
tail=0; head=0; ++tail; q[tail]=dui[S]; dis[dui[S]]=-zhi[dui[S]]; b[dui[S]]=false;
while(head<tail)
{
head=(head+1)%maxn;
int u=h1[q[head]];
while (u!=0)
{
if (dis[q[head]]-zhi[to1[u]]<dis[to1[u]])
{
dis[to1[u]]=dis[q[head]]-zhi[to1[u]];
if (b[to1[u]])
{
b[to1[u]]=false;
tail=(tail+1)%maxn; q[tail]=to1[u];
}
}
u=next1[u];
}
b[q[head]]=true;
}
int q=read();
int ans=-1;
for (int i=1;i<=q;++i)
{
int x=read();
if (dui[x]==dui[S]) ans=max(ans,zhi[S]);else
ans=max(ans,-dis[dui[x]]);
}
printf("%d",ans);
}