由有向图的邻接矩阵生成其可达矩阵
矩阵是研究图的性质的最有效工具之一,可运用矩阵运算求出图的路径、回路和其它性质。有些时候我们需要知道所给的图G中的某两个点之间是否存在边,为此前人引入了可达矩阵。定义不再赘述,在此给出一个由公式实
现的代码(P(G)= A(0) ∨A(1) ∨A(2) ∨ ··· ∨A(n-1))。
由于求解过程中需要多次使用矩阵相乘,故可将其单独拿出作为自定义的函数以减少总代码的篇幅。具体代码如下:
#include<stdio.h>
#define N 100
int mult(int a[N][N],int b[N][N],int n,int c[N][N])
{int i,j,k;for(i=0;i<n;i++){for(j=0;j<n;j++){c[i][j]=0;}}//每次接收矩阵乘积前先初始化为0 for(i=0;i<n;i++) { for(j=0;j<n;j++) { for(k=0;k<n;k++) { c[i][j]+=a[i][k]*b[k][j]; } } }
}//定义实现矩阵乘法的函数
int main()
{int n,a[100][100]={0},i,j,k,b[100][100]={0},c[100][100]={0},d[100][100]={0};scanf("%d",&n);for(i=0;i<n;i++){for(j=0;j<n;j++){scanf("%d",&a[i][j]);b[i][j]=a[i][j];d[i][j]=a[i][j];}}//输入邻接矩阵,并用另一个矩阵暂时存储 for(i=0;i<n-1;i++){mult(a,b,n,c);for(j=0;j<n;j++){for(k=0;k<n;k++){b[j][k]=c[j][k];//用c来接收前两个矩阵相乘的结果,并将其储存在b中以实现求解高次矩阵 d[j][k]+=c[j][k];}}} for(i=0;i<n;i++){for(j=0;j<n;j++){if(i==j) d[i][j]=1;else if(i!=j&&d[i][j]>0){d[i][j]=1;}}}//将矩阵d转化为可达矩阵 for(i=0;i<n;i++){for(j=0;j<n;j++){printf("%d ",d[i][j]);}printf("\n");}
}
也可以使用Warshall算法求解,不过由于Warshall算法的时间复杂度为O(n^3),故不适合计算大量数据(Warshall算法一般用于求解多源最短路径及传递闭包),若想进一步提高运算速度,可参看tarjan算法,对于普通的100*100的邻接矩阵,其速度大概能提高100倍左右。
由有向图的邻接矩阵生成其可达矩阵相关推荐
- 小世界网络邻接矩阵生成——python
小世界网络邻接矩阵生成--python 网络生成方法 python生成小世界网络最便捷的方式是导入networkx库,使用random_graphs.watts_strogatz_graph(n, k ...
- 有向图的邻接矩阵和邻接表画法
本例子展示有向图画邻接表和邻接矩阵,反过来也是如此 大三时自学了了数据结构中的图,发现图中有邻接矩阵和邻接表,作为笔记记录下来 如下图所示,给了一个有向图,求邻接矩阵 怎么做呢? 第一步:和无向图类似 ...
- 邻接矩阵中啥时候写0和无穷_带权有向图的邻接矩阵中0和无穷大
信息举报 时间:2020-12-23 本页为您甄选多篇描写带权有向图的邻接矩阵中0和无穷大,带权有向图的邻接矩阵中0和无穷大精选,带权有向图的邻接矩阵中0和无穷大大全,有议论,叙事 ,想象等形式.文章 ...
- matlab可达矩阵 结果,matlab求可达矩阵
matlab中可达矩阵的计算 2页 1下载券 邻接矩阵的可达矩阵及层... 2页 2下载券 邻接矩阵的一种计算方法 7页 免费 ISM算法(邻接矩阵求可达... 暂无...... JAVA矩阵算法代码 ...
- WarShall算法求传递闭包(可达矩阵)
最近在复习离散数学,顺便记录记录自己对warshall算法的理解. 1.传递闭包(可达矩阵) 传递闭包是有向图的一个重要性质,它指的是在有向图中从任意一个节点出发,可以到达的所有节点的集合.在某些应用 ...
- 用matlab做矩阵变换(理解+实践);假设随机发给某班所有同学任意数字,生成10*5的矩阵
问题 (1)班主任来查课,并带走 最后俩排同学谈话,显示剩余同学的矩阵A1. (2)在A中,使第i行的同学拿到的数字i显示该新矩阵A2. (3)在A中找到其中任意三行三列的同学,使他们的数字平方,显示 ...
- java 雷达图_POI生成原生Excel-雷达图
使用POI结合Excel的行列生成雷达图 一.先放最后输出的结果图: 雷达图 二.代码 @Test public void createRadar() throws IOException { fin ...
- 离散数学 —— 二元关系(图、零图与平凡图、度、握手定理、平行边、简单图与完全图、补图、子图与生成子图、同构、通路与回路、点与边割集、最短路线问题、强弱联通图、邻接矩阵与可达矩阵、欧拉图、平面图等)
5.1 关系及其表示: 补充一点: 简单路径:无重复边(点不作要求) 基本路径(初级路径):无重复点 一个路径是基本路径,那么必是简单路径 5.2 路与回路: 5.3 图的矩阵 ...
- 【有限马尔科夫链状态分解+Kosaraju 算法】基于Kosaraju 算法和可达矩阵的有限马尔科夫链状态分解
有限马尔科夫链状态分解+Kosaraju 算法 1 实验内容 2 理论基础 3 题目分析 4 按常返性和互通性对状态空间进行分解算法流程 4.1 强连通性和强连通分量 4.2 基于有向图 Kosara ...
- c++可达矩阵及连通性_3.9秒破百,矩阵式LED大灯加持,这台国产车可真香
比亚迪汉,是继唐.宋.秦.元之后,第五款以我国历史朝代命名的比亚迪量产车,汉朝是继秦朝之后的大一统王朝,享国四百零七年.将这一伟大王朝作为全新车型的名称,它不仅是对繁华盛世的记载,更代表着对未来的向往 ...
最新文章
- Alphabet量子公司横空出世!Sandbox将与谷歌、DeepMind成姊妹
- Matlab2013a安装简单教程以及遇到的问题解决(反复提示激活问题)
- mybatis_helloword(1)
- 《漫画算法》源码整理-4 大顶堆 小顶堆 优先队列
- Taro+react开发(18)--定义变量值
- redis-配置主从-实际操作-over
- matlab吧结果存数组,求助:如何将带有符号变量的运算结果储存到数组中
- Black Hat USA 2020 大会主议题大盘点(上)
- VC 项目支撑文件解释
- Maven3.5安装配置
- 【人脸质量评估】SDD-FIQA基于人脸相似度分布距离的无监督质量评估方法
- Java Date总结
- LM2586S 应用笔记
- 欧莱雅迎来入华25周年;万代南梦宫集团将启用新logo | 美通企业日报
- 截止频率计算公式wc_已知低通滤波器的传递函数是G(s)=(G0*Wc)/(s+Wc),截止频率不超过2HZ,怎么求?...
- C语言:str函数(一)
- 奥丁神叛虚拟机 台服登录不上游戏 账号无法登录游戏解决办法
- 让智能更落地,让中台更智能
- Linux“挂载”是什么意思
- 读取应用程序/类库配置文件(比如***.dll.config)的方法小结