pollard rho

这篇写的蛮好的：https://files-cdn.cnblogs.com/files/Doggu/Pollard-rho%E7%AE%97%E6%B3%95%E8%AF%A6%E8%A7%A3.pdf

miller rabin判断质数x：通过选取几个质数p，将x-1分解成2的次幂乘一个常数形式，然后算p的x-1次模x意义下的值，根据费马小定理得一定要是1才行，然后又跟据a^2 = 1 (mod x)，得x|(a-1)(a+1)，所以 a=1或x-1 (mod x),如果是一就继续，直到运行到了x-1或者所有的2用完了为止。选用多个质数p多次check基本上在一个较大范围内都不会有问题了。复杂度logx。

pollar rho：首先根据生日悖论，选多个值，他们的所有的差值（绝对值）是所求数的因数概率会比单纯的高很多(即选取根号级别就有约%50几率)，但是选了根号个算他们所有的差还是n级别的，因此考虑算他们的差与所求数最大公约数不是一，这样概率又被提高了一个级别，需选取根号（根号n）个数，但因这样会空间爆炸，所以并不随机生成 k 个数并两两进行比较，而是用函数?(?) = ( ? 2 + ? ) mod ?一个一个地生成并检查连续的两个数。但是这样rand可能会成环，因此用Floyd或brent判断出环时调整一波函数的初值和随机种子a即可。

代码（洛谷永远被卡一个点）：

#pragma GCC optimize(3,"inline","Ofast")
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const ll pr[]={2,3,5,7,13,11,17};
ll ans,seed;

ll Mul(ll a,ll b,ll m) {
    ll d=((long double)a/m*b+1e-8);
    ll r=a*b-d*m;
    return r<0?r+m:r;
}

ll qpow(ll x,ll y,ll mod)
{
    ll res=1;
    while(y)
    {
        if(y&1)res=Mul(res,x,mod);
        x=Mul(x,x,mod),y>>=1;
    }
    return res;
}

bool chk(ll x,ll p)
{
    ll y=x-1,tp;
    if(x==2||x==3||x==5||x==7||x==11||x==13||x==17)return 1;
    if(x%2==0)return 0;if(x%3==0)return 0;
    while(!(y&1))y>>=1;
    tp=qpow(p,y,x);
    while(y!=x-1&&tp!=x-1&&tp!=1)
    {
        tp=Mul(tp,tp,x);
        y<<=1;
    }
    return tp==x-1||(y&1)==1;
}

bool sol(ll x)
{
    bool isp=1;
    if(x<=1)return 0;
    for(int i=0;i<7;i++)
    {
        isp&=chk(x,pr[i]);
        if(!isp)break;
    }
    return isp;
}

ll gcd(ll x,ll y)
{return !y?x:gcd(y,x%y);}

ll f(ll nw,ll x)
{return (Mul(nw,nw,x)+seed)%x;}

ll pl(ll x)
{
    ll a,b,d;
    if(x%2==0)return 2;
    if(x%3==0)return 3;
    if(x%5==0)return 5;
    if(x%7==0)return 7;
    a=b=rand()%(x-1)+1;
    while(1)
    {
        d=gcd(abs(a-b),x);
        if(d!=1&&d!=x)return d;
        a=f(a,x),b=f(f(b,x),x);
        if(a==b)return x;
    }
}

void ask(ll x)
{
    if(sol(x)){ans=max(ans,x);return;}
    ll ys=x;seed=19260817;
    while(ys==x)--seed,ys=pl(x);
    ask(ys),ask(x/ys);
    return;
}

int main()
{
//	freopen("t.in","r",stdin);
    srand(19260817);
    int n;ll x;
    scanf("%d",&n);
    while(n--)
    {
        scanf("%lld",&x);
        ans=0,ask(x);
        if(ans==x)puts("Prime");
        else printf("%lld\n",ans);
    }
    return 0;
}