題目大意

如題，給你一個n行m列的矩陣,求它的秩。

樣例輸入

3 5
1 2 3 4 5
2 4 6 8 10
3 6 9 12 15

3 3
2 0 2
0 2 0
2 0 2

樣例輸出

1
2

代碼如下

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#include<map>
#include<vector>
#include<cmath>
#define CLR(x) memset(x,0,sizeof(x))
#define LL long long
using namespace std;

double mat[300][300];
int r,c;

void debug(){
    printf(" \n debug : \n");
    for(int i = 1 ; i <= r ; i++){
        for(int j = 1 ; j <= c ;j++){
            printf("%.0lf  ",mat[i][j]);
        }
        printf("\n");
    }
}
                                                      //保證一定得精度。
int cmp(double x,double y){
    double v = x - y;
    if(v > 1e-1) return 1;
    if(v < -1e-1) return -1;
    return 0;
}

void subrow(int r1 , int r2,double temp){
    for(int i = 1 ; i <= c ; i++){
        mat[r1][i] -= mat[r2][i]*temp;
    }
}

void swaprow(int r1 , int r2){
    for(int i = 1 ; i <= c ; i++){
        swap(mat[r1][i],mat[r2][i]);
    }
}

void solve(){
    for(int i = 1 ; i <= r ; i++){
        //debug();
        for(int cal = i ; cal <= c ; cal++){            //對這 i 行的每一列來找，要是這列本身和之下的列全是0，則找 i 行下一列。
            bool flag = true;
            
            if(cmp(mat[i][cal],0) == 0)
            {
                flag = false;                          //如果第 i 行 cal 列這個位置是 0 ，那找找這列下面的行有沒有
                for(int j = i+1 ; j <= r ; j++){       //一個不為 0 的數(shù) , 要是有 ， 就將它與 i 行交換。
                    int tmp = cmp(mat[j][cal],0);
                    if(tmp == 1 || tmp == -1){
                        flag = true;
                        swaprow(j,i);
                        break;
                    }
                }
            }
            
            
            if(!flag) continue;                        //如果這列全是 0 ， 到下一列。
            
            
            int v = i ;                                 //如果發(fā)現(xiàn)這列有值不為 0 并把它跟 i 行交換后，
            int maxn = mat[i][cal];                     //就再找這列有沒有比 i 行 這個位置的值更大的值，如果有，將那一行
            for(int j = i+1 ; j <= r ; j++){            //跟i行交換。（聽說是為了減小誤差）
                if(cmp(mat[j][cal],mat[i][cal]) == 1){
                    v = j; maxn = mat[j][cal];
                }
            }
            if(v!=i){
                swaprow(i,v);
                //debug();
            }
            
                                                        //相減。（全程確保精度）
            for(int j = 1 ; j <= r ; j++){
                if(j == i) continue;
                if(cmp(mat[i][cal],0) == 0) continue;
                double tmp = mat[j][cal] / mat[i][cal];
                subrow(j,i,tmp);
            }
            break;
        }
    }
}

int main()
{
    CLR(mat);
    scanf("%d%d",&r,&c);
    for(int i = 1 ; i <= r ; i++){
        for(int j = 1 ; j <= c ; j++){
            scanf("%lf",&mat[i][j]);
        }
    }
    solve();
    //debug();

    bool f = false;
    int ans = 0;
    for(int i = r ; i >= 1 ; i--){
        for(int j = 1 ; j <= c ; j++){
            int tmp = cmp(mat[i][j],0);
            if(tmp == 1 || tmp == -1){
                f = true;break;
            }
        }
        if(f){
            ans = i;break;
        }
    }

    printf("%d\n",ans);
}