代码语言
.
CSharp
.
JS
Java
Asp.Net
C
MSSQL
PHP
Css
PLSQL
Python
Shell
EBS
ASP
Perl
ObjC
VB.Net
VBS
MYSQL
GO
Delphi
AS
DB2
Domino
Rails
ActionScript
Scala
代码分类
文件
系统
字符串
数据库
网络相关
图形/GUI
多媒体
算法
游戏
Jquery
Extjs
Android
HTML5
菜单
网页交互
WinForm
控件
企业应用
安全与加密
脚本/批处理
开放平台
其它
【
C/C++
】
矩阵
作者:
睿智
/ 发布于
2016/4/21
/
690
define MAXN 100 #define fabs(x) ((x)>0?(x):-(x)) #define zero(x) (fabs(x)<1e-10) struct mat{ int n,m; double data[MAXN][MAXN]; }; int mul(mat& c,const mat& a,const mat& b){ int i,j,k; if (a.m!=b.n) return 0; c.n=a.n,c.m=b.m; for (i=0;i<c.n;i++) for (j=0;j<c.m;j++) for (c.data[i][j]=k=0;k<a.m;k++) c.data[i][j]+=a.data[i][k]*b.data[k][j]; return 1; } int inv(mat& a){ int i,j,k,is[MAXN],js[MAXN]; double t; if (a.n!=a.m) return 0; for (k=0;k<a.n;k++){ for (t=0,i=k;i<a.n;i++) for (j=k;j<a.n;j++) if (fabs(a.data[i][j])>t) t=fabs(a.data[is[k]=i][js[k]=j]); if (zero(t)) return 0; if (is[k]!=k) for (j=0;j<a.n;j++) t=a.data[k][j],a.data[k][j]=a.data[is[k]][j],a.data[is[k]][j]=t; if (js[k]!=k) for (i=0;i<a.n;i++) t=a.data[i][k],a.data[i][k]=a.data[i][js[k]],a.data[i][js[k]]=t; a.data[k][k]=1/a.data[k][k]; for (j=0;j<a.n;j++) if (j!=k) a.data[k][j]*=a.data[k][k]; for (i=0;i<a.n;i++) if (i!=k) for (j=0;j<a.n;j++) if (j!=k) a.data[i][j]-=a.data[i][k]*a.data[k][j]; for (i=0;i<a.n;i++) if (i!=k) a.data[i][k]*=-a.data[k][k]; } for (k=a.n-1;k>=0;k--){ for (j=0;j<a.n;j++) if (js[k]!=k) t=a.data[k][j],a.data[k][j]=a.data[js[k]][j],a.data[js[k]][j]=t; for (i=0;i<a.n;i++) if (is[k]!=k) t=a.data[i][k],a.data[i][k]=a.data[i][is[k]],a.data[i][is[k]]=t; } return 1; } double det(const mat& a){ int i,j,k,sign=0; double b[MAXN][MAXN],ret=1,t; if (a.n!=a.m) return 0; for (i=0;i<a.n;i++) for (j=0;j<a.m;j++) b[i][j]=a.data[i][j]; for (i=0;i<a.n;i++){ if (zero(b[i][i])){ for (j=i+1;j<a.n;j++) if (!zero(b[j][i])) break; if (j==a.n) return 0; for (k=i;k<a.n;k++) t=b[i][k],b[i][k]=b[j][k],b[j][k]=t; sign++; } ret*=b[i][i]; for (k=i+1;k<a.n;k++) b[i][k]/=b[i][i]; for (j=i+1;j<a.n;j++) for (k=i+1;k<a.n;k++) b[j][k]-=b[j][i]*b[i][k]; } if (sign&1) ret=-ret; return ret; }
试试其它关键字
矩阵
同语言下
.
C分鱼问题
.
链表
.
最大连续和
.
编码字符串
.
libiconv字符编码处理及判断字符串是否为utf8
.
一组数中两两二元组,差最大有几对,差最小呢?(数组
.
通过管道获取一个进程的执行状态
.
多关键字排序
.
字符串字典序排序
.
3元一次方程(牛顿迭代法求方程的根)
可能有用的
.
C#实现的html内容截取
.
List 切割成几份 工具类
.
SQL查询 多列合并成一行用逗号隔开
.
一行一行读取txt的内容
.
C#动态修改文件夹名称(FSO实现,不移动文件)
.
c# 移动文件或文件夹
.
c#图片添加水印
.
Java PDF转换成图片并输出给前台展示
.
网站后台修改图片尺寸代码
.
处理大图片在缩略图时的展示
睿智
贡献的其它代码
(
12
)
.
实现Oracle数据库备份
.
批量生成手机号码
.
根据当前文字选择返回被选中的文字
.
限制Textarea文本域字符个数
.
命令行CMD颜色设置
.
一个简单的Doctor类的封装
.
矩阵
.
获取网站缩略图
.
SQL查询一张表插入到另一张表
.
标签自适应高度和宽度
Copyright © 2004 - 2024 dezai.cn. All Rights Reserved
站长博客
粤ICP备13059550号-3