BZOJ2426 [HAOI2010]工厂选址

Description

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

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

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

Input

第1行： m b h n

第2行： a₁ a₂ … a_m（0<=a_i<=500, a₁+a₂+...+a_n>=b）

第3行： h₁ h₂ … h_n（0<=h_i<=100）

第4行： C₁₀ C₂₀ … C_m0（0<=Cij<=50）

第5行： C₁₁ C₂₁ … C_m1

… …

第n+4行：C_1n C_2n … C_mn

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

Source

Day2

正解：贪心

解题报告：

　　直接一遍贪心过去，跟上午考试的T1的apple几乎是一样的。

 //It is made by jump~

 #include <iostream>

 #include <cstdlib>

 #include <cstring>

 #include <cstdio>

 #include <cmath>

 #include <algorithm>

 #include <ctime>

 #include <vector>

 #include <queue>

 #include <map>

 #include <set>

 using namespace std;

 typedef long long LL;

 const int inf = (<<);

 const int MAXM = ;

 const int MAXN = ;

 int m,b,hh,n;

 int ans,jilu;

 int mei[MAXM],h[MAXN];

 int c[MAXM][MAXN];

 struct node{

     int x,y,z;

     int num;

 }a[MAXM];

 inline bool cmp(node q,node qq){ return q.z<qq.z; }

 inline int getint()

 {

     int w=,q=; char c=getchar();

     while((c<'' || c>'') && c!='-') c=getchar(); if(c=='-') q=,c=getchar();

     while (c>='' && c<='') w=w*+c-'', c=getchar(); return q ? -w : w;

 }

 inline void work(){

     m=getint(); b=getint(); hh=getint();  n=getint(); int now,bb,j; ans=inf;

     for(int i=;i<=m;i++) mei[i]=getint(); for(int i=;i<=n;i++) h[i]=getint();

     for(int i=;i<=n;i++) for(int j=;j<=m;j++) c[j][i]=getint();

     for(int i=;i<=n;i++) {

     now=h[i]+hh; bb=b;

     for(j=;j<=m;j++) a[j].num=mei[j];

     for(j=;j<=m;j++) a[j].y=c[j][i],a[j].x=c[j][],a[j].z=a[j].x-a[j].y;

     sort(a+,a+m+,cmp);

     for(j=;j<=m;j++) {

         if(a[j].num<=bb) {

         now+=a[j].num*a[j].x;

         bb-=a[j].num;

         }

         else { now+=bb*a[j].x,a[j].num-=bb; bb=; break; }

     }

     for(;j<=m;j++) now+=a[j].num*a[j].y;

     if(now<ans) ans=now,jilu=i;

     }

     printf("%d\n%d",jilu,ans);

 }

 int main()

 {

     work();

     return ;

 }

秒客网