【BZOJ2426】[HAOI2010]工廠選址（貪心）

題面

BZOJ
洛谷

題解

首先看懂題目到底在做什麽。
然而發現我們顯然可以對於每個備選位置跑一遍費用流，然後並不夠優秀。
不難發現所有的位置都要分配給兩個工廠，而其中一個工廠的用量是\(b\)。那麽我們先假裝把所有的全部分配給這一個工廠，這樣子我們每次把一些分給另外一個工廠的時候，對於答案的貢獻就已知了，那麽從貢獻小的往貢獻大的貪心即可。
還是記住這樣一句話，我們的貪心過程從某種意義上來說就是在模擬費用流的過程。
這樣子的復雜度就是\(O(nmlog)\)了。

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
#define ll long long
#define MAX 50500
inline int read()
{
    int x=0;bool t=false;char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')t=true,ch=getchar();
    while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
    return t?-x:x;
}
int m,b,H,n,a[MAX],h[MAX],c[MAX],c0[MAX];
int ans,sum=2e9;
int all,p[MAX],v[MAX];
bool cmp(int a,int b){return v[a]<v[b];}
int main()
{
    m=read();b=read();H=read();n=read();
    for(int i=1;i<=m;++i)a[i]=read();
    for(int i=1;i<=m;++i)all+=a[i];
    for(int i=1;i<=n;++i)h[i]=read();
    for(int i=1;i<=m;++i)c0[i]=read();
    for(int i=1;i<=m;++i)p[i]=i;
    for(int i=1;i<=n;++i)
    {
        for(int j=1;j<=m;++j)c[j]=read();
        int tot=H+h[i];
        for(int j=1;j<=m;++j)tot+=c0[j]*a[j];
        for(int j=1;j<=m;++j)v[j]=c[j]-c0[j];
        sort(&p[1],&p[m+1],cmp);
        int ht=all-b;
        for(int j=1;j<=m;++j)
            if(ht>a[p[j]])ht-=a[p[j]],tot+=v[p[j]]*a[p[j]];
            else{tot+=v[p[j]]*ht;break;}
        if(tot<sum)sum=tot,ans=i;
    }
    printf("%d\n%d\n",ans,sum);
    return 0;
}
 

【BZOJ2426】[HAOI2010]工廠選址（貪心）