HUD2012素数判定

Problem Description
对于表达式n^2+n+41，当n在（x,y）范围内取整数值时（包括x,y）(-39<=x<y<=50)，判定该表达式的值是否都为素数。

Input
输入数据有多组，每组占一行，由两个整数x，y组成，当x=0,y=0时，表示输入结束，该行不做处理。

Output
对于每个给定范围内的取值，如果表达式的值都为素数，则输出"OK",否则请输出“Sorry”,每组输出占一行。

Sample Input
0 1
0 0

Sample Output
OK

Author
lcy`

#include<iostream>
#include<math.h>
#include<string>
#include<time.h>
#include<cstring>
using namespace std;
const int Max=10000000;
bool primes[Max]={false};//假设全为质数
int prime[Max]={0};//假设在[0,n]的质数个数为0
void count_Prime(int n){
    primes[1]=true;
    prime[1]=1;//1为质数，在[0,1]内质数个数为1
    int i,j;
    for(i=2;i<=n;i++){
        prime[i]=prime[i-1];
        //[0,i-1]内的质数个数+第i个数是否为质数（是+1，不是+0）的总和为[0,i]内的质数个数
        if(primes[i]==false){//若第i个数为质数
            prime[i]++;//[0,i-1]内的质数个数+1
            for(j=i*2;j<=n;j+=i)//筛选出质数的相关合数
                primes[j]=true;
        }
    }
}
int main(){
    int x,y;
    count_Prime(2600);
    while(cin>>x>>y){
        if(x==0&&y==0)
            break;
        int p=1;
        for(int i=x+1;i<y;i++){
            int k=i*i+i+41;
            if(primes[k]==true)
                p=-1;
        }
        if(p==-1)
            cout<<"Sorry"<<endl;
        else
            cout<<"OK"<<endl;
    }
    return 0;
}