【bzoj3667】Rabin-Miller算法 素数测试

AC通道：http://www.lydsy.com/JudgeOnline/problem.php?id=3667

【吐槽】

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

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

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

/*************
  bzoj 3667
  by chty
  2016.11.7
*************/
#include
#include
#include
#include
#include
#include
#include
using namespace std;
typedef long long ll;
const ll prime[20]={2,3,5,7,11,13,17,19,23};
ll T,maxx;
inline ll read()
{
    ll x=0,f=1;  char ch=getchar();
    while(!isdigit(ch))  {if(ch=='-')  f=-1;  ch=getchar();}
    while(isdigit(ch))  {x=x*10+ch-'0';  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=1;for(;b;b>>=1,a=mul(a,a,mod))if(b&1)sum=mul(sum,a,mod);return sum;}
bool Rabin_Miller(ll p,ll a)
{
    if(p==2)  return 1;
    if((p&1==0)||p==1)  return 0;
    ll d=p-1;
    while(!(d&1))  d>>=1; 
    ll m=fast(a,d,p);
    if(m==1)  return 1;
    for(;d=0;m=mul(m,m,p),d<<=1)  {if(m==p-1)  return 1;}
    return 0;
}
bool isprime(ll x)
{
    for(ll i=0;i<9;i++)
    {
        if(prime[i]==x)  return 1;
        if(!Rabin_Miller(x,prime[i]))  return 0;
    }
    return 1;
}
void Pollord_Rho(ll x)
{
    if(isprime(x))  {maxx=max(maxx,x); return;}
    ll c=3;
    while(1)
    {
        ll x1(1),x2(1),i(1),k(2);
        while(1)
        {
            x1=(mul(x1,x1,x)+c)%x;
            ll d=gcd(abs(x1-x2),x);
            if(d>1&&d