#include<iostream> #include<cmath> #include<cstdio> #include<cstring> #include<algorithm> using namespace std; long long mul(long long a,long long n,long long mo){ ; while (n){ ) ans=(ans+a)%mo; a=(a+a)%mo; n/=; } return ans;

/* 直接费马小定理 */ #include<iostream> #include<cstdio> #include<cstdlib> #include<ctime> #define ll long long using namespace std; ll slow_mul(ll a,ll b,ll c) { ll ans=; a=a%c;b=b%c; while(b) { ) { b--; ans+=a; ans%=c; } a<<=;a%=c

How many prime numbers Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 12955    Accepted Submission(s): 4490 Problem Description   Give you a lot of positive integers, just to find out how many

Prime Test Time Limit: 6000MS   Memory Limit: 65536K Total Submissions: 29046   Accepted: 7342 Case Time Limit: 4000MS Description Given a big integer number, you are required to find out whether it's a prime number. Input The first line contains the

费马定理的逆定理几乎可以用来判断一个数是否为素数,但是有一些数是判断不出来的,因此,Miller_Rabin测试方法对费马的测试过程做了改进,克服其存在的问题. 推理过程如下(摘自维基百科): 摘自另一篇博文(手动滑稽): 原理明白了,就直接上代码了(KuangBin大神的板子): 代码思路是, Miller_Rabin()函数随机选取 s 个a,a用做“基底” check() 函数是用来判断x是否等于1,也就是判断a是否是n的凭证. Mul_mod()函数是 快速乘 ,求 a^t % n 之后

/** Miller_Rabin 算法进行素数测试 快速判断一个<2^63的数是不是素数,主要是根据费马小定理 */ #define ll __int128 ; ///随机化算法判定次数 ll MOD; ///计算ret=(a*b)%c a,b,c<2^63 ll mult_mod(ll a,ll b,ll c) { a%=c; b%=c; ll ret=; ll temp=a; while(b) { ) { ret+=temp; if(ret>c) ret-=c;//直接取模慢很多 }

//****************************************************************// Miller_Rabin 算法进行素数测试//速度快,而且可以判断 <2^63的数//****************************************************************const int S=20;//随机算法判定次数,S越大,判错概率越小 //计算 (a*b)%c.   a,b都是long long的数,直接相乘

Miller_Rabin素数测试     Miller_Rabin判断单个素数的方法运用了费马小定理,可以说非常之快了.     Miller_Rabin曾经被称作“黑科技”,但是根据费马小定理其实完全可以自己写出来大半. 其算法的运行过程如下: (1)对于奇数M,使得N=(2^r)*M+1 (2)选取随机数使得A<N (3)对于任意i(i<r),若(A^(2^i)) Mod N=N - 1,则N为素数 (4)或者,若(A^M) Mod N=1,则N通过随机数A的测试 若对素数N进行T次测试,

常用的hash函数是选一个数m取模(余数),这个数在课本中推荐m是素数,但是经常见到选择m=2^n,因为对2^n求余数更快,并认为在key分布均匀的情况下,key%m也是在[0,m-1]区间均匀分布的.但实际上,key%m的分布同m是有>关的. 证明如下: key%m = key - xm,即key减掉m的某个倍数x,剩下比m小的部分就是key除以m的余数.显然,x等于key/m的整数部分,以floor(key/m)表示.假设key和m有公约数g,即key=ag, m=bg, 则 key - x

Problem about GCD Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 470    Accepted Submission(s): 77 Problem Description Given integer m. Find multiplication of all 1<=a<=m such gcd(a, m)=1 (cop

#include<stdio.h> #include<string.h> #include<stdlib.h> #include<time.h> #include<iostream> #include<algorithm> using namespace std; //**************************************************************** // Miller_Rabin 算法进

集训队有人提到这个算法,就学习一下,如果用到可以直接贴模板,例题:POJ 1811 转自:http://www.cnblogs.com/kuangbin/archive/2012/08/19/2646396.html 传说中的随机算法. 效率极高. 可以对一个2^63的素数进行判断. 可以分解比较大的数的因子. #include<stdio.h> #include<string.h> #include<stdlib.h> #include<time.h> #

#include<iostream> #include<cstdio> #include<queue> #include<cstring> #include<algorithm> #include<cmath> #include<cstdlib> #include<ctime> #define lson l, m, rt<<1 #define rson m+1, r, rt<<1|1 #

给出一个64位的大数,如何快速判断其是否为素数 #include<algorithm> #include<cstdio> #include<cstring> #include<cmath> using namespace std; typedef long long LL; LL n,m; //**************************************************************** // Miller_Rabin 算法

Miller-rabin算法是一个用来快速判断一个正整数是否为素数的算法.它利用了费马小定理,即:如果p是质数,且a,p互质,那么a^(p-1) mod p恒等于1.也就是对于所有小于p的正整数a来说都应该复合a^(p-1) mod p恒等于1.那么根据逆否命题,对于一个p,我们只要举出一个a(a<p)不符合这个恒等式,则可判定p不是素数.Miller-rabin算法就是多次用不同的a来尝试p是否为素数. 但是每次尝试过程中还做了一个优化操作,以提高用少量的a检测出p不是素数的概率.这个优化叫做

题意:求n是否只有4个因子,如果是就输出除1外的所有因子. 模板题,就不排版了 #include<cstdio> #include<iostream> #include<algorithm> #include<cstring> #include<cmath> #include<queue> #include<map> #include<ctime> using namespace std; #define MO

题目链接:http://poj.org/problem?id=1811 题目解析:2<=n<2^54,如果n是素数直接输出,否则求N的最小质因数. 求大整数最小质因数的算法没看懂,不打算看了,直接贴代码,以后当模版用. 数据比较大,只能先用Miller_Rabin算法进行素数判断. 在用Pollard_rho分解因子.   #include <iostream> #include <stdio.h> #include <string.h> #include

Prime Test Time Limit: 6000MS   Memory Limit: 65536K Total Submissions: 27129   Accepted: 6713 Case Time Limit: 4000MS Description Given a big integer number, you are required to find out whether it's a prime number. Input The first line contains the

//**************************************************************** // Miller_Rabin 算法进行素数测试 //速度快,而且可以判断 <2^63的数 //**************************************************************** ;//随机算法判定次数,S越大,判错概率越小 LL mult_mod(LL a,LL b,LL mod) //(a*b)%c a,b,c<

大数因数分解Pollard_rho 算法 复杂度o^(1/4) #include <iostream> #include <cstdio> #include <algorithm> #include <cmath> #include <cstring> #include <map> using namespace std; ; ; map<long long, int>m; long long Random( long l

#include<stdio.h> #include<string.h> #include<stdlib.h> #include<time.h> #include<iostream> #include<algorithm> using namespace std; //**************************************************************** // Miller_Rabin 算法进