UVALive - 6185 Find the Outlier暴力填表+高斯消元+卡eps

https://cn.vjudge.net/problem/UVALive-6185

我真的是服了orz eps 1e5，1e6过不了 开1e2 1e1都能过

题意：给你一个d阶多项式f的f(0),f(1)...f(d+1),f(d+2) 有一个是错误的，问第几个是错的

题解：题目多给了两个方程（约束）。

想了一下如果只给一个，是找不出来的。

给两个的话，可以这么考虑：

　　先取出一个方程X，再取剩下的n个高斯消元一下，将解得的系数带入最后一个方程，if成立，说明X是错的，else再取另一个（说明错误的在n个 或最后一个）。

坑：下标无限搞错，对着样例1写了一下，虽然具体一些，但有个d+1写成了3orz

#include <iostream>
#include <vector>
#include <cstdlib>
#include <algorithm>
#include <cstring>
#include <cmath>
#include<cstdio>
#include<vector>
#define rep(i,t,n)  for(int i =(t);i<=(n);++i)
#define per(i,n,t)  for(int i =(n);i>=(t);--i)
#define mmm(a,b) memset(a,b,sizeof(a))
using namespace std;
typedef long long ll;
const int maxn= ;
double a[][],aa[][], del;
double eps = 5e-;
/*高斯消元，系数矩阵为a[i][j],i=1…n,j=1…n,常数为a[i][n+1],i=1…n，返回值为是否有唯一解，答案存在a[i][n+1]*/
bool gauss(int n) {
    for (int i = ; i <= n; i++) {
        int k = i;
        for (int j = i + ; j <= n; j++)if (fabs(a[j][i]) > fabs(a[k][i]))k = j;
        if (fabs(del = a[k][i]) < eps)return ;
        for (int j = i; j <= n + ; j++)swap(a[i][j], a[k][j]);
        for (int j = i; j <= n + ; j++)a[i][j] /= del;
        for (k = ; k <= n; k++)if (k != i) {
            del = a[k][i];
            for (int j = i; j <= n + ; j++)a[k][j] -= a[i][j] * del;
        }
    }
    return ;
}

int main() {
    int  d;
    while (cin >> d&&d) {
        mmm(a, );
        mmm(aa, );
        rep(i, , d + ) {
            //cin >> aa[i][d + 2];
            scanf("%lf", &aa[i][d + ]);
            //aa[i][d + 2] = aa[i][d + 1];
        }aa[][] = ;

        rep(i, , d + ) {
            rep(j, , d+) {
                aa[i+][j] =pow(1.0*i, j-);
            }
        }
        rep(k, , d + ) {
            int tot = ;
            rep(i, , d + )if (i != k&&tot<=d+) {  rep(j, , d+) a[tot][j] = aa[i][j]; tot++;}
            if (!gauss(d + ))continue;
            double temp = ;
            int last = d + ;
            if (k == d + )last--;
            rep(i, , d+)
                temp += a[i][d + ] * aa[last][i];
             if (abs(temp - aa[last][d + ]) < 0.5) {
                cout << k -  << endl; break;
            }

        }

    }
}
/*
1
1
2
2

4
42.4715310246
79.5420238202
28.0282396675
-30.3627807522
-49.8363481393
-25.5101480106
7.58575761381
*/
/*
2
1.0
4.0
12.0
16.0
25.0
1
-30.5893962764
5.76397083962
39.3853798058
74.3727663177
4
42.4715310246
79.5420238202
28.0282396675
-30.3627807522
-49.8363481393
-25.5101480106
7.58575761381
5
-21.9161699038
-48.469304271
-24.3188578417
-2.35085940324
-9.70239202086
-47.2709510623
-93.5066246072
-82.5073836498
0*/