大素数测试 求因子 poj 1811

抄别人的

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<stdlib.h>
#include<time.h>
#include<map>

using namespace std;
typedef long long ll;
map<ll,int>m1;

ll random(ll n)
{
    return ((double)rand()/RAND_MAX*n+0.5);
}

ll mult(ll a,ll b,ll c)
{
    ll ans=;
    a=a%c;
    while(b>)
    {
        if(b&)
            ans=(ans+a)%c;
        a=(a+a)%c;
        b=b>>;
    }
    return ans;
}
ll quick(ll a,ll b,ll c)
{
    ll ans=;
    a=a%c;
    while(b>)
    {
        if(b&)
            ans=mult(ans,a,c);
        a=mult(a,a,c);
        b=b>>;
    }
    return ans;
}
#define TIMES 10
bool witness(ll a,ll n)
{
    ll d =n-;
    while(!(d&))
        d>>=;
    ll t=quick(a,d,n);
    while(d!=n-&&t!=&&t!=n-)
    {
        t=mult(t,t,n);
        d<<=;
    }
    return t==n-||d&;

}
bool Miller_Rabin(ll n)
{
    if(n==)
        return ;
    if(n<||!(n&))
        return ;
    for(int i=;i<=TIMES;i++)
    {
        ll a=random(n-)+;
        if(!witness(a,n))
            return ;
    }
    return ;

}
ll gcd(ll a,ll b)
{
    return b?gcd(b,a%b):a;
}
ll pollard_rho(ll n,ll c)
{
    ll x,y,d,i=,k=;
    x=random(n-)+;
    y=x;
    while()
    {
        i++;
        x=(mult(x,x,n)+c)%n;
        d=gcd(y-x,n);
        if(d>&&d<n)
            return d;
        if(y==x)
            return n;
        if(i==k)
        {
            y=x;
            k<<=;
        }
    }

}
void find1(ll n,ll c)
{
    if(n==)
        return ;
    if(Miller_Rabin(n))
    {
        m1[n]++;
        return ;
    }
    ll p=n;
    while(p>=n)
        p=pollard_rho(p,c--);
    find1(p,c);
    find1(n/p,c);
}
int main()
{
    int t;
    scanf("%d",&t);

    while(t--)
    {
        ll a;
        scanf("%lld",&a);
        m1.clear();
        if(Miller_Rabin(a))
            printf("Prime\n");
        else
        {
            find1(a,);
            map<ll,int>::iterator i=m1.begin();
            printf("%lld\n",i->first);
        }
    }

    return ;
}