COJ 3018 求1~n之间的素数

求1~n之间的素数

难度级别：A； 运行时间限制：1000ms； 运行空间限制：256000KB； 代码长度限制：2000000B

试题描述

 素数是大于1，且除1和本身以外不能被其他整数所整除的数。要求输出1~n之间的素数。

输入

正整数n

输出

1~n之间的所有素数，包括n，数字之间用一个空格隔开，第一个数字前不能有空格。

输入示例

10

输出示例

2 3 5 7

其他说明

n是大于1且不大于100的正整数

题解：可以练习一下MR（都不会打了）

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<queue>
#include<cstring>
#include<ctime>
#define PAU putchar(' ')
#define ENT putchar('
')
using namespace std;
const int Safe=3;
int gcd(int a,int b){return !b?a:gcd(b,a%b);}
int mul(int a,int b,int p){
    int tmp=(a*b-(int)((double)a/p*b+1e-8)*p);
    return tmp<0?tmp+p:tmp;
}
int pow(int a,int b,int p){
    int ans=1;a%=p;
    for(int i=b;i;i>>=1,a=mul(a,a,p))
        if(i&1)ans=mul(ans,a,p);
    return ans;
}
bool check(int a,int n,int r,int s){
    int ans=pow(a,r,n),p=ans;
    for(int i=1;i<=s;i++){
        ans=mul(ans,ans,n);
        if(ans==1&&p!=1&&p!=n-1) return true;
        p=ans;
    } if(ans!=1)return true;return false;
}
bool MR(int n){
    if(n<=1) return false;
    if(n==2) return true;
    if(!(n&1)) return false;
    int r=n-1,s=0;
    while(!(r&1)) r>>=1,s++;
    for(int i=0;i<Safe;i++)
        if(check(rand()%(n-1)+1,n,r,s)) return false;
    return true;
}
inline int read(){
    int x=0,sig=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')sig=-1;ch=getchar();}
    while(isdigit(ch))x=10*x+ch-'0',ch=getchar();
    return x*=sig;
}
inline void write(int x){
    if(x==0){putchar('0');return;}if(x<0)putchar('-'),x=-x;
    int len=0,buf[15];while(x)buf[len++]=x%10,x/=10;
    for(int i=len-1;i>=0;i--)putchar(buf[i]+'0');return;
}
bool first=true;
void init(){
    srand(time(0));
    int n=read();
    for(int i=2;i<=n;i++){
        if(MR(i)){
            if(first) first=false;
            else PAU;
            printf("%d",i);
        }
    }
    return;
}
void work(){
    return;
}
void print(){
    return;
}
int main(){init();work();print();return 0;}