方格取数

方格取数

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

设有N*N的方格中，我们在其中的某些方格中填入正整数，而其他的方格中则放入数字0，如下表所示：
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

有多组输入数据，每组数据的第一行为一个整数N(N<=10)，表示N*N的方格图。
接下来的每行有三个整数，第一个为行号数，第二个为列号数，第三个为在该行、该列上所放的数。一行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

HYNU

// 二取方格数     多进程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;
}