您的位置:首页 > 其它

强连通分量(kosaraju)

2015-07-30 12:22 246 查看
强连通分量特点

1.有向图

2.强连通分量中任意2个点可以互达

3.在一个非强连通图中极大的强连通子图就是该图的强连通分量

4.一个非强连通图可以有多个强连通分量

kosaraju算法

事实:置换图(将图中的边的方向换向)和原始图有相同的强连通分量

描述:

1.首先使G为一个有向图,并且S是一个空栈
2.while S非空:任意选择一个不在栈S中的顶点。从这个顶点v出发运行一个DFS(深搜). 每次dfs扩展到一个顶点u,就将u压入栈S(按照深搜顺序将点编号)
3.反向图G中所有的边获取一个置换图
4. while S非空:取出栈顶元素v. 从v开始运行一个DFS. 访问的结点组成了一个包含v在内的强连通分量,记录,然后从图G和栈S中去掉这些顶点。同样的,BFS也可以来做。
注意:
第一次深搜的时候从搜索头开始编号,最后编号的为搜索点,第二次深搜的时候先增加点再开始搜索其余点
具体见
int ord[maxn];//正向搜索,顶点的编号
bool vis[maxn];//标记顶点是否遍历过
int out[maxn];//转化为DAG以后的每个缩点的出度
int belong[maxn];//当前顶点属于哪个集合,相当于染色,当前顶点被染成了什么颜色
int ans[maxn];//每种颜色包括多少顶点,也就是强联通分量点的个数
int color;//代表不同的颜色,即不同的联通分量
int no;//正向搜索排序的编号(1---n)
int n,m;//顶点数和边数
vector<int> g[maxn],gre[maxn];//存储正向图和逆图

void dfs1(int u)//从当前u顶点开始DFS
{
vis[u] = 1;
for(int i = 0; i < g[u].size(); i++)
{
int v = g[u][i];
if(!vis[v]) dfs1(v);
}
ord[no++]=u;//为每个顶点编号
}

void dfs2(int u)
{
vis[u] = 1;
belong[u] = color;//当前顶点u被染成了color
for(int i = 0; i < gre[u].size(); i++)
{
int v = gre[u][i];
if(!vis[v])
{
ans[color]++;
dfs2(v);
}
}
}

int kosaraju()//若仅有1个连通分量的出度为0,返回出度为0的那个连通分量的点的个数
{
color = 1,no = 1;
for(int i = 1; i <= n; i++)//初始化每个点为1个联通分量,即联通分量的点的个数为1
ans[i]=1;
memset(out,0,sizeof(out));//初始化每个联通分量的出度
memset(vis,0,sizeof(vis));
for(int i = 1; i <= n; i++)//搜索每个连通分量并对其进行正序搜索编号
if(!vis[i]) dfs1(i);
memset(vis,0,sizeof(vis));//一定要再次进行遍历初始化
for(int i = no-1; i >= 1; i--)//编好号以后,从排号最大的开始搜索
{
int v = ord[i];//最后入栈的点
if(!vis[v])
{
dfs2(v);
color++;//记录联通分量个数
}
}
//构造DAG
for(int i = 1; i <= n; i++)
{
for(int j = 0; j < g[i].size(); j++)
{
if(belong[i] == belong[g[i][j]]) //若2个点属于同一个联通分量
continue;
out[belong[i]]++;//若有边相连的2个点不是同一个联通分量,则该联通分量的出度+1
}
}
//查找每种颜色出度为0的点或缩点,out[i]代表第i种颜色的点或缩点的出度,出度为0唯一,输出该颜色集合中点的个数ans[i],否则输出0
int ret = 0;//记录出度为0的联通分量的点的个数,若有多于2个联通分量的出度为0,表示此图不连通
int cnt = 0;//记录出度为0的点个数,若出度为0的联通分量只有1个,即被其他点连通的点的个数为该联通分量的点数
for(int i = 1;i < color; i++)//联通分量的个数
{
if(!out[i])//若出度为0
{
cnt++;
ret = ans[i];
}
if(cnt == 2)//这两个if不能写颠倒,切记啊!
{
ret = 0;
break;
}
}
return ret;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: