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;
}