图论 —— 图的遍历 —— 哈密顿问题

2024-06-17 21:48
文章标签 遍历 问题 图论 哈密顿

本文主要是介绍图论 —— 图的遍历 —— 哈密顿问题,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

【基本概念】

  1. 哈密尔顿通路:经过图中每个结点且仅经过一次的通路。
  2. 哈密尔顿回路:经过图中每个结点且仅经过一次的回路。
  3. 哈密尔顿图:存在哈密尔顿回路的图。
  4. 竞赛图:每对顶点之间都有一条边相连的有向图,n 个顶点的竞赛图称为 n 阶竞赛图。
  5. 与欧拉回路的对比:欧拉回路是指不重复地走过所有路径的回路;哈密尔顿回路是指不重复地走过所有点并且最后回到起点的回路。

【判定】

1.哈密尔顿通路的判定

设一无向图有 n 个顶点,u、v 为图中任意不相邻的两点,deg(x) 代表 x 的度数

若 deg(u)+deg(v)\geqslant n-1  成立,则图中存在哈密尔顿通路

2.哈密尔顿回路的判定:Dirac 定理

设一无向图有 n 个顶点,u、v 为图中任意不相邻的两点,deg(x) 代表 x 的度数

deg(u)+deg(v)\geqslant n,则图中存在哈密尔顿回路

推论:对于 n \geqslant 3 的无向图,若其任一一点 u 的度数 deg(u)\geqslant \frac n2 ,则图中存在哈密尔顿回路

3.竞赛图

对于 n \geqslant 2 的竞赛图,一定存在哈密尔顿通路

【dfs 搜索求哈密尔顿回路】

以每个点为起点进行搜索,直到形成回路

#include <iostream>  
#include <cstring>  
#define N 101
using namespace std;  
int n,m;
int u,v;
int g[N][N];
int vis[N],appear[N];
int ans[N],num[N];
int length;void dfs(int last,int i)//last表示上次访问的点 
{  vis[i]=1;//标记为已经访问过  appear[i]=1;//标记为已在一张图中出现过  ans[length++]=i;//记录答案  for(int j=1;j<=num[i];j++)  {  if(g[i][j]==x&&g[i][j]!=last)//回到起点构成哈密顿环  {   ans[++length]=g[i][j];//存储答案for(int i=1;i<=length-1;i++) //找到了一个环,输出anscout<<ans[i]<<' ';  cout<<ans[length]<<endl;length--;//长度-1break;}  if(!vis[g[i][j]])//遍历与i相关联的所有未访问的点。  dfs(i,g[i][j]); }  length--;  vis[i]=0;//回溯
}  
int main()  
{  memset(vis,0,sizeof(vis));  memset(appear,0,sizeof(appear));  cin>>n>>m; //读入点数与边数for(int i=1;i<=m;i++){   cin>>u>>v; //读入两点g[u][++num[v]]=v;//记录u-v的边g[v][++num[u]]=u;//记录v-u的边}  for(x=1;x<=n;x++) //枚举每一个点,将其作为起点来尝试访问{  if(!appear[x])//如果点x不在之前曾经被访问过的图里  {  length=0;//记录答案的长度  dfs(0,x);  }  }  return 0;  
} 

【Dirac 定理下构造无向图的哈密顿回路】

1.过程

  1. 任意找两个相邻的节点 S 和 T,在其基础上扩展出一条尽量长的没有重复结点的路径,即如果 S 与结点 v 相邻,而且 v 不在路径 S -> T 上,则可以把该路径变成 v -> S -> T,然后 v 成为新的 S。从 S 和 T 分别向两头扩展,直到无法继续扩展为止,即所有与 S 或 T 相邻的节点都在路径 S -> T 上
  2. 若 S 与 T 相邻,则路径 S -> T 形成了一个回路
  3. 若 S 与 T 不相邻,可以构造出来一个回路。设路径 S -> T 上有 k+2 个节点,依次为 S, v1, v2, ..., vk, T。可以证明存在节点 vi(i属于[1, k]),满足 vi 与 T 相邻,且 vi+1 与 S 相邻,找到这个节点 vi,把原路径变成 S -> vi -> T -> vi+1 -> S,即形成了一个回路.
  4. 到此为止,已经构造出来了一个没有重复节点的的回路,如果其长度为 N,则哈密顿回路就找到了。如果回路的长度小于 N,由于整个图是连通的,所以在该回路上,一定存在一点与回路之外的点相邻。那么从该点处把回路断开,就变回了一条路径,同时还可以将与之相邻的点加入路径。再按照步骤 1 的方法尽量扩展路径,则一定有新的节点被加进来,接着回到路径 2

2.伪代码

