摘要： 定义: 设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有原根的充要 阅读全文

摘要： 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... 阅读全文

摘要： //矩阵类，支持矩阵的加减乘和幂运算 //------------------------------------------------------------------------------------ const int MAXN = 105; const int MAXM = 105; struct Martix { int n, m; ... 阅读全文

摘要： const double eps = 1e-15; //高斯消元模板 //----------------------------------------------------------------------------------- //把对应得系数矩阵化为对角矩阵，然后直接回代即可 const int maxn=100+5; const int maxm=100+5; //a为增广矩阵... 阅读全文

摘要： poj3420 阅读全文

摘要： 求单个数的欧拉函数 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; } 阅读全文

摘要： 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 阅读全文

摘要： 通过扩展欧几里得算法获取x或者y的最小整数解 阅读全文

摘要： //****************************************************************// Miller_Rabin 算法进行素数测试//速度快，而且可以判断 <2^63的数//************************************** 阅读全文