Miller-Rabin & Pollard-rho

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

利用费马小定理，但不能判断伪素数的情况

基于a的伪素数n：

\(a^{n-1} \equiv 1 \pmod n\)

如果对于所有与n互质的数都成立，则n为Carmichael数



定理：

对于质数\(p\)和\(e \ge 1\)

只有两个解\(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);

	}

}