二分小练

03:矩形分割

总时间限制: 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

//Geeksun 2018.1.3

#include <stdio.h>
#include <stdlib.h>
long long l,r,R,mid;
long long zuo,you,ans,N;
struct Rectangle
{
    long long x,y,w,h,xx,s;
} rec[10010];
long long sum(long mid)
{
    int i;
    long long s1 = 0,s2 = 0;
    for(i = 0; i < N; i++)
    {
        if(rec[i].xx <= mid)
        {
            s1 += rec[i].s;
        }
        else if(rec[i].x >= mid)
        {
            s2 += rec[i].s;
        }
        else
        {
            s1 += (mid - rec[i].x) * rec[i].h;
            s2 += (rec[i].xx - mid) * rec[i].h;
        }
    }
    return s1 - s2;
}

int main()
{
    int i;
    long long x,y,w,h,maxx = -1;
    scanf("%lld",&R);
    scanf("%lld",&N);
    for(i = 0; i < N; i++)
    {
        scanf("%lld %lld %lld %lld",&x,&y,&w,&h);
        rec[i].x = x;
        rec[i].y = y;
        rec[i].xx = x + w;
        rec[i].h = h;
        rec[i].s = w * h;
        if(x + w > maxx)
        {
            maxx = x + w;
        }
    }
    l = 0,r = R;
    while(l + 1 < r)
    {
        mid = (l + r) / 2;
        if(sum(mid) > 0)
        {
            r = mid;
        }
        else if(sum(mid < 0))
        {
            l = mid;
        }
        else
        {
            printf("%lld",mid);
            return 0;
        }
    }
    zuo = sum(l);
    you = sum(r);
    if(zuo > you||zuo < 0)
    {
        ans = r;
    }
    else
    {
        ans = l;
    }
    if(ans == maxx)
    {
        ans = R;
    }
    printf("%lld",ans);
    return 0;
}