HDU 1203 I NEED A OFFER!(DP, 01背包问题滚动数组优化)

Problem Description

Speakless很早就想出国，现在他已经考完了所有需要的考试，准备了所有要准备的材料，于是，便需要去申请学校了。要申请国外的任何大学，你都要交纳一定的申请费用，这可是很惊人的。Speakless没有多少钱，总共只攒了n万美元。他将在m个学校中选择若干的（当然要在他的经济承受范围内）。每个学校都有不同的申请费用a（万美元），并且Speakless估计了他得到这个学校offer的可能性b。不同学校之间是否得到offer不会互相影响。“I NEED A OFFER”，他大叫一声。帮帮这个可怜的人吧，帮助他计算一下，他可以收到至少一份offer的最大概率。（如果Speakless选择了多个学校，得到任意一个学校的offer都可以）。

Input

输入有若干组数据，每组数据的第一行有两个正整数n,m(0<=n<=10000,0<=m<=10000)
后面的m行，每行都有两个数据ai(整型),bi(实型)分别表示第i个学校的申请费用和可能拿到offer的概率。
输入的最后有两个0。

Output

每组数据都对应一个输出，表示Speakless可能得到至少一份offer的最大概率。用百分数表示，精确到小数点后一位。

Sample Input

10 3
4 0.1
4 0.2
5 0.3
0 0

Sample Output

44.0%

Hint

You should use printf("%%") to print a ‘%’.

这题题意一开始没读明白,最大概率都不知道怎么算出来的。后来想了想，至少获得一份的概率就是1 - (一份都不获得的概率)。这是一个关键思路！比如题中输出44.0%的那组数据是怎么得出的呢？就是只给第三个和第二个学校投申请，这样，两个都不获得通知书的概率为0.7 * 0.8 = 0.56。所以至少获得一个通知书的概率就是1 - 0.56 = 0.44！

那么实际上这一题就是一个01背包问题了，对于每一个学校，你都可以选择投申请书和不投申请书。01背包状态转移方程就很熟悉了，dp[j] = Min(dp[j], dp[j - a[i].cost] * a[i].chance); 这一题原本的max函数要换成min函数，dp[j]表示资金为j的情况下至少一个都不获得的最小概率，所以自然要用min函数。

01背包应该都是相当熟悉了，所以用二维数组可以轻松写出：

/*
01背包的变种，二维数组版 
*/
#include<iostream>
using namespace std;

const int maxn = 100;
struct Offer{
	int cost;
	double chance;				
};
Offer a[maxn];
int n, m;		//m表示行 
double dp[maxn][maxn];			//dp[i][j]表示不能得到一份offer的最小概率 

double Min(double x, double y)
{
	return x < y ? x : y;
}

double DP()
{
	for(int i = 1;i <= m;i++)
	{
		for(int j = 1;j <= n;j++)
		{
			if(i == 1)			//第一行是边界 
			{
				if(a[1].cost <= j)
					dp[i][j] = a[i].chance;
				else
					dp[i][j] = 1;
			}
			else
			{
				if(j - a[i].cost == 0)
				{
					if(a[i].chance < dp[i - 1][j])
						dp[i][j] = a[i].chance;
					else
						dp[i][j] = dp[i - 1][j];
				}
				else if(j - a[i].cost > 0)
					dp[i][j] = Min(dp[i - 1][j], dp[i - 1][j - a[i].cost] * a[i].chance);
				else
					dp[i][j] = dp[i - 1][j];
			}
		}
	}
	/*
	打印dp数组 
	for(int i = 1;i <= m;i++)
	{
		for(int j = 1;j <= n;j++)
		{
			cout << dp[i][j] << " ";
		}
		cout << endl;
	}
	*/
	return dp[m][n];
}

int main()
{
	while(scanf("%d%d", &n, &m) && m + n)
	{
		for(int i = 1;i <= m;i++)
		{
			double x;
			cin >> a[i].cost >> x;
			a[i].chance = 1.0 - x;
		}
		cout << DP() << endl;
	}
	return 0;
}

大概是这样的一个写法思路。但是，本题数据量的限制，我们不能使用二维数组来存放dp值，那么这个时候就涉及到01背包的优化写法，用一维的滚动数组优化。

说通俗点，因为我们在手写dp数表的时候就会发现，二维数组其实是多余的，因为我每次计算dp[i][j]的值时只需要用i - 1行的dp值就可以了，而i - 1行以上的值其实都没有必要存在数组中，因此诞生了这种优化空间的方法！但要注意，优化的话，遍历顺序就要有点改变，**列的遍历顺序就是从右往左了！**因为如果你从左往右遍历，你就会把后面需要用的dp值改掉，这就造成了错误！

贴上一维数组优化AC代码：

/*
01背包的变种，一维数组优化版 
*/
#include<iostream>
using namespace std;

const int maxn = 10005;
struct Offer{
	int cost;
	double chance;				
};
Offer a[maxn];
int n, m;		//m表示行 
double dp[maxn];	

double Min(double x, double y)
{
	return x < y ? x : y;
}

double DP()
{
	//重复利用dp数组，所以要从后面开始遍历了，因为利用的dp数组是前一行
	for(int i = 1;i <= m;i++)
	{
		for(int j = n;j >= a[i].cost;j--)
		{
			if(j - a[i].cost == 0)
				dp[j] = Min(dp[j], a[i].chance);
			else if(j - a[i].cost > 0)
				dp[j] = Min(dp[j], dp[j - a[i].cost] * a[i].chance);
		}	
	} 
	return dp[n];
}

void Clear()
{
	for(int i = 0;i <= n;i++)
		dp[i] = 1;
}

int main()
{
	while(scanf("%d%d", &n, &m) && m + n)
	{
		Clear();
		for(int i = 1;i <= m;i++)
		{
			double x;
			cin >> a[i].cost >> x;
			a[i].chance = 1.0 - x;
		}
		double res = 1.0 - DP();
		printf("%.1lf%%\n", res * 100);
	}
	return 0;
}

顺便说一句，hdu上面的评测是真的严格，每次思路感觉没问题，提交WA让人绝望。。。查错又花我不少时间。。。没办法