设 s 为哈密顿回路的起始点,t 为哈密顿回路中终点 s 之前的点,ans[] 为最终的哈密顿回路

  1. 初始化,令 s = 1,t 为 s 的任意一个邻接点.
  2. 如果 ans[] 中元素的个数小于 n,则从 t 开始向外扩展,如果有可扩展点 v,放入 ans[] 的尾部,并且 t=v,并继续扩展,如无法扩展进入步骤 3
  3. 将当前得到的 ans[] 倒置,s 和 t 互换,从 t 开始向外扩展,如果有可扩展点 v,放入 ans[] 尾部,并且 t=v,并继续扩展,如无法扩展进入步骤 4
  4. 如果当前 s 和 t 相邻,进入步骤 5,否则,遍历 ans[],寻找点 ans[i],使得 ans[i] 与 t 相连并且 ans[i +1] 与 s 相连,将从 ans[i+1] 到 t 部分的 ans[] 倒置,t=ans[i+1],进入步骤 5
  5. 如果当前 ans[] 中元素的个数等于 n,算法结束,ans[] 中保存了哈密顿回路(可看情况是否加入点 s),否则,如果 s 与 t 连通,但是 ans[] 中的元素的个数小于 n,则遍历 ans[],寻找点 ans[i],使得 ans[i] 与 ans[] 外的一点(j)相连,则令 s=ans[i-1],t=j,将 ans[] 中 s 到 ans[i-1] 部分的 ans[] 倒置,将 ans[] 中的 ans[i] 到 t 的部分倒置,将点 j 加入到 ans[] 的尾部,转步骤 2

3.时间复杂度

如果说每次到步骤 5 算一轮的话,那么由于每一轮当中至少有一个节点被加入到路径 S -> T 中,所以总的轮数肯定不超过 n 轮,所以时间复杂度为O(n^2)

空间上由于边数非常多,所以采用邻接矩阵来存储比较适合

4.实现

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<string>
#include<cstdlib>
#include<queue>
#include<set>
#include<map>
#include<stack>
#include<ctime>
#include<vector>
#define INF 0x3f3f3f3f
#define PI acos(-1.0)
#define N 1001
#define MOD 16007
#define E 1e-6
#define LL long long
using namespace std;
bool G[N][N];
bool vis[N];
int ans[N];void Reverse(int arv[N],int s,int t){//将数组anv从下标s到t的部分的顺序反向int temp;while(s<t){swap(arv[s],arv[t]);s++;t--;}
}void Hamilton(int n){int t;int s=1;//初始化取s为1号点for(int i=1;i<=n;i++)if(G[s][i]){t=i;//取任意邻接与s的点为tbreak;}memset(vis,false,sizeof(vis));vis[s]=true;vis[t]=true;ans[0]=s;ans[1]=t;int ansi=2;while(true){//从t向外扩展while(true){int i;for(i=1;i<=n;i++){if(G[t][i] && !vis[i]){ans[ansi++]=i;vis[i]=true;t=i;break;}}if(i>n)break;}//将当前得到的序列倒置Reverse(ans,0,ansi-1);//s和t互换swap(s,t);while(true){//从t继续扩展,相当于在原来的序列上从s向外扩展int i;for(i=1;i<=n;i++){if(G[t][i] && !vis[i]){ans[ansi++]=i;vis[i]=true;t=i;break;}}if(i>n)break;}//如果s和t不相邻,进行调整if(!G[s][t]){//取序列中的一点i,使得ans[i]与t相连,并且ans[i+1]与s相连int i;for(i=1;i<ansi-2;i++)if(G[ans[i]][t]&&G[s][ans[i+1]])break;i++;t=ans[i];Reverse(ans,i,ansi-1);//将从ans[i+1]到t部分的ans[]倒置}//此时s和t相连//如果当前序列包含n个元素,算法结束if(ansi==n)return;//当前序列中元素的个数小于n,寻找点ans[i],使得ans[i]与ans[]外的一个点相连int i,j;for(j=1;j<=n;j++){if(vis[j])continue;for(i=1;i<ansi-2;i++)if(G[ans[i]][j])break;if(G[ans[i]][j])break;}s=ans[i-1];t=j;//将新找到的点j赋给tReverse(ans,0,i-1);//将ans[]中s到ans[i-1]的部分倒置Reverse(ans,i,ansi-1);//将ans[]中ans[i]到t的部分倒置ans[ansi++]=j;//将点j加入到ans[]尾部vis[j]=true;}
}int main(){int n,m;while(scanf("%d%d",&n,&m)!=EOF&&(n||m)){n*=2;for(int i=0;i<=n;i++){for(int j=0;j<=n;j++){if(i==j){G[i][j]=false;G[j][i]=false;}else{G[i][j]=true;G[j][i]=true;}}}int ansi=0;memset(ans, 0, sizeof(ans));for(int i=1;i<=m;i++){int x,y;scanf("%d%d",&x,&y);G[y][x]=false;G[x][y]=false;}Hamilton(n);for(int i=0;i<n;i++)printf("%d ", ans[i]);printf("\n");}return 0;
}

【N 阶竞赛图下构造有向图的哈密顿通路】

含有N个顶点的有向图,且每对顶点之间都有一条边的图,一定存在哈密顿通路

