计蒜客 25985.Goldbach-米勒拉宾素数判定(大素数) (2018 ACM-ICPC 中国大学生程序设计竞赛线上赛 B)

若干年之前的一道题，当时能写出来还是超级开心的，虽然是个板子题。一直忘记写博客，备忘一下。

米勒拉判大素数，关于米勒拉宾是个什么东西，传送门了解一下:biubiubiu~

B. Goldbach

题目传送门

自己看题意吧，直接贴代码了。

代码：

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<cmath>
 #include<cstdlib>
 #include<algorithm>
 #include<queue>
 #include<vector>
 #include<stack>
 using namespace std;
 typedef unsigned long long ll;
 const int maxn=1e9+;
 const double eps=1e-;
 ll add_mod(ll a,ll b,ll mod)     //快乘法 基于快速幂的二分思想
 {
     ll ans=;                    //由于考虑到取模数很大 快速幂会溢出
     while(b)                     //必须使用该方法
     {
         if(b&)                    //我这里写的是非递归版
             ans=(ans+a)%mod;
         a=a*%mod;
         b>>=;
     }
     return ans;
 }

 ll pow_mod(ll a,ll n,ll mod)             //快速幂 递归版
 {
     if(n>)
     {
         ll tmp=pow_mod(a,n>>,mod)%mod;
         tmp=add_mod(tmp,tmp,mod);
         if(n&) tmp=add_mod(tmp,a,mod);
         return tmp;
     }
     return a;
 }

 bool Miller_Rabbin(ll n,ll a) //米勒拉宾素数判断函数主体
 {
     ll d=n-,s=,i;
     while(!(d&))             // 先把(2^s)*d 算出来
     {
         d>>=;
         s++;
     }
     ll t=pow_mod(a,d,n);    //a^d取一次余判断
     if(t== || t==-)        //一或负一则可以声明这可能是质数
         return ;
     for(i=; i<s; i++)               //不是的话继续乘上s个2
     {
         if(t==n-)            //(n-1)*(n-1)%n=1 这一步是优化
             return ;
         t=add_mod(t,t,n);    // 快乘
     }
     return ;
 }

 int is_prime(ll n)
 {
     ll i,tab[]= {,,,}; //本来应该取[1,n]内任意整数
     for(i=; i<; i++)               //但一般这几个数足以,不需要太多组测试
     {
         if(n==tab[i])
             return ;        //小判断小优化~
         if(!n%tab[i])
             return ;
         if(n>tab[i] && !Miller_Rabbin(n,tab[i]))
             return ;
     }
     return ;
 }
 ll max(ll a,ll b)
 {
     return a>b?a:b;
 }
 int main()
 {
     int t;
     scanf("%d",&t);
     while(t--)
     {
         ll n;
         scanf("%llu",&n);
         if(n==)
             printf("2 2\n");
         else
         {
             for(ll i=; i<=n/; ++i)
             {
                 if(is_prime(i)&&is_prime(n-i))
                 {
                     printf("%llu %llu\n",i,n-i);
                     break;
                 }
             }
         }
     }
 }

当时还是我队友帮我D的bug，他现在都不理我了((ಥ_ಥ))

难受