SPOJ Highways [矩阵树定理]

最新推荐文章于 2022-12-14 17:32:27 发布

转载最新推荐文章于 2022-12-14 17:32:27 发布 · 76 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/candy99/p/6421388.html

本文介绍了一个使用高斯消元法求解线性方程组的C++程序实现，该程序能够处理特定类型的矩阵方程，并通过消除法找到方程组的解。程序通过读取输入文件中的数据来构造矩阵，然后应用高斯消元法进行计算。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

裸题

注意：

1.消元时判断系数为0，退出

2.最后乘ans要用double....

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long ll;
const int N=15;
double eps=1e-8;
inline int read(){
    char c=getchar();int x=0;
    while(c<'0'||c>'9'){c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x;
}
int n,m,u,v;
double a[N][N];
void Gauss(){
    n--;
    for(int i=1;i<=n;i++){
        int r=i;
        for(int j=i+1;j<=n;j++) if(abs(a[j][i])>abs(a[r][i])) r=j;
        if(abs(a[r][i])<eps){puts("0");return;}
        if(r!=i) for(int k=1;k<=n;k++) swap(a[r][k],a[i][k]);
        for(int j=i+1;j<=n;j++){
            double t=a[j][i]/a[i][i];
            for(int k=i;k<=n;k++) a[j][k]-=t*a[i][k];
        }
    }
    double ans=1;
    for(int i=1;i<=n;i++) ans*=a[i][i];
    printf("%.0f\n",abs(ans));
}
int main(){
    freopen("in","r",stdin);
    int T=read();
    while(T--){
        n=read();m=read();
        memset(a,0,sizeof(a));
        for(int i=1;i<=m;i++){
            u=read();v=read();
            a[u][u]++;a[v][v]++;
            a[u][v]--;a[v][u]--;
        }
        Gauss();
    }
}