强连通分量一般用于缩点,可以将有环图变成DAG图,这时就可以使用树形dp
在有向图中,
- 按照1到n的顺序去dfs后序遍历原图,得到每个点的出栈序列q
- 按照q[n]………q[1]的顺序去dfs反图,依次得到每个强连通分量
(一共做两遍dfs)
模板:
int h[N][2],to[M][2],ne[M][2];
int cnt1,cnt0;
int add(int a,int b,int c){//1代表正向建图,0代表反向建图
if(c==1){
to[cnt1][c]=b;ne[cnt1][c]=h[a][c];h[a][c]=cnt1++;
}
else {
to[cnt0][c]=b;ne[cnt0][c]=h[a][c];h[a][c]=cnt0++;
}
}
int vis[N];
stack<int>rem;
void dfs1(int u){
vis[u]=1;
for( int i=h[u][1];i!=-1;i=ne[i][1]){
int j=to[i][1];
if(vis[j]==0) dfs1(j);
}
rem.push(u);
}
void dfs2(int u){
vis[u]=0;
for( int i=h[u][0];i!=-1;i=ne[i][0]){
int j=to[i][0];
if(vis[j]==1) {
f[j]=f[u];
dfs2(j);
这里可以维护强连通分量内点的某些性质例如最大值
}
}
return;
}
下面是重新建图的函数,注意初始化head
for( int i=1;i<=m;i++){
if(f[u[i]!=f[v[i]]])add(f[u[i]] ,f[v[i]],1);
}
经典例题:
点对统计(hdu)
给定一个有向图,统计有多少点对u,v(1≤u<v≤n)满足u可以到达v,且v可以到达u。
#include<bits/stdc++.h>
using namespace std;
#define N 100100
#define M 200100
typedef long long ll;
int h[N][2],to[M][2],ne[M][2];
int cnt1,cnt0;
int add(int a,int b,int c){//1代表正向建图,0代表反向建图
if(c==0) swap(a,b);
if(c==1){
to[cnt1][c]=b;ne[cnt1][c]=h[a][c];h[a][c]=cnt1++;
}
else {
to[cnt0][c]=b;ne[cnt0][c]=h[a][c];h[a][c]=cnt0++;
}
}
int vis[N];
stack<int>rem;
void dfs1(int u){
vis[u]=1;
for( int i=h[u][1];i!=-1;i=ne[i][1]){
int j=to[i][1];
if(vis[j]==0) dfs1(j);
}
rem.push(u);
}
int cnt=0;
void dfs2(int u){
vis[u]=0;cnt++;
for( int i=h[u][0];i!=-1;i=ne[i][0]){
int j=to[i][0];
if(vis[j]==1) dfs2(j);
}
return;
}
ll c ( ll a,ll b){
if(a%2==0){
return a/2*(a-1);
}
else return (a-1)/2*a;
}
int main(){
int t;
cin>>t;
while(t--){
int n,m;
cin>>n>>m;
memset(h,-1,sizeof(h));
cnt1=0;
cnt0=0;
for( int i=1;i<=m;i++){
int a,b;
scanf("%d%d",&a,&b);
add(a,b,1);
add(a,b,0);
}
if(n==1) {
cout<<0<<endl;
continue;
}
while(!rem.empty()) rem.pop();
memset(vis,0,sizeof(vis));
ll ans=0;
for( int i=1;i<=n;i++) if(vis[i]==0) {
dfs1(i);
}
while(!rem.empty()) {
int j=rem.top();rem.pop();
if(vis[j]==1) {
dfs2(j);
ans+=c((ll)cnt,2);cnt=0;
}
}
cout<<ans<<endl;
}
}
强连通分量与缩点
缩点的原理是利用并查集的方法选出一个点代表强连通分量中的所有点
int f[N];
void dfs1(int u){
vis[u]=1;
for( int i=h[u][0];i!=-1;i=ne[i][0]){
int j=to[i][0];
if(vis[j]==0) dfs1(j);
}
rem.push(u);
}
//这里初始化并查集
void ini(){
for( int i=1;i<=n;i++) f[i]=i;
}
void dfs2(int u){
vis[u]=0;
for( int i=h[u][1];i!=-1;i=ne[i][1]){
int j=to[i][1];
if(vis[j]==1) {
f[j]=f[u];
dfs2(j);
ri[u]=max(ri[u],ri[j]);
}
}
}
for( int i=1;i<=m;i++){
if(f[u[i]!=f[v[i]]])add(f[u[i]] ,f[v[i]],1);
}