AtCoder Grand Contest 012 E - Camel and Oases 状压dp

题意

平面上有n个点。初始有V的权值，每次可以从一个点走到与他距离不超过V的点，当V>0时也可以让V/2且到达任意一个点。问从每个点出发能否遍历所有点。
n,V<=200000

分析

显然只有$O(logV)$种不同的步数。
当权值$V$被固定后，整条链被分成若干个区间，每个区间内的点可以互达。
把每一种$V$对应的所有区间求出来后，不难发现两个区间要么是相离关系，要么是包含关系。
从某个点出发能否遍历所有点，就相当于把该点能通过权值$V$到达的点去掉，问能否用除了$V$以外的权值对应的区间去覆盖整条线段。
这个可以状压dp：
设$f[S]$表示用$S$里面的权值从左往右覆盖数轴，最大能到达的右端点。
$g[S]$表示最小能到达的左端点。
然后答案就很好求了。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define MAX(x,y) x=max(x,y)
#define MIN(x,y) x=min(x,y)
#define inf 0x3f3f3f3f
using namespace std;

const int N=400005;

int n,v,g[N*2],f[N*2],bin[25],tot,len[25],lef[25][N],rig[25][N],a[N],bel[N],p1[N],p2[N];
bool ans[N];

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;
}

void dp()
{
    bin[0]=1;
    for (int i=1;i<=tot;i++) bin[i]=bin[i-1]*2;
    for (int i=0;i<bin[tot];i++) g[i]=n+1;
    for (int i=0;i<bin[tot];i++)
        for (int j=0;j<=tot;j++)
            if (!(i&bin[j])) MAX(f[i+bin[j]],rig[j][f[i]+1]),MIN(g[i+bin[j]],lef[j][g[i]-1]);
    for (int i=0;i<bin[tot];i++)
    {
        int l=f[i],r=g[i^(bin[tot]-1)];
        if (l+1>=r)
        {
            for (int j=1;j<=n;j++) ans[j]=1;
            return;
        }
        if (p1[bel[l+1]]<=l+1&&p2[bel[l+1]]>=r-1) ans[bel[l+1]]=1;
    }
}

int main()
{
    n=read();v=read();
    for (int i=1;i<=n;i++) a[i]=read();
    tot=0;len[0]=0;
    while (v) len[++tot]=v,v>>=1;
    sort(len,len+tot+1);
    memset(lef,inf,sizeof(lef));
    for (int i=0;i<tot;i++)
        for (int j=1;j<=n;j++)
        {
            int l=j;
            while (l<n&&a[l+1]-a[l]<=len[i]) l++;
            rig[i][j]=l;lef[i][l]=j;j=l;
        }
    int id=0;
    for (int j=1;j<=n;j++)
    {
        int l=j;bel[j]=++id;
        while (l<n&&a[l+1]-a[l]<=len[tot]) l++,bel[l]=id;
        p1[id]=j;p2[id]=l;j=l;
    }
    dp();
    for (int i=1;i<=n;i++) puts(ans[bel[i]]?"Possible":"Impossible");
    return 0;
}