uvA Flooded!

Description

To enable homebuyers to estimate the cost of flood insurance, a real-estate firm provides clients with the elevation of each 10-meter by 10-meter square of land in regions where homes may be purchased. Water from rain, melting snow, and burst water mains will collect first in those squares with the lowest elevations, since water from squares of higher elevation will run downhill. For simplicity, we also assume that storm sewers enable water from high-elevation squares in valleys (completely enclosed by still higher elevation squares) to drain to lower elevation squares, and that water will not be absorbed by the land.  From weather data archives, we know the typical volume of water that collects in a region. As prospective homebuyers, we wish to know the elevation of the water after it has collected in low-lying squares, and also the percentage of the region's area that is completely submerged (that is, the percentage of 10-meter squares whose elevation is strictly less than the water level). You are to write the program that provides these results. 

Input

The input consists of a sequence of region descriptions. Each begins with a pair of integers, m and n, each less than 30, giving the dimensions of the rectangular region in 10-meter units. Immediately following are m lines of n integers giving the elevations of the squares in row-major order. Elevations are given in meters, with positive and negative numbers representing elevations above and below sea level, respectively. The final value in each region description is an integer that indicates the number of cubic meters of water that will collect in the region. A pair of zeroes follows the description of the last region.

 

Output

For each region, display the region number (1, 2, ...), the water level (in meters above or below sea level) and the percentage of the region's area under water, each on a separate line. The water level and percentage of the region's area under water are to be displayed accurate to two fractional digits. Follow the output for each region with a blank line.

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.

先对所有海拔排序，从小到大
然后，选取中间的海拔作为最终填满T1体积水后的海拔高度，判断其与T的关系，若T1>T，选取海拔较小的海拔，填满T2体积水。
直到得到两相邻海拔高度的T1,T2，使得（T1-T)*(T2-T)<0，则说明最后的水的高度在T1和T2之间。
考虑到数据个能比较多，采用二分法做该题。
但是，最后一直 结果错误，以下为采用二分法的错误代码
情况很多，细节也很多。

#include <iostream>
#include <algorithm>
#include <string.h>
#include <cmath>
#include <iomanip>
using namespace std;
const int maxn = *+;
int s[maxn], sum[maxn],m,n;    //s为各个网格的海拔高度，有正有负

void dataout(double a1,double a2,double ss,double numm,double T)
{
    //cout<<a1<<" "<<a2<< " "<<ss<< " "<<numm<<endl;
    if(numm == )numm=;
    double height = ss + (a1 - a2) / numm , per;
    if(T == )per=;
    else per =  * numm / (m*n) ;
    cout<<setiosflags(ios::fixed)<<setprecision();
    cout << "Water level is " << height << " meters."<<endl;
    cout << per << " percent of the region is under water."<<endl;
}

void deal2(double T)
{
    sort(s, s + m * n);
    int a = , b = n * m - ,mid;     //用二分法从中间向两边计算
    while(){
            mid = (a + b) / ;
            if(sum[mid]!=)sum[mid]=;
            for(int i = ;i < mid; i++){
                if(s[i] < s [mid]) {
                    sum[mid] += (s[mid] - s[i]);
                }
            }
            if (sum[mid] < T) a = mid;
            else b = mid;
            if(abs(a-b)<=)break;
    }
    if((sum[a-]-T)*(sum[a]-T) < &&a!=){
        if(sum[a-] > T) dataout(T,sum[a],s[a],a-,T);
        else dataout(T,sum[a-],s[a-],a,T);
    }
    else{
        if(sum[a+] > T) dataout(T,sum[a],s[a],a+,T);
        else dataout(T,sum[a+],s[a+],a,T);
    }
}

void deal1(double T)
{
    sort(s, s + m * n);
    int i;
    bool flag = false;
    sum[] = ;
    for(i=;i<n*m;i++){
        for(int j=;j<i;j++){
            sum[i]+=(s[i] - s[j]);
            /*cout<<"i "<<i<<" j "<<j<<endl;
            cout<<sum[i]<<" = "<<s[i]<<" - "<<s[j]<<endl;*/
        }
        if(sum[i] >= T){
            flag = true;
            break;
        }
    }
    /*cout<<i<<endl;
    cout<<"sum[i-1] "<<sum[0]<<"  sum[i] "<<sum[1]<<endl;*/
    if(flag){
        double height ,per;
        if(T == ){
            per=;
            height = s[];
        }
        else {
            per =  * i / m / n;
            height = s[i-] + (sum[i] - sum[i-]) / i;
        }
        cout<<setiosflags(ios::fixed)<<setprecision();
        cout << "Water level is " <<height << " meters."<<endl;
        cout << per<<" percent of the region is under water."<<endl;
    }
    else{
        i--;
        /*cout<<"false"<<endl;
        cout<<"i "<<i<<endl;
        cout<<" s[i-1] "<<s[i-1]<<" sum[i] "<<sum[i]<<endl;*/
        double height ,per;
        if(T == ){
            per = ;
            height = s[];
        }
        else {
            per = ;
            height = s[i] + (T - sum[i]) / (i + ) ;
        }
        cout<<setiosflags(ios::fixed)<<setprecision();
        cout << "Water level is " <<height << " meters."<<endl;
        cout << per<<" percent of the region is under water."<<endl;
    }
}

int main()
{
    int times=;
    double T;
    while( cin >> n >> m ){
        if(n == )break;
        for(int i = ;i < n*m;i++)cin >> s[i];
        cin >> T;
        cout<< "Region " << ++times << endl;
        memset(sum,,sizeof(sum));
        if(m == &&n == ){
                double height = T/ + s[],per;
                if(T == ) per = ;
                else per = ;
                cout << setiosflags(ios::fixed) << setprecision();
                cout << "Water level is " <<height << " meters."<<endl;
                cout << per << " percent of the region is under water."<<endl;
        }
        else if(m*n<)deal1(T/);
        else deal2(T / );
    }
    return ;
}

实在不知错在哪，就不用二分法试试。
一次AC
果然还是我想太多了，不用二分法也并没有超时。

#include <iostream>
#include <algorithm>
#include <iomanip>
using namespace std;
const int maxn = *+;
int s[maxn],m,n;    //s为各个网格的海拔高度，有正有负

void dataout(double a1,double a2,double ss,double numm,double T)
{
    if(numm == )numm=;
    double height = ss + (a1 - a2) / numm , per =  * numm / (m*n);
    cout<<setiosflags(ios::fixed)<<setprecision();
    cout << "Water level is " << height << " meters."<<endl;
    cout << per << " percent of the region is under water."<<endl;
}

void deal(double T)
{

    int a = ,t1 = ,t2 = ,i;
    for(i=;i<n*m;i++){
        t2 = t1;
        t1 = ;
        for(int j=;j<i;j++)t1+=(s[i] - s[j]);
        if(t2< T&& (t1 > T||t1 == T))break;

    }
    if(i == n*m) dataout(T,t1,s[i-],i,T);
    else dataout(T,t2,s[i-],i,T);
}

int main()
{
    int times=;
    double T;
    while( cin >> n >> m ){
        if(n == )break;
        for(int i = ;i < n*m;i++)cin >> s[i];
        cin >> T;
        cout<< "Region " << ++times << endl;
        sort(s, s + m * n);
        if(T == ){
                double height = s[], per = ;
                cout << setiosflags(ios::fixed) << setprecision();
                cout << "Water level is " <<height << " meters."<<endl;
                cout << per << " percent of the region is under water."<<endl;
        }
        else deal( T /  );
    }
    return ;
}