2018 ACM-ICPC 中国大学生程序设计竞赛线上赛 B(素数筛 + MR)

题意：给你一个合数，将其分解为两个质数，合数<2^63

题解：素数筛打表打到1e6，然后对每一个素数用MR判断n-这个素数是不是素数，如果是输出就好了。
（暴力大法好）
但是，所给的合数太大了，MR中很容易爆，所以要用unsigned longlong，输入输出用%llu。
一下午疯狂T，以为MR太慢。
（牛* 牛 *，本垃圾这厢有礼了。）

#include<cstdio>
#include<cstring>
#include<ctime>
#include<iostream>
#include<algorithm>
#include<bits/stdc++.h>
#define LL unsigned long long
using namespace std;
const int maxn=1e6+10;
int S;//规定随机多少次，次数越多正确率越高

LL quick_mod_multi(LL a,LL b,LL mod)//快速求a*b%mod
{
    LL t=0;
    a%=mod;
    while(b)
    {
        if(b&1) t=(t+a)%mod;
        a=(a<<1)%mod;
        b>>=1;
    }
    return t;
}
LL quick_mod_power(LL a,LL b,LL mod)//快速求a^b%mod
{
    LL res=1;
    a%=mod;
    while(b)
    {
        if(b&1)
            res=quick_mod_multi(res,a,mod);
        b>>=1;
        a=quick_mod_multi(a,a,mod);
    }
    return res;
}
bool miller_rabbin(LL n)//素性检验
{
    if(n==2) return true;
    if(n<2||!(n&1)) return false;
    int t=0;
    LL a,x,y,u=n-1;
    while((u&1)==0) t++,u>>=1;
    srand((unsigned)time(NULL));
    for(int i=0; i<S; i++) //随机多少次
    {
        a=rand()%(n-1)+1;
        x=quick_mod_power(a,u,n);
        for(int j=0; j<t; j++)
        {
            y=quick_mod_multi(x,x,n);
            if(y==1&&x!=1&&x!=n-1)
                return false;
            x=y;
        }
        if(x!=1) return false;
    }
    return true;
}
const LL N = 1e6;
LL prime[N] = {0},num_prime = 0;
int isNotPrime[N] = {1, 1};
void prime2()
{
    for(LL i = 2 ; i < N ; i ++)
    {
        if(! isNotPrime[i])
            prime[num_prime ++]=i;
        //关键处1
        for(LL j = 0 ; j < num_prime && i * prime[j] <  N ; j ++)
        {
            isNotPrime[i * prime[j]] = 1;
            if( !(i % prime[j] ) )  //关键处2
                break;
        }
    }
}
int main()
{
    S=8;
    int n;
    prime2();
    scanf("%d",&n);
    {
        int sum=0;
        for(int i=0; i<n; i++)
        {
            LL a;
            scanf("%llu",&a);
            int flag = 0;
            for (int j = 0; j < num_prime; j++)
            {
                if(miller_rabbin(a - prime[j]))
                {
                    flag = 1;
                    printf("%llu %llu\n",prime[j],a-prime[j]);
                    break;
                }
            }
        }
    }
}