c语言求素数的有关问题

c语言求素数的问题
用下面这个程序可以得出正确的结果:
#include <stdio.h>
#include <math.h>

int main(void)
{
int k;
int n, i;

printf("请输入您要求得的数：\n");
scanf("%d",&n);

k = sqrt(n + 1);

for (i = 2; i <= k; i++)
{
if (n%i == 0) break;
}
if (i>=k+1) printf("%d 是一个素数！\n",n);
else printf("%d 不是一个素数！\n",n);

return 0;
}

但是当k定义为double的时候，结果是错误的，比如输入17，就会显示不是素数，这是什么问题呢？


------解决方案--------------------
看错了，这个问题是发生在计算sqrt和比较时的错误。
1、计算sqrt时，由于k是double，则k的值可能会含有小数，比如5.00001。
2、i > k后，不一定i>k+1。因为k是double，在比较时i被提升为double。比如5>5.00001，但是5<5.00001+1。
------解决方案--------------------
如果n=17,则k=4.24
所以在出循环时i=5
则i<k+1所以显示为不是素数。

而4楼的数
如果n=17,则k=4.24
所以在出循环时i=5
则i>k所以显示为是素数。

所以i在等于17时不论k为sqrt(17)还是sqrt(18)民，n的值都在k与k+1之间。
------解决方案--------------------
本来在判断素数的时候，从2判断到sqrt(n)就行了，你判断到了sqrt(n+1)，然后最后在比较的时候，你又判断i≥k+1，两次都超出一点，最后完全有可能走到else部分。

C/C++ code

/*
 * 打印100以内的素数，每行显示5个。
 */

#include <math.h>
#include <stdio.h>

int isPrime(int i);

int main(void)
{
    int c = 0, i;

    for (i=1; i<=100; i++)
        if (isPrime(i))
        {
            printf("%3d", i);
            c++;

            if (c%5 == 0)
                printf("\n");

/*            if (c == 5)
            {
                printf("\n");
                c = 0;
            }
*/
        }

    return 0;
}

int isPrime(int i)
{
    int x;
    int y = (int)ceil(sqrt(i));  /* int y = (int)(sqrt(i + 1)); */

    if (i == 2)
        return 1;

    if (i%2==0 || i<2)  /* 过滤所有偶数及既不是素数又不是合数的情况 */
        return 0;

    for (x=3; x<=y; x++)
        if (i%x == 0)
            return 0;

    return 1;
}