BZOJ 4004 高斯消元

思路:
排个序 消元 完事~

但是!
坑爹精度毁我人生

我hhhh他一脸

红红火火恍恍惚惚

//By SiriusRen
#include <cmath>
#include <cstdio>
#include <algorithm>
using namespace std;
#define double long double
const int N=505;const double eps=1e-8;
int n,m,vis[N],ans2;double ans;
struct Node{double cost,eli[N];}node[N];
bool cmp(Node a,Node b){return a.cost<b.cost;}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            scanf("%Lf",&node[i].eli[j]);
    for(int i=1;i<=n;i++)scanf("%Lf",&node[i].cost);
    sort(node+1,node+1+n,cmp);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            if(abs(node[i].eli[j])>eps){
                if(!vis[j]){vis[j]=i,ans+=node[i].cost,ans2++;break;}
                else{
                    double temp=node[i].eli[j]/node[vis[j]].eli[j];
                    for(int k=j;k<=m;k++)
                        node[i].eli[k]-=temp*node[vis[j]].eli[k];
                }
            }
    printf("%d %d\n",ans2,(int)ans);
}

这里写图片描述

posted @ 2017-01-30 16:11  SiriusRen  阅读(113)  评论(0编辑  收藏  举报