SSL1644 取数字问题

---

原题链接

SSL1644 外网进不去

题目大意

给定$M\times N$的矩阵，其中的每个元素都是$-10$到$10$之间的整数。你的任务是从左上角$(1,1)$走到右下角$(M,N)$，每一步只能向右或向下，并且不能走出矩阵的范围。你所经过的方格里面的数字都必须被选取，请找出一条最合适的道路，使得在路上被选取的数字之和是尽可能小的正整数。

输入格式

第一行两个整数$M$，$N$，分别表示矩阵的行和列的数目。
接下来的$M$行，每行包括$N$个整数，就是矩阵中的每一行的$N$个元素。

输出格式

仅一行一个整数，表示所选道路上数字之和所能达到的最小的正整数。如果不能达到任何正整数就输出$-1$。

$\mathbf{S}\mathbf{a}\mathbf{m}\mathbf{p}\mathbf{l}\mathbf{e}$ $\mathbf{I}\mathbf{n}\mathbf{p}\mathbf{u}\mathbf{t}$

2 2
0 2
1 0

$\mathbf{S}\mathbf{a}\mathbf{m}\mathbf{p}\mathbf{l}\mathbf{e}$ $\mathbf{O}\mathbf{u}\mathbf{t}\mathbf{p}\mathbf{u}\mathbf{t}$

1

$\mathbf{H}\mathbf{i}\mathbf{n}\mathbf{t}\&\mathbf{E}\mathbf{x}\mathbf{p}\mathbf{l}\mathbf{a}\mathbf{i}\mathbf{n}$

0   2
v
1 > 0

数据范围

对于$100\%$的数据满足：$2\le N,M\le 10$。

解题思路

$dp$，还是$dp$。
很多人第一时间想到的就是从$(1,1)$开始递推，直接推到$(m,n)$，途中取最小值。但是题目要求的是求出最小的正整数，因此可以举出反例：

2 2
0 -2
1 0

此时如果按照上面的方法，路线应该是这样：

0 > -2
     v
1    0

输出-1。
但是正确的路线应该是这样的：

0   -2
v
1 >> 0

标准输出为1。
可见，上面的做法和标准答案有非常大的出入。

---

介绍完了错误的解法，来介绍一下正确的解法吧。
本题解将会从两个方面来解这道题。

---

1.$dfs$做法

我们可以用逆推，从终点$(m,n)$一直走到起点$(1,1)$。在$dfs$的途中，还要记录一个$state$变量，用来判断从当前位置推到起点，会不会有权值为$state$的路线。如果到$(1,1)$时$state$为$0$，那么就有权值为$state$的路线。
为什么可以这么做？
如果从起点走到终点，有一种解法是价值为$x$的，那么从终点倒退回来，应该价值为$0$。所以只要枚举到一个倒退回去价值为$0$的，就肯定存在这一个解法。
最后再顺便加上一个避免重复的剪枝，就可以快到飞起。
核心代码如下：

void dfs(int x,int y,int state)
{
    f[x][y][state]=1;
    if(x==1&&y==1)
        return;
    if(x>1&&!f[x-1][y][state-a[x-1][y]])
        dfs(x-1,y,state-a[x-1][y]);
    if(y>1&&!f[x][y-1][state-a[x][y-1]])
        dfs(x,y-1,state-a[x][y-1]);
    return;
}

---

2.$dp$做法

做完了$dfs$，$dp$也很容易，这里用的是顺推。
和$dfs$的做法差不多，都是看上一个位置有没有这个状态，然后再推过去。
状态转移方程：
设$f_{i,j,k}$代表从起点走到$(i,j)$有没有价值为$k$的方法，$a_{i,j}$为$(i,j)$上的数，则：
$$f_{i,j,k}=\{\begin{array}{ll}{f}_{i-1,j,k-a_{i,j}}& 1<i\le m\\ f_{i,j-1,k-a_{i,j}}& 1<j\le n\\ f_{i-1,j,k-ai,j}|f_{i-1,j,k-ai,j}& 同时满足上面两个条件\end{array}$$
注意要单独处理第一行和第一列的初始化。
核心代码：

f[1][1][base+mapp[1][1]]=true;
for(int i=2; i<=n; i++)
    for(int k=1; k<=2000; k++)
        if(f[i-1][1][k])
            f[i][1][k+mapp[i][1]]=true;
for(int i=2; i<=m; i++)
    for(int k=1; k<=2000; k++)
        if(f[1][i-1][k])
            f[1][i][k+mapp[1][i]]=true;
for(int i=2; i<=n; i++)
    for(int j=2; j<=m; j++)
        for(int k=1; k<=2000; k++)
        {
            if(f[i-1][j][k])
                f[i][j][k+mapp[i][j]]=true;
            if(f[i][j-1][k])
                f[i][j][k+mapp[i][j]]=true;
        }

---

最后，祝大家早日

上代码

代码$1(dfs)$

#include<iostream>

using namespace std;

const int     base=1000;
bool          f[20][20][2001];
int           a[20][20];
int           n,m;

void dfs(int x,int y,int state)
{
    f[x][y][state]=1;
    if(x==1&&y==1)
        return;
    if(x>1&&!f[x-1][y][state-a[x-1][y]])
        dfs(x-1,y,state-a[x-1][y]);
    if(y>1&&!f[x][y-1][state-a[x][y-1]])
        dfs(x,y-1,state-a[x][y-1]);
    return;
}

int main()
{
    cin>>n>>m;
    for(int i=1; i<=n; i++)
        for(int j=1; j<=m; j++)
            cin>>a[i][j];
    for(int i=1+base; i<=n*m*10+base; i++)
    {
        dfs(n,m,i-a[n][m]);
        if(f[1][1][base])
        {
            cout<<i-base<<endl;
            return 0;
        }
    }
    cout<<-1<<endl;
    return 0;
}

代码$2(dp)$

#include<iostream>

using namespace std;

const int      base=1000;
bool           f[20][20][2001];
int            n,m;
int            mapp[20][20];

int main()
{
    cin>>n>>m;
    for(int i=1; i<=n; i++)
        for(int j=1; j<=m; j++)
            cin>>mapp[i][j];
    f[1][1][base+mapp[1][1]]=true;
    for(int i=2; i<=n; i++)
        for(int k=1; k<=2000; k++)
            if(f[i-1][1][k])
                f[i][1][k+mapp[i][1]]=true;
    for(int i=2; i<=m; i++)
        for(int k=1; k<=2000; k++)
            if(f[1][i-1][k])
                f[1][i][k+mapp[1][i]]=true;
    for(int i=2; i<=n; i++)
        for(int j=2; j<=m; j++)
            for(int k=1; k<=2000; k++)
            {
                if(f[i-1][j][k])
                    f[i][j][k+mapp[i][j]]=true;
                if(f[i][j-1][k])
                    f[i][j][k+mapp[i][j]]=true;
            }
    for(int i=base+1; i<=base+1000; i++)
    {
        if(f[n][m][i])
        {
            cout<<i-base<<endl;
            return 0;
        }
    }
    cout<<-1<<endl;
    return 0;
}

完美切题$\sim$