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

Description

Input

Output

答案，保留4位小数

Sample Input

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

Sample Output

4.4286

Data Constraint

n*m<=100000

Solution

很简单的一道题
设sum[i][j]表示第i列向下挖j行的前缀和
假设第i列挖h[i]行
设二分的答案是m，答案如果合法，那么

$${{\sum_{i=1}^n sum[i][h[i]]}\over{\sum_{i=1}^n h[i]}}>=m$$
移项
$$\sum_{i=1}^n sum[i][h[i]]-h[i]*x>=0$$
那么二分后枚举得到 $\sum_{i=1}^n sum[i][h[i]]-h[i]*x$，判断是否大于0即可

Code

#include<cstdio>
#include<algorithm>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define N 100100
#define db double
using namespace std;
int n,m,q,tot=0;
db f[2][N],a[N];
int main()
{
    scanf("%d%d",&n,&m);
    fo(i,1,n*m) scanf("%lf",&a[i]);
    db l=0.0,r=1000000000.0;
    for(;r-l>0.00001;)
    {
        db x=(l+r)/2.0,jy=0;
        fo(i,1,n)
        {
            db mx=-21474836477777,s=a[(i-1)*m+1];
            fo(j,1,m) mx=max(mx,s-j*x),s+=a[(i-1)*m+j+1];
            jy+=mx;
        }
        if(jy<0) r=x;else l=x;
    }
    printf("%.4lf",l);
}