zoj 3540 Adding New Machine(map+离散化)

地址：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;
}