POJ 1163 数字三角形问题（DP）

http://poj.org/problem?id=1163

1、普通递归

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

#define __max(a,b)  (((a) > (b)) ? (a) : (b))
#define MAXNUM 101

int N;
int aMax[MAXNUM][MAXNUM];						//	aMax is memorandum
int matrix[MAXNUM][MAXNUM];

int Max(int i, int j)
{
	if (i == N)
		return matrix[i][j];
	return __max( Max(i + 1, j), Max(i + 1, j + 1) ) + matrix[i][j];
}

void Input(int _matrix[MAXNUM][MAXNUM])			//	the 2nd dimension must be given!
{
	for (int i = 1; i <= N; i++) {
		for (int j = 1; j <= i; j++)
			cin >> _matrix[i][j];
	}
}

int main(void)
{
	freopen("cin.txt", "r", stdin);	
	cin >> N;
	memset(aMax, -1, sizeof(aMax));				//	0xff
	memset(matrix, -1, sizeof(matrix));
	Input(matrix);
	cout << Max(1, 1) << endl;
	
	return 0;
}

2、记忆式搜索（动态规划）

int Max(int i, int j)
{
	if (i == N)
		return matrix[i][j];
	if (aMax[i + 1][j] == -1)
		aMax[i + 1][j] = Max(i + 1, j);
	if (aMax[i + 1][j + 1] == -1)
		aMax[i + 1][j + 1] = Max(i + 1, j + 1);
	return __max( aMax[i + 1][j], aMax[i + 1][j + 1] ) + matrix[i][j];
}

3、方法2的代码优化

int Max(int i, int j)
{
	if (i == N)
		return matrix[i][j];
	if (aMax[i][j] == -1)														//  在普通递归的程序中加上
		aMax[i][j] = __max( Max(i + 1, j), Max(i + 1, j + 1) ) + matrix[i][j];	//	改
	return aMax[i][j];															//  加
}