ZOJ 3645高斯消元

高斯消元是求方程组解的一般方法：

它在求矩阵的逆矩阵的时候也最方便的方法

View Code

#include <stdio.h>
#include <cstring>
#include <stdlib.h>
#include <math.h>
#include <iostream>
using namespace std;
#define max(a,b) ((a) > (b) ? (a) : (b))
#define min(a,b) ((a) < (b) ? (a) : (b))
#define eps 1e-6
const int Max = 19;double mat[Max][Max];
double c[Max][Max];
double ans[Max];int main()
{int t;scanf("%d",&t);while(t--){for(int i=0;i<12;i++){for(int j=0;j<12;j++){scanf("%lf",&mat[i][j]);}}for(int i=0;i<12;i++)c[i][11] = mat[i][11]*mat[i][11];for(int i=0;i<12;i++){for(int j=0;j<11;j++){c[i][j] = -2*mat[i][j];c[i][11] -= mat[i][j]*mat[i][j];}}for(int i=0;i<11;i++){for(int j=0;j<12;j++){c[i][j]-=c[i+1][j];}}for(int i=0;i<11;i++){if(c[i][i] == 0){for(int j=i+1;j<11;j++){if(c[j][i] != 0){for(int k=0;k<12;k++)swap(c[i][k],c[j][k]);break;}}}for(int j=i+1;j<11;j++){double be = c[j][i]/c[i][i];if(be == 0)continue;for(int k=i;k<12;k++){c[j][k] -= c[i][k]*be;}}}ans[10] = c[10][11]/c[10][10];for(int i=9;i>=0;i--){ans[i]=c[i][11];for(int j=i+1;j<11;j++){ans[i]-=c[i][j]*ans[j];}ans[i] /= c[i][i];}for(int i=0;i<10;i++){printf("%.2lf ",ans[i]+eps);}printf("%.2lf\n",ans[10]+eps);} return 0;
}

转载于:https://www.cnblogs.com/gray035/archive/2013/04/08/3008290.html

ZOJ 3645高斯消元相关推荐

高斯消元 zoj 3645 poj 1222/XOR消元
double a[maxn][maxn] , x[maxn] ; //a[i][j] 系数矩阵 , a[i][n+1] = y[i] , x解 int n ; //n个方程 void guass(){ ...
zoj 3353 Chess Board 高斯消元
高斯消元其本应用: 这题是要枚举自由元的: #include<iostream> #include<cstring> #include<cstdio> #defin ...
Rocksdb Ribbon Filter : 结合 XOR-filter 以及高斯消元算法实现的高效filter
文章目录前言 XOR-filter 实现原理 xor filter 的构造原理 xor filter 构造总结 XOR-filter 和 ADD-filter对比 XOR-filter 在计算上的优 ...
poj 1681 Painter#39;s Problem(高斯消元)
http://poj.org/problem? id=1681 求最少经过的步数使得输入的矩阵全变为y. 思路:高斯消元求出自由变元.然后枚举自由变元,求出最优值. 注意依据自由变元求其它解及求最优值 ...
AC自动机 + 概率dp + 高斯消元 --- HDU 5955 or 2016年沈阳icpc H [AC自动机 + 概率dp + 高斯消元]详解
题目链接题目大意: 就是有NNN个人,每个人都会猜一个长度为LLL的只包含{1,2,3,4,5,6}\{1,2,3,4,5,6\}{1,2,3,4,5,6}的序列,现在裁判开始投掷骰子,并且把每次的 ...
ICPC 2005 hangzhou Generator （UVA1358）KMP + 期望DP / 高斯消元
整理的算法模板合集: ACM模板点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 Generator Weblink https://www.luogu.com.cn/problem/ ...
2020 ACM / ICPC 济南 A Matrix Equation （高斯消元、乘法原理）
整理的算法模板合集: ACM模板点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划题目链接给你定义两种 010101 矩阵上的运算: Xi,j×Yi,j=(∑k=1NXi,kYk,j ...
luogu P4035 [JSOI2008]球形空间产生器（高斯消元 / 模拟退火）
整理的算法模板合集: ACM模板点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划数据范围只开到了10,而且是经典的力学结构,所以我们可以用模拟退火,可以做一下 nnn 维的正交分解h ...
BZOJ 2707: [SDOI2012]走迷宫 [高斯消元 scc缩点]
2707: [SDOI2012]走迷宫题意:求s走到t期望步数,\(n \le 10^4\),保证\(|SCC| \le 100\) 求scc缩点,每个scc高斯消元,scc之间直接DP 注意每次清 ...
[BZOJ 3143][Hnoi2013]游走（高斯消元+期望）
Description 一个无向连通图,顶点从1编号到N,边从1编号到M. 小Z在该图上进行随机游走,初始时小Z在1号顶点,每一步小Z以相等的概率随机选择当前顶点的某条边,沿着这条边走到下一个顶点, ...

ZOJ 3645高斯消元

ZOJ 3645高斯消元相关推荐

最新文章

热门文章