UVA - 815 Flooded（浮点数二分）

题目链接：https://vjudge.net/problem/UVA-815

题意：有一个n*m（1<=n,m<30）的网格，每个格子是边长为10米的正方形，网格四周是无限高的墙壁。输入每个格子的海拔高度，以及网格内雨水的总体积，输出水位的海拔高度以及有多少百分比的区域有水（即高度严格小于水平面）。

解题方案：二分水位高度进行check，记所有网格中高度最低的格子的高度为a[0]，当没有水的时候水位高度应该是a[0]，此为水位高度的下限；当所有雨水都囤积到最低的格子处时，水位高度为a[0]+V/100.0，此为水位高度的上限，接下来在区间[a[0] , a[0]+V/100.0]进行二分check即可。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <string>
#include <cmath>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>

using namespace std;

#define FOR(i,k,n) for(int i=k;i<n;i++)
#define FORR(i,k,n) for(int i=k;i<=n;i++)
#define scan(a) scanf("%d",&a)
#define scann(a,b) scanf("%d%d",&a,&b)
#define scannn(a,b,c) scanf("%d%d%d",&a,&b,&c)
#define mst(a,n)  memset(a,n,sizeof(a))
#define ll long long
#define N 1005
#define mod 1000000007
#define INF 0x3f3f3f3f

const double eps=1e-8;
const double pi=acos(-1.0);

int n,m;
double a[N];
double V;

bool check(double x)
{
    double sum=0;
    FOR(i,0,n*m)
    {
        if(a[i]<x)
            sum+=100.0*(x-a[i]);
    }
    //printf("%.5lf\n",sum);
    if(sum>=V) return true;
    else return false;
}

double found(double l,double r)
{
    while(fabs(l-r)>=eps)
    {
        double mid=(l+r)/2;
        if(check(mid))
            r=mid;
        else
            l=mid;
    }
    return l;
}

int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int Case=0;
    while(cin>>n>>m&&(n||m))
    {
        FOR(i,0,n*m) cin>>a[i];
        cin>>V;
        sort(a,a+n*m);
        double ans=found(a[0],a[0]+V/100.0);
        int cnt=0;
        FOR(i,0,n*m)
        {
            if(a[i]<ans) cnt++;
        }
        printf("Region %d\nWater level is %.2lf meters.\n%.2lf percent of the region is under water.\n\n"
               ,++Case,ans,(double)cnt*100.0/((double)n*m));
    }
    return 0;
}