【高斯消元】CDOJ1784 曜酱的线性代数课堂(二)
2017-10-05 01:30
281 查看
高斯消元求矩阵秩板子。
#include<cstdio> #include<cmath> #include<algorithm> #include<cstring> using namespace std; const double EPS=0.00000001; #define N 105 int n,m; double B ,A ; int guass_jordan() { memcpy(B,A,sizeof(A)); for(int i=1;i<=n;++i) { int pivot=i; for(int j=i+1;j<=n;++j) if(fabs(B[j][i])>fabs(B[pivot][i])) pivot=j; swap(B[i],B[pivot]); if(fabs(B[i][i])<EPS){ return i-1; } for(int j=i+1;j<=m;++j) B[i][j]/=B[i][i]; for(int j=1;j<=n;++j) if(i!=j) for(int k=i+1;k<=m;++k) B[j][k]-=B[j][i]*B[i][k]; } return n; } int main() { while(scanf("%d%d",&n,&m)!=EOF){ memset(A,0,sizeof(A)); for(int i=1;i<=n;++i){ for(int j=1;j<=m;++j){ scanf("%lf",&A[i][j]); } } printf("%d\n",guass_jordan()); } return 0; }
相关文章推荐
- 【高斯消元】CDOJ1783 曜酱的线性代数课堂(一)
- 【高斯消元】CDOJ1785 曜酱的线性代数课堂(三)
- CDOJ 1330 柱爷与远古法阵 概率DP,高斯消元,卡精度,极度恶心
- 关灯游戏 Lights out (三)(线性代数+高斯消元,搜索全部解)
- 【CF736D】Permutations 线性代数+高斯消元
- 【cdoj 1330】柱爷与远古法阵 高斯消元
- 【总结】线性基/(高斯消元 线性代数)
- CDOJ 1330 柱爷与远古法阵【高斯消元,卡精度】
- CDOJ 1330 柱爷与远古法阵(高斯消元)
- 线性代数(高斯消元):JSOI2008 球形空间产生器sphere
- [高斯消元 线性基 贪心 拟阵 逆元] BZOJ 4004 [JLOI2015]装备购买
- POJ 2345 Central heating 高斯消元
- [HDU 5755] Gambler Bo (高斯消元)
- POJ-1222-开关灯问题-高斯消元
- POJ 2947 高斯消元
- hdu 5755(高斯消元——模线性方程组模板)
- POJ 2947 Widget Factory(高斯消元)
- 混合颜料问题-网易的笔试题目(高斯消元找极大线性无关组)
- 高斯消元——混合
- 高斯消元