POJ 1877 Flooded!（二分）

Description
将一个区域分成r*c个方块，每个方块有有一个海拔（可正可负）。求当给区域注入指定容量的水时，水面的海拔是多少，以及被水淹没的方块占总方块数的百分比。每个方块的面积为100m^2，水的容量单位为立方米。
Input
多组用例，每组用例第一行为两个整数r和c表示区域行列数，然后是一个r*c的矩阵，矩阵元素为对应区域方块的海拔，最后是注入水的体积，以r=c=0结束输入
Output
对于每组用例，输出注入水后水面的海拔以及被水淹没方块占总方块数的百分比
Sample Input
3 3
25 37 45
51 12 34
94 83 27
10000
0 0
Sample Output
Region 1
Water level is 46.67 meters.
66.67 percent of the region is under water.
Solution
二分高度即可，注意精度和输出格式
Code

#include <iostream>
#include <cmath>
using namespace std;
#define MAXN 35
#define INF 9999999
#define esp 1e-12
int R,C,mCount=1,num;
double h[MAXN][MAXN],v,wl,per,low,high,mid,mv;
double min(double a,double b)//比较函数 
{
    return a<b?a:b;
}
int main()
{
    int r,c;
    while(scanf("%d%d",&R,&C)&&(R+C))
    {
        low=INF;//初始化 
        for(r=0;r<R;r++)
            for(c=0;c<C;c++)
            {
                scanf("%lf",&h[r][c]);
                low=min(low,h[r][c]);//更新二分下限 
            }
            scanf("%lf",&v);
            high=v/100.0+low;//二分上限 
            while(low+esp<high)//精度控制 
            {
                num=mv=0;//初始化,num为被水淹没方块数量,mv为按此高度淹没方块所需水量 
                mid=(low+high)/2.0;
                for(r=0;r<R&&mv+esp<v;r++)
                    for(c=0;c<C&&mv+esp<v;c++)
                        if(mid>h[r][c]+esp)
                        {
                            num++;
                            mv+=(mid-h[r][c])*100.0;
                        }
                if(mv+esp<v)//水量足够淹没方块,下限增 
                    low=mid;
                else//水量不够淹没方块,上限减 
                    high=mid;
            }           
            printf("Region %d\n",mCount++);//按格式输出 
            printf("Water level is %.2lf meters.\n",mid);//按格式输出 
            printf("%.2lf percent of the region is under water.\n\n",double(num)/double(R*C)*100.0);//按格式输出 
    }
    return 0;
}