洛谷P3265 装备购买

这个大毒瘤题....居然反向卡精度....

别的题eps要开小，这个毒瘤要开大...

我一开始是1e-12，挂的奇惨无比，50分......

然后改成1e-7，就70分了...

1e-5 90分 1e-4 AC

做法就是用高斯消元的变式，消出一个上三角形，然后把代价加起来就行了。

每次选择花费最小的主元来消。

如果这个位置是0，就横着挪一格。

#include <cstdio>
#include <algorithm>
#include <cmath>
const int N = 505;
const double eps = 1e-4;

double a[N][N];
int m, n;

inline void out() {
    for(int i = 1; i <= n; i++) {
        for(int j = 1; j <= m; j++) {
            printf("%lf ", a[i][j]);
        }
        printf("%.0lf 
", a[i][0]);
    }
    puts("");
    return;
}

inline void Gauss() {
    int i, t = 1;
    for(i = 1; i <= n; t++, i++) {
        if(t > m) {
            break;
        }
        for(int j = i; j <= n; j++) {
            if(fabs(a[i][t]) < eps && fabs(a[j][t]) > eps) {
                std::swap(a[i], a[j]);
            }
            else if(fabs(a[j][t]) > eps && a[j][0] < a[i][0]) {
                std::swap(a[i], a[j]);
            }
        }
        //out();
        if(fabs(a[i][t]) < eps) {
            i--;
            continue;
        }
        for(int j = i + 1; j <= n; j++) {
            double k = a[j][t] / a[i][t];
            for(int p = t; p <= m; p++) {
                a[j][p] -= k * a[i][p];
            }
        }
        //out();
    }
    i--;
    printf("%d ", i);
    double ans = 0;
    for(int j = 1; j <= i; j++) {
        ans += a[j][0];
    }
    printf("%.0lf", ans);
    return;
}

int main() {
    scanf("%d%d", &n, &m);
    for(int i = 1; i <= n; i++) {
        for(int j = 1; j <= m; j++) {
            scanf("%lf", &a[i][j]);
        }
    }
    for(int i = 1; i <= n; i++) {
        scanf("%lf", &a[i][0]);
    }
    //out();
    Gauss();
    return 0;
}