VIJOS 1052贾老二算算术 (高斯消元)
描述
贾老二是个品学兼优的好学生,但由于智商问题,算术学得不是很好,尤其是在解方程这个方面。虽然他解决 2x=2 这样的方程游刃有余,但是对于 {x+y=3 x-y=1} 这样的方程组就束手无策了。于是他要你来帮忙。前提是一次方程组且保证在integer的范围内可以处理所有问题。
格式
输入格式
第一行一个数字N(1≤N≤100)表示要求的未知数的个数,同时也是所给的方程个数。
第2到N+1行,每行N+1个数。前N个表示第1到N个未知数的系数。第N+1个数表示N个未知数乘以各自系数后的加和。(保证有唯一整数解)
输出格式
一行N个数,表示第1到N个未知数的值。
模板题!
1 #include<algorithm> 2 #include<cstdio> 3 #include<cstring> 4 #include<cmath> 5 #include<iostream> 6 double a[2005][2005]; 7 int n; 8 int gcd(int a,int b){ 9 if (b==0) return a; 10 else return gcd(b,a%b); 11 } 12 void gauss(){ 13 int now=1,to,ggcd; 14 double t; 15 for (int i=1;i<=n;i++){ 16 for (to=now;to<=n;to++) if (a[to][i]!=0) break; 17 if (to>n) continue; 18 if (to!=now) for (int j=1;j<=n+1;j++) std::swap(a[to][j],a[now][j]); 19 t=a[now][i]; 20 for (int j=1;j<=n+1;j++) a[now][j]/=t; 21 for (int j=1;j<=n;j++) 22 if (j!=now){ 23 t=a[j][i]; 24 for (int k=1;k<=n+1;k++) 25 a[j][k]-=t*a[now][k]; 26 } 27 now++; 28 } 29 } 30 int main(){ 31 scanf("%d",&n); 32 for (int i=1;i<=n;i++){ 33 for (int j=1;j<=n+1;j++) 34 scanf("%lf",&a[i][j]); 35 } 36 gauss(); 37 for (int i=1;i<=n;i++) 38 printf("%d ",(int)round(a[i][n+1])); 39 }
转载于:https://www.cnblogs.com/qzqzgfy/p/5543948.html
VIJOS 1052贾老二算算术 (高斯消元)相关推荐
- 【Vijos P1052】 贾老二算算术 高斯消元傻题
#include <stdio.h> int main() {puts("转载请注明出处谢谢");puts("http://blog.csdn.net/vmu ...
- 【高斯消元】兼 【期望dp】例题
[总览] 高斯消元基本思想是将方程式的系数和常数化为矩阵,通过将矩阵通过行变换成为阶梯状(三角形),然后从小往上逐一求解. 如:$3X_1 + 2X_2 + 1X_3 = 3$ $ ...
- hihoCoder #1195 : 高斯消元·一
时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 小Ho:喂不得了啦,那边便利店的薯片半价了! 小Hi:啥?! 小Ho:那边的便利店在打折促销啊. 小Hi:走走走,赶紧去看 ...
- 【BZOJ4820】[SDOI2017]硬币游戏(高斯消元)
[BZOJ4820][SDOI2017]硬币游戏(高斯消元) 题面 BZOJ 洛谷 题解 第一眼的感觉就是构\(AC\)自动机之后直接高斯消元算概率,这样子似乎就是\(BZOJ1444\)了.然而点数 ...
- 【SDOI2017】硬币游戏【KMP】【概率期望】【高斯消元】
题意:给 nnn 个长度为 mmm 的 01 串,一个 01 串初始为空,不断随机一个字符加在后面,当出现给定的 nnn 个串中的一个时停止.分别求在 nnn 个串处停止的概率. 考场思路历程: 显然 ...
- 线性代数五之高斯消元——[SDOI2010]外星千足虫,[HNOI2013]游走,[HNOI2011]XOR和路径,[hdu 4035]Maze
多类型高斯消元杂题 [SDOI2010]外星千足虫 description solution code [HNOI2013]游走 description solution code [HNOI2011 ...
- HDU 4870 Rating(高斯消元 )
HDU 4870 Rating 这是前几天多校的题目,高了好久突然听旁边的大神推出来说是可以用高斯消元,一直喊着赶快敲模板,对于从来没有接触过高斯消元的我来说根本就是一头雾水,无赖之下这几天做DP ...
- bzoj 2707: [SDOI2012]走迷宫(Trajan+高斯消元+Dp)
2707: [SDOI2012]走迷宫 Time Limit: 10 Sec Memory Limit: 128 MB Submit: 862 Solved: 328 [Submit][Statu ...
- 2021 ICPC Asia Jinan Regional Contest-J Determinant(取模高斯消元)
题面链接 https://pintia.cn/market/item/1459833348620926976 题面 题意 给你一个 n×nn\times nn×n 的矩阵,并给你一个长度为 10410 ...
最新文章
- dedecms /member/reg_new.php SQL Injection Vul
- CF570D Tree Requests(dsu on tree)
- epoll使用详解(精髓)
- 转的:SQL执行提高效率的多种方法
- 基于JAVA+SpringBoot+Vue+Mybatis+MYSQL的微信小程序便利店购物系统
- Python中使用random随机函数与Matplotlib库绘制随机漫步图
- 航天信息上传参数设置服务器设置,金税盘上传参数怎么设置?
- php获取客户端IP和php获取服务器端IP
- 树莓派学习-I2c通信
- CrystalDiskMark简介
- 基于线性回归对男性体脂率的预测
- 深恶痛绝的No mapping found for HTTP request with URI
- doctrine2 mysql_php – Doctrine2和MySQL分区
- 互联网装修O2O模式是否可行?
- 前端后端路径斜杆问题
- 一个比较全的C++农历算法
- 【蓝桥杯每日一练】—字符逆序
- shell之cut ---sed---awk--sort
- 跨考计算机面试英语自我介绍,2019考研复试面试英语自我介绍范文(1)
- Word文档很乱怎么办 杂乱的文章word排版教程