bzoj4802 欧拉函数

pollard-rho分解质因数。

#include<cstdio>#include<cstdlib>#include<algorithm>using namespace std;#define LL long longLL a[100];int tot;LL Rand(LL n){    return (((LL)rand()<<31)|rand())%(n-1)+1;}LL inc(LL x,LL y,LL p){    x+=y;    return x>=p?x-p:x;}LL dec(LL x,LL y,LL p){    x-=y;    return x<0?x+p:x;}LL mul(LL x,LL y,LL p){    LL ret=0;    for (;y;y>>=1,x=inc(x,x,p))        if (y&1) ret=inc(ret,x,p);    return ret;}LL pow(LL x,LL y,LL p){    LL ret=1;    for (;y;y>>=1,x=mul(x,x,p))        if (y&1) ret=mul(ret,x,p);    return ret;}LL gcd(LL x,LL y){    return y?gcd(y,x%y):x;}int check(LL n){    LL x=n-1,y,z;    int k=0;    while (!(x&1))    {        x>>=1;        k++;    }    for (int i=1;i<=5;i++)    {        y=pow(Rand(n),x,n);        for (int j=1;j<=k;j++)        {            z=mul(y,y,n);            if (z==1&&y!=1&&y!=n-1) return 0;            y=z;        }        if (y!=1) return 0;    }    return 1;}void solve(LL n){    if (check(n))    {        a[++tot]=n;        return;    }    LL x,y,c,d;    x=y=Rand(n);    c=Rand(n);    for (int i=1,k=1;;i++)    {        x=inc(mul(x,x,n),c,n);        d=gcd(dec(x,y,n),n);        if (d==n)        {            solve(n);            return;        }        if (d>1)        {            solve(d);            solve(n/d);            return;        }        if (i==k)        {            y=x;            k<<=1;        }    }}int main(){    LL n;    scanf("%lld",&n);    if (n==1)    {        printf("1\n");        return 0;    }    solve(n);    sort(a+1,a+tot+1);    for (int i=1;i<=tot;i++)        if (a[i]!=a[i-1]) n=n/a[i]*(a[i]-1);    printf("%lld\n",n);}