[动态规划][usaco3.1.6]stamps

Stamps 邮票

译 by Felicia Crazy

---

[编辑] 描述

已知一个 N 枚邮票的面值集合（如，{1 分，3 分}）和一个上限 K —— 表示信封上能够贴 K 张邮票。计算从 1 到 M 的最大连续可贴出的邮资。

例如，假设有 1 分和 3 分的邮票；你最多可以贴 5 张邮票。很容易贴出 1 到 5 分的邮资（用 1 分邮票贴就行了），接下来的邮资也不难：

6 = 3 + 3 
7 = 3 + 3 + 1 
8 = 3 + 3 + 1 + 1 
9 = 3 + 3 + 3 
10 = 3 + 3 + 3 + 1 
11 = 3 + 3 + 3 + 1 + 1 
12 = 3 + 3 + 3 + 3 
13 = 3 + 3 + 3 + 3 + 1

然而，使用 5 枚 1 分或者 3 分的邮票根本不可能贴出 14 分的邮资。因此，对于这两种邮票的集合和上限 K=5，答案是 M=13。

[规模最大的一个点的时限是3s]

小提示:因为14贴不出来,所以最高上限是13而不是15

[编辑] 格式

PROGRAM NAME: stamps

INPUT FORMAT:

(file stamps.in)

第 1 行： 两个整数，K 和 N。K（1 <= K <= 200）是可用的邮票总数。N（1 <= N <= 50）是邮票面值的数量。

第 2 行 .. 文件末： N 个整数，每行 15 个，列出所有的 N 个邮票的面值，每张邮票的面值不超过 10000。

OUTPUT FORMAT:

(file stamps.out)

第 1 行：一个整数，从 1 分开始连续的可用集合中不多于 K 张邮票贴出的邮资数。

[编辑] SAMPLE INPUT

5 2
1 3

[编辑] SAMPLE OUTPUT

13

比较入门的动态规划题目，但是看了题解后才做出来。

原始的方程可以是：

f ( i , j , k ) = Or f ( i-1 , j2 , k-a[i]*(j-j2) )

f ( i , j , k ) 表示前i种邮票贴j张，能不能表示出k。

Or 表示对于所有的j2，将 f( i-1 , j2 , k-a[i]*(j-j2) )求或。

容易想到：第一维i，可以省略以节省空间，但是不能少这一层循环。

接下来我们可能想到用求最优解来代替判断可行解，以减少枚举量。

即有两种方案：

1、前i种邮票用掉j张，能表示出的从1开始的连续的最大邮资。

2、前i种邮票表示出的邮资为j，最少需要的邮票数。

第2种应该算是反向的思维，然而很明显第2种描述更简单也更易实现。

我们最大的瓶颈在于，k的取值不确定，最大可能是10000*200=2000000。这里不处理就不可能实现。

但是我们注意到：我们需要从1到m的一段连续区间。因此我们猜想一旦不符合条件，就可以立即退出。

如果按照原来的方案，因为和与每一种邮票都相关，所以无法实现。

正确的方法通过调整枚举顺序而完美地解决了这个问题：

外层枚举邮资j，内层枚举邮票种类i。这样便消除了后效性：用1~j-1推出了j，每一种邮资都考虑到了所有邮票种类，因此与j+1及之后的状态都无关。

/*
ID: wuyihao1
LANG: C++
TASK: stamps
*/
#include <cstdio>
#include <cstring>

int f[3000000];
int a[80];

int min(int a,int b)
{
	if (a < b)
		return a;
	return b;
}

int main()
{
	freopen("stamps.in","r",stdin);
	freopen("stamps.out","w",stdout);

	int k,n;
	
	scanf("%d%d",&k,&n);
	memset(f,0x3f,sizeof f);
	f[0] = 0;
	for (int i=1;i<n+1;i++)
	{
		scanf("%d",a+i);
	}
	
	int i = 0;
	while (1)
	{
		i ++;
		for (int j=1;j<n+1;j++)
		{
			for (int k=0;i-a[j]*k>=0;k++)
			{
				f[i] = min(f[i],f[i-a[j]*k] + k);
			}
		}
		if (f[i] > k)
			break;
	}
	printf("%d\n",i-1);

	return 0;
}