[hihoCoder] 高斯消元·一 [TPLY]
2024-10-10 20:08:29
高斯消元一
题目链接 : http://hihocoder.com/problemset/problem/1195?sid=1269842
很"好aoaoaoaoaoaoa"的高斯消元模板题
题意
多个方程组,要求:输出解、判无解、判多解
保证方程解非负
做法
第一点注意
首先要会高斯消元(废话)
然后还要卡精度
所以一定要用eps卡精度
这是第一点
第二点注意
然后就是最恶心的:判无解多解
我们先明确一点:无解优先级大于多解
什么意思呢?
对于一个方程,我们要先检查无解,再检查多解
明白这一点之后,咱们继续。
第三点注意
无解怎么判断?
如果 方程系数都等于0并且结果大于 0 则无解(因为题目保证解是非负数)
(形如 x × 0 = y , x × -1 = y | x>0 && y>0,肯定无解)
第四点注意
如何判断多解?
如果在消元过程中,某一列都被消成了0,并且保证该方程有解,那么这个方程是多解的。
为什么呢?
因为如果一个未知数上的系数都是0,那么这个未知数有无穷多种取法,所以方程就有多组解了。
第五点注意
如果你发现有多解,但是不确定是不是无解怎么办?
如果在最后用倒三角求未知数的值时
我们求到一个未知数的系数为0
但是它的值不为0的时候
那么它就是无解的
反之就是多解的
总结
综上所述 这道题是"不折不扣"的"模板题"
代码
#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <queue>
#include <cmath>
#define rg register int
#define ll long long
#define RG register
#define il inline
using namespace std;
il int gi()
{
rg x=0,o=0;RG char ch=getchar();
while(ch!='-'&&(ch<'0'||'9'<ch)) ch=getchar();
if(ch=='-') o=1,ch=getchar();
while('0'<=ch&&ch<='9') x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
return o?-x:x;
}
int n,m;
#define db double
db x[1001],fc[1001][501];
il int gauss()
{
RG bool flg=0;
for(rg now,k=1; k<=n; ++k)
{
now=k;
for(rg i=k+1; i<=m; ++i)
if(fabs(fc[i][k])>fabs(fc[now][k]))
now=i;
if(now==k && fabs(fc[k][k])<1e-7)
{
flg=1;
continue;
}
if(now!=k) swap(fc[now],fc[k]);
for(rg i=k+1; i<=n+1; ++i) fc[k][i]/=fc[k][k];fc[k][k]=1;
for(rg i=k+1; i<=m; ++i)
{
for(rg j=k+1; j<=n+1; ++j)
fc[i][j]-=fc[i][k]*fc[k][j];
fc[i][k]=0;
}
}
for(rg j,i=1; i<=m; ++i)
{
for(j=1; j<=n; ++j)
if(fabs(fc[i][j])>1e-6)
break;
if(j==n+1 && fabs(fc[i][n+1])>1e-6) return 0; // 如果 方程系数小于0并且结果大于 0 则无解
}
for(rg i=n; i>=1; --i)
{
for(rg j=i+1; j<=n; ++j) fc[i][n+1]-=fc[i][j]*fc[j][n+1];
if(fc[i][n+1] && !fc[i][i] ) return 0;
}
if(flg) return -1;
return 1;
}
int main()
{
n=gi(),m=gi();
for(rg i=1; i<=m; ++i)
for(rg j=1; j<=n+1; ++j)
scanf("%lf",&fc[i][j]);
rg ans=gauss();
if(ans==-1) puts("Many solutions");
else if(ans==0) puts("No solutions");
else for(rg i=1; i<=n; ++i) printf("%d\n",(int)(fc[i][n+1]+0.5));
return 0;
}
最新文章
- 当忘记mysql数据库密码时如何进行修改
- Java提高篇——设计模式
- Mob.com 短信验证的简单使用
- Java常用正则表达式验证工具类RegexUtils.java
- centOS6.6升级gcc4.8
- TCP经受时延的ACK
- Gradle Tips#1-tasks
- jquery用div模拟一个下拉列表框
- openstack操作之一 命令行
- 使用xUnit为.net core程序进行单元测试(4)
- TCP发送源码学习(3)--tcp_transmit_skb
- 学习-xlsxwriter模块
- 控制结构(7): 程序计数器(PC)
- 【原创】大叔问题定位分享(15)spark写parquet数据报错ParquetEncodingException: empty fields are illegal, the field should be ommited completely instead
- Chrome 调试技巧
- 网站开发,推荐使用SuperSlide 插件-Tab标签切换,图片滚动,无缝滚动,焦点图
- ERP系统知识笔记
- 【转】系统去掉 Android 4.4.2 的StatusBar和NavigationBar
- Java中-classpath和路径的使用
- java获取屏幕密度
热门文章
- OpenLayer3调用天地图示例
- BCDEdit命令添加WinPE启动项
- devexpress entity framework 与 asp.net mvc的坑
- Python自动化--语言基础2--运算符、格式化输出、条件语句、循环语句、列表、元组
- $.ajax的一些坑啊
- mac idea中的Application Server was not connected before run configuration stop, reason: Unable to ping server at localhost:1099问题
- applicationContext.xml最基本配置文件
- 初识vue——语法初解
- iOS应用如何得知用户有新拍的图片?
- 未找到与命令“dotnet-ef”匹配的可执行文件