[BZOJ4082][Wf2014]Surveillance[倍增]
题意
给你一个长度为 \(len\) 的环,以及 \(n\) 个区间,要你选择尽量少的区间,使得它们完全覆盖整个环。问最少要多少个区间。
\(len,n\leq 10^6\) .
分析
-
考虑普通的区间覆盖的贪心做法,这里只需要倍增一下。
-
如果区间 \(r< l\) 把 \(r\) 设置成 \(r+len\) 。
-
每个解都一定存在一个区间满足其标号最小,且可以通过后面的区间跳回自己的右边,只需要以这类区间开头即可。
-
总时间复杂度为 \(O(nlogn)\) .
处理区间覆盖的技巧:记录**后缀最小值 **\(mn\),如果\(mn\leq l+1\) 就让指针++.这样就可以找到最靠后的满足要求的位置。
倍增的时候注意:如果不存在这么长的转移是不能够转移的,否则可能答案(步数)会算大。
代码
#include<bits/stdc++.h>
using namespace std;
#define go(u) for(int i=head[u],v=e[i].to;i;i=e[i].last,v=e[i].to)
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define pb push_back
typedef long long LL;
inline int gi(){
int x=0,f=1;char ch=getchar();
while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-48;ch=getchar();}
return x*f;
}
template<typename T>inline bool Max(T &a,T b){return a<b?a=b,1:0;}
template<typename T>inline bool Min(T &a,T b){return b<a?a=b,1:0;}
const int N=1e6 + 7,inf=0x3f3f3f3f;
int m,len,ans=inf;
int fa[N][21],mn[N];
struct data{
int l,r;
bool operator <(const data &rhs)const{
if(r!=rhs.r) return r<rhs.r;
return l<rhs.l;
}
}t[N];
int main(){
len=gi(),m=gi();
rep(i,1,m) {
t[i].l=gi(),t[i].r=gi();
if(t[i].r<t[i].l) t[i].r+=len;
}
sort(t+1,t+1+m);
mn[m+1]=inf;
for(int i=m;i;i--) mn[i]=min(t[i].l,mn[i+1]);
for(int i=1,j=1;i<=m;++i){
while(j+1<=m&&mn[j+1]-1<=t[i].r)++j;
if(i!=j) fa[i][0]=j;
}
for(int i=m;i;--i)
for(int j=1;j<=20;++j) fa[i][j]=fa[fa[i][j-1]][j-1];
for(int i=1;i<=m;++i){
int s=1,x=i;
for(int j=20;~j;--j) if(fa[x][j]&&t[fa[x][j]].r-t[i].l+1<len) x=fa[x][j],s+=1<<j;
if(fa[x][0]&&t[x].r-t[i].l+1<len) ++s,x=fa[x][0];
if(t[x].r-t[i].l+1>=len) Min(ans,s);
}
if(ans==inf) puts("impossible");
else printf("%d\n",ans);
return 0;
}