【图论】连通性问题

已于 2024-04-19 12:39:25 修改
阅读量359 收藏 10
点赞数 5
文章标签： 图论 深度优先 算法
于 2024-04-06 09:05:15 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/2302_80811345/article/details/136678180
版权
本文介绍了如何使用倍增法求解LCA问题，以及涉及图论中的强连通分量的tarjan算法和割点的检测。还提及了与PopularCows相关的问题和可达性的概念。
摘要由CSDN通过智能技术生成

lca

1.1倍增法求lca

【manim | 算法】7分钟学会倍增法求解LCA_哔哩哔哩_bilibili

递归初始化

void dfs(long u,long father){
dep[u]=dep[father]+1;//只在这里初始化dep

for(long i=1;(1<<i)<=dep[u];i++)
    fa[u][i]=fa[fa[u][i-1]][i-1];//只这里用的倍增

for(long i=head[u];~i;i=edge[i].next){
    long v=edge[i].to;
    if(v==father)continue;
    fa[v][0]=u;
    dfs(v,u);
}

}

lca主体 



long lca(long x,long y){
if(dep[x]<dep[y])swap(x,y);

for(int i=20;i>=0;i--){//跳到同一个深度
    if(dep[fa[x][i]]>=dep[y])x=fa[x][i];
    if(x==y)return x;
}

for(int i=20;i>=0;i--){
    if(fa[x][i]!=fa[y][i]){//一起跳
        x=fa[x][i];
        y=fa[y][i];
    }
}
return fa[x][0];
}

开始交了几发都不对，以为必须用tarjan等快方法

今天看到输出少了一个，n-1条边，s的父是0，ac了

#include<iostream>
using namespace std;

long n,q,s;
struct EDGE{
long  next;
long to;

}edge[1000001];
long head[500001];
long tot=0;

void addedge(long u,long v){
edge[++tot].next=head[u];
edge[tot].to=v;

head[u]=tot;
}
long dep[500001];
long fa[500001][21];//最多跳20次
void dfs(long u,long father){
dep[u]=dep[father]+1;

for(long i=1;(1<<i)<=dep[u];i++)
    fa[u][i]=fa[fa[u][i-1]][i-1];

for(long i=head[u];~i;i=edge[i].next){
    long v=edge[i].to;
    if(v==father)continue;
    fa[v][0]=u;
    dfs(v,u);
}

}

long lca(long x,long y){
if(dep[x]<dep[y])swap(x,y);

for(int i=20;i>=0;i--){
    if(dep[fa[x][i]]>=dep[y])x=fa[x][i];
    if(x==y)return x;
}

for(int i=20;i>=0;i--){
    if(fa[x][i]!=fa[y][i]){
        x=fa[x][i];
        y=fa[y][i];
    }
}
return fa[x][0];
}


int main(){

cin>>n>>q>>s;
for(long i=1;i<=n;i++)head[i]=-1;

for(long i=1;i<=n-1;i++){
    long u,v;cin>>u>>v;
    addedge(u,v);addedge(v,u);
}
dfs(s,0);
while(q--){
    long a,b;cin>>a>>b;
    if(a==b)cout<<a<<endl;
    else
    cout<<lca(a,b)<<endl;
}
}

1.2Design the city

Design the city - ZOJ 3195 - Virtual Judge (vjudge.net.cn)

先mark着

#include<iostream>
using namespace std;
//连接三个点最短路径
//无向无环
//必然有一个深度最小

//就当图来做，这三个点的最小生成树

//两两找最近公共祖先

-------------------------------

2.1强连通分量

B3609 [图论与代数结构 701] 强连通分量 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

#include<iostream>
#include<vector>
#include<stack>
using namespace std;

stack<int> s;
int low[10001];int dfn[10001];
bool vis[10001];
int scc[10001];
int block_cnt;
int scc_cnt;
vector<int> scc_zong[10001];
vector<int> G[10001];
int n,m;

