OpenJudge：百练 4136:矩形分割

OpenJudge：百练 4136:矩形分割

描述

平面上有一个大矩形，其左下角坐标（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 5
2
1 1 2 1
5 1 2 1
题中给的整数值较大，程序中使用long long，否则可能 Wrong answer

#include<iostream>
using namespace std;
const int MAX=1e5;
long long R,N,L[MAX],sz=0,mid,sy=0,T[MAX],W[MAX],H[MAX];

long long S(long long k){//计算左边减右边矩形面积
	sz=0;sy=0;
	for(long long i=1;i<=N;i++){//全在左边
		if(L[i]+W[i]<=k)
			sz+=W[i]*H[i];
		else if(L[i]>=k)//全在右边
			sy+=W[i]*H[i];
		else {//左右都有
			sz+=(k-L[i])*H[i];
			sy+=(L[i]+W[i]-k)*H[i];
		}
	}
	return sz-sy;//返回左边减右边的值
}

int main(){
	cin>>R;
	cin>>N;
	for(long long  i=1;i<=N;i++)
		cin>>L[i]>>T[i]>>W[i]>>H[i];
	long long j=0,m=R;

	while(j<m){//二分法求左边减右边面积面积之差最小
		mid=(j+m)/2;
		if(S(mid)>=0)
			m=mid;
		else
			j=mid+1;
	}
	
	
	while(S(j)==S(j+1)&&j+1<=R){//在左边减右边面积差不变的情况下，让题中k，也就是程序中的j尽可能地大
		j++;
	}
	cout<<j;
	return 0;
	
	
}