luogu P4035 [JSOI2008]球形空间产生器（高斯消元 / 模拟退火）

整理的算法模板合集： ACM模板

点我看算法全家桶系列！！！

实际上是一个全新的精炼模板整合计划

数据范围只开到了10，而且是经典的力学结构，所以我们可以用模拟退火，可以做一下 nnn 维的正交分解hhh，经典物理题。

然后实际上本题就是一个高斯消元模板题

apple pencil用习惯了，普通的笔写起来，字丑见谅

我用的是普通的高斯消元模板，因为听说高斯 - 约旦消元法有可能会被卡掉，所以学了，但是一直不敢用QWQ

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <map>
#include <queue>
using namespace std;
typedef long long ll;
typedef int itn;
typedef pair<int, int>PII;
const int N = 2007, mod = 1e9 + 7, INF = 2.1e9;
const double eps = 1e-8;int n, m;
double b[N][N];//增广矩阵
double a[N][N];int guass(double a[][N])
{int c;//当前最左边的一列int r;//当前最上面的一行for(c = 1, r = 1; c <= n; ++ c){int t = r;for(int i = r + 1; i <= n; ++ i){if(fabs(a[i][c]) > fabs(a[t][c]))t = i;}if(fabs(a[t][c]) < eps)//等于0就下一个continue;//换到第一行（最上面的一行）for(int i = c; i <= n + 1; ++ i)swap(a[t][i], a[r][i]);//把该行第一个数变成 1（该行每一列都/=第c列的值）for(int i = n + 1; i >= c; -- i)a[r][i] /= a[r][c];//中间所有列（所有项）全部消掉for(int i = r + 1; i <= n; ++ i)if(fabs(a[i][c]) > eps)for(int j = n + 1; j >= c; -- j)a[i][j] -= a[r][j] * a[i][c];r ++ ;}if(r <= n){for(int i = r; i <= n; ++ i)if(fabs(a[i][n + 1]) > eps)//非零，无解return 2;return 1;//无穷多组解}//回代//行最简形矩阵for(int i = n; i >= 1; -- i)for(int j = i + 1; j <= n + 1; ++ j)a[i][n + 1] -= a[j][n + 1] * a[i][j];    return 0;
}int main()
{scanf("%d", &n);for(int i = 1; i <= n + 1; ++ i)for(int j = 1; j <= n; ++ j)scanf("%lf", &a[i][j]);for(int i = 1; i <= n; ++ i)for(int j = 1; j <= n; ++ j){b[i][j] = 2 * (a[i][j] - a[i + 1][j]);b[i][n + 1] += a[i][j] * a[i][j] - a[i + 1][j] * a[i + 1][j];}guass(b);for(int i = 1; i < n; ++ i)printf("%.3f ", b[i][n + 1]);printf("%.3f\n", b[n][n + 1]);return 0;
}

luogu P4035 [JSOI2008]球形空间产生器（高斯消元 / 模拟退火）相关推荐

BZOJ.1013.[JSOI2008]球形空间产生器(高斯消元)
题目链接 HDU3571 //824kb 40ms //HDU3571弱化版跟那个一比这个太水了,练模板吧. //列出$n+1$个二次方程后两两相减,就都是一次方程了. #include <c ...
高斯消元简单线性代数线性基学习记录
线性代数,唉高斯消元 P4035 [JSOI2008]球形空间产生器题目描述有一个球形空间产生器能够在n维空间中产生一个坚硬的球体.现在,你被困在了这个n维球体中,你只知道球面上n+1个点的坐标 ...
BZOJ 1013: [JSOI2008]球形空间产生器sphere 高斯消元
1013: [JSOI2008]球形空间产生器sphere Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://www.lydsy.com/Judg ...
【高斯消元】球形空间产生器（luogu 4035/金牌导航高斯消元-1）
球形空间产生器 luogu 4035 金牌导航高斯消元-1 题目大意给出n+1个n维的点,让你求一个点,使该点到所有点欧几里得距离相等输入样例 2 0.0 0.0 -1.0 1.0 1.0 0. ...
BZOJ 1013: [JSOI2008]球形空间产生器sphere( 高斯消元 )
可以得到N条N元一次方程, 高斯消元就OK了..时间复杂度O(N3) ----------------------------------------------------------------- ...
bzoj 1013: [JSOI2008]球形空间产生器sphere（高斯消元）
1013: [JSOI2008]球形空间产生器sphere Time Limit: 1 Sec Memory Limit: 162 MB Submit: 5598 Solved: 2943 [Su ...
【高斯消元】[JSOI2008]球形空间产生器sphere
[JSOI2008]球形空间产生器sphere 省选难度的板子 + 普及思维 = 紫题 (没毛病日推天天推紫题再见了同学们我自请退群为什么要我一个三维生物想象n维空间呢以二维举例(糊一波我的手稿 ...
【BZOJ】1013: [JSOI2008]球形空间产生器sphere（高斯消元）
http://www.lydsy.com/JudgeOnline/problem.php?id=1013 只要列出方程组就能套高斯来解了. 显然距离相等,所以开不开平方都无所谓. b表示圆心,可列 s ...
BZOJ 1013 JSOI2008 球形空间产生器sphere 高斯消元
题目大意:给定n维空间下的n+1个点,求这n个点所在的球面的球心以前尝试了非常久的模拟退火0.0 至今仍未AC 0.0 今天挖粪涂墙怒学了高斯消元-- 我们设球心为X(x1,x2,...,xn) 如 ...

luogu P4035 [JSOI2008]球形空间产生器（高斯消元 / 模拟退火）

luogu P4035 [JSOI2008]球形空间产生器（高斯消元 / 模拟退火）相关推荐

最新文章

热门文章