求素数p的原根
摘要：
定义: 设m>1,gcd(a,m)=1,使得成立的最小正整数d为a对模m的阶，记为δm(a) 如果δm(a)=φ(m),则称a是模m的原根 定理：设m>1,gcd(a,m)=1,那么正整数x是同于方程的一个根当且仅当δm(a) | x 定理：由欧拉定理得 gcd(a,n)=1 定理：模m有原根的充要 阅读全文
posted @ 2018-12-06 11:30 欣崽 阅读(3397) 评论(0) 推荐(0) 编辑
单变元模线性方程模板
摘要：
template void exgcd(T a,T b,T &d,T &x,T &y){ if(!b) {d=a;x=1;y=0;} else {exgcd(b,a%b,d,y,x);y-=x*(a/b);} } //求解a*x=b mod n //如果有解那么解的个数有d=gcd(a,n)个 vector line_mod_equation(long long a,long l... 阅读全文
posted @ 2018-12-04 18:56 欣崽 阅读(243) 评论(0) 推荐(0) 编辑
常系数线性齐次递推
摘要：
//矩阵类，支持矩阵的加减乘和幂运算 //------------------------------------------------------------------------------------ const int MAXN = 105; const int MAXM = 105; struct Martix { int n, m; ... 阅读全文
posted @ 2018-12-03 22:46 欣崽 阅读(150) 评论(0) 推荐(0) 编辑
Gauss消元模板
摘要：
const double eps = 1e-15; //高斯消元模板 //----------------------------------------------------------------------------------- //把对应得系数矩阵化为对角矩阵，然后直接回代即可 const int maxn=100+5; const int maxm=100+5; //a为增广矩阵... 阅读全文
posted @ 2018-12-03 22:25 欣崽 阅读(195) 评论(0) 推荐(0) 编辑
欧拉函数
摘要：
求单个数的欧拉函数 int euler(int n){ //返回euler(n) int res=n,a=n; for(int i=2;i*i<=a;i++){ if(a%i==0){ res=res/i*(i-1);//先进行除法是为了防止中间数据的溢出 while(a%i==0) a/=i; } 阅读全文
posted @ 2018-12-01 00:33 欣崽 阅读(211) 评论(0) 推荐(0) 编辑
快速幂
摘要：
template<class T,class T1> T fast_mod(T a,T b,T1 Mod){ a%=Mod; if(b==0) return 1; T ans=1,base=a; while(b!=0){ if(b&1)ans=(ans*base)%Mod; base=(base*b 阅读全文
posted @ 2018-11-29 23:39 欣崽 阅读(137) 评论(0) 推荐(0) 编辑
Miller_Rabin()算法素数判定 +ollard_rho 算法进行质因数分解
摘要：
//****************************************************************// Miller_Rabin 算法进行素数测试//速度快，而且可以判断 <2^63的数//************************************** 阅读全文
posted @ 2018-11-27 23:12 欣崽 阅读(322) 评论(0) 推荐(0) 编辑
