专题讲座4 图的连通性和Tarjan算法 学习心得

题目都好抽象，但又有点套路

题目地址：

ZJNU-2022暑期专题-图的连通性和tarjan - Virtual Judge (vjudge.net)

目录

概念：

板子例题：

强连通分量

拓扑排序

割点（点双）

欧拉回路 

最小树形图

实战例题分析：

---

概念：

强连通分量---有向边

极大连通子图------无向边（还有边双和点双的概念）

板子例题：

强连通分量

struct node{
	int to,nxt;
}e[N];
int head[N];
int cnt=1,tot=0,num=0;
stack<int> s;
int dfn[N],low[N],scc[N];
void add(int u,int v){
	e[cnt].to=v;
	e[cnt].nxt=head[u];
	head[u]=cnt++;
}
void tarjan(int x){
	if (dfn[x]) return;
	dfn[x]=low[x]=++num;
	s.push(x);
	for (int i=head[x];i;i=e[i].nxt){
		int v=e[i].to;
		if (!dfn[v]){
			tarjan(v);
			low[x]=min(low[x],low[v]);
		}
		else if (!scc[v]){
			low[x]=min(low[x],dfn[v]);
		}
	}
	if (dfn[x]==low[x]){
		tot++;
		while(true){
			int t=s.top();
			s.pop();
			scc[t]=tot;
			if (t==x){
				break;
			}
		}
	}
}

拓扑排序

queue<int> q;
int tuopu(){
	rep(i,1,n){
		if (scc[i]==i&&!ein[i]){
			q.push(i);
			dis[i]=p[i];
		}
	}
	while (!q.empty()){
		int u=q.front();q.pop();	
		for (int i=h[u];i;i=e[i].nxt){
			int v=e[i].to;
			dis[v]=max(dis[v],dis[u]+p[v]);
			ein[v]--;
			if (ein[v]==0) q.push(v);
		}
	}
	int ma=0;
	rep(i,1,n){
		ma=max(ma,dis[i]);
	}
	return ma;
}

割点（点双）

struct node{
	int to,nxt;
}e[N];
int head[N];
int cnt=1,tot=0,num=0;
int dfn[N],low[N],scc[N];
int cut[N];
void add(int u,int v){
	e[cnt].to=v;
	e[cnt].nxt=head[u];
	head[u]=cnt++;
}
void tarjan(int x,int mr){//mr是本次搜索树的根 
	dfn[x]=low[x]=++num;
	int child=0;
	for (int i=head[x];i;i=e[i].nxt){
		int v=e[i].to;
		if (!dfn[v]){
			tarjan(v,mr);
			low[x]=min(low[x],low[v]);
			if (low[v]>=dfn[x]&&x!=mr) cut[x]=true;
			if (x==mr) child++;
		}
		low[x]=min(low[x],dfn[v]);
	}
	if (x==mr&&child>=2) cut[mr]=true;
}
void work(){
	int n,m;
	cin>>n>>m;
	int x,y;
	rep(i,1,m){
		cin>>x>>y;
		add(x,y);
		add(y,x);
	}
	rep(i,1,n){
		if (!dfn[i]) tarjan(i,i);
	}
	int ans=0;
	rep(i,1,n){
		if (cut[i]) ans++;
	}
	cout<<ans<<endl;
	rep(i,1,n){
		if (cut[i]) cout<<i<<" ";
	}
}

欧拉回路 

char s[5];
multiset<int> mp[N];
int du[N];
char ans[N];
int n;
void findoula(int u){
	while(mp[u].size()){
		int v=*mp[u].begin();
		mp[u].erase(mp[u].begin());
		mp[v].erase(mp[v].find(u));
		findoula(v);
	}
	ans[n--]=u;
}
void work(){
	cin>>n;
	rep(i,1,n){
		cin>>s;
		mp[s[0]].insert(s[1]);
		mp[s[1]].insert(s[0]);
		du[s[0]]++;
		du[s[1]]++;
	}
	int cnt=0,head=0;
	rep(i,64,125){
		if (du[i]&1){
			cnt++;
			if (head==0) head=i;
		}
	}
	if (cnt&&cnt!=2){
		cout<<"No Solution"<<endl;
		return;
	}
	if (head==0){
		rep(i,64,125){
			if (du[i]){
				head=i;
				break;
			}
		}
	}
	findoula(head);
	cout<<ans<<endl;
}

