2018宁夏网络赛 B Goldbach （米勒拉宾素数测试）

2018宁夏网络赛 B Goldbach （米勒拉宾素数测试）

题目链接

题目大意：

给你一个偶数n (2<n<=1e18)
让你把n分解成两个素数的和。（如果有多个输出任意一个）

解题思路：

米勒拉宾素数测试板子
这个板子可以在 log 级的时间内判断一个数是不是素数。 用于判断大数是不是素数。

AC代码:

#include<bits/stdc++.h>
#define ll unsigned long long int
using namespace std;
 
ll ModMul(ll a,ll b,ll n){//¿ìËÙ»ýÈ¡Ä£ a*b%n
    ll ans=0;
    while(b){
        if(b&1)
        ans=(ans+a)%n;
        a=(a+a)%n;
        b>>=1;
    }
    return ans;
}
ll ModExp(ll a,ll b,ll n){  //¿ìËÙÃÝÈ¡Ä£ a^b%n
    ll ans=1;
    while(b){
        if(b&1)
        ans=ModMul(ans,a,n);
        a=ModMul(a,a,n);
        b>>=1;
    }
    return ans;
}
bool miller_rabin(ll n){//Miller-RabinËØÊý¼ì²âËã·¨
    ll i,j,a,x,y,t,u,s=10;
    if(n==2)
        return true;
    if(n<2||!(n&1))
        return false;
    for(t=0,u=n-1;!(u&1);t++,u>>=1); //n-1=u*2^t
    for(i=0;i<s;i++){
        a=rand()%(n-1)+1;
        x=ModExp(a,u,n);
        for(j=0;j<t;j++){
            y=ModMul(x,x,n);
            if(y==1&&x!=1&&x!=n-1)
                return false;
            x=y;
        }
        if(x!=1)
            return false;
    }
    return true;
}
int main()
{
	int t;
	ll n;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%lld",&n);
		if(n==4)
		{
			printf("2 2\n");
			continue;
		}
		for(ll i=3;i<=n/2;i+=2)
		{
			if(miller_rabin(i)&&miller_rabin(n-i))
			{
				printf("%lld %lld\n",i,n-i);
				break;
			}
		}
	}
    return 0;
}