出现free(): invalid next size (fast)错误，疑惑transform_coords函数内存分配问题

我最近碰到个棘手的问题：多次调用transform_coords函数后，程序直接抛出free(): invalid next size (fast)错误，完全搞不懂这个函数里为啥会有内存分配相关的Bug。先贴一下函数代码：

template <typename T, typename U> // x -> ax+b
void transform_coords(double a, double b, int n, const T* oldc, U* newc) {
    int i, j, q0, q;
    if (n<2) return;
    int *p = new int[n]; // Pascal's triangle
    p[0] = 1;
    for (i=1; i<n; ++i) p[i] = 0;
    for (i=0; i<n; ++i) newc[i] = 0;
    for (j=0; j<n; ++j) {
        for (i=0; i<n; ++i) {
            if (p[i]==0) break;
            double m = p[i];
            if (i!=j) m *= std::pow(a,j-i);
            if (i... // 代码此处被截断
        }
    }
}

问题根源排查 #

首先最扎眼的问题：你用new int[n]分配了堆内存，但从头到尾没有释放！每次调用这个函数都会泄漏n个int大小的内存，次数多了堆内存会被过度占用、碎片化，甚至被非法修改，最终触发glibc堆管理器的检测机制，抛出free(): invalid next size (fast)这类堆损坏报错。

另外从截断的代码推测，你后续逻辑里可能还存在数组越界访问的情况：比如循环中i或j不小心超出了0~n-1的范围，写入p数组之外的内存区域，直接破坏堆的元数据，这也会触发同样的错误。

修复建议 #

• 补上内存释放逻辑：在函数所有可能的退出路径前，加上delete[] p;。更稳妥的方式是用智能指针替代手动内存管理，彻底避免泄漏风险：

  // 替换原来的int* p = new int[n];
  std::unique_ptr<int[]> p(new int[n]);
  

  智能指针会在函数结束时自动释放内存，不用手动写delete。

• 检查数组访问边界：仔细核对所有涉及p、oldc、newc的循环索引，确保不会出现i >= n或j >= n的情况，比如内层循环的break条件是否能提前终止非法访问。

• 额外逻辑校验：注释标注的是线性变换x -> ax+b，但代码里用到了std::pow(a,j-i)，当j < i时会计算a的负次幂，这可能和预期逻辑不符，虽然和内存错误无关，但可以顺便排查下。

内容的提问来源于stack exchange，提问作者SU3