strassen算法java_使用java写的矩阵乘法实例(Strassen算法)_Java_软件编程
Strassen算法于1969年由德国数学家Strassen提出,该方法引入七个中间变量,每个中间变量都只需要进行一次乘法运算。而朴素算法却需要进行8次乘法运算。
原理
Strassen算法的原理如下所示,使用sympy验证Strassen算法的正确性
import sympy as s
A = s.Symbol("A")
B = s.Symbol("B")
C = s.Symbol("C")
D = s.Symbol("D")
E = s.Symbol("E")
F = s.Symbol("F")
G = s.Symbol("G")
H = s.Symbol("H")
p1 = A * (F - H)
p2 = (A + B) * H
p3 = (C + D) * E
p4 = D * (G - E)
p5 = (A + D) * (E + H)
p6 = (B - D) * (G + H)
p7 = (A - C) * (E + F)
print(A * E + B * G, (p5 + p4 - p2 + p6).simplify())
print(A * F + B * H, (p1 + p2).simplify())
print(C * E + D * G, (p3 + p4).simplify())
print(C * F + D * H, (p1 + p5 - p3 - p7).simplify())
复杂度分析
$$f(N)=7\times f(\frac{N}{2})=7^2\times f(\frac{N}{4})=...=7^k\times f(\frac{N}{2^k})$$
最终复杂度为$7^{log_2 N}=N^{log_2 7}$
java矩阵乘法(Strassen算法)
代码如下,可以看看数据结构的定义,时间换空间。
public class Matrix {
private final Matrix[] _matrixArray;
private final int n;
private int element;
public Matrix(int n) {
this.n = n;
if (n != 1) {
this._matrixArray = new Matrix[4];
for (int i = 0; i < 4; i++) {
this._matrixArray[i] = new Matrix(n / 2);
}
} else {
this._matrixArray = null;
}
}
private Matrix(int n, boolean needInit) {
this.n = n;
if (n != 1) {
this._matrixArray = new Matrix[4];
} else {
this._matrixArray = null;
}
}
public void set(int i, int j, int a) {
if (n == 1) {
element = a;
} else {
int size = n / 2;
this._matrixArray[(i / size) * 2 + (j / size)].set(i % size, j % size, a);
}
}
public Matrix multi(Matrix m) {
Matrix result = null;
if (n == 1) {
result = new Matrix(1);
result.set(0, 0, (element * m.element));
} else {
result = new Matrix(n, false);
result._matrixArray[0] = P5(m).add(P4(m)).minus(P2(m)).add(P6(m));
result._matrixArray[1] = P1(m).add(P2(m));
result._matrixArray[2] = P3(m).add(P4(m));
result._matrixArray[3] = P5(m).add(P1(m)).minus(P3(m)).minus(P7(m));
}
return result;
}
public Matrix add(Matrix m) {
Matrix result = null;
if (n == 1) {
result = new Matrix(1);
result.set(0, 0, (element + m.element));
} else {
result = new Matrix(n, false);
result._matrixArray[0] = this._matrixArray[0].add(m._matrixArray[0]);
result._matrixArray[1] = this._matrixArray[1].add(m._matrixArray[1]);
result._matrixArray[2] = this._matrixArray[2].add(m._matrixArray[2]);
result._matrixArray[3] = this._matrixArray[3].add(m._matrixArray[3]);;
}
return result;
}
public Matrix minus(Matrix m) {
Matrix result = null;
if (n == 1) {
result = new Matrix(1);
result.set(0, 0, (element - m.element));
} else {
result = new Matrix(n, false);
result._matrixArray[0] = this._matrixArray[0].minus(m._matrixArray[0]);
result._matrixArray[1] = this._matrixArray[1].minus(m._matrixArray[1]);
result._matrixArray[2] = this._matrixArray[2].minus(m._matrixArray[2]);
result._matrixArray[3] = this._matrixArray[3].minus(m._matrixArray[3]);;
}
return result;
}
protected Matrix P1(Matrix m) {
return _matrixArray[0].multi(m._matrixArray[1]).minus(_matrixArray[0].multi(m._matrixArray[3]));
}
protected Matrix P2(Matrix m) {
return _matrixArray[0].multi(m._matrixArray[3]).add(_matrixArray[1].multi(m._matrixArray[3]));
}
protected Matrix P3(Matrix m) {
return _matrixArray[2].multi(m._matrixArray[0]).add(_matrixArray[3].multi(m._matrixArray[0]));
}
protected Matrix P4(Matrix m) {
return _matrixArray[3].multi(m._matrixArray[2]).minus(_matrixArray[3].multi(m._matrixArray[0]));
}
protected Matrix P5(Matrix m) {
return (_matrixArray[0].add(_matrixArray[3])).multi(m._matrixArray[0].add(m._matrixArray[3]));
}
protected Matrix P6(Matrix m) {
return (_matrixArray[1].minus(_matrixArray[3])).multi(m._matrixArray[2].add(m._matrixArray[3]));
}
protected Matrix P7(Matrix m) {
return (_matrixArray[0].minus(_matrixArray[2])).multi(m._matrixArray[0].add(m._matrixArray[1]));
}
public int get(int i, int j) {
if (n == 1) {
return element;
} else {
int size = n / 2;
return this._matrixArray[(i / size) * 2 + (j / size)].get(i % size, j % size);
}
}
public void display() {
for (int i = 0; i < n; i++) {
for (int j = 0; j < n; j++) {
System.out.print(get(i, j));
System.out.print(" ");
}
System.out.println();
}
}
public static void main(String[] args) {
Matrix m = new Matrix(2);
Matrix n = new Matrix(2);
m.set(0, 0, 1);
m.set(0, 1, 3);
m.set(1, 0, 5);
m.set(1, 1, 7);
n.set(0, 0, 8);
n.set(0, 1, 4);
n.set(1, 0, 6);
n.set(1, 1, 2);
Matrix res = m.multi(n);
res.display();
}
}
总结
到此这篇关于使用java写的矩阵乘法的文章就介绍到这了,更多相关java矩阵乘法(Strassen算法)内容请搜索我们以前的文章或继续浏览下面的相关文章希望大家以后多多支持我们!
郑重声明:本文版权归原作者所有,转载文章仅为传播更多信息之目的,如作者信息标记有误,请第一时间联系我们修改或删除,多谢。
strassen算法java_使用java写的矩阵乘法实例(Strassen算法)_Java_软件编程相关推荐
- 矩阵生发java_使用java写的矩阵乘法实例(Strassen算法)
Strassen算法于1969年由德国数学家Strassen提出,该方法引入七个中间变量,每个中间变量都只需要进行一次乘法运算.而朴素算法却需要进行8次乘法运算. 原理 Strassen算法的原理如下 ...
- java矩阵连乘算法_使用java写的矩阵乘法实例(Strassen算法)
Strassen算法于1969年由德国数学家Strassen提出,该方法引入七个中间变量,每个中间变量都只需要进行一次乘法运算.而朴素算法却需要进行8次乘法运算. 原理 Strassen算法的原理如下 ...
- 《算法导论》——矩阵乘法的Strassen算法
前言: 很多朋友看到我写的<算法导论>系列,可能会觉得云里雾里,不知所云.这里我再次说明,本系列博文时配合<算法导论>一书,给出该书涉及的算法的c++实现.请结合<算法导 ...
- 矩阵相乘的strassen算法_矩阵乘法的Strassen算法+动态规划算法(矩阵链相乘和硬币问题)...
矩阵乘法的Strassen 这个算法就是在矩阵乘法中采用分治法,能够有效的提高算法的效率. 先来看看咱们在高等代数中学的普通矩阵的乘法 两个矩阵相乘 上边这种普通求解方法的复杂度为: O(n3) 也称 ...
- 矩阵相乘的strassen算法_4-2.矩阵乘法的Strassen算法详解
题目描述 请编程实现矩阵乘法,并考虑当矩阵规模较大时的优化方法. 思路分析 根据wikipedia上的介绍:两个矩阵的乘法仅当第一个矩阵B的列数和另一个矩阵A的行数相等时才能定义.如A是m×n矩阵和B ...
- 循环取矩阵的某行_1.2 震惊! 某大二本科生写的矩阵乘法吊打Mathematica-线性代数库BLAS-矩阵 (上)...
本文是 1. 线性代数库BLASzhuanlan.zhihu.com 系列的第二篇, 将讲述矩阵类的结构和矩阵基础运算的AVX2加速算法. 1. 矩阵类的结构 在讲述矩阵各种算法之前很有必要详解一下 ...
- 稀疏矩阵加法运算_1.2 震惊! 某大二本科生写的矩阵乘法吊打Mathematica-线性代数库BLAS-矩阵 (上)...
本文是 1. 线性代数库BLASzhuanlan.zhihu.com 系列的第二篇, 将讲述矩阵类的结构和矩阵基础运算的AVX2加速算法. 1. 矩阵类的结构 在讲述矩阵各种算法之前很有必要详解一下 ...
- java编写螺旋矩阵讲解_Java如何实现螺旋矩阵 Java实现螺旋矩阵代码实例
本篇文章小编给大家分享一下Java实现螺旋矩阵代码实例,小编觉得挺不错的,现在分享给大家供大家参考,有需要的小伙伴们可以来看看. 给定一个包含 m x n 个元素的矩阵(m 行, n 列),请按照顺时 ...
- python矩阵乘法分治算法_矩阵乘法的Strassen算法详解 --(算法导论分治法求矩阵)...
1 题目描述 2 思路分析 3 解法 4 小结 1 题目描述 请编程实现矩阵乘法,并考虑当矩阵规模较大时的优化方法. 2 思路分析 根据wikipedia上的介绍:两个矩阵的乘法仅当第一个矩阵B的列数 ...
最新文章
- java 读取流的字符编码格式_如何使用Java代码获取文件、文件流或字符串的编码方式...
- Java中“==”的使用,以及“==”和equal的比较
- c#.net课程设计:ZCMU通讯录(待更新)
- java B2B2C电子商务平台分析之十五-----EureKa服务注册与发现
- Hibernate的配置详解
- c# 无法加载oraops.dll_Robotstudio软件二次开发:基于C#语言的Smart组件开发基础
- 利用JOrgChart只需2分钟即可配置简单组织机构图
- VPP /什么是VPP?读这篇文章就够了
- 手机通讯录备份代码实现二
- 从零开始学前端:CSS复合选择器 --- 今天你学习了吗?(CSS:Day10)
- java tlv协议_看懂通信协议:自定义通信协议设计之TLV编码应用
- Deep Learning Notes: Chapter 1 Introduction
- ykhmi是什么触摸屏软件_YKHMI 中达优控7寸一体机兼容台达ES2
- java读取word目录
- 提高DevExpress winform 控件加载速度
- 第三方登录/分享最佳实践
- 挑战大数据 金仓助力三农自助支付系统
- 工欲善其事,必先利其器 - 修炼Android Studio操作技巧
- 如何利用大数据软件进行家政服务的选址要素收集与分析
- Charles 的简单使用
热门文章
- oracle使用nbu的必要,nbu恢复oracle,使用与源备份端不同用户问题解决
- html实现验证码效果,js实现验证码功能
- 解决pip无法安装bayes-opt报错:ERROR: Could not find a version that satisfies the requirement bayes-opt
- jetbrains(IDEA/Pycharm)修改断点的颜色与调试时debugger的颜色
- linux网卡主备,linux网卡bounding的主备模式下上层路由端需要什么设置?
- @RequestParam用法与@PathVariable用法的区别
- 测绘专业学生对测绘软件的一点认识
- ActiveMQ消息队列的使用
- java源码阅读Object
- Python3连接MySQL