传送门:UVA-815
二分法,水位最低在最低的格子,最高不超过maxn+v/(n*m)/100+1,其中maxn表示最高的格子的高度,v代表水体积。最后保留两位小数,精度eps设为1e-4足够了
AC代码:
#include<iostream>
#include<cmath>
using namespace std;
const double eps=1e-4;
int main()
{
int v,n,m,pq[40][40],num=1;//v表示能留住的水的总体积,num代表数据集的次序
while(cin>>m>>n){
int minn=(1<<30),maxn=0;
if(!m) break;
for(int i=0;i<m;++i)
for(int j=0;j<n;++j){
cin>>pq[i][j];
if(pq[i][j]<minn) minn=pq[i][j];
if(pq[i][j]>maxn) maxn=pq[i][j];//这里不能写成else if,否则只有一个格子时会出错
}
cin>>v;
double i=minn,j=maxn+v/(n*m)/100+1,pos;//二分法,i代表左边界,j代表右边界,pos代表水位
while(j>i&&fabs(i-j)>eps){
double sum=0;
pos=(i+j)/2.0;
for(int p=0;p<m;++p)
for(int q=0;q<n;++q)
if(pos>pq[p][q]) sum+=pos-pq[p][q];
sum*=100;
if(sum>v) j=pos;
else i=pos;
}
int cnt=0;
for(int p=0;p<m;++p)//统计在水平面下的格子
for(int q=0;q<n;++q)
if(pq[p][q]<pos) ++cnt;
double tmp=100.0*cnt/(n*m);
printf("Region %d\nWater level is %.2lf meters.\n",num++,pos);
printf("%.2lf percent of the region is under water.\n\n",tmp);
}
return 0;
}