Mathematics
Wearry
这个作者很懒,什么都没留下…
展开
-
拓展Lucas定理求C(n, m) % p
#include<bits/stdc++.h>using namespace std;const int maxn = 100000 + 10;typedef long long LL;LL Pow(LL n, LL m, LL mod) { LL ans = 1; while(m > 0) { if(m & 1) ans = (LL)ans * n % mod原创 2016-08-11 15:47:48 · 832 阅读 · 0 评论 -
浅谈基于随机性算法的素数判定
引:考虑这样的一个问题判断一个long long范围内的数是否是素数。引理1-1: ap≡a(modp)|(a,p)=1 a^p \equiv a \pmod p | (a, p) = 1 其实这就是熟悉的Fermat小定理了~~引理1-2: 对于 a2≡1(modn)|n为素数a^2 \equiv 1 \pmod n | n 为素数 的解仅有 a = 1 或 a = n - 1原创 2016-09-12 22:13:39 · 645 阅读 · 2 评论 -
Radix-2 Recursive-FFT
#include<bits/stdc++.h>using namespace std;const int INF = 0x3f3f3f3f;const int maxl = 300000 + 10;const double PI = acos(-1.0);int base, c[maxl];complex<double> a[maxl], b[maxl], p[maxl];int read原创 2016-08-18 14:14:41 · 580 阅读 · 0 评论 -
Radix-2 Iterative-FFT
#include<bits/stdc++.h>using namespace std;const double PI = acos(-1.0);const int maxn = 300000 + 10;int n, m, base;int c[maxn], rev[maxn];complex<double> a[maxn], b[maxn], p[maxn];int read() {原创 2016-08-24 17:25:42 · 684 阅读 · 1 评论