【洛谷P3389 【模板】高斯消元法】

这是个版子题,当然本蒟蒻也是看了好几天才明白

对于这样的线性方程组,我们可以看成是一个矩阵

对于百度百科给的定义(我感到很迷)赶脚和行列式有的一拼

但我们要注意的是:

行列式是一个确切的值(有关行列式求值,下片博文再说),可以看做是一个数值;

而矩阵不同,矩阵是一个数表,无特殊的值。当然,性质也略有不同(这里就先不细讲了)

高斯消元的核心就是把这个n*n的矩阵消成一个对角矩阵(就是除了a【i】【i】为1以外,其余全为零但是不要乱划a【i】【n+1】这是我们要输出的东西)然后就ok了

其中一些性质可以通过度娘或下面的(李昊的)过程代码验证一下

//这不是AC代码,只是一个检测的(黈力的小心一点啊QAQ)
#include<iostream> #include<algorithm> #include<cstdio> #include<cstring> #include<cmath> #include<cstdlib> using namespace std; typedef long long ll; typedef long double ld; typedef pair<int,int> pr; const double pi=acos(-1); #define rep(i,a,n) for(int i=a;i<=n;i++) #define per(i,n,a) for(int i=n;i>=a;i--) #define Rep(i,u) for(int i=head[u];i;i=Next[i]) #define clr(a) memset(a,0,sizeof a) #define pb push_back #define mp make_pair #define fi first #define sc second ld eps=1e-9; ll pp=1000000007; ll mo(ll a,ll pp){if(a>=0 && a<pp)return a;a%=pp;if(a<0)a+=pp;return a;} ll powmod(ll a,ll b,ll pp){ll ans=1;for(;b;b>>=1,a=mo(a*a,pp))if(b&1)ans=mo(ans*a,pp);return ans;} ll read(){ ll ans=0; char last=' ',ch=getchar(); while(ch<'0' || ch>'9')last=ch,ch=getchar(); while(ch>='0' && ch<='9')ans=ans*10+ch-'0',ch=getchar(); if(last=='-')ans=-ans; return ans; } //head int n,m; double a[100][100]; bool check(int k){ if(fabs(a[k][n+1])<eps)return 1; rep(i,1,n) if(fabs(a[k][i])>eps)return 1; return 0; } int main(){ freopen("1.txt","r",stdin); n=read();m=read(); // a_i,1 a_i,2 ... a_i,n a_i,n+1 rep(i,1,m) rep(j,1,n+1)a[i][j]=read(); rep(j,1,m){ rep(k,1,n+1)cout<<a[j][k]<<" "; puts(""); } int flag=0; rep(i,1,n){ int t=i; while(a[t][i]==0 && t<=n)t+=1; if(t==n+1){ flag=1; continue; } rep(j,1,n+1)swap(a[i][j],a[t][j]); double kk=a[i][i]; rep(j,1,n+1)a[i][j]/=kk; rep(j,1,m) if(i!=j){ double kk=a[j][i]; rep(k,1,n+1) a[j][k]-=kk*a[i][k]; } puts("------------"); rep(j,1,m){ rep(k,1,n+1)cout<<a[j][k]<<" "; puts(""); } } if(flag){ rep(i,1,m) if(!check(i)){ printf("No solution\n"); return 0; } printf("So many solutions\n"); } }

通过这个代码,可以看到化简的每一步过程

下面才是真正的AC代码:

//真正的AC代码
#include<iostream> #include<cstdio> #include<cmath> using namespace std; int n,s; double a[1001][1001]; int main() { cin>>n; for(int i=1;i<=n;i++)//输入矩阵 for(int j=1;j<=n+1;j++)//别忘记最后那一列   cin>>a[i][j]; for(int i=1;i<=n;i++) { s=i; while(a[s][i]==0&&s<=n) s++; if(s==n+1) { cout<<"No Solution";//判是否有解(若有一行或一列全为0,则无解) return 0; } for(int j=1;j<=n+1;j++) swap(a[i][j],a[s][j]); double k=a[i][i]; for(int j=1;j<=n+1;j++) a[i][j]=a[i][j]/k; //为防止a【i】【i】为零,可以换一下 for(int j=1;j<=n;j++) { if(i!=j) { double ki=a[j][i]; for(int m=1;m<=n+1;m++) a[j][m]=a[j][m]-ki*a[i][m];//这就保证可以将除了a【i】【i】的全部为0 } } } for(int i=1;i<=n;i++) printf("%.2lf\n",a[i][n+1]);//输出ouo return 0; }

这其中用到了一些神奇的性质,还是提前掌握为妙

转载于:https://www.cnblogs.com/gongcheng456/p/10705246.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值