BZOJ 2426: [HAOI2010]工厂选址

2426: [HAOI2010]工厂选址

Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 364 Solved: 248
[Submit][Status][Discuss]

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

[Submit][Status][Discuss]

枚举一下选择哪个新工厂，然后贪心的分配每个煤矿的分配即可得到当前的最有答案，更新全局最优解。

 #include <bits/stdc++.h>

 #define rep(i,x,y) for(int i=x;i<=y;++i)

 template <class T>

 inline void read(T &x)

 {

     static char c;

     x = ;

     c = getchar();

     while (c < )c = getchar();

     while (c > )x = x *  + c - '', c = getchar();

 }

 const int mxn = ;

 const int mxm = ;

 const int inf = 1E9 + ;

 int m, b, h, n;

 int A[mxm];

 int H[mxn];

 int C[mxn][mxm];

 int ans = inf, id = ;

 int seq[mxm], sub[mxm];

 inline bool cmp(int a, int b)

 {

     return sub[a] < sub[b];

 }

 inline void calculate(int t)

 {

     int cost = h + H[t], tmp = b;

     rep (i, , m)cost += A[i] * C[t][i];

     rep (i, , m)seq[i] = i, sub[i] = C[][i] - C[t][i];

     std::sort(seq + , seq + m + , cmp);

     rep (i, , m)if (tmp)

     {

         int q = seq[i];

         if (tmp >= A[q])

             cost += sub[q] * A[q], tmp -= A[q];

         else

             cost += sub[q] * tmp, tmp = ;

     }

     if (ans > cost)

         ans = cost, id = t;

 }

 signed main(void)

 {

     read(m);

     read(b);

     read(h);

     read(n);

     rep (i, , m)read(A[i]);

     rep (i, , n)read(H[i]);

     rep (i, , n)rep (j, , m)read(C[i][j]);

     rep (i, , n)calculate(i);

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

 }

@Author: YouSiki