BestCoder Round #63 (div.1) B.matrix

matrix

 

 Accepts: 124

 

 Submissions: 171

 Time Limit: 6000/3000 MS (Java/Others)

 

 Memory Limit: 65536/65536 K (Java/Others)

问题描述

给定$n\ast m$($n+m$为奇数)的矩阵，从$(1,1)$走到$(n,m)$且只能往右往下走，设经过的数为a_1, a_2 ... a_{2k}a​1​​,a​2​​...a​2k​​,贡献为a_1*a_2+a_3*a_4+...+a_{2k-1}*a_{2k}a​1​​∗a​2​​+a​3​​∗a​4​​+...+a​2k−1​​∗a​2k​​，求最小贡献。

输入描述

若干组数据(大概$5$组)。
每组数据第一行两个数$n,m(1\le n,m\le 1000$且$n+m$为奇数)。
接下来$n$行每行$m$个数a_i, j(1 \leq ai,j \leq 100)a​i​​,j(1≤ai,j≤100)描述这个矩阵。

输出描述

对于每组数据，输出一行表示答案。

输入样例

2 3
1 2 3
2 2 1
2 3
2 2 1
1 2 4

输出样例

4
8

求这种题马上想到dp。

我用的是记忆化搜索。

重(1, 1)到(n, m)的最小贡献度

(n, m)只能从(n-1, m)和(n, m-1)来

我们就可以用这种思想求

#include
     
     
      
      
#include
      
      
       
       
#include
       
       
        
        
#include
        
        
         
         
using namespace std;

const int N = 1001;
long long a[N][N];
long long dp[N][N];

void debug(int n, int m)
{
    for(int i = 1; i <= n; i++)
    {
        for(int j = 1; j <= m; j++)
            printf("%I64d ", dp[i][j]);
        printf("\n");
    }
}

long long dfs(int n, int m)
{
    if(dp[n][m]) return dp[n][m];
    dp[n][m] = 1234567890;
    if(n == 1 && m == 1)
        return dp[1][1] = 0;
    if(n > 1)
    {
        if( (n+m+1) % 2 == 0)
            dp[n][m] = min(dp[n][m], dfs(n-1, m)+a[n-1][m]*a[n][m]);
        else
            dp[n][m] = min(dp[n][m], dfs(n-1, m));
    }
    if(m > 1)
    {
        if( (n+m+1) % 2 == 0)
            dp[n][m] = min(dp[n][m], dfs(n, m-1)+a[n][m-1]*a[n][m]);
        else
            dp[n][m] = min(dp[n][m], dfs(n, m-1));
    }
    return dp[n][m];
}

int main(void)
{
    int n, m;
    while(~scanf("%d%d", &n, &m))
    {
        int i, j;
        memset(dp, 0, sizeof(dp));
        for(i = 1; i <= n; i++)
            for(j = 1; j <= m; j++)
                scanf("%I64d", &a[i][j]);
        printf("%I64d\n", dfs(n, m));
//        debug(n, m);
    }
    return 0;
}