地址:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=4495
题意:给你w*h个格子,并且一些矩形块占用着一些格子,问剩下的空间,能装下1*m的矩形的方案数
分析:我们可以枚举一个下标,然后维护另一个下标构成的线段,由于题目的特殊性,线段不会重叠,所以变得简单了许多,对于每次添加一条线段l~r,它的左边第一条线段的右边界是l',右边第一条线段的左边界是r',当前的答案是now,那么新的答案是now-(r'-l'-1 -m +1) + (l -l'-1 -m +1) +(r'-r-1-m+1),当然得注意下于0的情况,也就是每次只需要快速得到l'和r',所以用map就能够维护这样一个东西,具体看代码吧
PS:这题wa了一天多才发现错误,但是居然是个诡异的错误,当然大家也懒得知道这样的错误,为了测处这个错我交了近50次T_T
还有,以后尽量都写构造函数 ,尽量用stl来实现吧,至少对于比赛是这样的.......
现在已经搞得没体力值了= =
代码:
#include <cstdio>
#include <algorithm>
#include <iostream>
#include <map>
#include <vector>
using namespace std;
#define REP(i,n) for(int i=0;i<(n);++i)
#define N 1000010
typedef long long LL;
int g[N][4];
struct que
{
int l,r,a,t;
que(int _a, int _l, int _r, int _t):l(_l),r(_r),a(_a),t(_t){}
};
bool cmp(const que &p, const que &q)
{
return p.a<q.a||(p.a==q.a&&p.t>q.t);
}
vector<que> q;
map<int, int> sgt;
map<int, int>::iterator it;
LL solve(int w, int h, int n, int len)
{
q.clear();
REP(i,n)
{
swap(g[i][0],g[i][1]),swap(g[i][2],g[i][3]);
q.push_back(que(g[i][0],g[i][1],g[i][3],1));
q.push_back(que(g[i][2],g[i][1],g[i][3],-1));
}
sort(q.begin(),q.end(),cmp);
sgt.clear();
sgt[0]=0,sgt[h+1]=h+1;
LL ret=0,now=max(h-len+1,0);
int i,l,r,low=0;
for(i=0;i<q.size();)
{
ret+=(q[i].a-1-low)*now;
low=q[i].a;
for(;i<q.size();++i)
if(q[i].a==low&&q[i].t>0)
{
sgt[q[i].l]=q[i].r;
it=sgt.find(q[i].l);
--it;
l=it->second;
++it,++it;
r=it->first;
if(r-l>len)
now+=max(q[i].l-l-len,0)+max(r-q[i].r-len,0)-(r-l-len);
}else break;
ret+=now;
for(;i<q.size();++i)
if(q[i].a==low)
{
it=sgt.find(q[i].l);
--it;
l=it->second;
++it,++it;
r=it->first;
--it;
sgt.erase(it);
if(r-l>len)
now-=max(q[i].l-l-len,0)+max(r-q[i].r-len,0)-(r-l-len);
}else break;
}
ret+=now*(w-low);
return ret;
}
int main()
{
int w,h,n,m;
LL ans;
while(~scanf("%d%d%d%d",&w,&h,&n,&m))
{
REP(i,n)REP(j,4)scanf("%d",&g[i][j]);
ans=solve(h,w,n,m);
if(m!=1)ans+=solve(w,h,n,m);
printf("%lld\n",ans);
}
return 0;
}