@[rainygame](/user/804607) 读入问题,题目给了 $n+1$ 个点的坐标,不是按方程的格式给的。
by Flanksy @ 2023-04-12 23:52:47
另外造增广矩阵的代码写错了。
```cpp
for (int i=1; i<=n; i++){
for (int j=1; j<=n; j++) a[i][j] = 2*(b[i][1]-b[i+1][1]);//here
for (int j=1; j<=n; j++) a[i][n+1] += pow(a[i][j], 2)-pow(a[i+1][j], 2);
}
```
by Flanksy @ 2023-04-12 23:59:51
最后是上边的代码片段第三行,应该用原有信息更新增广矩阵最右边那一列的值,而不是增广矩阵自身的信息。
改完这三个就能过样例了。
by Flanksy @ 2023-04-13 00:02:57
@[Flanksy](/user/662295) 谢谢!已关
by rainygame @ 2023-04-13 06:21:14
输出难好像是用long double 的原因吧,我也是
by sjzez__chess @ 2023-05-03 19:58:27