int ans[105];
int map[105][105];
void Insert(int arv[], int &len, int index, int key){if(index>len)index=len;len++;for(int i=len-1; i>=0; i--){if(i!=index && i)arv[i]=arv[i-1];else{arv[i]=key;return;}}
}
void Hamilton(int n){int ansi = 1;ans[ansi++] = 1;for(int i=2; i<=n; i++){//第一种情况,直接把当前点添加到序列末尾if(map[i][ans[ansi-1]]==1)ans[ansi++]=i;else{int flag=0;//当前序列从后往前找到第一个满足条件的点j,使得存在<Vj,Vi>且<Vi,Vj+1>.for(int j=ansi-2; j>0; j--){if(map[i][ans[j]]==1){//找到后把该点插入到序列的第j+1个点前.flag=1;Insert(ans,ansi,j+1,i);break;}}if(!flag)//否则说明所有点都邻接自点i,则把该点直接插入到序列首端.Insert(ans,ansi,1,i);}}
}
int main(){int n,m;scanf("%d", &n);m=n*(n-1)/2;for(int i=0;i<m;i++){int u,v;scanf("%d%d",&u,&v);if(u<v)map[v][u]=1;}Hamilton(n);for(int i=1;i<=n;i++)printf(i==1? "%d":" %d",ans[i]);printf("\n");return 0;
}

这篇关于图论 —— 图的遍历 —— 哈密顿问题的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1070561

相关文章

springboot循环依赖问题案例代码及解决办法

《springboot循环依赖问题案例代码及解决办法》在SpringBoot中,如果两个或多个Bean之间存在循环依赖(即BeanA依赖BeanB,而BeanB又依赖BeanA),会导致Spring的... 目录1. 什么是循环依赖?2. 循环依赖的场景案例3. 解决循环依赖的常见方法方法 1:使用 @La

SpringBoot启动报错的11个高频问题排查与解决终极指南

《SpringBoot启动报错的11个高频问题排查与解决终极指南》这篇文章主要为大家详细介绍了SpringBoot启动报错的11个高频问题的排查与解决,文中的示例代码讲解详细,感兴趣的小伙伴可以了解一... 目录1. 依赖冲突:NoSuchMethodError 的终极解法2. Bean注入失败:No qu

MySQL新增字段后Java实体未更新的潜在问题与解决方案

《MySQL新增字段后Java实体未更新的潜在问题与解决方案》在Java+MySQL的开发中,我们通常使用ORM框架来映射数据库表与Java对象,但有时候,数据库表结构变更(如新增字段)后,开发人员可... 目录引言1. 问题背景:数据库与 Java 实体不同步1.1 常见场景1.2 示例代码2. 不同操作

如何解决mysql出现Incorrect string value for column ‘表项‘ at row 1错误问题

《如何解决mysql出现Incorrectstringvalueforcolumn‘表项‘atrow1错误问题》:本文主要介绍如何解决mysql出现Incorrectstringv... 目录mysql出现Incorrect string value for column ‘表项‘ at row 1错误报错

如何解决Spring MVC中响应乱码问题

《如何解决SpringMVC中响应乱码问题》:本文主要介绍如何解决SpringMVC中响应乱码问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录Spring MVC最新响应中乱码解决方式以前的解决办法这是比较通用的一种方法总结Spring MVC最新响应中乱码解

pip无法安装osgeo失败的问题解决

《pip无法安装osgeo失败的问题解决》本文主要介绍了pip无法安装osgeo失败的问题解决,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一... 进入官方提供的扩展包下载网站寻找版本适配的whl文件注意:要选择cp(python版本)和你py

解决Java中基于GeoTools的Shapefile读取乱码的问题

《解决Java中基于GeoTools的Shapefile读取乱码的问题》本文主要讨论了在使用Java编程语言进行地理信息数据解析时遇到的Shapefile属性信息乱码问题,以及根据不同的编码设置进行属... 目录前言1、Shapefile属性字段编码的情况:一、Shp文件常见的字符集编码1、System编码

Spring MVC使用视图解析的问题解读

《SpringMVC使用视图解析的问题解读》:本文主要介绍SpringMVC使用视图解析的问题解读,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录Spring MVC使用视图解析1. 会使用视图解析的情况2. 不会使用视图解析的情况总结Spring MVC使用视图

Redis解决缓存击穿问题的两种方法

《Redis解决缓存击穿问题的两种方法》缓存击穿问题也叫热点Key问题,就是⼀个被高并发访问并且缓存重建业务较复杂的key突然失效了,无数的请求访问会在瞬间给数据库带来巨大的冲击,本文给大家介绍了Re... 目录引言解决办法互斥锁(强一致,性能差)逻辑过期(高可用,性能优)设计逻辑过期时间引言缓存击穿:给

Java程序运行时出现乱码问题的排查与解决方法

《Java程序运行时出现乱码问题的排查与解决方法》本文主要介绍了Java程序运行时出现乱码问题的排查与解决方法,包括检查Java源文件编码、检查编译时的编码设置、检查运行时的编码设置、检查命令提示符的... 目录一、检查 Java 源文件编码二、检查编译时的编码设置三、检查运行时的编码设置四、检查命令提示符