100分求矩阵运算的C++库

2004-12-30 08:45:59 391 阅读 AI 写同款 · GEO 优化AI 选题数据监测一键分发 ›

c++求最大子矩阵和问题最大子矩阵和问题。给定m行n列的整数矩阵A，求矩阵A的一个子矩阵，使其元素之和最大。输入格式: 第一行输入矩阵行数m和列数n（1≤m≤100，1≤n≤100），再依次输入m×n个整数。输出格式: 输出第一行为最大子矩阵各元素之和，第二行为子矩阵在整个矩阵中行序号范围与列序号范围。输入样例1: 5 6 60 3 -65 -92 32 -70 -41 14 -38 54 2 29 69 88 ... 阅读详情

哪里有免费的C++矩阵运算的库，可以在VC6.0下编译通过的？

...全文

矩阵的c++实现在大学数学课程《线性代数》中，就有矩阵和行列式的出现，这篇文章主要讲矩阵在c++中的实现和一些用途。此篇文章只写c++的内容，不具体写到数学中矩阵的一些公式、性质。阅读详情

oldboy1234

C++ 语言 65,214 成员 250,514 内容

+ 关注 + 加入社区

切换为时间正序

zhengwei1984222 2004-12-31

打赏
举报

我想知道为什么不用matlab
可以生成c代码的

oldboy1234 2004-12-30

打赏
举报

刚才查了Blitz的官方网站，不支持VC6.0啊
Visual C++ 6.0. Does not support member templates outside of class declarations, and crashes (CRASHES!) if you try to do so. Definitely not recommended.

xuzheng318 2004-12-30

打赏
举报

http://www.osl.iu.edu/research/mtl/
The Matrix Template Library, 比Fortran 90还快的矩阵运算程序库

http://www.netlib.org

http://math.nist.gov/tnt

goodluckyxl 2004-12-30

打赏
举报

矩阵的计算用STL可以很好的实现
效率可以很高
但是如果使用编译期进行求值
可能效率会很差很差

求特征根，秩，逆这些涉及的相对复杂
自己实现考虑异常不小的工作量
-_-#

oldboy1234 2004-12-30

打赏
举报

Blitz库我有了，不过感觉好复杂，就是想找一个简单点的库，支持基本运算即可

oldboy1234 2004-12-30

打赏
举报

谢谢楼上的，不过我需要一些稍微复杂的运算，比如求矩阵的特征根，求矩阵的逆这些。效率上不需要考虑多少，计算量不是很大。我本来找了个metaware，是linux下的，可惜库中用了模版的偏特化，而VC6.0又不支持，郁闷。

BluntBlade 2004-12-30

打赏
举报

使用Blitz库吧。

goodluckyxl 2004-12-30

打赏
举报

//以前写的部分矩阵，对于矩阵运算应该采用延迟策略来的比较好一些
//自己考虑考虑
#include <stdio.h>
#include <string.h>
#include <iostream.h>
#include <assert.h>
class Matrix
{
public:
Matrix();
Matrix( int RowRecord, int ColRecord );
~Matrix();

Matrix& operator = ( Matrix& OneMatrix );
Matrix& operator +=( Matrix& OneMatrix );
friend ostream& operator << ( ostream& os, Matrix& TheMatrix );

private:
void InitMatrix();
void DestoryMatrix();
private:
int RowLength;
int ColLength;
int** DataStore;
};

Matrix::Matrix()
{
RowLength = 0;
ColLength = 0;
DataStore = 0;
}

Matrix::Matrix( int RowRecord, int ColRecord )
{
assert( RowRecord != 0 && ColRecord != 0 );
RowLength = RowRecord;
ColLength = ColRecord;
InitMatrix();
}

void Matrix::InitMatrix()
{
DataStore = new int*[RowLength];
for( int Num = 0; Num < RowLength; Num++ )
{
try
{
*(DataStore+Num) = new int[ColLength];
memcpy( DataStore+Num, 0, ColLength*sizeof(int) );
}
catch (...)
{
for( int i = Num-1; i >= 0; i--)
delete[] *(DataStore+i);

delete[] DataStore;
}
}

}
void Matrix::DestoryMatrix()
{
if( ColLength*RowLength != 0 )
{

for( int i = ColLength-1; i >= 0; i-- )
{
delete[] *(DataStore+i);
}
delete[] DataStore;

}
}

Matrix::~Matrix()
{
DestoryMatrix();
}

Matrix& Matrix::operator += (Matrix& OneMatrix )
{
if( RowLength != OneMatrix.RowLength || ColLength != OneMatrix.ColLength )
{
return *this;
}
else
{
for( int i = 0; i < RowLength; i++ )
{
for( int j = 0; j < ColLength; j++ )
{
*(*(DataStore+i)+j) += *(*(OneMatrix.DataStore+i)+j);
}
}
return *this;
}
}

Matrix& Matrix::operator = ( Matrix& OneMatrix )
{
DestoryMatrix();
RowLength = OneMatrix.RowLength;
ColLength = OneMatrix.ColLength;

InitMatrix();
*this += OneMatrix;
return *this;
}

ostream& operator << ( ostream& os, Matrix& TheMatrix )
{
for( int i = 0; i < TheMatrix.RowLength; i++ )
{
for( int j = 0; j < TheMatrix.ColLength; j++ )
{
os<<*(*(TheMatrix.DataStore+i)+j);
}
os<<endl;
}
return os;
}

100分求矩阵运算的C++库

相关推荐