poj 1595

先用筛选法生成大范围内的素数，关键要注意输出的要求，具体的看题目；

http://poj.org/problem?id=1595

#include<iostream>
#include<cmath>
const int maxn=1222; 
using namespace std;
int L,U; 
int p[maxn],sz; 
int flag[maxn];
void sieve(int L,int U)
{
        int d=U-L+1,i,limit=(int)sqrt((double)U);
        
                  for(i=0;i<d;i++)      
                          flag[i]=true;
                  for(i=L&0x1;i<d;i+=2) 
                          flag[i]=false;
                  for(i=3;i<=limit;i+=2)
                  {
                          if(i>L&&flag[i-L]==false)
                                  continue;
                          
                          int j=L/i*i;
                          if(j<L)   j+=i;
                          if(j==i)  j+=i; j=j-L;
                          for(;j<d;j+=i)  
                                  flag[j]=false;
                  }
                  if(L<=1)  flag[1-L]=false;
                  if(L<=2)  flag[2-L]=true;
                  sz=0;
                  for(i=1;i<d;i++)
                  {       p[0]=1;  if(flag[i]==true) p[++sz]=i+L;}
}

int main()
{
        
        int i,n,j,c;
        
        while(cin>>n>>c)
        {
                sieve(1,n);
                sz++;
                cout<<n<<" "<<c<<": ";
                if(n==2&&c==1) { cout<<"2"<<endl<<endl;continue;}
        if(c>n||2*c>n||2*c>=sz)
                {
                        for(i=0;i<sz-1;i++)
                                cout<<p[i]<<" ";
                        cout<<p[sz-1]<<endl<<endl;
                }
                else
                {
                        
                        if(sz%2==0)
                        {
                                
                                int t=(sz-2*c)/2;
                                for(i=0,j=0;j<2*c-1;j++,i++)
                                        cout<<p[i+t]<<" ";
                                cout<<p[t+2*c-1]<<endl<<endl;
                        }
                        else
                        {
                                int t=(sz-2*c+1)/2;
                                for(i=0,j=0;j<2*c-2;j++,i++)
                                        cout<<p[i+t]<<" ";
                                cout<<p[t+2*c-2]<<endl<<endl; 
                        }
                }
        }
        
        return 0;
}