1. 程式人生 > >【bzoj4082】[Wf2014]Surveillance 倍增

【bzoj4082】[Wf2014]Surveillance 倍增

Description

給你一個長度為len的環,以及n個區間,要你選擇儘量少的區間,使得它們完全覆蓋整個環。問最少要多少個區間。

Input

輸入資料的第一行是兩個整數len和n,代表環的長度以及區間個數。之後n行描述的是n個區間,每個區間分別用一對數字(a,b)表示,若a≤b則表示這個區間覆蓋的是[a,b]部分,否則表示這個區間覆蓋的是除掉[a+1,b-1]以外的其他部分。

Output

輸出只有一行,一個整數,代表覆蓋整個環所需要的最少區間個數。

Sample Input

100 7

1 50

50 70

70 90

90 40

20 60

60 80

80 20

Sample Output

3

程式碼

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#define ll long long
#define mod 10000007
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if
(ch=='-') f=-1;ch=getchar();} while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} return x*f; } struct node{int l,r;}a[1000005]; int ans=1000000000,len,n; int mn[1000005],f[1000005][21]; bool cmp(node a,node b){return a.r<b.r;} int main() { len=read();n=read(); for (int i=1;i<=n;i++) { a[i].l=read();a[i].r=read(); if
(a[i].r<a[i].l) a[i].r+=len; } sort(a+1,a+n+1,cmp); mn[n]=a[n].l;for (int i=n-1;i;i--) mn[i]=min(a[i].l,mn[i+1]); for (int i=1,j=1;i<=n;i++) { while (j<n&&mn[j+1]-1<=a[i].r) j++; if (i!=j)f[i][0]=j; } for (int k=1;k<=20;k++) for (int i=1;i<=n;i++) f[i][k]=f[f[i][k-1]][k-1]; for (int i=1;i<=n;i++) { int sum=1,now=i; for (int j=20;j>=0;j--) if (f[now][j]&&a[f[now][j]].r<a[i].l+len) now=f[now][j],sum+=1<<j; if (a[now].r<a[i].l+len-1&&f[now][0]) now=f[now][0],sum++; if (a[now].r>=a[i].l+len-1) ans=min(ans,sum); } if (ans==1000000000) printf("impossible");else cout<<ans; return 0; }