洛谷P3265 装备购买

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

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

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

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

1e-5 90分 1e-4 AC

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

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

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

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

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

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

 inline void Gauss() {
     int i, t = ;
     for(i = ; 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][] < a[i][]) {
                 std::swap(a[i], a[j]);
             }
         }
         //out();
         if(fabs(a[i][t]) < eps) {
             i--;
             continue;
         }
         for(int j = i + ; 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 = ;
     for(int j = ; j <= i; j++) {
         ans += a[j][];
     }
     printf("%.0lf", ans);
     return;
 }

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

AC代码