void tarjan(int x){
low[x]=dfn[x]=++block_cnt;
s.push(x);
vis[x]=true;

for(int i=0;i<G[x].size();i++){

    int v=G[x][i];
    if(!vis[v]){
            //为什么是low[v]
            //如果v是没走过，当场走的话，可以追溯v更新后的
            //如果v走过，但不属于联通分量，v就没有追溯的，现在收入
            //觉得一样，一回试试low[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]){
        scc_cnt++;
while(1){
int temp=s.top();
    s.pop();
    scc[temp]=scc_cnt;
    scc_zong[scc_cnt].push_back(temp);
    if(temp==x)break;
}

}
}

void solu(){

for(int i=0;i<=n;i++)vis[i]=0;

cout<<scc_cnt-1;
for(int i=1;i<=n;i++){
        if(!vis[i]){
                cout<<endl;
                vis[i]=1;
                cout<<i;
                    for(int j=i+1;j<=n;j++){
                    if(scc[i]==scc[j]){
                        vis[j]=1;
                    cout<<" "<<j;
                    }
            }


        }

}


//想欧拉筛
/*for(int i=1;i<=n;i++){
        if(!vis[i]){
            int gai_scc=scc[i];

            for(int j=0;i<gai_scc[])
        }

}*/


/*for(int i=scc_cnt;i>=1;i--){
    cout<<endl;
    for(int j=scc_zong[i].size()-1;j>=0;j--){
            cout<<scc_zong[i][j];
        if(j!=0)cout<<" ";
    }

}
*/

/*for(int i=1;i<=scc_cnt;i++){if(i!=1)cout<<endl;
    for(int j=0;j<scc_zong[i].size();j++){
        if(j!=0)cout<<" ";
        cout<<scc_zong[i][j];
    }*/



}


int main(){

cin>>n>>m;
for(int i=1;i<=n;i++)G[0].push_back(i);
//这样应该就能按顺序遍历1，2，3，点所在的强连通分量了
//啊啊错了
while(m--){
    int u,v;
    cin>>u>>v;
    G[u].push_back(v);
}

tarjan(0);//0不会只追溯到自己的
solu();


}

这里是加了一个万能源点0，从后面模板割点学到了tarjan整图不连通的另一个解决方案

for(int i=1;i<=n;i++)
if(!dfn[x])tarjan(x);

2.2牛的舞会

P2863 [USACO06JAN] The Cow Prom S - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

2.3popular cows

2186 -- Popular Cows (poj.org)

P2341 [USACO03FALL / HAOI2006] 受欢迎的牛 G - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

2.4可达性

可达性 (nowcoder.com)

每个入度为0的联通块取一个最小点


#include<iostream>
#include<vector>
#include<stack>
using namespace std;

bool ans[100001];
stack<long> s;
long low[100001];long dfn[100001];
bool vis[100001];
long color[100001];
long block_cnt;
long scc_cnt;
vector<long> scc[100001];
vector<long> G[100001];
long n,m;

void tarjan(long x){
low[x]=dfn[x]=++block_cnt;
s.push(x);
vis[x]=true;

for(long i=0;i<G[x].size();i++){

    long v=G[x][i];
    if(!vis[v]){

        tarjan(v);
        low[x]=min(low[x],low[v]);
    }
    else if(!color[v]){
        low[x]=min(low[x],dfn[v]);
    }
}

//这样走过一通后，如果是最栈底的点
if(dfn[x]==low[x]){
        scc_cnt++;
while(1){
long temp=s.top();
    s.pop();
    color[temp]=scc_cnt;
    scc[scc_cnt].push_back(temp);
    if(temp==x)break;
}

}
}

//有那种能到任意点的，说明暗暗给出图是弱联通
void solve(){
//如果u和v不是一个联通酷爱
//v记录入度

//入度为0的快加上

for(long i=1;i<=scc_cnt;i++){
    ans[i]=true;
}

for(long i=1;i<=n;i++){
    for(long j=0;j<G[i].size();j++){
        long v=G[i][j];
        if(color[i]!=color[v]){
            ans[color[v]]=0;

                  }
    }
}
long out=0;

for(long i=1;i<=scc_cnt;i++){//i是联通块编号
    if(ans[i]){out++;}}

    cout<<out<<endl;
//入度为0的块，每个取一个最小的点

bool flag=1;
for(long i=1;i<=scc_cnt;i++){//i是联通块编号
    if(ans[i]){
            long temp=n;
        for(long j=0;j<scc[i].size();j++){
            temp=min(temp,scc[i][j]);
        }
        if(flag==1)flag=0;
        else cout<<" ";
    cout<<temp;
    }
}
}

int main(){

cin>>n>>m;

while(m--){
    long u,v;
    cin>>u>>v;
    G[u].push_back(v);
}
for(long i=1;i<=n;i++){


if(!dfn[i]){tarjan(i);}
}


solve();


}

------------------------------- 

 模板缩点

 P3387 【模板】缩点 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

stl不知道怎么用，先mark

更新：stl会用了，dp不会了

#include<iostream>
#include<vector>
#include<stack>
using namespace std;

bool vis[10001];
int dfn[10001];
int low[10001];
int block_cnt;
stack<int> s;
vector<int> G[10001];

void tarjan(int x){
dfn[x]=low[x]=++block_cnt;
s.push(x);
vis[x]=1;


for(int i=head[x];~i;i=edge[i].next){
    int v=edge[i].to;
    if(!vis[v]){
        tarjan(v);
        low[x]=min(low[x],low[v]);

    }
    else if(!scc[v]){
        low[x]=min(low[x],dfn[v]);
    }
}
if(low[x]==dfn[x]){
    scc_cnt++;
    while(1){
        int temp=s.top();
        s.pop();
        scc[temp]=scc_cnt;
        quanscc[scc_cnt]+=quanpoint[temp];
        if(temp==x)break;
    }
}

}


void suo(){
for(int i=1;i<=n;i++){
    for(int j=0;j<G[u].size();j++){
        if(scc[i]!=scc[G[i][j]]){
            addedge(scc[i],scc[G[i][j]]);
        }
    }
}

//dp
//状压
}


int main(){
cin>>n>>m;
for(int i=1;i<=n;i++)cin>>quanpoint[i];
while(m--){
    cin>>u>>v;
    G[u].push_back(v);
}


}

-------------------------------

模板割点

P3388 【模板】割点（割顶） - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

#include<iostream>
#include<vector>
#include<stack>
using namespace std;


int n;long m;
stack<int> s;
int block_cnt;
int dfn[20001];
int low[20001];
int root;
bool ge[20001];
vector<int> G[20001];
//割点ge[x]=0
//非root,有儿子，low[v]>=dfn[x]
//root,没有连线

//传参root，就对了

void tarjan(int x,int root){
dfn[x]=low[x]=++block_cnt;
s.push(x);
int ch=0;

for(int i=0;i<G[x].size();i++){
    int v=G[x][i];

    if(!dfn[v]){
        ch++;//如果遍历过，说明是两个子树中间有边的情况
        //不是一个v一条边
        tarjan(v,root);//ch在tarjan(x)内定义，不会到v内
        low[x]=min(low[x],low[v]);
        if(x!=root&&low[v]>=dfn[x]){ge[x]=1;}
//无论如何，认为所有v的low都必须》=dfn[x],
            //只一个就决定了怎么能行

    }
    else {
        low[x]=min(low[x],dfn[v]);
    }

}

if(ch>=2&&root==x){
    ge[x]=1;
}

}

int main(){
cin>>n>>m;

while(m--){
    int u,v;
    cin>>u>>v;
    G[u].push_back(v);
    G[v].push_back(u);
}
for(int i=1;i<=n;i++){
    if(!dfn[i]){

            //root=i;
    //root是一个滚动的根
        tarjan(i,i);
    }
}

int ans=0;
for(int i=1;i<=n;i++){
    if(ge[i])ans++;
}
cout<<ans<<endl;

bool flag=1;
for(int i=1;i<=n;i++)

if(ge[i]){
    if(flag==1){flag=0;}
    else {cout<<" ";}
        cout<<i;}

}

arin876
关注
  • 5
    点赞
  • 踩
  • 10
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
图论——连通性
Albert_Jw的博客
08-01 1559
割点： 1.无向图 2.删去这个点及其所连边后，图不再联通 点双连通图： 1.无向图 2.没有割点（删去任意一个点图仍联通） 点双联通分量： 无向图G中所有子图G’ 如果G’ 1.是点双联通子图 2.不是其他点双联通子图的真子集， 则G’是G的极大点双联通子图，也称点双联通分量。 桥（割边）： 1.无向图 2.删此边（不删其连着的点），剩下的图不再联通 边双连通图： 1.无向图 2.删任意一边，剩下的图仍联通 边双联通分量： 无向图G中所有子图G’ 如果G’ 1.是边双联通子图 2.不是其他边双联通
【图论】连通性
weixin_44413191的博客
08-19 829
本文为图论复习，总结连通性相关知识点。欢迎交流
图论连通性问题总结
wwxy1995的博客
02-19 239
静态图的连通性问题不是必须用并查集，可以通过DFS/BFS解决。
图论的连通性问题
2301_80018137的博客
12-12 494
这两种算法可以用于判断图的连通性，即从一个节点出发，看是否能够访问到图中的所有节点。大佬: 不客气，连通性问题是图论中的一个基础概念，希望你在学习中能够更深入地理解它。大佬: 是的，如果两个节点是连通的，那么它们之间一定存在至少一条路径。连通性问题是图论中一个重要的主题，它关注的是图中节点之间的连接情况。小白: 明白了，连通性问题看起来很有用，可以用来判断图中节点之间的关系。小白: 你好，我听说图论可以解决一些连通性的问题，你能给我讲讲吗？小白: 那么，在图中如果有两个节点是连通的，它们之间肯定有路径吗？
图论 —— 图的连通性
热门推荐
Alex_McAvoy的博客
10-20 1万+
【基本概念】 1.连通图与连通分量 1）连通图：无向图 G 中，若对任意两点，从顶点 Vi 到顶点 Vj 有路径，则称 Vi 和 Vj 是连通的，图 G 是一连通图 2）连通分量：无向图 G 的连通子图称为 G 的连通分量 任何连通图的连通分量只有一个，即其自身，而非连通的无向图有多个连通分量 以上图为例，总共有四个连通分量，分别是：ABCD、E、FG、HI。...
图论中的连通性问题
500g的博客
02-24 427
无向图连通性问题： 有向图连通性问题： 并查集： 一种树型数据结构，处理一些不相交集合的合并及查询问题。通常在使用中以森林来表示。 因为我们不关心树的具体结构，所以不需要记录节点的儿子，只要记录节点的父亲即可。 初始时，每个节点自己就是树根。在合并两个集合时，找到这两个集合的根，把其中一个根的父亲设为另一个集合的根。 int find(int x){ if(!fa[x]) return x; return find(fa[x]); } void merge(int x,int y){ int rtx=
图论.连通性专题
weixin_30596343的博客
10-26 72
开一个关于图的连通性的专题,内容包括: 强连通分量 双连通分量 最大团 极大团 全局最小割 最短路相关 其他的想到了再补 转载于:https://www.cnblogs.com/Patt/p/6001434.html...
图的连通性
ZhuRanCheng的博客
06-03 2402
图的连通性总结
C++图论之强连通图
y6123236的专栏
12-26 1148
强连通分量算法还有算法。有兴趣者可自行了解。
图论------有向图的连通性问题
weixin_43743711的博客
09-22 6227
前言： 一些概念需要我们理解一下，以便更好地进行下面的内容。 在有向图G中，如果两个顶点间至少存在一条互相可达路径，称两个顶点强连通(strongly connected)。 如果有向图G的每两个顶点都强连通，称G是一个强连通图。 非强连通图有向图的极大强连通子图，称为强连通分量(strongly connected components)。 很明显强连通分量出现在有向有环图中，每一个连通的分量都可以被当作是强连通分量。如左下图每一个虚线框内都是一个强连通分量（ＳＣＣ），右下图则将强连通分量当作一个结点
图论中的连通性问题2222
05-06
图论中的连通性问题是图论中的一个重要问题，包括有向图的连通性问题和无向图的连通性问题。在有向图中，连通性问题可以分为强连通图和弱连通图两种。强连通图是指从图中的任意两个点A和B，存在一条从A到B的路径和一...
图论中的连通性问题2.111111
05-06
图论中的连通性问题2.111111 图论中连通性问题是研究图的连接性和分离性的一个重要部分。连通性问题的研究可以分为无向图和有向图两大类。在无向图中，我们可以使用深度优先搜索（Depth-First Search，DFS）来遍历...
图论 连通性 讲义 acm
07-12
总之，理解图的连通性及其相关的概念，如割点、割边、连通度和边连通度，是深入学习图论和解决ACM竞赛中涉及图的问题的基础。这些理论工具对于分析和设计算法至关重要，因为它们帮助我们识别和操纵图结构，从而找到...
matlab判断图的连通性
05-16
2. 调用`concom.m`函数处理图的连通性问题。 3. 显示结果，包括图的连通性状态（连通或非连通）以及每个顶点所在的连通块。 这两个脚本的组合提供了全面的图连通性分析工具，对于学习图论、图算法和MATLAB编程具有...
新建文本文档_计算连通性_
10-03
"计算连通性"是网络分析中的核心概念，它涉及到图论、数据结构和算法等多个方面。在此程序中，主要目标是计算并存储复杂网络图的连通性矩阵，这在社交网络分析、计算机网络、生物网络等领域都有广泛应用。 首先，...
图论·Day01
2301_80132162的博客
07-06 231
【代码】图论·题解1。
代码随想录算法训练营第64天：图论2[1]
2302_79934729的博客
07-03 667
当然广搜还有很多细节需要注意的地方，后面我会针对广搜的题目还做针对性的讲解。因为在理论篇讲太多细节，可能会让刚学广搜的录友们越看越懵，所以细节方面针对具体题目在做讲解。本篇我们重点讲解了广搜的使用场景，广搜的过程以及广搜的代码框架。其实在二叉树章节的层序遍历​**** 中，我们也讲过一次广搜，相当于是广搜在二叉树这种数据结构上的应用。这次则从图论的角度上再详细讲解一次广度优先遍历。相信看完本篇，大家会对广搜有一个基础性的认识，后面再来做对应的题目就会得心应手一些。
【leetcode52-55图论、56-63回溯】
qq_46277517的博客
07-03 138
图论 回溯
算法训练 | 图论Part4 | 107. 寻找存在的路径
最新发布
qq_48896570的博客
07-06 83
【代码】算法训练 | 图论Part4 | 107. 寻找存在的路径。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值