uva11983扫描线k次覆盖

自己做的是从下往上扫描的，一直wa，不知道坑在哪里。。但是作为模板。我还是找了份不错的ac代码

/*
被覆盖不低于k次的点
每个点对应了一个单位面积，本题把点转面积即是被覆盖不低于k次的面积
可以当做求k次面积覆盖模板的题目！
*/
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<map>
using namespace std;
#define maxn 60005
#define lson l,m,rt<<1
#define rson m,r,rt<<1|1
#define ll long long
int n,k;
struct Seg{
    int l,r,h,c;
    Seg(){}
    Seg(int l,int r,int h,int c):l(l),r(r),h(h),c(c){}
    bool operator<(const Seg& a)const {
        return h<a.h;
    }
}segs[maxn];
int tot,totx,x[maxn];//线段树建立在x轴上
int len[maxn<<][],flag[maxn<<];
map<int,int>mp;

inline void pushup(int rt,int l,int r){
    if(flag[rt]>=k){
        for(int i=;i<=k;i++) len[rt][i]=;
        len[rt][k]=x[r]-x[l];

    }
    else if(flag[rt]>){
        int cur=flag[rt];//当前区间覆盖次数
        for(int i=;i<=k;i++) len[rt][i]=;//先把当前区间所有覆盖情况置零
        len[rt][cur]=x[r]-x[l];
        if(l+==r) return;

        for(int i=;i<=k;i++){
            if(i+cur>=k) len[rt][k]+=len[rt<<][i]+len[rt<<|][i];
            else len[rt][i+cur]+=len[rt<<][i]+len[rt<<|][i];
        }
        for(int i=cur+;i<=k;i++)
            len[rt][cur]-=len[rt][i];
    }
    else {
        for(int i=;i<=k;i++) len[rt][i]=;
        if(l+==r) return;
        for(int i=;i<=k;i++)
            len[rt][i]=len[rt<<][i]+len[rt<<|][i];
    }
}
void update(int L,int R,int c,int l,int r,int rt){
    if(L<=l && R>=r){
        flag[rt]+=c;;
        pushup(rt,l,r);
        return;
    }
    int m=l+r>>;
    if(L<m) update(L,R,c,lson);
    if(R>m) update(L,R,c,rson);
    pushup(rt,l,r);
}

void init(){
    tot=totx=;
    mp.clear();
    memset(len,,sizeof len);
    memset(flag,,sizeof flag);
}
int main(){
    int T,a,b,c,d;
    cin >> T;
    for(int tt=;tt<=T;tt++){
        init();
        scanf("%d%d",&n,&k);
        for(int i=;i<n;i++){
            scanf("%d%d%d%d",&a,&b,&c,&d);
            c++;d++;
            segs[tot++]=Seg(a,c,b,);
            segs[tot++]=Seg(a,c,d,-);
            x[totx++]=a;x[totx++]=c;
        }
        sort(x,x+totx);
        sort(segs,segs+tot);
        totx=unique(x,x+totx)-x;
        for(int i=;i<totx;i++) mp[x[i]]=i;

        ll res=;
        for(int i=;i<tot;i++){
            if(i!=)
                res+=(segs[i].h-segs[i-].h)*len[][k];
            //  cout<<segs[i].h-segs[i-1].h<< " "<<len[1][k]<<'\n';
            update(mp[segs[i].l],mp[segs[i].r],segs[i].c,,totx-,);

        }
        printf("Case %d: %lld\n",tt,res);
    }
    return ;
}