lightoj 1036 - A Refining Company （dp）

题意:在一个m x n的地图里有金矿和银矿，当你在某个点（x, y）时， 你可以建造一个只能向西或者向北挖取矿石的履带车，金矿石只能被向北移动的履带车挖取，银矿石只能被向西移动的履带车挖取，并且履带车不能交叉。给出你每个地点的金矿和银矿的数量为你用履带车挖取的最大矿石数量是多少.（履带车会一直走到边界）

思路:

当你位于（x, y）时， 你只有两个操作：

1）建立一个向北运输的履带车，挖一路上的金矿石。

2）建立一个想西运输的履带车，挖一路上的银矿石。

定义dp[i][j] 从(i, j) 出发可以得到的最大矿石数量

于是有状态方程:dp[i][j] = max(dp[i-1][j] + a[0][i][j], dp[i][j-1] + a[1][i][j])

a[0][i][j]表示从(i, j)向西出发的履带车可以挖取的银矿石数量。

a[1][i][j]表示从(i, j)向北出发的履带车可以挖取的金矿石数量。

/***********************************************
 * Author: fisty
 * Created Time: 2015-08-19 19:48:34
 * File Name   : 1036.cpp
 *********************************************** */
#include <iostream>
#include <cstring>
#include <deque>
#include <cmath>
#include <queue>
#include <stack>
#include <list>
#include <map>
#include <set>
#include <string>
#include <vector>
#include <cstdio>
#include <bitset>
#include <algorithm>
using namespace std;
#define Debug(x) cout << #x << " " << x <<endl
#define Memset(x, a) memset(x, a, sizeof(x))
const int INF = 0x3f3f3f3f;
typedef long long LL;
typedef pair<int, int> P;
#define FOR(i, a, b) for(int i = a;i < b; i++)
#define lson l, m, k<<1
#define rson m+1, r, k<<1|1
#define MAX_N 550
int G[2][MAX_N][MAX_N];
int a[2][MAX_N][MAX_N];
int t;
int m, n;
int dp[MAX_N][MAX_N];
void init(){
    for(int i = 1;i <= m; i++){
        for(int j = 1;j <= n; j++){
            a[0][i][j] = a[0][i][j-1] + G[0][i][j];
        }
    }
    for(int i = 1;i <= m; i++){
        for(int j = 1;j <= n; j++){
            a[1][i][j] = a[1][i-1][j] + G[1][i][j];
        }
    }
}
int main() {
    //freopen("in.cpp", "r", stdin);
    //cin.tie(0);
    //ios::sync_with_stdio(false);
    scanf("%d", &t);
    int cnt = 1;
    while(t--){
        scanf("%d%d", &m, &n);
        for(int z = 0;z < 2; z++){
            for(int i = 1;i <= m; i++){
                for(int j = 1;j <= n; j++){
                    scanf("%d", &G[z][i][j]);
                }
            }
        }
        init();
        for(int i = 1;i <= m; i++){
            for(int j = 1;j <= n; j++){
                dp[i][j] = max(dp[i][j-1] + a[1][i][j], dp[i-1][j] + a[0][i][j]);
            }
        }
        printf("Case %d: %d\n",cnt++, dp[m][n]);
    }
    return 0;
}