洛谷3387：【模板】缩点——题解

https://www.luogu.org/problemnew/show/P3387

有环我们很不好判断，但是显然在一个强连通分量里面的点我们都可以到达。

所以缩起来之后就是一个DAG，之后随便你怎么算都行（bfs/dfs）

（我觉得并不算是dp emmm）

#include<stack>
#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=1e4+5;
const int M=1e5+5;
inline int read(){
int x=0,w=1;char ch=0;
while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*w;
}
struct node{
int w,to,nxt;
}e[2][M];
int cnt,h[2]
;
inline void add(int u,int v,int k){
e[k][++cnt].to=v;e[k][cnt].nxt=h[k][u];h[k][u]=cnt;
}
int dfn
,low
,val
,to
,sum
,f
,indeg
,t,l;
stack<int>q;
queue<int>p;
bool inq
,vis
;
void tarjan(int u){
int v;
dfn[u]=low[u]=++t;
q.push(u);inq[u]=1;
for(int i=h[0][u];i;i=e[0][i].nxt){
v=e[0][i].to;
if(!dfn[v]){
tarjan(v);
low[u]=min(low[u],low[v]);
}else if(inq[v]){
low[u]=min(low[u],dfn[v]);
}
}
if(dfn[u]==low[u]){
l++;
do{
v=q.top();
q.pop();
inq[v]=0;
to[v]=l;
sum[l]+=val[v];
}while(v!=u);
}
}
void bfs(int u){
vis[u]=1;f[u]=sum[u];
p.push(u);
while(!p.empty()){
u=p.front();p.pop();vis[u]=0;
for(int i=h[1][u];i;i=e[1][i].nxt){
int v=e[1][i].to;
if(f[v]<f[u]+sum[v]){
f[v]=f[u]+sum[v];
if(!vis[v]){
vis[v]=1;
p.push(v);
}
}
}
}
}
int main(){
int n=read(),m=read();
for(int i=1;i<=n;i++)val[i]=read();
for(int i=1;i<=m;i++){
int u=read(),v=read();
add(u,v,0);
}
for(int i=1;i<=n;i++){
if(!dfn[i])tarjan(i);
}
for(int u=1;u<=n;u++){
for(int j=h[0][u];j;j=e[0][j].nxt){
int v=to[e[0][j].to];
bool link=0;
for(int k=h[1][to[u]];k&&!link;k=e[1][k].nxt){
int w=e[1][k].to;
if(v==w)link=1;
}
if(!link&&to[u]!=v){
add(to[u],v,1);
indeg[v]++;
}
}
}
memset(f,-127,sizeof(f));
for(int i=1;i<=l;i++){
if(!indeg[i]){
bfs(i);
}
}
int ans=0;
for(int i=1;i<=l;i++)ans=max(ans,f[i]);
printf("%d\n",ans);
return 0;
}

+++++++++++++++++++++++++++++++++++++++++++

+本文作者：luyouqi233。 　　　　　　　　　　　　　　+

+欢迎访问我的博客：http://www.cnblogs.com/luyouqi233/+

+++++++++++++++++++++++++++++++++++++++++++