习题5-5：矩形分割(POJ 4136)（***）

问题链接：http://bailian.openjudge.cn/practice/4136/

Note:

1.只判断左右矩阵面积是否相等是不够的，需要返回差值

2.求左右矩阵面积差值的函数返回值要用long long 型，因为（1<=R<=1000000）

描述

平面上有一个大矩形，其左下角坐标（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<cstdio>
#include<iostream>
using namespace std;
int R, N;
struct node{
	int L, T;//左上角坐标 
	int W;//宽度
	int H;//高度 
};
node A[10001];

long long area(int mid){//求左右矩阵面积的差值 
	long long sumleft = 0, sumright = 0;
	for(int i = 0; i < N; i++){
		int a = A[i].L, b = A[i].L + A[i].W - 1;//小矩阵的左右边 
	    if(b < mid)
	        sumleft += A[i].W * A[i].H;
	    else if(a >= mid)
	        sumright += A[i].W * A[i].H;
	    else{
	    	sumleft += (mid - A[i].L) * A[i].H;
			sumright += (A[i].L + A[i].W - mid) * A[i].H;
		}
	}	
	return sumleft - sumright;
}

int compare(int left, int right){
	int mid;
	while(left <= right)
	{
		mid = (left + right)/2;
		if(area(mid) < 0)//如果左边的面积小于右边的面积
		    left = mid + 1;
		else if(area(mid) > 0)//如果左边的面积大于右边的面积 
		    right = mid - 1;
		else
		    return mid;     
	}
	return left;
}

int main(){
	cin >> R >> N;
	
	for(int i = 0; i < N; i++)
	{
		cin >> A[i].L >> A[i].T >> A[i].W>> A[i].H;
	}
	
	if(N == 1 && A[N-1].W == 1){
		cout << R << endl;
		return 0;
	}
	     
	
	int x = compare(0, R);
	
	while(area(x) == area(x + 1) && x < R)
	    x++;
	cout << x;

	return 0;
}