NYOJ 712 探寻宝藏 &NYOJ 61 传纸条 双线程DP

 NYOJ 61  传纸条

状态:dp(A,B)

转移方程：dp(A,B)=max{  dp(A1,B1),dp(A2,B2),dp(A1,B2),dp(A2,B1) }+map(A)+map(B);

如图所示：

依据此可以写出.

提交运行结果如下:

限定A在B下

#include <stdio.h>
#include <string.h>
int dp[53][53][53][53];
int map[53][53];
int max(int a,int b)
{
  return a>b? a:b;   
}
 
int main()
{
  int T,row,col,i,j,ii,jj,k,t1,t2,ans;
  
  scanf("%d",&T);
   
    do
    {
       scanf("%d%d",&row,&col);         
        
       for(i = 1 ; i <= row ; ++i)         
       for(j = 1 ; j <= col ; ++j)
       scanf("%d",&map[i][j]);  
       
       for(i = 1 ; i <=row ; ++i)
       for(j = 1 ; j <= row ; ++j)
       for(ii = 1 ; ii <= col ; ++ii)
       for(jj = 1 ; jj <= col ; ++jj)
       if(i > j &&ii <jj)
       {
             t1 = max(dp[i][ii-1][j][jj-1],dp[i][ii-1][j-1][jj]);    
             t2 = max(dp[i-1][ii][j][jj-1],dp[i-1][ii][j-1][jj]); 
             dp[i][ii][j][jj] = max(t1,t2) + map[i][ii] + map[j][jj];
       }     
              
       ans = max(dp[row][col-1][row-1][col],dp[row-1][col][row][col-1]);      
       printf("%d\n",ans);      
                  
    }while(--T);
      
 // system("pause");
  return 0;   
}

注意到题中条件
ii+i=jj+j=k;
这样可以进行优化，将复杂度为O(n^3)化为O(n^2)
A(i,k-i),B(j,k-j)
若A在B上,则i>j (k-i>=0 ,k-j>=0)

#include <stdio.h>
#include <string.h>
int dp[105][53][53];
int map[53][53];
int max(int a,int b)
{
  return a>b? a:b;   
}
 
int main()
{
  int T,row,col,i,j,all,k,t1,t2,ans;
  scanf("%d",&T);
   
    do
    {
       scanf("%d%d",&row,&col);         
        
       for(i = 1 ; i <= row ; ++i)         
       for(j = 1 ; j <= col ; ++j)
       scanf("%d",&map[i][j]);  
        
       all = row + col ;
       for(k = 2 ; k <= all ; ++k)
       for(i = 1 ; i <=row ; ++i)
       for(j = 1 ; j <= row ; ++j)
       if(i >j && k >= i && k >= j)
       {
             t1 = max(dp[k-1][i][j],dp[k-1][i][j-1]);    
             t2 = max(dp[k-1][i-1][j-1],dp[k-1][i-1][j]); 
             dp[k][i][j] = max(t1,t2) + map[i][k-i] + map[j][k-j];
       }     
              
       ans = max(dp[all-1][row-1][row],dp[all-1][row][row-1]);   
           
       printf("%d\n",ans);      
                  
    }while(--T);
    
  return 0;   
}

NYOJ  712探 寻 宝 藏

    题意：给一个矩阵，求两条不相交的线从左上角到右下角经过的元素的最大和

 

    双线dp - 即同时考虑两条不相交的线，使其线上的和最大

    显然我们需要记录每一步时两个线同时往前走的位置(这样可以容易的控制其不相交).

    状态：dp[k,(x1,y1),(x2,y2)] 在第k步，双线里一线在(x1,y1) 二线在(x2,y2) 的最大和

    转移：dp[k,(x1,y1),(x2,y2)] =max(dp[k-1,(x1-1,y1),(x2-1,y2)],

                                     dp[k-1,(x1,y1-1),(x2-1,y2)],

                                      dp[k-1,(x1-1,y1),(x2,y2-1)],

                                     dp[k-1,(x1,y1-1),(x2,y2-1)]) + map[x1][y1] + map[x2][y2];

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;

template <typename _T>
inline _T Max(_T a , _T b){    return (a<b)?(b):(a);    }
template <typename _T>
inline _T Max(_T a , _T b , _T c , _T d){
    return (Max(a,b) < Max(c,d)) ? (Max(c,d)) : (Max(a,b));
}

const int M = 52;
int map[M][M];
int dp[M+M][M][M]; // dp[k][i][j]: 第k步时，双线里 一线在i位置（横坐标为i，纵坐标为k-i） 另一线在j位置（横坐标为j，纵坐标为k-j）的最大和

int main()
{
    freopen("in.txt","r",stdin);
    int T;cin >> T;
    while(T--){
        int row , col;
        cin >> row >> col;
        memset(dp,0,sizeof(dp));
        for(int i = 1 ; i <= row ; i++)
            for(int j = 1 ; j <= col ; j++)
                scanf("%d" , &map[i][j]);
		int beg = map[1][1];
		int end = map[row][col];
		map[1][1] = map[row][col] = 0;
        int bound = col + row;
        for(int k = 1 ; k <= bound ; k++){ // 起点和终点为唯一交叉处，特判
            for(int x1 = 1 ; x1 <= row ; x1++){
                for(int x2 = 1 ; x2 <= row ; x2++){
                    if (x1 == x2 && !(k==bound && x1==row))continue;
                    if ( !(1<=x1 && x1<=row) || !(1<=x2 && x2<=row))continue;
                    int y1 = k - x1;
                    int y2 = k - x2;
                    if ( !(1<=y1 && y1<=col) || !(1<=y2 && y2<=col))continue;
                    dp[k][x1][x2] = Max(dp[k-1][x1-1][x2-1],
                                        dp[k-1][ x1 ][x2-1],
                                        dp[k-1][x1-1][ x2 ],
                                        dp[k-1][ x1 ][ x2 ]);
                    dp[k][x1][x2] += map[x1][y1] + map[x2][y2];
                }
            }
        }
        cout << beg + dp[bound][row][row] + end << endl;
    }
    return 0;
}