p3627&bzoj1179 抢掠计划(ATM)
题目
Input
Output
输出一个整数,表示Banditji从市中心开始到某个酒吧结束所能抢劫的最多的现金总数。
Sample Input
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
分析
Tarjan缩点的好题,先用Tarjan缩点,求出每个环的总钱数和这个环中是否有酒吧,然后重新建图跑spfa求最大点权和,但由于要在酒吧结束,所以最终答案只能在有酒吧的环中取。
插曲
作为一个最短路只会写folyd和dijstar的蒟蒻,我自己yy了一个spfa,然而它居然对了(QAQ)
代码
#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<cctype>
#include<cmath>
#include<cstdlib>
#include<queue>
#include<ctime>
#include<vector>
#include<set>
#include<map>
#include<stack>
using namespace std;
int bar[510000],mey[510000],belong[510000],cnt,tmy[510000],tbr[510000];
int dfn[510000],low[510000],vis[510000],sum;
int d[510000],ist[510000],maxn,iqu[510000];
vector<int>v[500010];
vector<int>nv[500010];
stack<int>a;
queue<pair<int,int> >q;
int n,m,s,p;
inline void read(int &x){
int f=1;x=0;
char s=getchar();
while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
while(s>='0'&&s<='9'){x=x*10+(s-'0');s=getchar();}
x*=f;
}
inline void tarjan(int x){
low[x]=dfn[x]=++cnt;
a.push(x);
ist[x]=1;
int i,j,k;
for(i=0;i<v[x].size();i++)
if(!dfn[v[x][i]]){
tarjan(v[x][i]);
low[x]=min(low[x],low[v[x][i]]);
}else if(ist[v[x][i]]){
low[x]=min(low[x],dfn[v[x][i]]);
}
if(dfn[x]==low[x]){
sum++;
int ok=0;
while(1){
int u=a.top();
ist[u]=0;
if(bar[u])ok=1;
belong[u]=sum;
tmy[sum]+=mey[u];
a.pop();
if(u==x)break;
}
if(ok)tbr[sum]=1;
}
return;
}
int main()
{ int i,j,k,x,y;
read(n),read(m);
for(i=1;i<=m;i++){
read(x),read(y);
v[x].push_back(y);
}
for(i=1;i<=n;i++){
read(mey[i]);
}
scanf("%d%d",&s,&p);
for(i=1;i<=p;i++){
read(x);
bar[x]=1;
}
for(i=1;i<=n;i++)
if(!dfn[i])
tarjan(i);
for(i=1;i<=n;i++)
for(j=0;j<v[i].size();j++)
if(belong[i]!=belong[v[i][j]]){
nv[belong[i]].push_back(belong[v[i][j]]);
}
d[belong[s]]=tmy[belong[s]];
q.push(make_pair(belong[s],tmy[belong[s]]));
iqu[belong[s]]=1;
while(!q.empty()){
int u=q.front().first;
q.pop();
iqu[u]=0;
if(tbr[u])maxn=max(maxn,d[u]);
for(i=0;i<nv[u].size();i++){
d[nv[u][i]]=max(d[nv[u][i]],d[u]+tmy[nv[u][i]]);
if(!iqu[nv[u][i]]){
q.push(make_pair(nv[u][i],d[nv[u][i]]));
iqu[nv[u][i]]=1;
}
}
}
cout<<maxn<<endl;
return 0;
}