UVA

来源:互联网 发布:unity3d虚拟现实源码 编辑:程序博客网 时间:2024/06/16 02:43

题目链接: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 0x3f3f3f3fconst 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;}


0 0