BZOJ 1013: [JSOI2008]球形空间产生器sphere

/* ???BZOJ 1013: [JSOI2008]球形空间产生器sphere ???高斯消元 ???QAQ SB的我也能终于能秒题了啊 ???设球心的坐标为(x,y,z...) ???那么就可以列n+1个方程，化化式子高斯消元即可*/#include <cstdio>#include <iostream>#include <cstring>#define rg register#define Max 50typedef double flo; flo v[Max][Max], s[Max];inline flo P (flo x) { return x * x; }inline void swap (flo &a, flo &b) { flo c = a; a = b, b = c; }int main (int argc, char *argv[]){ ???int N, p; scanf ("%d", &N); rg int i, j, k; flo x; ???for (i = 1; i <= N; ++ i) scanf ("%lf", &s[i]); ???for (i = 1; i <= N; ++ i) ???????for (j = 1; j <= N; ++ j) ???????????scanf ("%lf", &x), v[i][j] = (x - s[j]) * 2.0, v[i][N + 1] += P (x) - P (s[j]); ???memset (s, 0, sizeof s); ???for (i = 1; i <= N; ++ i) ???{ ???????for (p = i, j = i + 1; j <= N; ++ j) if (v[j][i] > v[p][i]) p = i; ???????if (p != i) for (j = i; j <= N + 1; ++ j) swap (v[p][j], v[i][j]); ???????for (j = i + 1; j <= N; ++ j) ???????????for (k = i + 1, x = v[i][i] / v[j][i]; k <= N + 1; ++ k) ???????????????v[j][k] = v[i][k] - x * v[j][k]; ???} ???for (i = N; i >= 1; -- i) ???{ ???????for (j = i + 1; j <= N; ++ j) v[i][N + 1] -= s[j] * v[i][j]; ???????s[i] = v[i][N + 1] / v[i][i]; ???????} ???for (i = 1; i < N; ++ i) printf ("%.3lf ", s[i]); printf ("%.3lf", s[i]); return 0;}