矩阵相乘的strassen算法_矩阵乘法的Strassen算法+动态规划算法(矩阵链相乘和硬币问题)...

矩阵乘法的Strassen

这个算法就是在矩阵乘法中采用分治法,能够有效的提高算法的效率。

先来看看咱们在高等代数中学的普通矩阵的乘法

两个矩阵相乘

上边这种普通求解方法的复杂度为: O(n3)

也称之为暴力求解或者朴素求解

这是暴力求解的代码,三重循环,显然复杂度是O(n3)

voidMul(int** matrixA,int** matrixB,int** matrixC)

{

for(inti = 0; i < 2; ++i)

{

for(intj = 0; j < 2; ++j)

{

matrixC[i][j] = 0;

for(intk = 0; k < 2; ++k)

{

matrixC[i][j] += matrixA[i][k] * matrixB[k][j];

}

}

}

由于不满足上边这种复杂度太大,德国某位牛人开始找另外的解法了

先分析一下下边的

将一个矩阵分成四块

如上图,A和B矩阵都被分成了四块,该算法复杂度依然是n3,于是上边那位老哥不服,他觉得这不是最优的解,还有更优的,于是他分析了上边是四个等式,四个等式中有八个乘法,四个加法

矩阵乘法的复杂度主要就是体现在相乘上,而多一两次的加法并不会让复杂度上升太多。故此,老哥思考,是否可以让矩阵乘法的运算过程中乘法的运算次数减少,从而达到降低矩阵乘法的复杂度,我们都知道,想要获取时间上的效率,很多时候都是以空间换时间,于是老哥定义了七个变量

这七个变量均是矩阵,ABCDEFGH原来两个相乘矩阵里边划分好的八个小矩阵

图三

或者看这个图,总之七个矩阵变量是要求的(PPT上和这差不多,只是变量顺序换了)

图四

求出则七个矩阵,就能求出A*B的值

这个图就是A*B的值,至于为什么能求出来,归功于牛人构造的七个巧妙的式子,利用七个式子之间的关系就求出了下边四个变量,也就是解

图五

最后那老哥证明了,这个复杂度是这个

图六

顺带复习一下PPT上这个

如对于上边图六那个公式,a=7,k=2,b=2  显然7>2^2,所以套第三个T(n)

动态规划算法

动态规划和分治法相似,都是通过组合字问题来求解原问题,不同之处在于分治法的子问题互不干涉、互不交叉,而动态规划相反,它会利用已经求解的子问题进而求解新的子问题

先举个简单的例子感受一蛤什么是动态规划

钱币问题——用面值1元、3元、5元的硬币,如何用最少的硬币凑到11块钱?

第一步要想的就是,怎么把一个大问题变小问题

既然要求最少的硬币凑到11块钱,这里用c[i]=表示凑到i元最小要j个硬币

那我先求最少的硬币凑到0块钱,显然需要0个硬币,所以才c[0]=0

接下来求最少的硬币凑到1块钱,现在只有面值1块的能用,我就用一个,用完之后还需凑0元,这时才c[0]=0已知,所以才c[1]=1+c[0]=1

接下来求最少的硬币凑到2块钱,现在只有面值1块的能用,我也先用一个,用完之后还需凑1元,这时才c[1]=1已知,所以c[2]=1+c[1]=2

接下来求最少的硬币凑到3块钱,现在有面值1块的和三块的,如果我先用一个3块的,用完之后还需凑0元,这时才c[0]=0已知,所以c[3]=1+c[0]=1;如果我先用一个1块的,用完之后还需凑2元,这时才c[2]=2已知,所以c[3]=1+c[2]=3,取这两种中的最小的那种情况

后边的以此类推....

矩阵链乘法

如果要求n个给定序列的矩阵相乘的乘积(比如ABCDEFG),矩阵具有结合律,所以计算的步骤有很多种选择,但如果结合律用的不好会产生比较大的代价

在了解这个咱们要研究算法是干啥的之前,先了解几个概念

1、矩阵相容:也就是两个矩阵要能够相乘,即A的列数等于B的行数

2、标量乘法:若A是p*q,B是 q*r,则A*B的代价就是其标量乘法,也就是pqr

所以要求n个给定序列的矩阵相乘的乘积,我们要研究使得该成绩代价最小,也就是其标量乘法次数之和最少(这块最好参照一下算法导论211页很详细),说白了,就是在乘法式子中如何打括号

官方的话就不说了,直接上一串矩阵,你应该干什么和怎么干,哈哈,怎么干

图中给出了6个矩阵相乘,你应该做的就是给它大括号,决定计算顺序,使得计算代价最小

这个就是m[ ][ ]的算法    int t= m[i][k]+ m[k+1][j]+p[i-1]*p[k]*p[j]          、

现在来解释一下上边的这个算法,我只能说老师的PPT略傻逼,都没怎么解释

m[i][j]表示矩阵从第i个矩阵乘到第j个矩阵的最小代价

int t= m[i][k]+ m[k+1][j]+p[i-1]*p[k]*p[j] :

上边这个算法的意思是,第i个矩阵到第k个矩阵相乘的代价+第k个矩阵到第j个矩阵相乘的代价,加上这两个乘好了的前后两个矩阵相乘的代价

然后理解了怎么算,从小到大算就OK了,按照斜线的顺序算,i和j挨着越近越好算,先算对角线,全是0,再算m[1][2],m[2][3],m[3][4]...以此类推,因为后边计算的斜线,会用到上一条斜线上那些数

比如算m[1][3]会用到m[1][1],m[1][2],m[2][3],m[3][3]

最后解释一下怎么找分解点,也就是在哪打括号,下边图的矩阵是标记矩阵,也就是在动态规划的过程中,你每次的最优解是在哪划分的它会记录下来,如果要求A[1][6]怎么打括号,找到s[1][6]=3,然后在A[1][3]和A[3][6]里重复上边的步骤,每个找到的点都是分界点..比如第一个找到的3

【6】】

来自moonsmile的祝福~

明天将推出贪心算法

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/349113.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

计算机c盘哪些东西可以清理,细说电脑c盘哪些文件可以删除

有些网友反映&#xff0c;自己看C盘里的文件太多了&#xff0c;电脑又太卡&#xff0c;情急之下就把里面的东西删掉了&#xff0c;现在系统都不能用了。为了避免大家再入这个坑&#xff0c;我给大家讲一下哪些是C盘里的无用文件&#xff0c;并且删除后不会影响系统使用C盘是指电…

springboot 注解动态赋值_SpringBoot 使用 @Value 注解读取配置文件给静态变量赋值

1、application.properties 配置文件CxU免费资源网mail.usernamexue163.commail.passwordxuemail.hostsmtp.163.commail.smtp.authtrue2、给普通变量赋值&#xff0c;直接在变量上添加 Value 注解CxU免费资源网import org.springframework.beans.factory.annotation.Value;publ…

软件测试度量计算方法有哪些,软件测试度量(三)

进度差异趋势6.4.3 范围变化(SC)这个指标指出如何固定测试范围。下面总范围 以前的范围 新范围&#xff0c;如果范围扩大的话总范围 以前的范围 - 新范围&#xff0c;如果范围缩小的话一个发布版本范围变化趋势7、结论度量是评估的重要组成部分以及任何业务改进的基础。是应…

ceph与hdfs的比较_分布式存储中HDFS与Ceph两者的区别是什么,各有什么优势?

过去两年&#xff0c;我的主要工作都在Hadoop这个技术栈中&#xff0c;而最近有幸接触到了Ceph。我觉得这是一件很幸运的事&#xff0c;让我有机会体验另一种大型分布式存储解决方案&#xff0c;可以对比出HDFS与Ceph这两种几乎完全不同的存储系统分别有哪些优缺点、适合哪些场…

使用带有OAuth的Spring Security保护资源

1.简介 在本教程中&#xff0c;我们将研究如何使用Spring Security和OAuth来基于路径模式&#xff08; / api / ** &#xff09;保护服务器上的管理资源。 我们配置的另一个路径模式&#xff08; / oauth / token &#xff09;将帮助已配置的授权服务器生成访问令牌。 请注意&a…

2080ti服务器支持什么系统,2080ti深度学习性能

AI开发平台ModelArtsModelArts是面向开发者的一站式AI开发平台&#xff0c;为机器学习与深度学习提供海量数据预处理及半自动化标注、大规模分布式Training、自动化模型生成&#xff0c;及端-边-云模型按需部署能力&#xff0c;帮助用户快速创建和部署模型&#xff0c;管理全周…

pythonset是什么类型的游戏_Python集合(set)类型的操作

原文详见&#xff1a;http://blog.csdn.net/business122/article/details/7541486python的set和其他语言类似, 是一个无序不重复元素集, 基本功能包括关系测试和消除重复元素. 集合对象还支持union(联合), intersection(交), difference(差)和sysmmetric difference(对称差集)等…

服务器选购seo优化规则,网站seo优化注意事项1—域名和服务器选择

原标题&#xff1a;网站seo优化注意事项1—域名和服务器选择一个网站的域名和服务器是很重要的&#xff0c;这关乎网站的优化前提&#xff0c;作为一个seo优化人员必须要对自己网站的域名和服务器有所了解。网站seo优化注意事项1&#xff1a;域名的关联性域名怎么选择呢&#x…

as cast float server sql_Sql Server中Float格式转换字符串varchar方法(转)

SELECT CONVERT(varchar(100), CAST(testFloat AS decimal(38,2)))SELECT STR(testFloat, 38, 2)从Excel中导入到sql2000&#xff0c;有一列“联系方式”变成了float类型&#xff0c;我想转换成nvarchar类型&#xff0c;用下面的语句select convert(nvarchar(30),convert(int,联…

openjpa_OpenJPA:内存泄漏案例研究

openjpa本文将提供完整的根本原因分析详细信息以及解决影响Oracle Weblogic Server 10.0生产环境的Java堆内存泄漏&#xff08;Apache OpenJPA泄漏&#xff09;的方法。 这篇文章还将演示在管理javax.persistence.EntityManagerFactory生命周期时遵循Java Persistence API最佳实…

美国凯斯西储大学计算机硕士专业怎么样,在凯斯西储大学读硕士大约需要多少花费?...

凯斯西储大学是美国著名大学之一&#xff0c;始建于1826年&#xff0c;坐落于俄亥俄州的克里夫兰&#xff0c;是一所以独立研究闻名的世界顶级私立大学&#xff0c;美国一级国家级大学。美国作为当今世界留学费用最高的国家之一&#xff0c;费用问题是所有赴美留学的学生都非常…

ffmpeg 缩放算法_图像尺寸调整算法介绍并手动实现近邻算法

课程更新完毕&#xff0c;后面会根据同学的意见维护此课程课程讲师&#xff1a;【夏曹俊】课程基于opencv**版3.2讲解&#xff0c;全部使用**的c接口&#xff0c;旧的接口会逐步被opencv抛弃&#xff0c;所以我们学习还是尽量学习新的接口&#xff0c;QT使用**的5.8版本&#x…

win7如何修改dns服务器地址,Win7系统DNS怎么设置?Win7系统DNS设置方法

Win7系统DNS怎么设置?众所周知&#xff0c;DNS地址是一个域名服务器地址&#xff0c;它可以把用户的网站地址解析成IP地址。如果这个服务器出现问题&#xff0c;可能就上不了网了。我们在使用Win7系统的时候&#xff0c;就是因为域名解析服务器不能将要访问的域名解析为正确的…

gradle web_简单的Gradle Web应用程序

gradle webGradle通过“ war”和“ jetty”插件轻松支持Java Web应用程序。 本教程将向您展示如何快速创建一个简单的Java Servlet Web应用程序。 我们的基本项目结构将是&#xff1a; basic-websrcmainjavacomcodetutrHelloWorldServlet.javawebappWEB-INFweb.xml首先&#x…

密钥文件登录服务器,密钥文件登录云服务器

密钥文件登录云服务器 内容精选换一换远程桌面协议(Remote Desktop Protocol&#xff0c;RDP)&#xff0c;是微软提供的多通道的远程登录协议。本节为您介绍如何使用RDP文件远程登录Windows弹性云服务器。从管理控制台下载的RDP文件对应唯一的云服务器&#xff0c;当前RDP文件命…

centos格式化优盘命令_u盘格式化工具fat32u盘安装centos6.5

3、声卡驱动不完全匹配&#xff0c;重装下声卡驱动(这种情况很少)朋友电脑开不了机&#xff0c;找我帮他看看&#xff0c;雨林木风win10家庭版&#xff0c;他电脑系统是有做过备份的现在做64位系统引起的■灵活设置背景2、电脑主板电池怎么换新的主板电池不贵的&#xff0c;才不…

maven 安装本地jar_使用Maven安装本地jar

maven 安装本地jar1.问题与选择 Maven是一种非常通用的工具&#xff0c;其可用的公共存储库首屈一指。 但是&#xff0c;总会有一个不在任何地方托管的工件&#xff0c;或者托管该工件的存储库有依赖的风险&#xff0c;因为当您需要它时&#xff0c;它可能不会启动。 发生这种情…

一个网站服务器有多少个ip,一个服务器可以有多少个ip地址

一个服务器可以有多少个ip地址 内容精选换一换华为云帮助中心&#xff0c;为用户提供产品简介、价格说明、购买指南、用户指南、API参考、最佳实践、常见问题、视频帮助等技术文档&#xff0c;帮助您快速上手使用华为云服务。会话保持&#xff0c;指负载均衡器可以识别客户与服…

redis集群scan_Redis scan命令的一次坑

Redis作为当前服务架构不可或缺的Cache&#xff0c;其支持丰富多样的数据结构&#xff0c;Redis在使用中其实也有很多坑&#xff0c;本次博主遇到的坑或许说是Java程序员会遇到的多一点&#xff0c;下面就听博主详细道来。线上服务堵塞String key keyOf(appid);int retryCount…

jbpm小项目测试_尝试使用jBPM Console NG(测试版)

jbpm小项目测试大家好&#xff01; 这是有关jBPM Console NG的另一篇文章。 经过6个月的辛苦工作&#xff0c;我很高兴为开发人员社区撰写这篇文章&#xff0c;以进行尝试。 在这篇文章中&#xff0c;我将解释如何从源代码构建应用程序。 这背后的主要思想是知道如何在测试过程…