JZOJ 4809. 挖金矿

最新推荐文章于 2019-07-08 16:44:12 发布

Luckfort

最新推荐文章于 2019-07-08 16:44:12 发布

阅读量439

点赞数 3

本文链接：https://blog.csdn.net/huangjingyuan107/article/details/52738737

版权

Problem

Description

这里写图片描述

Input

这里写图片描述

Output

这里写图片描述

Sample Input

4 3
4 3 3
5 1 6
2 6 1
3 2 9

Sample Output

4.4286

Data Constraint

这里写图片描述

Hint

这里写图片描述

Solution

这题显然二分。
令 $sum[i][...]$ 表示第i行的前缀和，第i行挖 $h[i]$ 层，二分答案为 $mid$ ，合法情况当且仅当

Σ n i = 1 s u m [ i ] [ h [ i ] ] Σ n i = 1 h [ i ] - x 非 负

$\frac{\Sigma_{i=1}^n sum[i][h[i]]}{\Sigma_{i=1}^n h[i]}-x非负$
那么我们很容易得到一条式子：

Σ n i = 1 (s u m [i] [h [i]] - h [i] \cdot x) \geq 0

$\Sigma_{i=1}^n (sum[i][h[i]]-h[i]·x)≥0$
那么我们只要算一下

Σ n i = 1 m a x (s u m [i] [j] - j x) 是 否 \geq 0 (1 \leq j \leq h)

$\Sigma_{i=1}^n max(sum[i][j]-jx)是否≥0 (1≤j≤h)$ 就OK了。
由于n，h不确定，所以二维的sum数组要改成一维。

Code

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#define N 100010
#define eps 0.00001
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
int i,j,n,h,tot;
double l,r,mid,x,ans,sum,mx,a1,a2,a[N],s[N];
int main()
{
    freopen("4809_1.in","r",stdin);
    scanf("%d%d",&n,&h);
    fo(i,1,n*h)
    {
        scanf("%lf",&a[i]);
        sum+=a[i];
        s[i]=sum;
        if (i%h==0) sum=0;
    }
    l=0;r=1000000000;
    while (l+eps<r)
    {
        mid=(l+r)/2;
        sum=0;
        fo(i,1,n)
        {
            mx=-2147483647;
            fo(j,1,h)
                mx=max(mx,s[(i-1)*h+j]-j*mid);
            sum+=mx;
        }
        if (sum>=0) l=mid;else r=mid;
    }
    printf("%.4lf",l);
}