CF1523G - Try Booking——暴力、结论、树套树

最新推荐文章于 2022-06-11 12:03:17 发布

偶耶XJX

最新推荐文章于 2022-06-11 12:03:17 发布

阅读量169

点赞数

分类专栏：信息竞赛解题文章标签： Codeforces 算法数据结构

本文链接：https://blog.csdn.net/weixin_43960287/article/details/117484221

版权

信息竞赛解题专栏收录该内容

195 篇文章 5 订阅

订阅专栏

G - Try Booking

题目描述

给定 $m$ 个区间，按如下方式操作：

只操作区间大小不小于 $x$ 的区间，按标号从小到大依次操作；
若该区间内没有任何地方被涂色，则将该区间内每处都涂色。

对于 $[1, n]$ 内的每个 $x$ ，求操作完成后数轴上被涂色的总长。

数据范围与提示

$1\le n\le 5\cdot 10^4,1\le m\le 10^5,1\le l_i\le r_i\le n$ 。

前言

$n+\frac{n}{2}+\frac{n}{3}+...+\frac{n}{n-1}+1 \approx n*\ln n$ ，请熟记这个结论。

思路

有一个暴力的思路，假设我们处理了那些区间该加入，那么我们可以按标号从小到大覆盖上去。假设我们找到了一个标号最小的区间 $[l, r]$ ，由于后面加入的区间不能和 $[l, r]$ 有交，所以我们可以把问题化成 $[1, l - 1]$ 、 $[r + 1, n]$ 内的两个子问题求解，最后覆盖上的区间数量不超过 $\frac{n}{x}$ 个。

这样一来，我们把每个 $x$ 对应的答案加起来发现它不超过 $n+\frac{n}{2}+\frac{n}{3}+...+\frac{n}{n-1}+1 \approx n*\ln n$ ，说明只要查找标号最小区间的速度够快，此暴力是可行的！

于是考虑经典做法，用树套树维护区间内的区间信息。我们通过把 $x$ 倒着求，可以把问题简化为仅加入区间的问题，所以用树状数组套线段树即可。时间复杂度 $O(n\log^3n+m\log^2n)$ 。

代码

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<ctime>
#include<map>
#define ll long long
#define MAXN 100005
#define uns unsigned
#define INF 0x3f3f3f3f
#define MOD 1000000007ll
#define lowbit(x) ((x)&(-(x)))
using namespace std;
inline ll read(){
	ll x=0;bool f=1;char s=getchar();
	while((s<'0'||s>'9')&&s>0){if(s=='-')f^=1;s=getchar();}
	while(s>='0'&&s<='9')x=(x<<1)+(x<<3)+s-'0',s=getchar();
	return f?x:-x;
}
int n,m;
int L[MAXN],R[MAXN],ans[MAXN];
vector<int>ad[MAXN];//↓线段树
struct itn{
	int ls,rs,a;itn(){ls=rs=0,a=INF;}
}t[MAXN<<6];
int IN;
inline void add(int x,int l,int r,int z,int d){
	t[0]=itn();
	if(l==r){t[x].a=min(t[x].a,d);return;}
	int mid=(l+r)>>1;
	if(z<=mid){
		if(!t[x].ls)t[x].ls=++IN,t[IN]=itn();
		add(t[x].ls,l,mid,z,d);
	}
	else{
		if(!t[x].rs)t[x].rs=++IN,t[IN]=itn();
		add(t[x].rs,mid+1,r,z,d);
	}
	t[x].a=min(t[t[x].ls].a,t[t[x].rs].a);
}
inline int query(int x,int l,int r,int a,int b){
	if(!x||a>b)return INF;
	if(l==a&&r==b)return t[x].a;
	int mid=(l+r)>>1,res=INF;
	if(a<=mid)res=min(res,query(t[x].ls,l,mid,a,min(mid,b)));
	if(b>mid)res=min(res,query(t[x].rs,mid+1,r,max(a,mid+1),b));
	return res;
}//↓树状数组
int rt[MAXN];
inline void add_(int x,int l,int r){
	for(;r<=n;r+=lowbit(r)){
		if(!rt[r])rt[r]=++IN,t[IN]=itn();
		add(rt[r],1,n,l,x);
	}
}
inline int query_(int l,int r){
	int res=INF;
	for(;r>0;r-=lowbit(r))res=min(res,query(rt[r],1,n,l,n));
	return res;
}
inline int solve(int l,int r){
	if(l>r)return 0;
	int id=query_(l,r);
    if(id>m)return 0;
	return solve(l,L[id]-1)+solve(R[id]+1,r)+R[id]-L[id]+1;
}
signed main()
{
	n=read(),m=read();
	for(int i=1;i<=m;i++){
		L[i]=read(),R[i]=read();
		ad[R[i]-L[i]+1].push_back(i);
	}
	for(int X=n;X>0;X--){
		for(uns i=0;i<ad[X].size();i++)add_(ad[X][i],L[ad[X][i]],R[ad[X][i]]);
		ans[X]=solve(1,n);
	}
	for(int i=1;i<=n;i++)printf("%d\n",ans[i]);
	return 0;
}