被粉碎的线段树seg

被粉碎的线段树(seg)
[题目自行脑补]...


又是一种进阶版,由于题目图片无法复制,我在这里概括一下题意:
咳咳,这道T3部分分给的有点水啊,60分很好拿(其实能拿85分的只是本蒟蒻只拿了60分)
回到正题,题意是...
原本线段树区间[l,r]是划分为[l,(l+r)/2]和[(l+r)/2+1,r]的。
然而没这么简单,所以小S自己给这棵线段树每个区间定了个mid。
构成了一种全新的线段树,现在给定m个区间,求最少要多少个线段树上的区间能够将这个区间完全不重叠的覆盖。




此题60分的数据范围是n<=30000 m<=1000 显然可以用nm的暴力来做(严谨点是(2n-1)*m)构造出每一个区间一共(2n-1)个,按深度从左到右存储,然后用dfs从树根开始往下查找,如果当前区间与给定要覆盖的区间没有任何交集,则不用往下找,如果完全覆盖,也不必往下找,答案数加一,部分覆盖就继续往下找。其实每一层只会有两个点会继续往下找,所以其实这个算法的时间复杂度是O(m*树的深度)。然后这样其实就变成80分了啊!赞,再加上一个mid=r-1的部分分就有85分。




上述方法与AC此题毫无关系hhh
现在讲正解叭。。。结论:区间为[l,r],则答案其实就是(r-l+1)*2-完全包含在这个区间内的区间的个数。
为什么呢?..想想,这个区间肯定是分散的包含在几棵子树中的,那么这几棵子树的叶子结点的个数就是区间的长度。
叶子结点的个数*2-1就是这些树的节点总个数。辣么减掉被完全覆盖的点的个数不就是答案了嘛~


现在的问题就是:给定区间[l,r],求有多少个区间满足x>=l && y<=r 。
这其实是一个经典的树状数组题,按前一个排序,后端点加入树状数组即可。
我的做法略有些劣,我是将询问和原树的那些区间放在一起排序,遇到询问的就做sum操作
否则就修改。关键字1个与2个3个的差别就是20分40分和100分啊啊啊qwq
这真的是郁闷死,由于我这个做法,所以先按左端点排序,再按右端点排序,最后一个关键字是询问或者原有区间。

真心郁闷啊...



#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#define N 100000+1000
using namespace std;
int k,x,n,m,t[N],le[N+N+N],ri[N+N+N],son[N+N][2],res[N],q[N],q1[N],num,num1,l1[N],r1[N];
struct node{
	int l,r,id,belong;
}a[N+N+N+N];
bool cmp(node x,node y){if(x.l!=y.l) return x.l>y.l;else if(x.r!=y.r) return x.r<y.r;else return x.id<y.id;}
int read(){
	char ch='+',s=getchar();
	while(s<'0'||s>'9'){if (s=='-') ch='-';s=getchar();}
  	int ans=0;
  	while(s>='0'&&s<='9'){ans=ans*10+(s-'0');s=getchar();}
  	if (ch=='-')ans=-ans;
  	return ans;
}
void dfs(int t){
	if(le[t]==ri[t]) return;
	if(t>1){
		son[t][1]=++k;
		son[t][2]=++k;
	}
	x=read();
	le[son[t][1]]=le[t];ri[son[t][1]]=x;
	le[son[t][2]]=x+1;ri[son[t][2]]=ri[t];
	dfs(son[t][1]);dfs(son[t][2]);
}
void modify(int x){
	while(x<=n){
		t[x]++;
		x+=x&-x;
	}
}
int sum(int x){
	int ans=0;
	while(x>0){
		ans+=t[x];
		x-=x&-x;
	}
	return ans;
}
int main()
{
	freopen("seg.in","r",stdin);
	freopen("seg.out","w",stdout);
	n=read(),m=read();
	le[1]=1;ri[1]=n;
	son[1][1]=2;son[1][2]=3;k=3;
	dfs(1);
	for(int i=1;i<=k;i++) a[i].l=le[i],a[i].r=ri[i],a[i].id=1,a[i].belong=0;
	for(int i=1;i<=m;i++){
		k++;
		a[k].l=read(),a[k].r=read();
		a[k].id=2;a[k].belong=i;
		l1[i]=a[k].l;r1[i]=a[k].r;
	}
	sort(a+1,a+1+k,cmp);
	for(int i=1;i<=k;i++){
		if(a[i].id==1) modify(a[i].r);
		else res[a[i].belong]=sum(a[i].r);
	}
	for(int i=1;i<=m;i++) printf("%d\n",(r1[i]-l1[i]+1)*2-res[i]);
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值