方格取数
Time Limit : 3000/1000ms (Java/Other) Memory Limit : 65535/32768K (Java/Other)
Total Submission(s) : 8 Accepted Submission(s) : 5
Font: Times New Roman | Verdana | Georgia
Font Size: ← →
Problem Description
0 0 0 0 0 0 0 0
0 0 13 0 0 6 0 0
0 0 0 0 7 0 0 0
0 0 0 14 0 0 0 0
0 21 0 0 0 4 0 0
0 0 15 0 0 0 0 0
0 14 0 0 0 0 0 0
0 0 0 0 0 0 0 0
某人从图中的左上角出发,可以向下行走,也可以向右行走,直到到达右下角。在走过的路上,他可以取走方格中的数(取走后的方格中将变为数字0)。
此人从左上角到右下角共走了两次,试找出两条这样的路径,使得取得的数字和为最大。
Input
接下来的每行有三个整数,第一个为行号数,第二个为列号数,第三个为在该行、该列上所放的数。一行0 0 0表示结束。
Output
Sample Input
8 2 3 13 2 6 6 3 5 7 4 4 14 5 2 21 5 6 4 6 3 15 7 2 14 0 0 0
Sample Output
67
Author
// 二取方格数 多进程dp
我们令 dp[x1,y1,x2,y2]表示第一条路径走到了(x1,y1),第二条路径走到了(x2,y2)是的最优解,可以很自然的得到状态转移方程:
设p=max(sum[i1-1,j1,i2-1,j2],sum[i1-1,j1,i2,j2-1],sum[i1,j1-1,i2-1,j2],sum[i1,j1-1,i2,j2-1])
dp[x1][y1][x2][y2]= 0 (x1||y1||x2||y2==0)
= p+data[x1][y1] (x1==x2&&y1==y2)
= p+data[x1][y1]+data[x2][y2] (x1!=x2||y1!=y2)
#include<cstdio>
#include<cstring>
int dp[15][15][15][15],p[15][15];
int main()
{
int n,i1,j1,i2,j2,x,y,z;
while(scanf("%d",&n)!=EOF)
{
memset(dp,0,sizeof(dp));
memset(p,0,sizeof(p));
while(1)
{
scanf("%d%d%d",&x,&y,&z);
if(x+y+z==0) break;
p[x][y]=z;
}
for(i1=1;i1<=n;i1++)
for(j1=1;j1<=n;j1++)
for(i2=1;i2<=n;i2++)
for(j2=1;j2<=n;j2++)
{
if(dp[i1-1][j1][i2-1][j2]>dp[i1][j1][i2][j2])
dp[i1][j1][i2][j2]=dp[i1-1][j1][i2-1][j2];
if(dp[i1-1][j1][i2][j2-1]>dp[i1][j1][i2][j2])
dp[i1][j1][i2][j2]=dp[i1-1][j1][i2][j2-1];
if(dp[i1][j1-1][i2-1][j2]>dp[i1][j1][i2][j2])
dp[i1][j1][i2][j2]=dp[i1][j1-1][i2-1][j2];
if(dp[i1][j1-1][i2][j2-1]>dp[i1][j1][i2][j2])
dp[i1][j1][i2][j2]=dp[i1][j1-1][i2][j2-1];
dp[i1][j1][i2][j2]+=p[i1][j1];
if(i1!=i2||j1!=j2)
dp[i1][j1][i2][j2]+=p[i2][j2];
}
printf("%d\n",dp[n][n][n][n]);
}
return 0;
}
上面是四维的解法,当n很大时可能会超出内存,所以为了解决这个问题可以把四维压缩成3维节约空间,并可以用滚动数组进一步优化。
因为两条路径同时从(1,1)出发,那么当两条路径都已经包含k个方格时,两条路径的末端必同在矩阵的第k条对角线上,
因为对于每一条路径,向右延伸的格子数+向下延伸的格子数=k,也就是末端两个格子的纵横坐标之和=k
所以我们只需要知道两路径末端所在的行编号x1,x2以及两末端所在对角线编号k,就可以确定末端节点的位置(x1,k+1-x1),(x2,k+1-x2)。很显然,路径每延伸一个格子所能得到的最优值,只与上次延伸的最优值有关,那么我们就可以以对角线划分阶段进行dp。至此,我们可以建立新的状态转移方程:
dp[x1,x2,k]表示第一个路径末端在第x1行,第二路径末端在第x2行,且两末端同在第k条对角线上时的最优解,所以
令 p=max{dp[x1-1,x2,k-1],dp[x1,x2-1,k-1],
dp[x1-1,x2-1,k-1],dp[x1,x2,k-1]}
dp[x1,x2,k]={ p+map[x1,k+1-x1]; (x1=x2);
p+map[x1,k+1-x1]+map[x2,k+1-x2]; else
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int INF = 999999999;
int dp[22][12][12],p[12][12];
int main()
{
//freopen("a.txt","r",stdin);
int n,i,j,k,a,b,c,x;
while(scanf("%d",&n)!=EOF)
{
memset(dp,0,sizeof(dp));
memset(p,0,sizeof(p));
while(1)
{
scanf("%d%d%d",&a,&b,&c);
if(a==0&&b==0&&c==0)break;
p[a][b]=c;
}
//t=0; //dp[1][1][0]=p[1][1];
for(k=1;k<=2*n;k++)
{
//t=1-t;
for(i=1;i<=k;i++)
for(j=1;j<=k;j++)
{
x=-INF;
x=max(x,dp[k-1][i][j-1]);
x=max(x,dp[k-1][i-1][j]);
x=max(x,dp[k-1][i-1][j-1]);
x=max(x,dp[k-1][i][j]);
if(i==j)
dp[k][i][j]=x+p[i][k+1-i];
else dp[k][i][j]=x+p[i][k+1-i]+p[j][k+1-j];
}
}
printf("%d\n",dp[2*n][n][n]);
}
return 0;
}
下面用滚动数组:
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int INF = 999999999;
int dp[2][12][12],p[12][12];
int main()
{
//freopen("a.txt","r",stdin);
int n,i,j,k,a,b,c,x;
while(scanf("%d",&n)!=EOF)
{
memset(dp,0,sizeof(dp));
memset(p,0,sizeof(p));
while(1)
{
scanf("%d%d%d",&a,&b,&c);
if(a==0&&b==0&&c==0)break;
p[a][b]=c;
}
int t=0;
for(k=1;k<=2*n;k++)
{
t=1-t;
for(i=1;i<=k;i++)
for(j=1;j<=k;j++)
{
x=-INF;
x=max(x,dp[1-t][i][j-1]);
x=max(x,dp[1-t][i-1][j]);
x=max(x,dp[1-t][i-1][j-1]);
x=max(x,dp[1-t][i][j]);
if(i==j)
dp[t][i][j]=x+p[i][k+1-i];
else dp[t][i][j]=x+p[i][k+1-i]+p[j][k+1-j];
}
}
printf("%d\n",dp[t][n][n]);
}
return 0;
}