HDU 6981.Rise in Price

最新推荐文章于 2021-08-27 11:47:31 发布

「已注销」

最新推荐文章于 2021-08-27 11:47:31 发布

阅读量247

点赞数 4

分类专栏： enroll

本文链接：https://blog.csdn.net/qq_46907547/article/details/119360570

版权

enroll 专栏收录该内容

35 篇文章 0 订阅

订阅专栏

Issue:

在这里插入图片描述

meaning:

给出一个n $\times$ n的矩阵，每个位置都有两个权值 $a_{ij}$ ， $b_{ij}$ ，要求从（1,1）出发，每次只能向右或者向下走，问：走到（n,n）的所有路径中， $\sum{a_{ij}}\times\sum{b_{ij}}$ 的最大值为多少。

Thinking:

n的数据范围是100；
用暴力把每一条路径都走一遍，找最优解；
考虑剪枝函数；
假设当前已经拿到的钻石 $\sum{a_{ij}}=S_a$ ，当前已经提高的价格 $\sum{b_{ij}}=S_b$ ，还没有取到的钻石 $\sum{a_{ij}}=A$ ，还没有提高的价格 $\sum{b_{ij}}=B$ ；
答案表示为 $(S_a+A)\times(S_b+B)=Sa \times S_b +S_a \times B+A \times S_b+A\times B$ ；
如果我即将要走的点（向下或者向右）的估价值不能够更新答案的话，就剪枝掉；
为了写估价函数，令未来某条路径的权值之和为 $C_i=A_i+B_i$ ；
用动态规划写出C数组，使每个位置都对应一个 $C_{max}$ ；
$C_{max} = A + B$ ；
已经走过的路径 $S_a$ 和 $S_b$ 已经是最优解不能更改了，能改变的只有未来怎么走，即 $A$ 和 $B$ ;
那么我们就要找到一个合适的 $A$ 和 $B$ 使得函数 $(S_a+A)\times(S_b+B)$ 最大，就是整体最优解;
将 $B = C_{max}-A$ 代入到 $(S_a+A)\times(S_b+B)$ 中，得到二次函数 $F = -A^2 + (S_b - S_a + C_{max})A + S_aS_b + S_aC_{cmax}$ ，要使F最大，A就要取 $-\frac{b}{2a}$ ;
$\frac{S_a - S_b + C_{max}}{2}, B = \frac{S_b - S_a + C_{max}}{2}$ ；
将A和B回代到函数F中就是一条路径的答案，如果这个答案不能更新ans，则剪枝;
代入得到 $\frac{1}{4}S_a^2 + \frac{1}{4}S_b^2 + \frac{1}{4}C_{max}^2 + \frac{1}{2}S_aS_b + \frac{1}{2}S_aC_{max} + \frac{1}{2}S_bC_{max}$ ;
运用小学完全平方公式 $a+b+c)^2 = a^2 +b^2+c^2+2ab+2ac+2bc$ ;
得到 $(\frac{S_a+S_b+C_{max}}{2})^2$
下面深搜的过程看一下代码吧

Algorithm:

启发式搜索

Code:

#include <bits/stdc++.h>
#define int long long

using namespace std;	

const int N = 105;

int a[N][N], b[N][N], c[N][N];
int ans, n;

int get(int Sa, int Sb, int Cmax)
{
	return Sa + Sb + Cmax >> 1;
} 

void dfs(int x, int y, int Sa, int Sb)
{
	if (x == n && y == n)
	{
		ans = max(ans, Sa*Sb);
		return;
	}
	
	int t = get(Sa, Sb, c[x + 1][y]);
	if (x < n && t * t > ans) 
	{
		dfs(x + 1, y, Sa + a[x + 1][y], Sb + b[x + 1][y]);
	}

	t = get(Sa, Sb, c[x][y + 1]);
	if (y < n && t * t > ans) 
	{
		dfs(x, y + 1, Sa + a[x][y + 1], Sb + b[x][y + 1]);
	}
}

signed main()
{
	int T;
	cin >> T;
	while (T -- )
	{
		cin >> n;
		for (int i = 1; i <= n; i ++ )
			for (int j = 1; j <= n; j ++ ) cin >> a[i][j];
		for (int i = 1; i <= n; i ++ )
			for (int j = 1; j <= n; j ++ ) cin >> b[i][j];
		
		for (int i = n; i >= 1; i -- )
			for (int j = n; j >= 1; j -- )
			{
				c[i][j] = a[i][j] + b[i][j];
				if (i < n) c[i][j] = max(c[i][j], c[i + 1][j] + a[i][j] + b[i][j]);
				if (j < n) c[i][j] = max(c[i][j], c[i][j + 1] + a[i][j] + b[i][j]);
			}
			
		ans = 0;
		dfs(1, 1, a[1][1], b[1][1]);
		cout << ans << endl;
	}
	return 0;
}