【NOIP2016提高A组五校联考1】挖金矿

最新推荐文章于 2019-07-03 20:02:24 发布

无尽的蓝黄

最新推荐文章于 2019-07-03 20:02:24 发布

阅读量288

点赞数 1

分类专栏：题解二分

本文链接：https://blog.csdn.net/chen1352/article/details/52735022

版权

题解同时被 2 个专栏收录

177 篇文章 0 订阅

订阅专栏

二分

14 篇文章 0 订阅

订阅专栏

题目

这里写图片描述

分析

我们二分答案
设 $sum_{i,j}$ 表示的i列前个数的和，
假设当前出的二分答案为x，第i列挖了 $h_j$ 层，则

\sum n i = 1 s u m i , h i \sum n i = 1 h i > = x

$\dfrac{\sum_{i=1}^{n}sum_{i,h_i}}{\sum_{i=1}^{n}h_i}>=x$
转移得

∑ni=1sumi,hi−∑ni=1hix>=0 $\sum_{i=1}^{n}sum_{i,h_i}-\sum_{i=1}^{n}h_ix>=0$ 。
那么对于每一列，把最大的

sumi,h−hx $sum_{i,h}-hx$ 求出了，加在一起就可以了。

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const long long maxlongint=214748364700000;
const int mo=1000000007;
const int N=50005;
using namespace std;
int n,m,p,g,a[N*2];
long long sum[N*2];
int rf(double l,double r)
{
    while(l+0.000001<r)
    {
        double mid=(l+r)/2,num=0,p;
        for(int i=1;i<=n;i++)
        {
            p=-maxlongint;
            for(int j=1;j<=m;j++)
            {
                if(sum[(i-1)*m+j]-j*mid>p)
                {
                    p=sum[(i-1)*m+j]-j*mid;
                }
            }
            num+=p;
        }
        if(num<0)
            r=mid;
        else
            l=mid;
    }
    printf("%.4lf",l);
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            scanf("%d",&a[(i-1)*m+j]);
            if(j==1)
                sum[(i-1)*m+j]=a[(i-1)*m+j];
            else
                sum[(i-1)*m+j]=sum[(i-1)*m+j-1]+a[(i-1)*m+j];
        }
    }
    rf(1,maxlongint);
}