【筛法打表】poj1595 Prime Cuts

传送门：http://poj.org/problem?id=1595

题目大意：给定一个n，求出1..n的所有质数（包括1），然后输出中间的2*c个（质数的个数为偶数）/2*c-1个（质数的个数为奇数）质数。

思路：由于数据范围极小，随便怎么搞都行。这里先打个表，然后找出小于n的质数，然后输出规定个数的质数。

注意：要求的个数超过实际个数时，直接输出实际的；打表的话不要只打1000以内的，多打几个，方便判断小于n的质数（不小心在这里wa了一发）。

#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std;

int prime[1000],cnt,boo[1100+10],n,c;

int main(){
    memset(boo,0,sizeof(boo));
    cnt=1;
    prime[cnt]=1;
    for (int i=2;i<=1000+100;i++){
        if (!boo[i]){
            cnt++;
            prime[cnt]=i;
            for (int j=i+i;j<=1000+100;j+=i)
                boo[j]=1;
        }
    }
    //for (int i=1;i<=100;i++) cout<<boo[i]<<" ";

    while (scanf("%d%d",&n,&c)!=EOF){
        int pos=0;

        cout<<n<<" "<<c<<":";

        for (int i=1;i<=cnt;i++){
            if (prime[i]>n) {
                pos=i-1;
                break;
            }
        }
//算法问题，需要对prime的计算大于1000
        if (pos%2==0){
            if (c*2>=pos){
                for (int i=1;i<=pos;i++) cout<<" "<<prime[i];
                cout<<endl<<endl;
            }
            else {
                for (int i=pos/2-c+1;i<=pos/2+c;i++) cout<<" "<<prime[i];
                cout<<endl<<endl;
            }
        }
        else {
            if (c*2-1>=pos){
                for (int i=1;i<=pos;i++) cout<<" "<<prime[i];
                cout<<endl<<endl;
            }
            else {
                for (int i=(pos+1)/2-c+1;i<=(pos+1)/2+c-1;i++) cout<<" "<<prime[i];
                cout<<endl<<endl;
            }
        }
    }
    return 0;
}