Summary ------ Tarjan 连通分量

强连通分量

Tarjan算法

模板 牛的舞会
求图中点数大于1的强连通分量个数 即为答案
做法： tarjan判强连通分量，每个强连通分量统计结点个数，

Tarjan算法：
dfn[] 第i个结点第几个背访问到
low[] 第i个结点可以走到的 dfn[] 最小的结点 的 dfn[] 的值
stack[] 栈
vis[] 第i个结点是否在栈里
color[] 第i个强连通分量的相关信息

void tarjan(int u){
	dfn[u]=++deep;  
	low[u]=dfn[u];      //第一次访问结点， low即为它本身dfn
	vis[u]=1;
	stack[++top]=u;  //压入栈
	for(int i=head[u];i;i=e[i].nxt){    //枚举它可以走到的点
		int v=e[i].v;
		if(!dfn[v]){    //  没搜过就递归搜索
			tarjan(v);
			low[u]=min(low[u],low[v]);   // 用 low[v]更新low[u]
		}
		else{  //已经搜过
			if(vis[v]){
				low[u]=min(low[u],dfn[v]);     //用dfn[v] 更新 low[u] 
			}
		}
	}
	if(dfn[u]==low[u]){    // Tarjan核心   若一点dfn == low 则它所能走到的点和它构成一个强连通分量
		color[++sum]++;
		vis[u]=0;
		while(stack[top]!=u){    //将这个强连通分量的结点弹出栈内
			color[sum]++;
			vis[stack[top--]]=0;
		}
		top--;
	}
}

//牛的舞会AC code

#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long ll;
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=x*10+ch-'0';ch=getchar();}
	return x*f;
}
int n,m,cnt;
int head[10001];
struct node{
	int nxt;
	int v;
}e[50001];
void add(int u,int v){
	e[++cnt].nxt=head[u];
	e[cnt].v=v;
	head[u]=cnt;
}
int deep;
int dfn[10001];
int low[10001];
bool vis[10001];
int stack[10001];
int color[10001];
int top,sum;
void tarjan(int u){
	dfn[u]=++deep;
	low[u]=dfn[u];
	vis[u]=1;
	stack[++top]=u;
	for(int i=head[u];i;i=e[i].nxt){
		int v=e[i].v;
		if(!dfn[v]){
			tarjan(v);
			low[u]=min(low[u],low[v]);
		}
		else{
			if(vis[v]){
				low[u]=min(low[u],dfn[v]);
			}
		}
	}
	if(dfn[u]==low[u]){
		color[++sum]++;
		vis[u]=0;
		while(stack[top]!=u){
			color[sum]++;
			vis[stack[top--]]=0;
		}
		top--;
	}
}
int main(){
	n=read(),m=read();
	for(int i=1;i<=m;i++){
		int u,v;
		u=read(),v=read();
		add(u,v);
	}
	for(int i=1;i<=n;i++){
		if(!dfn[i]) tarjan(i);
	}
	int ans=0;
	for(int i=1;i<=sum;i++){
		if(color[i]>1) ans++;
	}
	cout<<ans<<endl;
	return 0;
}

缩点

还是Tarjan
找出强连通分量，一个强连通分量里的点点权加和，用连接不同强连通分量的边建新图
图成为DAG 可以dp等操作 强连通分量里的点可以特殊处理

模板：缩点

// luogu-judger-enable-o2
#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long ll;
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=x*10+ch-'0';ch=getchar();}
	return x*f;
}

int n,m,cnt1,cnt2,top,sum,deep;
int val[10001],head1[10001],belong[10001],head2[10001],in[10001],topo[10001],value[10001];
int dfn[10001],low[10001],stack[10001],f[10001];
bool vis[10001];

struct node{
	int nxt,v;
}e1[100001],e2[100001];


void add1(int u,int v){
	e1[++cnt1].nxt=head1[u];
	e1[cnt1].v=v;
	head1[u]=cnt1;
}

void add2(int u,int v){
	e2[++cnt2].nxt=head2[u];
	e2[cnt2].v=v;
	head2[u]=cnt2;
}


void tarjan(int u){
	dfn[u]=++deep;
	low[u]=deep;
	vis[u]=1;
	stack[++top]=u;
	for(int i=head1[u];i;i=e1[i].nxt){
		int v=e1[i].v;
		if(!dfn[v]){
			tarjan(v);
			low[u]=min(low[u],low[v]);
		}
		else{
			if(vis[v]){
				low[u]=min(low[u],dfn[v]);
			}
		}
	}
	if(dfn[u]==low[u]){
		belong[u]=++sum;
		vis[u]=0;
		while(stack[top]!=u){
			belong[stack[top]]=sum;
			vis[stack[top--]]=0;
		}
		top--;
	}         
}


int main(){
	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();
		add1(u,v);
	}
	for(int i=1;i<=n;i++){
		if(!dfn[i]){
			tarjan(i);
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=head1[i];j;j=e1[j].nxt){
			int v=e1[j].v;
			if(belong[i]!=belong[v]){
				add2(belong[i],belong[v]);
				in[belong[v]]++;
			}
		}
	}
	
	int cnt=0;
	for(int i=1;i<=sum;i++){
		if(in[i]==0){
			topo[++cnt]=i;
		}
	}
	
	for(int i=1;i<=sum;i++){
		int u=topo[i];
		for(int j=head2[u];j;j=e2[j].nxt){
			int v=e2[j].v;
			in[v]--;
			if(in[v]==0){
				topo[++cnt]=v;
			}
		}
	}
	
	for(int i=1;i<=n;i++){
		value[belong[i]]+=val[i];
	}
	for(int i=1;i<=sum;i++){
		f[i]=value[i];
	}

	for(int i=1;i<=sum;i++){
		int u=topo[i];
		for(int j=head2[u];j;j=e2[j].nxt){
			int v=e2[j].v;
			f[v]=max(f[v],f[u]+value[v]);
		}
	}
	int ans=0;
	for(int i=1;i<=sum;i++){
		ans=max(ans,f[i]);
	}
	cout<<ans<<endl;
	return 0;
}