Analysis
考虑两条来回的不重合的路径等价于两条从左上角走到右下角不重合的路径
如果只求一条从左上走到右下的路径,,使其路径权值和最大,这就是很基础的dp了
f
[
i
]
[
j
]
=
m
a
x
(
f
[
i
−
1
]
[
j
]
,
f
[
i
]
[
j
−
1
]
)
+
a
[
i
]
[
j
]
f[i][j]=max(f[i-1][j],f[i][j-1])+a[i][j]
f[i][j]=max(f[i−1][j],f[i][j−1])+a[i][j]
现在要求两条,考虑定义四维数组,一条路走到(x1,y1),另一条走到(x2,y2)的可获得的最大权值
f
[
x
1
]
[
y
1
]
[
x
2
]
[
y
2
]
=
m
a
x
(
f
[
x
1
]
[
y
1
−
1
]
[
x
2
]
[
y
2
−
1
]
,
f[x1][y1][x2][y2]=max(f[x1][y1-1][x2][y2-1],
f[x1][y1][x2][y2]=max(f[x1][y1−1][x2][y2−1],
f
[
x
1
]
[
y
1
−
1
]
[
x
2
−
1
]
[
y
2
]
,
f[x1][y1-1][x2-1][y2],
f[x1][y1−1][x2−1][y2],
f
[
x
1
−
1
]
[
y
1
]
[
x
2
]
[
y
2
−
1
]
,
f[x1-1][y1][x2][y2-1],
f[x1−1][y1][x2][y2−1],
f
[
x
1
−
1
]
[
y
1
]
[
x
2
−
1
]
[
y
2
]
)
f[x1-1][y1][x2-1][y2])
f[x1−1][y1][x2−1][y2])
+
a
[
x
1
]
[
y
1
]
+
a
[
x
2
]
[
y
2
]
+a[x1][y1]+a[x2][y2]
+a[x1][y1]+a[x2][y2]
现在
n
,
m
<
=
50
n,m<=50
n,m<=50,可以水过去。比较慢。但思路简单
需要注意x2的枚举要从x1+1开始,这样就可以避免走重复的情况了
最后输出的时候答案是存在
f
[
n
−
1
]
[
m
]
[
n
]
[
m
−
1
]
f[n-1][m][n][m-1]
f[n−1][m][n][m−1]里面的,不要写错了(因为这个贡献了一发WA)
再考虑另外一种解法
两张纸条同时传的话,同一时刻它们穿过的人数是相同的,根据这点可以进行降维。
令f[i][x1][x2]表示第i步两张纸条的x坐标分别是x1和x2,则y1=i-x1+2,y2=i-x2+2。
这样就可以得出dp方程
f
[
i
]
[
x
1
]
[
x
2
]
=
m
a
x
(
f
[
i
−
1
]
[
x
1
]
[
x
2
]
,
f
[
i
−
1
]
[
x
1
−
1
]
[
x
2
]
,
f
[
i
−
1
]
[
x
1
]
[
x
2
−
1
]
,
f
[
i
−
1
]
[
x
1
−
1
]
[
x
2
−
1
]
)
+
a
[
x
1
]
[
i
−
x
1
+
2
]
+
a
[
x
2
]
[
i
−
x
2
+
2
]
f[i][x1][x2] = max( f[i-1][x1][x2], f[i-1][x1-1][x2], f[i-1][x1][x2-1], f[i-1][x1-1][x2-1] )+ a[x1][i-x1+2] + a[x2][i-x2+2]
f[i][x1][x2]=max(f[i−1][x1][x2],f[i−1][x1−1][x2],f[i−1][x1][x2−1],f[i−1][x1−1][x2−1])+a[x1][i−x1+2]+a[x2][i−x2+2]
注意答案在 f [ n + m − 3 ] [ n − 1 ] [ n ] f[n+m-3][n-1][n] f[n+m−3][n−1][n]里,因为 f [ n + m − 2 ] [ n ] [ n ] f[n+m-2][n][n] f[n+m−2][n][n]在我们计算的时候是算做不合法情况,没有被纳入考虑
Code
优秀的方法二:
#include<bits/stdc++.h>
using namespace std;
int n,m,a[55][55],f[105][55][55];
int main(){
scanf("%d%d",&n,&m);
int i,j,x1,y1,x2,y2;
for(i=1;i<=n;++i) for(j=1;j<=m;++j) scanf("%d",&a[i][j]);
for(i=1;i<=n+m-2;++i){//枚举走的步数
for(x1=max(1,i+2-m);x1<=min(i+1,n);++x1){
y1=i-x1+2;
for(x2=max(x1+1,i+2-m);x2<=min(i+1,n);++x2){
y2=i-x2+2;
f[i][x1][x2]=max(max(f[i-1][x1][x2],f[i-1][x1-1][x2-1]),max(f[i-1][x1-1][x2],f[i-1][x1][x2-1]))+a[x1][y1]+a[x2][y2];
}
}
}
cout<<f[n+m-3][n-1][n];
return 0;
}
简单的方法一:
#include<bits/stdc++.h>
using namespace std;
int n,m,a[55][55],dp[55][55][55][55];
inline int Max(int a,int b,int c,int d){
return max(max(a,b),max(c,d));
}
int main(){
scanf("%d%d",&n,&m);
int i,j,x1,y1,x2,y2;
for(i=1;i<=n;++i) for(j=1;j<=m;++j) scanf("%d",&a[i][j]);
for(x1=1;x1<=n;++x1){
for(y1=1;y1<=m;++y1)
for(x2=x1+1;x2<=n;++x2)
for(y2=1;y2<=m;++y2)
dp[x1][y1][x2][y2] = Max(dp[x1][y1-1][x2][y2-1],dp[x1-1][y1][x2-1][y2],dp[x1-1][y1][x2][y2-1],dp[x1][y1-1][x2-1][y2])+ a[x1][y1] + a[x2][y2];
}
cout<<dp[n-1][m][n][m-1];
return 0;
}