2018北大暑校acm算法训练课程 矩形分割 二分

总时间限制: 1000ms 内存限制: 65536kB
描述
平面上有一个大矩形，其左下角坐标（0，0），右上角坐标（R,R)。大矩形内部包含一些小矩形，小矩形都平行于坐标轴且互不重叠。所有矩形的顶点都是整点。要求画一根平行于y轴的直线x=k（k是整数) ，使得这些小矩形落在直线左边的面积必须大于等于落在右边的面积，且两边面积之差最小。并且，要使得大矩形在直线左边的的面积尽可能大。注意：若直线穿过一个小矩形，将会把它切成两个部分，分属左右两侧。

输入
第一行是整数R，表示大矩形的右上角坐标是(R,R) (1 <= R <= 1,000,000)。
接下来的一行是整数N,表示一共有N个小矩形(0 < N <= 10000)。
再接下来有N 行。每行有4个整数，L,T, W 和 H, 表示有一个小矩形的左上角坐标是(L,T),宽度是W，高度是H (0<=L,T <= R, 0 < W,H <= R). 小矩形不会有位于大矩形之外的部分。
输出
输出整数n，表示答案应该是直线 x=n。 如果必要的话，x=R也可以是答案。
样例输入
1000
2
1 1 2 1
5 1 2 1
样例输出
5

分析:二分横线位置,判断答案(不多说了,这个题和前面的太类似了)

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<stack>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<sstream>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int MAX_N=100000;
const int INF = 0x3f3f3f3f;
#define eps 1e-7
//#define Kongxiangzhouye
#define _for(i,a,b) for( int i=(a); i<(b); ++i)
#define _rep(i,a,b) for( int i=(a); i<=(b); ++i)
ll x1[10005];
ll x2[10005];
ll h[10005];
ll s[10005];   
ll R,N,ans;
ll C(ll x){
    ll sum=0;
    for(ll i=1;i<=N;i++){  
        if(x2[i]<=x)
            sum+=s[i];
        else if(x1[i]>=x)
            sum-=s[i];
        else
            sum=sum+(x-x1[i])*h[i]-(x2[i]-x)*h[i];
    }
    return sum;
}
int main(){
    ios::sync_with_stdio(false);
    #ifdef Kongxiangzhouye
        freopen("in.txt","r",stdin);
        freopen("out.txt","w",stdout);
    #endif
    cin>>R>>N;
    int te;
    for(ll i=1;i<=N;i++){
        cin>>x1[i]>>te>>w[i]>>h[i];
        x2[i]=x1[i]+w[i];
        s[i]=w[i]*h[i];
    }
    ll l=0;
    ll r=R;
    while(l<r){
        ll mid=(l+r)/2;
        if(C(mid)>=0)
            r=mid;  
        else
            l=mid+1;
    }
    if(l==r)
        ans=l;
    while(f(ans)==f(ans+1)&&ans<R)//就是因为这儿,唉 
        ans++;
    cout<<ans<<endl;
    return 0;
}