随笔分类 -  ACM数学模板

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

摘要：const int maxn=30000+5; int prime[maxn]; void marktable(int n){ memset(prime,0,sizeof(prime)); for(int i=2;i<=n;i++){ if(!prime[i]) prime[++prime[0]]=i; for(int j=1;j<=prim... 阅读全文

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

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