c - SIGSEGV报道使用pow()

时间:2013-01-29 09:03:06

标签: c debugging sigsegv

编译器在tmp = pow(p[i],j);处投了一个SIGSEGV而p[i]j是两个整数,p[i]是数组p的有效元素,我真的不喜欢不知道为什么......

原始代码在此处:http://pastebin.com/DYhGeHxm

2 个答案:

答案 0 :(得分:6)

你是否想过这样:

int i,j,p[2000], a[5000000],num,count,tmp;

可能会推动你离你很近,或者直接 over 你的堆栈空间的边缘?那是

4 + 4 + 8000 + 20000000 + 4 + 4 + 4

<强>字节

即。你有一个 19.08兆字节堆栈空间声明。考虑至少动态分配a。当我把它改为:

int *a = malloc(5000000 * sizeof(*a));

并重新编写代码,它已经超越了你所拥有的seg-fault。不幸的是,它在这个地方死了:

count = 0;
for(i = 0; i < num; i++) {
    for(j = 2; ;j++) {
        tmp = pow(p[i],j);
        if(tmp > 5000000) break;
        a[count++] = tmp; // <=== faulted here, count was 5000193
    }
}

当达到分配的最大大小a[]时,两个循环都应该中断。我做了以下。位于main()的顶部:

static const int a_max = 5000000;
int *a = malloc(a_max*sizeof(*a));

在循环中:

count = 0;
for(i = 0; i < num && count < a_max; i++)
{
    for(j = 2; count < a_max; j++)
    {
        tmp = pow(p[i],j);
        if(tmp > 5000000)
            break;
        a[count++] = tmp;
    }
}

这会让你通过所有设置。最后一件事是快速排序算法本身,它也显得很破旧。我高度建议从较小的数据大小开始调试。

祝你好运。


编辑如果你需要一个参考快速排序算法,我有一个坐在我的一个垃圾文件夹的源文件中。不保证它甚至是正确的(非常确定它是,并且跳过长度-1也排序),但我知道它不会挂起,所以它有它为它= P

// quicksort for ints
static void quicksort_ints(int *arr, int left, int right)
{
    int p = (left+right)/2;    // as good as any
    int l = left, r = right;   // movable indicies

    while (l <= r)
    {
        while (arr[l] < arr[p])
            ++l;
        while (arr[r] > arr[p])
            --r;
        if (l <= r)
        {
            int tmp = arr[l];
            arr[l] = arr[r];
            arr[r] = tmp;
            ++l;
            --r;
        }
    }

    if (left < r)
        quicksort_ints(arr, left, r);
    if (l < right)
        quicksort_ints(arr, l, right);
}

答案 1 :(得分:0)

检查i的值。对p数组来说可能太大了。您可以使用调试器执行此操作。