HDU 2012 素数判定

#include<bits/stdc++.h>
using namespace std;
int prime[3000];
void IsPrime()
{
    memset(prime,1,sizeof(prime));
    for(int i=2;i<sqrt(3000);i++)
    {
        if(prime[i])
        {
            for(int j=i*i;j<3000;j+=i)
                prime[j]=0;
        }
    }
}
int main()
{
    int x,y;
    IsPrime();
    while(cin>>x>>y)
    {
        if(x==0&&y==0)break;
        int flag=1;
        for(int i=x;i<=y;i++)
        {
            int id=i*i+i+41;//cout<<id<<endl;
            if(prime[id]==0)
            {
                flag=0;
                break;
            }
        }
        if(flag==0) cout<<"Sorry"<<endl;
        else cout<<"OK"<<endl;
    }
    return 0;
}