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