spoj Prime Generator

题意：判断ll-rr范围内的质数。

一个个用miller-rabin算法判断

//#pragma comment(linker,"/STACK:1024000000,1024000000")
#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<vector>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<algorithm>
#include <stack>
using namespace std;
typedef long long lon;
const lon SZ=,INF=0x7FFFFFFF;

lon pow(lon x,lon n,lon mod)
{
    lon res=,ele=x;
    for(;n;)
    {
        if(n&)res=ele*res%mod;
        ele=ele*ele%mod;
        n/=;
    }
    return res;
}

bool isp(lon x)
{
    if(x==)return ;
    int sz=;
    lon arr[sz]={,,};
    for(lon i=;i<sz;++i)
    {
        if(x==arr[i])return ;
        lon id=x-;
        for(;!(id&)&&id;id/=);
        lon m=pow(arr[i],id,x);
        if(m==)continue;
        for(;id!=x-&&m!=;)
        {
            //if(x==5)cout<<m<<endl;
            if(m==x-)break;
            m=m*m%x;
            //if(m!=1&&m!=x-1)return 0;
            id*=;
        }
        if(m!=x-)return ;
    }
    return ;
}

int main()
{
    std::ios::sync_with_stdio();
    //freopen("d:\\1.txt","w",stdout);
    lon casenum;
    cin>>casenum;
    for(lon time=;time<=casenum;++time)
    //for(;scanf("%d",&n)!=EOF;)
    {
        lon ll,rr;
        cin>>ll>>rr;
        for(lon i=ll;i<=rr;++i)
        {
            if(isp(i))cout<<i<<endl;
        }
        cout<<endl;
    }
    return ;
}