[Miller-Rabin & Pollard-rho]【学习笔记】

Miller-Rabin & Pollard-rho

很久之前就学过了...今天重学一遍

---

利用费马小定理，但不能判断伪素数的情况
基于a的伪素数n：
\(a^{n-1} \equiv 1 \pmod n\)
如果对于所有与n互质的数都成立，则n为Carmichael数

定理：
对于质数\(p\)和\(e \ge 1\)
\[ x^2 \equiv 1 \pmod p^e \]
只有两个解\(x=1,\ x=-1\)

分解\(n=u*2^t\),反复平方的时候如果存在非平凡平方根则不是质数
可以证明Carmicheal数一定不是\(p^e\)

Pollard-rho启发式因子分解期望\(O(\sqrt{p})\)找到一个为p的质因子

快速乘要用long double黑科技

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
inline ll read(){
    char c=getchar();ll x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}
ll n;
ll Mul(ll a, ll b, ll P) {
    ll t = (a*b - (ll)((long double)a/P*b+1e-8)*P); 
    return t<0 ? t+P : t;
}
ll Pow(ll a, ll b, ll P) {
    ll ans=1; a%=P;
    for(; b; b>>=1, a=Mul(a, a, P))
        if(b&1) ans=Mul(ans, a, P);
    return ans;
}
bool witness(ll a, ll n, ll u, int t) { 
    ll x=Pow(a, u, n), y=x;
    for(int i=1; i<=t; i++) { 
        x=Mul(x, x, n);
        if(x==1 && y!=1 && y!=n-1) return true;
        y=x;
    }
    return x!=1;
}
bool MillerRabin(ll n) { 
    if(n==2) return true;
    if(n<=1 || !(n&1)) return false;
    ll u=n-1, t=0;
    while(!(u&1)) u>>=1, t++;
    for(int i=1; i<=10; i++) 
        if(witness(rand()%(n-1)+1, n, u, t)) return false;
    return true;
}
ll gcd(ll a, ll b) {return b==0?a:gcd(b, a%b);}
ll rho(ll n, ll c) {
    int k=2; ll x=rand()%n, y=x, d=1;
    for(int i=1; d==1; i++) {
        x=(Mul(x,x,n)+c)%n;
        d=gcd(n, y>x?y-x:x-y);
        if(i==k) y=x, k<<=1;
    }
    return d;
}
ll Max;
void solve(ll n) { 
    if(n==1) return;
    if(MillerRabin(n)) {Max=max(Max, n); return;}
    ll t=n;
    while(t==n) t=rho(n, rand()%(n-1)+1);
    solve(t); solve(n/t);
}

int main() {
    freopen("in","r",stdin);
    srand(317);
    int T=read();
    while(T--) {
        n=read();
        Max=0;
        solve(n);
        if(Max==n) puts("Prime");
        else printf("%lld\n",Max);
    }
}

转载于:https://www.cnblogs.com/candy99/p/6624643.html