BZOJ.3693.圆桌会议(Hall定理 线段树)

题目链接

先考虑链。题目相当于求是否存在完备匹配。那么由Hall定理,对于任意一个区间[L,R],都要满足[li,ri]完全在[L,R]中的ai之和sum小于等于总位置数,即R-L+1。(其实用不到Hall定理,显然)
为什么不是子集呢,因为区间并和子集等价,所有区间并都是要验证的。
而且可以发现,只有当R为某个r[i],L为某个l[j]时,[L,R]才有必要验证。
所以我们将区间按r[]排序,枚举每个r[i]作为R。限制条件为\(sum<=R-L+1\)\(sum+L-1<=R\),对于前面所有位置,即\(L∈[1,R]\),我们可以直接用线段树维护每个位置的sum+L-1,查询最大值就可以了。
每次枚举到i,会使\(L∈[1,l[i]]\)时+ai,这个区间加就行了。

至于环,拆环为链就可以了。
但是有个问题,比如m=5,([1,3],2),([3,2],3),会分成区间[1,3],[3,7],[6,8],[3,2]实际是包含[1,3]的,但是[1,3]被分成了[3,3]和[6,7],它的a就不会就加上了。
因为r=l-1时,长为m,会把另一个[l,r]与复制出来的[l+m,r+m]分割开,导致都不包含在[r,l-1]内。
因为这些区间长度都为m,所以只需要先判一下整个区间即可。

r[]不需要离散化,否则注意是与ref[r]比较。。

不难 又浪费很长时间。。

//10496kb   888ms(Rank1!)
#include <cstdio>
#include <cctype>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 300000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
const int N=2e5+5;

int n,m,ref[N];
char IN[MAXIN],*SS=IN,*TT=IN;
struct Opt
{
    int l,r,a;
    Opt() {}
    Opt(int l,int r,int a):l(l),r(r),a(a) {}
    bool operator <(const Opt &x)const{
        return r<x.r;
    }
}q[N];
struct Segment_Tree
{
    #define lson rt<<1
    #define rson rt<<1|1
    #define ToL l,m,rt<<1
    #define ToR m+1,r,rt<<1|1
    int mx[N<<2],tag[N<<2];

    #define Update(rt) mx[rt]=std::max(mx[lson],mx[rson])
//  #define Add(rt,v) mx[rt]+=v, tag[rt]+=v
//  #define PushDown(rt) Add(lson,tag[rt]), Add(rson,tag[rt]), tag[rt]=0
    inline void Add(int rt,int v){
        mx[rt]+=v, tag[rt]+=v;
    }
    void PushDown(int rt){
        Add(lson,tag[rt]), Add(rson,tag[rt]), tag[rt]=0;
    }
    void Build(int l,int r,int rt)
    {
        tag[rt]=0;
        if(l==r) mx[rt]=ref[l]-1;//ref!
        else Build(l,l+r>>1,lson), Build((l+r>>1)+1,r,rson), Update(rt);
    }
    void Modify(int l,int r,int rt,int L,int R,int v)
    {
        if(L<=l && r<=R) {Add(rt,v); return;}
        if(tag[rt]) PushDown(rt);
        int m=l+r>>1;
        if(L<=m) Modify(ToL,L,R,v);
        if(m<R) Modify(ToR,L,R,v);
        Update(rt);
    }
    int Query(int l,int r,int rt,int L,int R)
    {
        if(L<=l && r<=R) return mx[rt];
        if(tag[rt]) PushDown(rt);
        int m=l+r>>1;
        if(L<=m)
            if(m<R) return std::max(Query(ToL,L,R),Query(ToR,L,R));
            else return Query(ToL,L,R);
        return Query(ToR,L,R);
    }
}T;

inline int read()
{
    int now=0;register char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);now=now*10+c-'0',c=gc());
    return now;
}
inline int Find(int x,int r)
{
    int l=1,mid;
    while(l<r)
        if(ref[mid=l+r>>1]<x) l=mid+1;
        else r=mid;
    return l;
}

int main()
{
    int Case=read();
    while(Case--)
    {
        n=read(),m=read(); int Q=0,t=0; long long sum=0;
        for(int i=1,l,r,a; i<=n; ++i)
        {
            ref[++t]=l=read(),r=read(),sum+=(a=read());
            if(l<=r) q[++Q]=Opt(l,r,a), q[++Q]=Opt(l+m,r+m,a), ref[++t]=l+m/*!*/;
            else q[++Q]=Opt(l,r+m,a);
        }
        if(sum>m) {puts("No"); continue;}
        std::sort(q+1,q+1+Q);
        std::sort(ref+1,ref+1+t); int cnt=1;
        for(int i=2; i<=t; ++i) if(ref[i]!=ref[i-1]) ref[++cnt]=ref[i];
        for(int i=1; i<=Q; ++i) q[i].l=Find(q[i].l,cnt);//, q[i].r=Find(q[i].r,cnt);
        n=cnt, T.Build(1,n,1);
        bool f=1;
        for(int i=1; i<=Q; ++i)
        {
            T.Modify(1,n,1,1,q[i].l,q[i].a);
            if(T.Query(1,n,1,1,q[i].l/*l对?用r还得离散化...*/)>q[i].r) {f=0; break;}//ref[r]
        }
        puts(f?"Yes":"No");
    }
    return 0;
}

猜你喜欢

转载自www.cnblogs.com/SovietPower/p/9315295.html