[codeforces232E]Quick Tortoise

题目大意

给定一张n*m的地图，每个位置(i,j)可能是空地和障碍物。有q个询问，每次询问给出四个数x1,y1,x2,y2，问从(x1,y1)是否能走到(x2,y2)（只能向右或向下走到空地，保证两个位置都是空地）

数据范围

1≤n,m≤500
1≤q≤600,000

分析

点比较多，无法储存每对点之间的连通性，而且似乎用图论的知识难以下手。
考虑离线去做。
假如现在确定了一个数y，然后处理所有跨过第y列的询问。
从第y列的某个点出发，求出哪些点可以到达它、它可以到达哪些点，然后对于每个询问，判断一下能否经过这个点到达目的地。
这样做是很慢的，所以考虑用二进制压缩第y列的所有点的状态，然后做一次dp。如果是第y列左边的点，则f[i][j]=f[i+1][j] or f[i][j+1]，如果是右边的点，则f[i][j]=f[i-1][j] or f[i][j-1]，就可以处理它能到达第y列哪些点（或哪些点到达它）。
然后对于每个询问，把f数组and起来，判断是否为0即可。

#include <cstdio>
#include <cstring>
#include <algorithm>

#define X1 q[e[i]].x1
#define Y1 q[e[i]].y1
#define X2 q[e[i]].x2
#define Y2 q[e[i]].y2

using namespace std;

const int maxn=505,maxm=600005;

typedef long long LL;

int n,m,Q,tot,f[maxn][maxn][18],N,h[maxn<<2],e[maxm],next[maxm],tmp[18];

bool ans[maxm];

char map[maxn][maxn];

struct query
{
    int x1,y1,x2,y2;
}q[maxm];

void add(int x,int y)
{
    e[++tot]=y; next[tot]=h[x]; h[x]=tot;
}

int init(int l,int r,int x,int y1,int y2)
{
    int mid=l+r>>1;
    if (y1<=mid && mid<=y2) return x;
    if (y2<mid) init(l,mid-1,x<<1,y1,y2);
    else init(mid+1,r,(x<<1)+1,y1,y2);
}

void solve(int l,int r,int x)
{
    int mid=l+r>>1;
    if (l<mid) solve(l,mid-1,x<<1);
    if (mid<r) solve(mid+1,r,(x<<1)+1);
    int i,j,k;
    for (i=1;i<=n;i++)
        for (j=l;j<=r;j++)
            for (k=0;k<=N;k++) f[i][j][k]=0;
    for (i=1;i<=n;i=j+1)
    {
        memset(tmp,0,sizeof(tmp));
        for (j=i;j<=n && map[j][mid]=='.';j++)
        {
            tmp[j/30]+=(1<<(j%30));
            memcpy(f[j][mid],tmp,sizeof(tmp));
        }
    }
    for (j=mid+1;j<=r;j++)
        for (i=1;i<=n;i++) if (map[i][j]=='.')
            for (k=0;k<=N;k++) f[i][j][k]=f[i-1][j][k]|f[i][j-1][k];
    for (i=1;i<=n;i=j+1)
    {
        memset(tmp,0,sizeof(tmp));
        for (j=i;j<=n && map[j][mid]=='.';j++);
        for (k=j-1;k>=i;k--)
        {
            tmp[k/30]+=(1<<(k%30));
            memcpy(f[k][mid],tmp,sizeof(tmp));
        }
    }
    for (j=mid-1;j>=l;j--)
        for (i=n;i;i--) if (map[i][j]=='.')
            for (k=0;k<=N;k++) f[i][j][k]=f[i+1][j][k]|f[i][j+1][k];
    for (i=h[x];i;i=next[i])
    {
        if (Y1==mid || Y2==mid)
        {
            if (Y1==mid && (f[X2][Y2][X1/30]&(1<<(X1%30)))>0 ||
            Y2==mid && (f[X1][Y1][X2/30]&(1<<(X2%30)))>0) ans[e[i]]=1;
            continue;
        }
        for (j=0;j<=N;j++)
        if ((f[X1][Y1][j]&f[X2][Y2][j])>0)
        {
            ans[e[i]]=1; break;
        }
    }
}

int main()
{
    scanf("%d%d",&n,&m); N=n/30;
    for (int i=1;i<=n;i++) scanf("%s",map[i]+1);
    scanf("%d",&Q);
    for (int i=0;i<Q;i++)
    {
        scanf("%d%d%d%d",&q[i].x1,&q[i].y1,&q[i].x2,&q[i].y2);
        add(init(1,m,1,q[i].y1,q[i].y2),i);
    }
    solve(1,m,1);
    for (int i=0;i<Q;i++)
        if (ans[i]) printf("Yes\n");else printf("No\n");
    return 0;
}