OJ, VJ素数判定

对于表达式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

    # include<stdio.h>
    #include<math.h>
    int main()
    {
        int x, y, temp, i, j, n, flag = 1;
        while(scanf("%d%d", &x, &y) != EOF )
          {
            if(x == 0 && y == 0)
                break;
            if(x > y)
            {
                temp = x;
                x = y;
                y = temp;
            }
            flag = 1;
           for(i = x; i <= y; i++)
               {
                   for(j = 2; j <= sqrt(i*i+i+41); j++)
                    {
                        if((i*i+i+41) % j == 0)//写这里的时候i*i+i+41这个式子没加括号。相当于41和j在取余。结果代码看了半天，我吐了。
                        {
                            flag = 0;
                            printf("Sorry\n");
                            break;
                        }
                    }
                    if(flag == 0)
                        break;
                }
            if(i == y + 1)
            printf("OK\n");
          }
        return 0;
    }