[贪心] BZOJ 2811 [Apio2012]Guard

好恶心的贪心

大概就是先把不可能有的去掉

然后把区间合并 有包含的大区间去掉

如果只剩K个空 输出 结束

剩下区间扫一遍 如果一个区间内只剩一个空 删掉这个空 加入答案

然后对于每个位置 如果让它不放 覆盖剩余区间所需的个数如果大于K 这个空必选

*对于覆盖前t个区间至少要几个可以用DP

#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;

inline char nc()
{
	static char buf[100000],*p1=buf,*p2=buf;
	if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; }
	return *p1++;
}

inline void read(int &x)
{
	char c=nc(),b=1;
	for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
	for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}

const int N=100005;

int n,m,K;
int del[N],sum[N],last[N];

struct Seg{
	int l,r;
	Seg(int _=0,int __=0):l(_),r(__) { }
	bool operator < (const Seg &B) const{
		return l==B.l?r<B.r:l<B.l;
	}
}a[N],b[N],tmp[N];
int acnt,bcnt;

int tot,ans[N];

int f[N];

int main()
{
	int x,y,z,t,icnt;
	freopen("t.in","r",stdin);
	freopen("t.out","w",stdout);
	read(n); read(K); read(m);
	for (int i=1;i<=m;i++)
		read(x),read(y),read(z),z?a[++acnt]=Seg(x,y):b[++bcnt]=Seg(x,y);
	sort(b+1,b+bcnt+1);
	x=1; y=0;
	for (int i=1;i<=bcnt;i++)
		if (b[i].l<=y)
			y=max(y,b[i].r);
		else
		{
			for (int j=x;j<=y;j++) del[j]=1; 
			x=b[i].l; y=b[i].r;
		}
	for (int j=x;j<=y;j++) del[j]=1;
	x=0;
	for (int i=1;i<=n;i++)
		if (!del[i])
			sum[i]=sum[i-1]+1,last[i]=x=i;
		else
			sum[i]=sum[i-1],last[i]=x;
	if (sum[n]==K){
		for (int i=1;i<=n;i++) if (!del[i]) printf("%d\n",i); return 0;
	}
	sort(a+1,a+acnt+1);
	x=n+1; m=0;
	for (int i=acnt;i;i--)
		if (a[i].r<x)
           tmp[++m]=a[i],x=a[i].r;
	for (int i=1;i<=m;i++) a[i]=tmp[i];
	for (int i=1;i<=m;i++)
		if (sum[a[i].r]-sum[a[i].l-1]==1)
		{
			x=last[a[i].r];
			if (!del[x]) ans[++tot]=x,del[x]=1;
		}
	x=0;
	for (int i=1;i<=n;i++)
		if (!del[i])
			sum[i]=sum[i-1]+1,last[i]=x=i;
		else
			sum[i]=sum[i-1],last[i]=x;
	sort(ans+1,ans+tot+1);
	acnt=0;
	for (int i=m;i;i--)
	{
		x=lower_bound(ans+1,ans+tot+1,a[i].l)-ans;
		y=lower_bound(ans+1,ans+tot+1,a[i].r+1)-ans;
		if (x==y) tmp[++acnt]=a[i];
	}
	m=acnt;
	for (int i=1;i<=acnt;i++) a[i]=tmp[i];
	for (int i=m;i;i--)
	{
		x=last[a[i].r];
       	y=lower_bound(a+1,a+m+1,Seg(x+1,0))-a;
       	f[i]=f[y]+1;
	}
	K-=tot;
   	z=0; icnt=0;
	for (int i=1;i<=m;i++)
	    if (a[i].l>z)
		{
	        x=last[a[i].r]; y=last[x-1];
	        t=lower_bound(a+1,a+m+1,Seg(y+1,0))-a;
	        if (icnt+f[t]>=K) ans[++tot]=x;
	        z=x; icnt++;
	    }
	sort(ans+1,ans+tot+1);
	if (!tot) 
		printf("-1\n");
   	else
       	for (int i=1;i<=tot;i++) printf("%d\n",ans[i]);
	return 0;
}