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);
}

这里写图片描述

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值