2015.4.7
来源:互联网 发布:淘宝店网址手机打开 编辑:程序博客网 时间:2024/06/05 21:11
今天看了算法导论上的strassen算法,懂了个大概。下面是该算法的java实现:
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();
}
}
感觉程序还是没有读太懂,过几天再看看。
- 2015.4.7
- 【2015.4.7】心情贴
- OC笔记 - NSThread(2015.4.7)
- 2015.4.7-4.8 widget中使用Listview (上)分析说明
- 2015.4.1
- 2015.4感悟
- 2015.4.16
- 开博首文(2015.4.4)
- 2015.4.9学习mybatis
- 2015.4.10学习mybatis
- 2015.4.15友元类时间
- 2015.4.18博客说明
- 移动设计-2015.4.21
- 2015.4.23问题集
- 2015.4-5笔记
- 2015.4.8 虚拟存储器
- 2015.4.18 Bernstein条件
- Xilinx petalinux 2015.4 试用
- JSP包含其他页面的三种方式
- TLD源码解析-tldGenerateFeatures
- 手机卫士-05
- 如何去除MFC窗口的自动记忆功能 (如去除工具栏,菜单栏上次运行程序自动保存的内容)
- openfire安装 Unknown character set: 'utf8mb4'错误
- 2015.4.7
- 小马哥----高仿三星s5 主板E118型号 6572芯片拆机刷机图与开机界面图
- web 锚点定位
- 致:同年代的童真童鞋们
- 程序员要学会问问题!
- 去除JSP页面自动生成的空行
- 【转发】关于Java性能的9个谬论
- NYOJ 104 最大和(DP)
- 如何开启webview的LBS功能