【筛法打表】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;
}