Miller_Rabbin大素数测试

传送门：点击打开链接

题意：验证哥德巴赫猜想，任何一个偶数都可以表示为两个质数之和，输出任意一对质数。（这题有点坑，主要复习Miller_Rabbin）

思路：这题主要考察大素数的判断，利用Miller_Rabbin素数测试。

 

引理：（费马小定理） 
a^p−1≡1(modp),(a,p)=1
快速判定一个数是否为素数的方法： 
如果存在一个整数a，满足a^p−1≡1(modp)的话，那么n称为基于a的伪素数。 
当有多个满足关系的a时，n为素数的概率趋向于1。 
所以多取几个a测试即可。

 

代码：

 

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
using namespace std;
#define ll unsigned long long

ll random(ll n){///产生大随机数 0-n 
    return (ll)((double)rand()/RAND_MAX*n+0.5);
}

ll qml(ll a,ll b,ll m){///乘法换加法取余防止溢出
    ll s=0;
    while(b){
        if(b&1) s=(s+a)%m;
        a=(a+a)%m,b/=2;///用ull 防止溢出
    }
    return s;
}

ll mypow(ll a,ll b,ll m){
    ll s=1;
    while(b){
        if(b&1) s=qml(s,a,m);
        a=qml(a,a,m),b/=2;
    }
    return s;
}

bool Miller_Rabbin(ll x){
    if(x==2) return true;        ///2要直接判断
    for(int i=1;i<=11;++i){
        ll a=rand()%(x-2)+2;
        ///ll a=random(x-2)+1;居然会超时
        if(mypow(a,x-1,x)!=1)
            return false;
    }
    return true;
}

int main(){
    ll n,t;
    scanf("%llu",&t);
    while(t--){
        scanf("%llu",&n);
        for(ll i=2;;i++){
            if(Miller_Rabbin(i)&&Miller_Rabbin(n-i)){
                printf("%llu %llu\n",i,n-i);break;
            }
        }
    }
    return 0;
}