[bzoj2426][贪心]工厂选址

最新推荐文章于 2019-03-13 21:45:06 发布

原创最新推荐文章于 2019-03-13 21:45:06 发布 · 349 阅读

0 ·

CC 4.0 BY-SA版权

bzoj 同时被 2 个专栏收录

482 篇文章

订阅专栏

贪心

22 篇文章

订阅专栏

本文探讨了一个涉及煤炭产量分配的问题，目标是最小化发电厂运行成本和煤炭运输成本的总和。通过对不同发电厂选址和煤炭分配策略进行分析，提出了一种贪心算法来解决这一优化问题。

Description

某地区有m座煤矿，其中第i号矿每年产量为ai吨，现有火力发电厂一个，每年需用煤b吨，每年运行的固定费用（包括折旧费，不包括煤的运费）为h元，每吨原煤从第i号矿运到原有发电厂的运费为Ci0（i=1，2，…，m）。

现规划新建一个发电厂，m座煤矿每年开采的原煤将全部供给这两座发电厂。现有n个备选的厂址。若在第j号备选厂址建新厂，每年运行的固定费用为hj元。每吨原煤从第i号矿运到j号备选厂址的运费为Cij（i=1，2，…，m；j=1，2，…，n）。

试问：应把新厂厂址选取在何处？m座煤矿开采的原煤应如何分配给两个发电厂，才能使每年的总费用（发电厂运行费用与原煤运费之和）为最小。

Input

第1行： m b h n

第2行： a1 a2 … am （0<=ai<=500, a1+a2+…+an>=b）

第3行： h1 h2 … hn （0<=hi<=100）

第4行： C10 C20 … Cm0 （0<=Cij<=50）

第5行： C11 C21 … Cm1
第n+4行：C1n C2n … Cmn

Output

第1行：新厂址编号，如果有多个编号满足要求，输出最小的。

第2行：总费用

Sample Input

4 2 7 9

3 1 10 3

6 3 7 1 10 2 7 4 9

1 2 4 3

6 6 8 2

4 10 8 4

10 2 9 2

7 6 6 2

9 3 7 1

2 1 6 9

3 1 10 9

4 2 1 8

2 1 3 4

Sample Output

8

49

HINT

对于所有数据, n<=50, m<=50000, b<=10000

题解

这题题意超级坑
我本来以为，是两个发电厂都是b吨矿的。。然后开心地敲了个费用流然后样例都过不去
结果第一个发电厂只要b吨其他全部运去第二个发电厂？？？
~~开心的TLE~~
那这样还打啥子费用流。。直接贪心啊，反正你要么运去第一个要么运去第二个
枚举发电厂位置，然后存一下每个矿到第二个厂与到第一个厂的差，按差从大到小排个序。差越大当然越要全部运去第一个厂啦。
运不来第一个就只能去第二个咯，贪个心就行了
告诫自已一定要细看题意

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
struct node
{
    int x,y,pos,g;
}a[51000];
bool cmp(node n1,node n2){return n1.pos>n2.pos;}
int m,b,h,n;
int u[51000],v[51000];
int mp[51][51000];
int mcost,mout;
int main()
{
    scanf("%d%d%d%d",&m,&b,&h,&n);
    for(int i=1;i<=m;i++)scanf("%d",&u[i]);
    for(int i=1;i<=n;i++)scanf("%d",&v[i]);
    for(int i=1;i<=n+1;i++)for(int j=1;j<=m;j++)scanf("%d",&mp[i][j]);
    mcost=999999999;mout=0;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)a[j].x=mp[1][j],a[j].y=mp[i+1][j],a[j].pos=a[j].y-a[j].x,a[j].g=u[j];
        sort(a+1,a+1+m,cmp);
        int sum=0,A=0;bool bk=false;
        for(int j=1;j<=m;j++)
        {
            if(A+a[j].g<=b)
            {
                A+=a[j].g;
                sum+=a[j].x*a[j].g;
            }
            else
            {
                int lea=b-A;
                A+=lea;sum+=a[j].x*lea;
                sum+=a[j].y*(a[j].g-lea);
            }
        }
        sum+=h;sum+=v[i];
        if(sum<mcost){mcost=sum,mout=i;}
    }
    printf("%d\n%d\n",mout,mcost);
    return 0;
}