【bzoj2811】APIO2012 守卫 guard

       题目，点击就送

           这道题竟然不是差分？!（虽然看着很像，但是注意他有一个K的限制，而且差分求出的只是一个合理的方案，也不满足这里的必须。。）       在网上搜到了一个说用差分A了的人，然而他说为了不带坏小朋友，就不说这种骗分的做法了（%%%）。

           标准解法是贪心，首先删除那些一定没有的区间，zro hzwer orz 是把区间重新编号，我懒得编，就每次去二分位置，然后去掉那些区间长度为1的区间，这是必定有的。对于剩下的，这里要判一下是否剩下的和可能点个数相等，若相等，则全都是啦。首先在区间右边端点贪心可以得到解，然后我们考虑由右端点向左走，如果这时算出来的答案已经大于了我点的限制，那么前面那个一定是一个解了，这里可以用动归处理，细节处理解释在代码里面有。。

           

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#include<queue>
#define N 200005
using namespace std;
int n,m,u,tot=0;
struct eee{
   int a,b,c;
   eee(int _a=0,int _b=0,int _c=0):a(_a),b(_b),c(_c){}
   bool operator<(const eee &o)const{
      return a<o.a||(a==o.a&&b>o.b);
   }
}a[N],b[N];
int c[N],res[N],sum[N],sum2[N],f[N],rp[N],la[N];
void add(int x,int v){
	while(x<=n){
		c[x]+=v;
		x+=(x&(-x));
	}
}
int ask(int x){
	int ret=0;
	while(x){
		ret+=c[x];
		x-=(x&(-x));
	}
	return ret;
}
int ans[N];
int main(){
	//freopen("","r",stdin);
	//freopen("","w",stdout);
	int i,j;
	scanf("%d%d%d",&n,&u,&m);
	for(i=1;i<=m;i++){
		scanf("%d%d%d",&a[i].a,&a[i].b,&a[i].c);
		if(!a[i].c){add(a[i].a,1);add(a[i].b+1,-1);}
	}
	for(i=1;i<=n;i++){
		if(ask(i))res[i]=0;//一定没有标记为0 
		else res[i]=1;//反之 
		sum[i]=sum[i-1]+res[i];
	}
	for(i=1;i<=m;i++)
	    if(a[i].c&&sum[a[i].b]==sum[a[i].a-1]+1){
	    	int wei=upper_bound(sum+a[i].a,sum+a[i].b+1,sum[a[i].a-1])-sum;
	    	if(res[wei]!=2){//不确定是否一定有人 
	    		ans[++ans[0]]=wei;
	    		res[wei]=2;
	    		u--;
			}
		}
		else if(a[i].c&&!(sum[a[i].b]-sum[a[i].a-1])){puts("-1");return 0;}
	for(i=1;i<=n;i++){
		if(res[i]==2)sum2[i]=sum2[i-1]+1;
		else sum2[i]=sum2[i-1];
	}
	for(i=n;i>=1;i--){
		if(res[i]==1)sum[i]=sum[i+1]+1;
		else sum[i]=sum[i+1];
	}
	if(u==sum[1]&&u){
		for(i=1;i<=n;i++)
		  if(res[i])printf("%d\n",i);
		return 0;
	}
	else if(sum[1]<u){
		puts("-1");
		return 0;
	}
//	for(i=1;i<=ans[0];i++)cout<<"()()()()"<<ans[i]<<endl;
	sort(a+1,a+m+1);
//	for(i=1;i<=m;i++)cout<<a[i].a<<" "<< a[i].b<<" "<<a[i].c<<endl;
	for(i=1;i<=m;i++)
	   if(a[i].c){
		 if(sum2[a[i].b]>sum2[a[i].a-1])continue;//此区间已满足， 不要了。 
		 while(tot&&a[i].b<=b[tot].b)tot--;//若b把a完全包含，则不要b，因为若a满足了b一定满足 
		 b[++tot]=a[i];
	}
//	cout<<"***"<<tot<<endl;
	f[tot+1]=0;
	for(i=tot;i>=1;i--){
		int ll=b[i].a,rr=b[i].b;
		int L=ll,R=rr;
		while(L<R){
			int mid=(L+R+1)/2;
			if(sum[mid]-sum[rr+1]>0)L=mid;
			else R=mid-1;
		}
		rp[i]=L;//找到此区间真正的右端点 
		la[i]=upper_bound(b+i+1,b+tot+1,eee(rp[i],-1))-b;//找到左端点在我这里的那个 
		f[i]=f[la[i]]+1;//DP 
	}
	int xz=1;
	while(xz<=tot){
		int ll=b[xz].a,rr=b[xz].b;
		int L=ll,R=rr;
		while(L<R){
			int mid=(L+R+1)/2;
			if(sum[mid]-sum[rr+1]>0)L=mid;
			else R=mid-1;
		}
		int wei=L;//找到真正的右端点 
		L=ll;R=rr;
		while(L<R){
			int mid=(L+R+1)/2;
			if(sum[mid]-sum[rr+1]>1)L=mid;
			else R=mid-1;
		}
		int wei2=L;//找到真正的左端点，即wei~wei2中间必定有一个  
		j=upper_bound(b+xz+1,b+tot+1,eee(wei2,-1))-b;//找到左端点在wei2的区间
		if(f[j]+1>u)ans[++ans[0]]=wei;//一定要放一个更优 
		u--;
		xz=upper_bound(b+xz+1,b+tot+1,eee(wei,-1))-b;
	}
	if(!ans[0])puts("-1");
	else{
		sort(ans+1,ans+ans[0]+1);
		for(i=1;i<=ans[0];i++)printf("%d\n",ans[i]);
	}
	//fclose(stdin);
	//fclose(stdout);
    return 0;
}

这题怎么会写这么长我去QAQ