P1650 田忌赛马,贪心,线性dp

P1650 田忌赛马 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

题目描述

我国历史上有个著名的故事: 那是在 2300 年以前。齐国的大将军田忌喜欢赛马。他经常和齐王赛马。他和齐王都有三匹马:常规马,上级马,超级马。一共赛三局,每局的胜者可以从负者这里取得 200 银币。每匹马只能用一次。齐王的马好,同等级的马,齐王的总是比田忌的要好一点。于是每次和齐王赛马,田忌总会输 600 银币。

田忌很沮丧,直到他遇到了著名的军师――孙膑。田忌采用了孙膑的计策之后,三场比赛下来,轻松而优雅地赢了齐王 200 银币。这实在是个很简单的计策。由于齐王总是先出最好的马,再出次好的,所以田忌用常规马对齐王的超级马,用自己的超级马对齐王的上级马,用自己的上级马对齐王的常规马,以两胜一负的战绩赢得 200 银币。实在很简单。

如果不止三匹马怎么办?这个问题很显然可以转化成一个二分图最佳匹配的问题。把田忌的马放左边,把齐王的马放右边。田忌的马 A 和齐王的 B 之间,如果田忌的马胜,则连一条权为 200 的边;如果平局,则连一条权为 0 的边;如果输,则连一条权为 −200 的边……如果你不会求最佳匹配,用最小费用最大流也可以啊。 然而,赛马问题是一种特殊的二分图最佳匹配的问题,上面的算法过于先进了,简直是杀鸡用牛刀。现在,就请你设计一个简单的算法解决这个问题。

输入格式

第一行一个整数 n ,表示他们各有几匹马(两人拥有的马的数目相同)。第二行 n 个整数,每个整数都代表田忌的某匹马的速度值(0≤ 速度值 ≤100)。第三行 n 个整数,描述齐王的马的速度值。两马相遇,根据速度值的大小就可以知道哪匹马会胜出。如果速度值相同,则和局,谁也不拿钱。

输出格式

仅一行,一个整数,表示田忌最大能得到多少银币。

输入输出样例

输入 #1复制

3
92 83 71
95 87 74

输出 #1复制

200

说明/提示

数据规模与约定

  • 对于 20%的数据,1≤N≤65;
  • 对于 40% 的数据,1≤N≤250;
  • 对于 100%的数据,11≤N≤2000。

 解析:

贪心:
性质:先从大到小排序,
 若 a 中剩下的最快的马比 b 中剩下最快的马快,则两只马去比;
 若 a 中剩下的最快的马比 b 中剩下的最快的马慢,则用 a 中剩下的最慢的马去输;
 若 a 中剩下的最快的马与 b 中剩下最快的马相等,则:
 1.a 中剩下的最慢的马与 b 中剩下最慢的马快,则两只最慢的马比
 2.a 中剩下的最慢的马不比 b 中剩下最慢的马快,且比 b 中最快的马慢,则用 a 中最慢的马输给 b 中最快的马,否则说明剩下的马都能两两打平

 那么最后一种情况存不存在两只最慢的马相等呢?若存在为什么不让这两组相等的马打平手呢?
 其实如果存在这种情况,那么 a 中倒数第二慢的马一定能赢 b 中最慢的马(因为性质的最后一句),且 a 中最快的马一定能赢 b 中第二快的马

 dp:
 根据上述贪心的性质我们已经有了推理的方向,因此我们可以据此将集合划分为两个不重不漏的子集
 f[i][j] :表示第 i 从比赛时从 a 的前面选了 j 匹马所赢的最大金币数量
 状态转移方程为:f[i][j]=max(f[i-1][j]+g[n-i+j+1][i],f[i-1][j-1]+g[j][i])
 其中 g[j][i] 表示 a 中的第 j 匹马与 b 中的第 i 匹马比是能赢200,-200,0

贪心代码:

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;
typedef long long LL;
const int N = 2e3 + 3, INF = 2e8;
int n;
int a[N], b[N];

int cmp(const int& a, const int& b) {
	return a > b;
}

int main() {
	cin >> n;
	for (int i = 1; i <= n; i++) {
		scanf("%d", &a[i]);
	}
	for (int j = 1; j <= n; j++) {
		scanf("%d", &b[j]);
	}
	sort(a + 1, a + 1 + n, cmp);
	sort(b + 1, b + 1 + n, cmp);
	int ans = 0,la,lb,ra,rb;
	la = lb = 1;
	ra = rb = n;
	for (int i = 1; i <= n; i++) {
		if (a[la] > b[lb]) {
			ans+=200;
			la++, lb++;
		}
		else if (a[la] < b[lb]) {
			ans -= 200;
			ra--;
			lb++;
		}
		else if (a[ra] > b[rb]) {
			ans += 200;
			ra--;
			rb--;
		}
		else {
			if (a[ra] < b[lb])ans -= 200;
			ra--;
			lb++;
		}
	}
	cout << ans << endl;
	return 0;
}

dp:

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;
typedef long long LL;
const int N = 2e3 + 3, INF = 2e8;
int n;
int a[N], b[N], g[N][N], f[N][N];

int cmp(const int& a, const int& b) {
	return a > b;
}

int main() {
	cin >> n;
	for (int i = 1; i <= n; i++) {
		scanf("%d", &a[i]);
	}
	for (int j = 1; j <= n; j++) {
		scanf("%d", &b[j]);
	}
	sort(a + 1, a + 1 + n, cmp);
	sort(b + 1, b + 1 + n, cmp);
	for (int i = 1; i <= n; i++) {
		for (int j = 1; j <= n; j++) {
			if (a[i] > b[j])g[i][j] = 200;
			else if (a[i] < b[j])g[i][j] = -200;
			else g[i][j] = 0;
			f[i][j] = -INF;
		}
	}
	for (int i = 1; i <= n; i++) {

		f[i][i] = f[i - 1][i - 1] + g[i][i];
		f[i][0] = f[i - 1][0] + g[n - i + 1][i];
		for (int j = 1; j < i; j++) {
			f[i][j] = max(f[i - 1][j] + g[n - i + j + 1][i], f[i - 1][j - 1] + g[j][i]);
		}
	}
	int ans = f[n][1];
	for ( int i = 1; i <= n; i++) {
		ans = max(ans, f[n][i]);
	}
	cout << ans << endl;
	return 0;
}

你可能感兴趣的:(#,线性dp,贪心,动态规划,算法)