高斯消元法 可以用于求解线性方程组，即n元1次方程组。利用矩阵，大致思路与普通解方程方法类似。只是更具一般性。将系数与右侧的常数存成一个矩阵，然后每次用第i行消去下面每行的第i个系数，最后就会得到一个一元方程，然后从后到前依次代回即可。 然后就是精度的问题，因为计算机中没有分数，所以只能用doubl ...
分类：
其他好文 时间：
2018-05-19 15:40:12
阅读次数：
189
题目：https://www.lydsy.com/JudgeOnline/problem.php?id=2337 异或就一位一位考虑； x为到n的概率，解方程组即可； 考虑了n就各种蜜汁错误，所以索性不管n了，这样的题好像不管n比较方便。 代码如下： ...
分类：
其他好文 时间：
2018-05-17 19:57:45
阅读次数：
187
Description 有一个球形空间产生器能够在 n 维空间中产生一个坚硬的球体。现在，你被困在了这个 n 维球体中，你只知道球面上 n+1 个点的坐标，你需要以最快的速度确定这个 n 维球体的球心坐标，以便于摧毁这个球形空间产生器。 Input & Output Input 第一行是一个整数 n ...
分类：
Web程序 时间：
2018-05-16 20:54:33
阅读次数：
262
解析解思路一—矩阵变换 抛物线属于圆锥曲线，圆锥曲线有六个参数，若令F为1，则还剩下五个自由度。 圆锥曲线的投影变换，可以用参数矩阵乘以投影矩阵来实现。矩阵乘法的运算结果是投影后的参数矩阵代表的圆锥曲线。 解析解思路二—高斯消元 上图所示，圆锥曲线有五个自由度。在投影前的抛物线上选出五个点，乘以投影 ...
分类：
其他好文 时间：
2018-05-15 22:57:27
阅读次数：
399
刚才想着把算法复杂度降为O(n^2)，但其实对于第i+1行~第n行加上第k行*一个系数，第k行的各个数的值不是相等的，是我想多了…………………… 程序： AX=B，已知A,B，求X A(i,j) i>=j 都不用处理 1.化为上三角矩阵 (实际上是A(i,j) [i>=j] 被忽略) 2.根据上三角 ...
分类：
其他好文 时间：
2018-05-14 23:01:57
阅读次数：
227
逐次下降法的定义： 下面给出Python实现 # -*- coding: utf-8 -*-import numpy as np from numpy import *from common_libs import *import matplotlib.pyplot as plt #消元发求解方程组 ...
分类：
编程语言 时间：
2018-05-13 19:52:08
阅读次数：
295
这题有人用毒瘤的高斯消元做，也有人DFS，BFS，还有人证明性质然后直接计算，而我选择了最暴力的方法： 2^16枚举！ 一算复杂度发现最坏情况是65536×16×4，无压力...... 秒A 1 #include <cstdio> 2 using namespace std; 3 4 bool a[ ...
分类：
其他好文 时间：
2018-05-10 14:20:25
阅读次数：
159
```cpp include const int N=104; double a[N][N]; int n; double fabs(double x) {return x 0?x: x;} void swap(int i,int j) { double tmp; for(int k=i;k=1;i ...
分类：
其他好文 时间：
2018-05-06 22:26:49
阅读次数：
183
这里仅作为自我检查用，模板代码请移步其他博文 标+的表示已学完，标?的表示需要进一步学习，标-的表示有计划但未开始学习，标*的表示暂时没有计划学习 数学 *BSGS ?FFT&NTT *Lucas&扩展Lucas *差分约束 ?多项式处理 -高斯消元 -线性基 ?矩阵快速幂 -卡特兰数 +扩展欧几里 ...
分类：
编程语言 时间：
2018-05-06 13:29:29
阅读次数：
182
