数学,数论
文章平均质量分 73
star-city
蒟蒻OIer一枚
展开
-
高斯消元
高斯消元其实就是个模拟……想想初中怎么学的加减消元,模拟这个过程就好了。 主要看看代码里的一些技巧。 Code #include <iostream> #include <cstring> #include <cstdio> #include <cmath> const int maxn = 107; const double eps =...原创 2018-07-21 13:37:25 · 103 阅读 · 0 评论 -
扩展欧几里得算法
这个算法是用来解决这样的问题的: 给定a,b,ca,b,ca, b, c,求 ax+by=cax+by=cax + by = c 这个不定方程的最小非负整数解。 扩欧求的实际上是 ax+by=gcd(a,b)ax+by=gcd(a,b)ax + by = gcd(a, b) 这个方程的解,所以只有当 gcd(a,b)|cgcd(a,b)|cgcd(a, b)|c 的时候才有整数解。 ...原创 2018-07-21 14:31:48 · 173 阅读 · 0 评论 -
Miller-rabbin算法
这个算法用于快速判断一个数是否是质数。 那怎么判呢? 首先,既然要快速,我们就不能用朴素的O(n‾√)O(n)O(\sqrt n)的试除法。我们会联想到质数的一些性质: ap−1≡1(modp)(0&lt;a&lt;p)ap−1≡1(modp)(0&lt;a&lt;p)a^{p-1} \equiv 1 (mod p)(0 < a < p) 也就是费马小定理了,我们可以由此想到一个简单的方法...原创 2018-07-21 18:28:17 · 635 阅读 · 0 评论 -
BSGS大法——有趣的分块打表
对于三个整数A,B,CA,B,CA, B, C,满足A,B<C<109A,B<C<109A, B < C < 10^9且A,CA,CA, C互质,求一个最小的x使得Ax≡B (mod C)Ax≡B (mod C)A^x \equiv B~(mod~C),怎么做? 分块打表! 令t=⌈C‾‾√⌉,x=i∗t−j (i,j&...原创 2018-07-31 23:56:55 · 299 阅读 · 0 评论