目录 [TOC] 前言 sympy不仅在符号运算方面强大，在解方程方面也是很强大。 本章节学习对应官网的：Solvers "官方教程" https://docs.sympy.org/latest/tutorial/solvers.html （一）求解多元一次方程 solve() 1.说明： 解多元一 ...
分类：
编程语言 时间：
2019-03-17 23:24:22
阅读次数：
473
二次剩余 $Cipolla$ 算法 概述 大概就是在模 $p$ 意义下开根号,如求解方程$x^2\equiv n(mod\ p)$. 这里只考虑 $p$ 为素数的情况.若 $p=2$ ,则$x=0\ when\ n=0,x=1\ when\ n=1$. 若 $p?$ 为奇素数,定义勒让德符号: $$ ...
分类：
其他好文 时间：
2019-03-11 13:39:22
阅读次数：
187
"洛咕" 题意:求关于x的同余方程$ax\equiv1\pmod{b}$的最小正整数解. 方程$ax\equiv1\pmod{b}$有解当且仅当$gcd(a,b)=1$.所以方程可写为$a x+b y=1$,用扩展欧几里得算法求出一组特解$x_0,y_0$,通解是所有模b与$x_0$同余的整数,题目 ...
分类：
其他好文 时间：
2019-03-09 14:16:17
阅读次数：
233
一次同余方程 前面已经提到，剩余类可以看作一个特殊的“数”，剩余类环可以看作定义了剩余类加法和乘法的“数集”.类似于实数集情形，我们也可以在剩余类环中解方程或方程组。 例如，在模6的剩余类环中解方程[5][x]=3，这里[x]是模6的剩余类环中的未知剩余类，注意到 $$[5][x] = [3]\Le ...
分类：
其他好文 时间：
2019-02-27 20:20:17
阅读次数：
192
2019-02-13 P3704 [SDOI2017]数字表格：莫比乌斯反演 P3702 [SDOI2017]序列计数：快速幂+多项式 2019-02-14 P3703 [SDOI2017]树点涂色：LCT+线段树 Yahoo Programming Contest 2019 D - Ears：dp ...
分类：
其他好文 时间：
2019-02-15 22:28:11
阅读次数：
209
求解方程 求解一个给定的方程，将x以字符串"x=#value"的形式返回。该方程仅包含'+'，' - '操作，变量 x 和其对应系数。 如果方程没有解，请返回"No solution"。 如果方程有无限解，则返回"Infinite solutions"。 如果方程中只有一个解，要保证返回值 x 是一 ...
分类：
其他好文 时间：
2019-02-15 13:39:56
阅读次数：
365
//弦截法求解方程的根 //要求：输入左右两个端点x值 //结果：在一定精度范围内求解出方程的根 //难点：1)推导出x处的横坐标的求解公式 2）迭代掉原来的左端点或者右端点 #include "pch.h" #include #include #include using namespace st... ...
分类：
其他好文 时间：
2019-01-23 10:35:55
阅读次数：
191
PART1 <求解方程> 1，二分法 def bisect(f,a,b,TOL=0.000004): u_a = a u_b = b while(u_b-u_a)/2.0 > TOL: c = (u_a+u_b)/2.0 if f(c) == 0: break if f(u_a)*f(c) < 0: ...
分类：
其他好文 时间：
2019-01-23 00:20:26
阅读次数：
175
这里用到的还是最小二乘方法，和上一次这篇文章原理差不多。 就是首先构造最小二乘函数，然后对每一个系数计算偏导，构造矩阵乘法形式，最后解方程组。 比如有一个二次曲面：z=ax^2+by^2+cxy+dx+ey+f 首先构造最小二乘函数，然后计算系数偏导（我直接手写了）： 解方程组（下图中A矩阵后面求和 ...
分类：
其他好文 时间：
2019-01-21 13:55:20
阅读次数：
192
【训练题】高斯消元[1] 描述 贾老二是个品学兼优的好学生，但由于智商问题，算术学得不是很好，尤其是在解方程这个方面。虽然他解决 2x=2 这样的方程游刃有余，但是对于 {x+y=3 x-y=1} 这样的方程组就束手无策了。于是他要你来帮忙。前提是一次方程组且保证在integer的范围内可以处理所有 ...
分类：
其他好文 时间：
2019-01-13 17:04:17
阅读次数：
212
