C语言中数组排序时丢失数据问题求助

嘿，我太懂这种挫败感了——自己写的shell_sort和标准库qsort都掉链子，肯定不是排序算法本身的锅，得从float数组的细节操作里找问题。咱们一步步来排查：

最可能的几个原因 #

• qsort比较函数的致命错误：很多人写float的比较函数时会直接返回*(float*)a - *(float*)b，但这个写法有大坑——当两个float的差值小于1时，会被截断成int类型的0，导致qsort误以为这两个值相等，排序逻辑直接混乱，看起来就像数据丢了。
• 自定义shell_sort的数组越界：如果你的shell_sort里步长计算、循环边界写错了，很容易访问到数组外面的内存，把原本的有效数据给冲掉，这也是排序后数据丢失的常见原因。
• float精度不够：float只有6-7位有效数字，如果你的变量边界数值差异极大（比如一个是1e6，一个是1e-6），精度丢失会让相近的数值被“合并”，看起来像数据丢了。

针对你的代码的排查步骤 #

先假设你的代码大概是这样（如果和实际有出入，你可以对应调整）：

自定义shell_sort示例（可能存在的问题） #

void shell_sort(float arr[], int n) {
    int gap, i, j;
    float temp;
    for (gap = n/2; gap > 0; gap /= 2) {
        for (i = gap; i < n; i++) {
            temp = arr[i];
            // 这里要注意j >= gap的条件，会不会出现j - gap为负？
            for (j = i; j >= gap && arr[j - gap] > temp; j -= gap) {
                arr[j] = arr[j - gap];
            }
            arr[j] = temp;
        }
    }
}

qsort比较函数的错误写法 #

// 错误！差值小于1时会被截断为0，导致排序逻辑错误
int compare_float(const void *a, const void *b) {
    return *(float*)a - *(float*)b;
}

那你可以先做这几件事：

1. 修正qsort的比较函数：必须明确判断大小关系，返回正确的int值：

   int compare_float(const void *a, const void *b) {
       const float fa = *(const float*)a;
       const float fb = *(const float*)b;
       if (fa < fb) return -1;
       if (fa > fb) return 1;
       return 0;
   }
   

   这是qsort处理float时最容易踩的坑，先把这个改了试试。

2. 检查shell_sort的循环边界：在排序过程中加个打印，或者用断言验证j - gap不会小于0，比如在循环里加assert(j - gap >= 0);（记得包含<assert.h>），如果触发断言，说明越界了。

3. 验证排序前的数据正确性：写个辅助函数打印数组所有元素，确认排序前数据都是对的——有时候数据丢失根本不是排序的问题，是之前赋值的时候就出错了（比如把int直接转float没处理，或者读取数据时的错误）：

   void print_float_arr(float arr[], int len) {
       for (int i = 0; i < len; i++) {
           printf("arr[%d] = %.10f\n", i, arr[i]);
       }
       puts("---");
   }
   

   排序前调用一次，排序后再调用一次，对比每一个元素的变化，就能精准定位是哪个环节出问题了。

4. 测试double类型是否正常：如果你的数值精度要求高，换成double数组试试，如果排序正常，那就是float的精度不足以存储你的变量边界数据，得换类型。

内容的提问来源于stack exchange，提问作者Alan Garduño III