AcWing 1018. 最低通行费（线性dp 数字三角形模型）

一个商人穿过一个 N×N 的正方形的网格，去参加一个非常重要的商务活动。

他要从网格的左上角进，右下角出。

每穿越中间 1 个小方格，都要花费 1 个单位时间。

商人必须在 (2N−1) 个单位时间穿越出去。

而在经过中间的每个小方格时，都需要缴纳一定的费用。

这个商人期望在规定时间内用最少费用穿越出去。

请问至少需要多少费用？

注意：不能对角穿越各个小方格（即，只能向上下左右四个方向移动且不能离开网格）。

输入格式
第一行是一个整数，表示正方形的宽度 N。

后面 N 行，每行 N 个不大于 100 的正整数，为网格上每个小方格的费用。

输出格式
输出一个整数，表示至少需要的费用。

数据范围
1≤N≤100

输入样例：
5
1 4 6 8 10
2 5 7 15 17
6 8 9 18 20
10 11 12 19 21
20 23 25 29 33

输出样例：
109

样例解释
样例中，最小值为 109=1+2+5+7+9+12+19+21+33。

思路：与摘花生 、数字三角形思路几乎一模一样，只不过这里我们要的是最小值

为什么这题的思路和那两道例题一致？

注意题目给了一个比较隐晦的条件：

商人必须在 (2N−1) 个单位时间穿越出去。

对于起点的(1, 1)到终点(n, n)，它们之间的曼哈顿距离是2n−2，而本题又要求我们在2n−1 的时间内，从起点走到终点。
因此，得出结论，我们的走的路线不是完全随机的，而是遵循最短路的原则走的。也就是说，每次移动，至少要使曼哈顿距离缩短1。于是，规定了我们每次在不越过边界的情况下只能向右或向下移动。

这就顺利的将本题思路和摘花生 、数字三角形联系到一起了。

状态表示和状态计算参见：摘花生 、数字三角形。

写法一：标准写法

#include<bits/stdc++.h>

using namespace std;
const int N = 110;
int w[N][N];
int dp[N][N];
int n;

int main()
{
    cin>>n;
    for(int i=1;i<=n;++i)
    {
        for(int j=1;j<=n;++j)
        {
            cin>>w[i][j];
        }
    }
        
    memset(dp,0x3f,sizeof dp);
    
    for(int i=1;i<=n;++i)
    {
        for(int j=1;j<=n;++j)
        {
            //如果dp[i][i] = a[1][1]写在外面的话会被dp[0][1]=dp[1][0]=inf更新，应该在循环内部特判一下
            if(i==1&&j==1) dp[i][j]=w[i][j];
            else dp[i][j]=min(dp[i-1][j]+w[i][j],dp[i][j-1]+w[i][j]);
        }
    }
        
    cout<<dp[n][n]<<endl;
    
    return 0;    
}

写法二：写法一稍作变形

#include<bits/stdc++.h>

using namespace std;

const int N = 100 + 10;
int n;
int w[N][N];
int dp[N][N];

int main()
{
	ios_base::sync_with_stdio(false), cin.tie(0);
	cin >> n;
	for (int i = 1; i <= n; ++i)
	{
		for (int j = 1; j <= n; ++j)
		{
			cin >> w[i][j];
		}
	}
	memset(dp, 0x3f, sizeof dp);//由于要求最小值则应当初始化为正无穷，这就相对于下面的写法减少了边界判断
	dp[1][1] = w[1][1];
	for (int i = 1; i <= n; ++i)
	{
		for (int j = 1; j <= n; ++j)
		{
			int& x = dp[i][j];//引用，相当于给变量取了个别名，x与dp[i][j]共同进退
			
			//dp[1][1] = w[1][1]语句如果不像写法一一样写循环里面特判，下面的用min更新状态就应该分步进行
			x = min(x, dp[i - 1][j] + w[i][j]);
			x = min(x, dp[i][j - 1] + w[i][j]);
		}
	}
	cout << dp[n][n] << endl;
	return 0;
}

写法三：主要是特殊判断 yls写法

#include<bits/stdc++.h>

using namespace std;

const int N = 110, INF = 1e9;

int n;
int w[N][N];
int f[N][N];

int main()
{
    scanf("%d", &n);

    for (int i = 1; i <= n; i ++ )
        for (int j = 1; j <= n; j ++ )
            scanf("%d", &w[i][j]);

    for (int i = 1; i <= n; i ++ )
        for (int j = 1; j <= n; j ++ )
            if (i == 1 && j == 1) f[i][j] = w[i][j];    // 特判左上角
            else
            {
                f[i][j] = INF;
                if (i > 1) f[i][j] = min(f[i][j], f[i - 1][j] + w[i][j]);   // 只有不在第一行的时候，才可以从上面过来
                if (j > 1) f[i][j] = min(f[i][j], f[i][j - 1] + w[i][j]);   // 只有不在第一列的时候，才可以从左边过来
            }

    printf("%d\n", f[n][n]);

    return 0;
}

写法四：记忆化搜索

#include<bits/stdc++.h>

using namespace std;

const int N = 110;
#define inf 0x3f3f3f3f;

int n;
int w[N][N];
int f[N][N];

int dp(int x, int y)
{
    if (f[x][y] >= 0) return f[x][y];
    if (x == 1 && y == 1) return w[x][y];
    if (x < 1 || y < 1) return inf;

    int t = inf;
    t = min(t, dp(x - 1, y));
    t = min(t, dp(x, y - 1));
    return f[x][y] = t + w[x][y];
}
int main()
{
    cin >> n;
    for (int i = 1; i <= n; ++ i)
    {
        for (int j = 1; j <= n; ++ j)
        {
            cin >> w[i][j];
        }
    }

    
    memset(f, -1, sizeof f);
    cout << dp(n, n) << endl;

    return 0;
}