【极值问题】【CF1063B】 Labyrinth

传送门

Description

给你一个$n~\times~m$的矩阵，一开始你在第$r$行第$c$列。你的上下移动不受限制，向左最多移动$x$次，向右最多移动$y$次。求你最多能到多少个点。包括起始点。

Input

第一行是$n$和$m$，代表矩阵规模。

第二行是$r$和$c$，代表你的位置

第三行是$x$和$y$，代表移动限制

下面$n$行每行$m$个字符，有且仅有'.'和''两种。如果第$i$行第$j$列是''代表你不能经过这个点。

Output

输出一行一个数代表能到的最多点数

Sample Input

4 5

3 2

1 2

.....

.***.

...**

*....

Sample Output

Hint

$For~All:$

$0~\leq~n,m,r,c~\leq~2000$,$0~\leq~x,y~\leq~10^9$

Solution

朴素的$bfs$显然是对的，可以状态太多存不下。

考虑如果从$(sx,sy)$点到一个点$(x,y)$时，假设共向右走了$r$步，向左走了$l$步，显然$r-l$是一个定值。具体的，$r-l~=~y-sy$。于是，对于任意一个目标$(x,y)$，发现$l$事实上与$r$线性正相关。对于一个点，显然到该点的$r$越小越好，同时由于$l$和$r$线性正相关，所以最小化$r$的同时，$l$已经被最小化了。于是可以直接建图跑最短路，所有向右的边权为1，其他边权为0。跑完后扫描整个地图就可以判断合法性了。

这里的一个新姿势是$0/1bfs$。当边权只有$0/1$时，可以使用双端队列进行bfs，具体的，当当前边的权值时$0$时，将终点压入队首，否则压入队尾。考虑这么做的正确性：易证任意一时刻队列中的点dist差值不超过1。于是正确性显然。$0/1bfs$的复杂度为$O(V+E)$。相比dij少了一个log。

Code

#include<queue>

#include<cstdio>

#include<cstring>

#define rg register

#define ci const int

#define cl const long long int

typedef long long int ll;

namespace IO{

    char buf[110];

}

template <typename T>

inline void qr(T &x) {

    rg char ch=getchar(),lst=' ';

    while((ch > '9') || (ch < '0')) lst=ch,ch=getchar();

    while((ch >= '0') && (ch <= '9')) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();

    if(lst == '-') x=-x;

}

template <typename T>

inline void qw(T x,const char aft,const bool pt) {

    if(x < 0) {putchar('-');x=-x;}

    rg int top=0;

    do {

        IO::buf[++top]=x%10+'0';

    } while(x/=10);

    while(top) putchar(IO::buf[top--]);

    if(pt) putchar(aft);

}

template <typename T>

inline T mmax(const T a,const T b) {return a > b ? a : b;}

template <typename T>

inline T mmin(const T a,const T b) {return a < b ? a : b;}

template <typename T>

inline T mabs(const T a) {return a < 0 ? -a : a;}

template <typename T>

inline void mswap(T &a,T &b) {

    T _temp=a;a=b;b=_temp;

}

const int maxn = 2010;

const int fx[]={0,-1,0,1};

const int fy[]={1,0,-1,0};

const int fv[]={1,0,0,0};

int n,m,sx,sy,x,y;

int MU[maxn][maxn];

char mp[maxn][maxn];

std::deque<int>Qx,Qy;

int main() {

    qr(n);qr(m);qr(sx);qr(sy);qr(x);qr(y);

    for(rg int i=1;i<=n;++i) scanf("%s",mp[i]+1);

    memset(MU,0x3f,sizeof MU);

    MU[sx][sy]=0; Qx.push_front(sx);Qy.push_front(sy);

    while(!Qx.empty()) {

        int hx=Qx.front(),hy=Qy.front();Qx.pop_front();Qy.pop_front();

        for(rg int i=0;i<4;++i) {

            int dx=hx+fx[i],dy=hy+fy[i];

            if((dx > n) || (dy > m) || (!dx) || (!dy) || (mp[dx][dy] == '*') || (MU[dx][dy] <= MU[hx][hy]+fv[i])) continue;

            MU[dx][dy]=MU[hx][hy]+fv[i];

            if(i) {Qx.push_front(dx);Qy.push_front(dy);}

            else {Qx.push_back(dx);Qy.push_back(dy);}

        }

    }

    rg int _ans=0;

    for(rg int i=1;i<=n;++i) {

        for(rg int j=1;j<=m;++j) if(mp[i][j] != '*') {

            if((MU[i][j] <= y) && ((MU[i][j]-j+sy) <= x)) ++_ans;

        }

    }

    qw(_ans,'\n',true);

    return 0;

}

Summary

当题设需要最小化多个变量时，不妨考虑变量间的相关关系，从此转化成单变量的极值问题。

当边权只有$0$和$1$的时候，可以考虑使用$0/1bfs$,省去dij的log。