【bzoj3667】Rabin-Miller算法

3667: Rabin-Miller算法

Time Limit: 60 Sec  Memory Limit: 512 MB
Submit: 1200  Solved: 363
[Submit][Status][Discuss]

Input

第一行：CAS,代表数据组数（不大于350），以下CAS行，每行一个数字，保证在64位长整形范围内，并且没有负数。你需要对于每个数字：第一，检验是否是质数，是质数就输出Prime 
第二，如果不是质数，输出它最大的质因子是哪个。

Output

第一行CAS(CAS<=350，代表测试数据的组数) 
以下CAS行：每行一个数字，保证是在64位长整形范围内的正数。 
对于每组测试数据：输出Prime，代表它是质数，或者输出它最大的质因子，代表它是和数

Sample Input

6
2
13
134
8897
1234567654321
1000000000000

Sample Output

Prime
Prime
67
41
4649
5

 

 

 

【吐槽】

这是一道模板题，但是出了很玄学的错误，交到bzoj上一直wa。

然后要到了数据，用cena评测，然后发现并没有错误。。。

哪位大神知道这种玄学错误的话，欢迎指正，感激不尽。

 /*************
   bzoj 3667
   by chty
   2016.11.7
 *************/
 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<cstdlib>
 #include<ctime>
 #include<cmath>
 #include<algorithm>
 using namespace std;
 typedef long long ll;
 const ll prime[]={,,,,,,,,};
 ll T,maxx;
 inline ll read()
 {
     ll x=,f=;  char ch=getchar();
     while(!isdigit(ch))  {if(ch=='-')  f=-;  ch=getchar();}
     while(isdigit(ch))  {x=x*+ch-'';  ch=getchar();}
     return x*f;
 }
 ll gcd(ll a,ll b) {return !b?a:gcd(b,a%b);}
 ll mul(ll x,ll y,ll mod) {return ((x*y-(ll)(((long double)x*y+0.5)/mod)*mod)%mod+mod)%mod;}//一行快速乘
 ll fast(ll a,ll b,ll mod) {ll sum=;for(;b;b>>=,a=mul(a,a,mod))if(b&)sum=mul(sum,a,mod);return sum;}
 bool Rabin_Miller(ll p,ll a)
 {
     if(p==)  return ;
     if((p&==)||p==)  return ;
     ll d=p-;
     while(!(d&))  d>>=;
     ll m=fast(a,d,p);
     if(m==)  return ;
     for(;d<p&&d>=;m=mul(m,m,p),d<<=)  {if(m==p-)  return ;}
     return ;
 }
 bool isprime(ll x)
 {
     for(ll i=;i<;i++)
     {
         if(prime[i]==x)  return ;
         if(!Rabin_Miller(x,prime[i]))  return ;
     }
     return ;
 }
 void Pollord_Rho(ll x)
 {
     if(isprime(x))  {maxx=max(maxx,x); return;}
     ll c=;
     while()
     {
         ll x1(),x2(),i(),k();
         while()
         {
             x1=(mul(x1,x1,x)+c)%x;
             ll d=gcd(abs(x1-x2),x);
             if(d>&&d<x)
             {
                 Pollord_Rho(d);
                 Pollord_Rho(x/d);
                 return;
             }
             if(x1==x2)  break;
             if(++i==k)  k<<=,x2=x1;
         }
         c++;
     }
 }
 void solve(ll n)
 {
     if(isprime(n))  {puts("Prime"); return;}
     maxx=;
     Pollord_Rho(n);
     printf("%lld\n",maxx);
 }
 int main()
 {
     freopen("cin.in","r",stdin);
     freopen("cout.out","w",stdout);
     T=read();
     for(ll i=;i<=T;i++) {ll n=read();solve(n);}
     return ;
 }