HDOJ 2012 素数判定

#include<iostream>
#include<math.h>

using namespace std;

int judge(int n)
{
    int a=n*n+n+41;
    if((a%6!=5) && (a%6!=1))
    {
        return  0;
    }
    float b=sqrtf((float)a);
    for(int i=5;i<=b;i+=6)
    {
        if(a%i==0||a%(i+2)==0)
        {

            return 0;
        }

        else
            continue;
    }
    return 1;
}

int main()
{
    int x,y;
    while(cin>>x>>y&&x||y)
    {
        int cnt=1;
        for(int i=x;i<=y;i++)
        {
            cnt=judge(i);
            if(cnt==0)
                break;
        }

        if(cnt==1)
            cout<<"OK"<<endl;
        else
            cout<<"Sorry"<<endl;
    }
}