BZOJ 4004 [JLOI2015]装备购买 拟阵+极大线性无关组

题意:链接

方法:拟阵+极大线性无关组

解析:题解同上上

这题就是消元的部分变了一下而已..

原来是01串现在变消矩阵了.

代码:

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 510
#define eps 1e-5
using namespace std;
int n,m;
struct node
{
    double z[N];
    double cost;
}a[N];
int b[N];
int cmp(node x,node y)
{
    return x.cost<y.cost;
}
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].z[j]);
        }
    }
    for(int i=1;i<=n;i++)scanf("%lf",&a[i].cost);
    sort(a+1,a+n+1,cmp);
    int cnt=0;
    double ans=0;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            if(fabs(a[i].z[j])>eps)
            {
                if(!b[j])
                {
                    b[j]=i;
                    break;
                }else
                {
                    double rate=a[i].z[j]/a[b[j]].z[j];
                    for(int k=1;k<=m;k++)
                    {
                        a[i].z[k]-=rate*a[b[j]].z[k];   
                    }
                }
            }
        }
        int flag=0;
        for(int j=1;j<=m;j++)
        {
            if(fabs(a[i].z[j])>eps){flag=1;break;}
        }
        if(flag)
            cnt++,ans+=a[i].cost;
    }
    printf("%d %.0lf\n",cnt,ans);
}