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('\n')
 using namespace std;
 const int Safe=;
 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-)*p);
     return tmp<?tmp+p:tmp;
 }
 int pow(int a,int b,int p){
     int ans=;a%=p;
     for(int i=b;i;i>>=,a=mul(a,a,p))
         if(i&)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=;i<=s;i++){
         ans=mul(ans,ans,n);
         if(ans==&&p!=&&p!=n-) return true;
         p=ans;
     } if(ans!=)return true;return false;
 }
 bool MR(int n){
     if(n<=) return false;
     if(n==) return true;
     if(!(n&)) return false;
     int r=n-,s=;
     while(!(r&)) r>>=,s++;
     for(int i=;i<Safe;i++)
         if(check(rand()%(n-)+,n,r,s)) return false;
     return true;
 }
 inline int read(){
     int x=,sig=;char ch=getchar();
     while(!isdigit(ch)){if(ch=='-')sig=-;ch=getchar();}
     while(isdigit(ch))x=*x+ch-'',ch=getchar();
     return x*=sig;
 }
 inline void write(int x){
     if(x==){putchar('');return;}if(x<)putchar('-'),x=-x;
     int len=,buf[];while(x)buf[len++]=x%,x/=;
     for(int i=len-;i>=;i--)putchar(buf[i]+'');return;
 }
 bool first=true;
 void init(){
     srand(time());
     int n=read();
     for(int i=;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 ;}