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;
}