BZOJ2811: [Apio2012]Guard

一开始以为是差分约束yy了个结论WA了qwq，虽然至今还不知道为什么那个结论不对

先将序列中不可能有忍者的位置标出来
如果剩下的位置恰好有K个，要特判直接全部输出

然后将每个区间缩掉两端的0
如果存在某两个区间A,B，A包含B，那么A是不必要的
处理完后剩下的区间两两互不包含，按右端点排序后，从左到右左右端点递增
易知只有每个区间的右端点有可能一定有忍者，因为贪心的放肯定是放右端点更优

用f[i]表示保证了1~i区间的限制最少需要几个忍者，g[i]表示保证了i~n的区间限制至少需要几个忍者，F[i]表示i位置一定有忍者，保证所有限制，至少需要几个忍者
f[i]和g[i]可以简单的做两个贪心得到，F[i]用两个指针找i向左第一个不能覆盖的区间，向右第一个不能覆盖的区间，O(n)扫过去得到

然后从左到右扫区间，维护一个单调队列，对于每个区间$[l,r]$，维护$[l,r)$（左闭右开）内F[i]的最小值，如果这个最小值<=K，就代表右端点不是必须取的

标记，缩区间可以用并查集做到O(n)
dp，判断都可以做到O(n)
所以除去排序的复杂度是O(n)的，排序可以用桶排做到总复杂度O(n)
（然而跑的和nlogn的线段树一样快并没有什么用qwq

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

const int maxn = 210000;
const int maxm = 210000;

int n,m,K;
struct node { int l,r; }a[maxm]; int cnt;
inline bool cmp(const node x,const node y){return x.r==y.r?x.l>y.l:x.r<y.r;}
int fal[maxn],far[maxn],s[maxn];
int findfa(int fa[],const int x){return fa[x]==x?x:fa[x]=findfa(fa,fa[x]);}

int f[maxn],g[maxn],F[maxn];
int q[maxn],head,tail;
int ans[maxn],tp;

int main()
{
    scanf("%d%d%d",&n,&K,&m);
    for(int i=1;i<=n;i++) fal[i]=far[i]=i,s[i]=1; far[n+1]=n+1;
    for(int i=1;i<=m;i++)
    {
        int l,r,c; scanf("%d%d%d",&l,&r,&c);
        if(!c)
        {
            int j=l;
            while(j<=r)
            {
                j=findfa(far,j);
                if(j<=r) s[j]=0,far[j]=j+1;
            }
            j=r;
            while(j>=l)
            {
                j=findfa(fal,j);
                if(j>=l) s[j]=0,fal[j]=j-1;
            }
        }
        else a[++cnt]=(node){l,r};
    }
    int oth=0; for(int i=1;i<=n;i++) oth+=s[i];
    if(oth==K)
    {
        for(int i=1;i<=n;i++) if(s[i]) printf("%d\n",i);
        return 0;
    }

    for(int i=1;i<=cnt;i++)
        a[i].l=findfa(far,a[i].l),a[i].r=findfa(fal,a[i].r);
    sort(a+1,a+cnt+1,cmp);
    int temp=0;
    for(int i=1;i<=cnt;i++)
        if(a[i].l>a[temp].l) a[++temp]=a[i];
    cnt=temp;

    int r=a[1].r; f[1]=1;
    for(int i=2;i<=cnt;i++)
    {
        if(a[i].l<=r) f[i]=f[i-1];
        else f[i]=f[i-1]+1,r=a[i].r;
    }
    int l=a[cnt].l; g[cnt]=1;
    for(int i=cnt-1;i>=1;i--)
    {
        if(a[i].r>=l) g[i]=g[i+1];
        else g[i]=g[i+1]+1,l=a[i].l;
    }
    l=0,r=2; oth=1;
    for(;oth<=n;oth++)
    {
        while(a[l+1].r<oth&&l<cnt) l++;
        while(r<=cnt&&oth>=a[r].l) r++;
        F[oth]=f[l]+g[r]+1;
    }
    head=1,tail=0; r=0;
    for(int i=1;i<=cnt;i++)
    {
        while(r+1<a[i].r)
        {
            ++r;
            if(s[r])
            {
                while(head<=tail&&F[q[tail]]>=F[r]) tail--;
                q[++tail]=r;
            }
        }
        while(head<=tail&&q[head]<a[i].l) head++;
        if(head>tail||F[q[head]]>K) ans[++tp]=a[i].r;
    }
    if(!tp) puts("-1");
    //for(int i=1;i<=cnt;i++) printf("%d %d\n",a[i].l,a[i].r); puts("");
    for(int i=1;i<=tp;i++) printf("%d\n",ans[i]);

    return 0;
}