最小树形图

const int N=2e3+5;
const int M=2e5+5;
const int inf=2e9;
struct Edge{
	int u,v,w;
} e[M];
int n,m,root;
int pre[N],ine[N];
int vis[N],id[N];
int zhuliu(){
	int ans=0;
	while(1){
		rep(i,1,n) ine[i]=inf; 
		rep(i,1,m){
			int u=e[i].u,v=e[i].v;
			if(u!=v&&e[i].w<ine[v])  
				ine[v]=e[i].w,pre[v]=u;
		}
		rep(i,1,n)  
			if(i!=root&&ine[i]==inf) return -1;
		int cnt=0;
		rep(i,1,n) vis[i]=id[i]=0;
		rep(i,1,n){
			if(i==root) continue;
			ans+=ine[i];
			int v=i;
			while(vis[v]!=i&&!id[v]&&v!=root) {
				vis[v]=i;
				v=pre[v];
			}
			if(!id[v]&&v!=root){
				id[v]=++cnt; 
				for(int u=pre[v];u!=v;u=pre[u])
					id[u]=cnt;
			}
		}
		if(cnt==0) break; 
		rep(i,1,n)
			if(!id[i]) id[i]=++cnt;
		rep(i,1,m){
			int u=e[i].u,v=e[i].v;
			e[i].u=id[u],e[i].v=id[v];
			if(id[u]!=id[v]) e[i].w-=ine[v];
		}
		root=id[root];
		n=cnt;
	}
	return ans;
}
void work(){
    cin>>n>>m>>root;
    rep(i,1,m){
        cin>>e[i].u>>e[i].v>>e[i].w;
        if (e[i].u==e[i].v){
        	e[i].w=inf;
		}
    }
    cout<<zhuliu()<<endl;
}

实战例题分析：

#include <bits/stdc++.h>
#define rep(i,l,r) for (int i=l;i<=r;i++)
#define nep(i,l,r) for (int i=r;i>=l;i--)
#define pii pair<int,int>
#define int long long
#define CIO std::ios::sync_with_stdio(false)
using namespace std;
const int N=2e5+5;
const int MAX=1e9+7;
struct node{
	int to,nxt;
}e[N];
int head[N];
int cnt=1,tot=0,num=0;//保持链式前向星 从一开始的习惯
stack<int> s;
int dfn[N],low[N],scc[N],rd[N],sum[N];
int val[N];
void add(int u,int v){
	e[cnt].to=v;
	e[cnt].nxt=head[u];
	head[u]=cnt++;
}
void tarjan(int x){
	dfn[x]=low[x]=++num;
	s.push(x);
	for (int i=head[x];i;i=e[i].nxt){
		int v=e[i].to;
		if (!dfn[v]){//没有被dfs遍历到
			tarjan(v);
			low[x]=min(low[x],low[v]);
		}
		else if (!scc[v]){//遍历到后的三种非法边
			low[x]=min(low[x],dfn[v]);
		}
	}
	if (dfn[x]==low[x]){
		tot++;
		while(true){
			int t=s.top();
			s.pop();
			scc[t]=tot;//进行缩点，把环压缩成一个点，放入scc中
			sum[tot]++;//sum一般是根据题目意思变换
			if (t==x){
				break;
			}
		}
	}
}
void work(){
	int n,m;
	cin>>n>>m;
	int x,y;
	rep(i,1,m){
		cin>>x>>y;
		add(x,y);
	}
	rep(i,1,n){
		if (!dfn[i]){
			tarjan(i);
		}
	}
	rep(i,1,n){
		for(int j=head[i];j;j=e[j].nxt){//两层循环遍历，用缩点后的无环图寻找出入度求解
			int v=e[j].to;
			if (scc[i]!=scc[v]){
				rd[scc[i]]++;
			}
		}
	}
	int ans=0;
	rep(i,1,tot){
		if (!rd[i]){
			if (ans){
				cout<<'0'<<endl;
				return;
			}
			ans=i;
		}
	}
	cout<<sum[ans]<<endl;
}
signed main(){
	CIO;
	{
		work();
	}
	return 0;
}