2023-12-23 poj3040 Allowance 数学证明 贪心算法 代码实现

一、算法伪码

allowance_greedy_distribution(coins, c):
// coins记录了硬币数量和面值，c是每天至少发放的津贴
	将coins按照coins.value从大到小排序
	day=0
	while true:
	    today=0
	    for coin in coins:
		    尽可能选择面值大的硬币但不使其溢出
		if 还有硬币:
			使用还有的最小一种硬币，并使其数量-1
		if 没有硬币:
			break

二、引理

在硬币数目不受限制的情况下，

引理$P_n$：一组货币{$B_1$,$B_2$,…,$B_n$}满足整除性质，具有价值$$W_n=\sum _{i=1}^n{k}_i{B}_i$$，可以硬币数不增的将{$B_1$,$B_2$,…,$B_{n-1}$}尽数替换为$B_n$，并且这个过程结束后，满足$$W_{n-1}^{{}^{\prime }}<B_n$$详细证明请见2023-12-23 硬币选择问题 硬币找零问题 poj3040 贪心选择性质 数学证明-CSDN博客的引理部分。

三、贪心替换引理

硬币数量受限的贪心性质：至少存在一个从大到小选择硬币且满足硬币剩余数量都大于0的替换

引理：在硬币数量受限制时，存在一个满足贪心性质的替换。

证明：采用构造解法证明。
根据题意，至少有一组不一定满足贪心性质的解{${\lambda }_1,{\lambda }_2,...,{\lambda }_n$}使得$$W_n=\sum _{i=1}^n{\lambda }_i{B}_i=W$$
先假设硬币数量不受限制，由$P_n$，一定有一组满足贪心性质的替换{${\lambda }_1^{{}^{\prime }},{\lambda }_2^{{}^{\prime }},...,{\lambda }_n^{{}^{\prime }}$}，且满足$${\lambda }_i<={\lambda }_i^{{}^{\prime }}，i>=2$$以及$${\lambda }_i<=coins[i].numbers，i>=1$$
再按照以下方法构造一个满足贪心性质的替换：
从${\lambda }_n^{{}^{\prime }}到{\lambda }_2^{{}^{\prime }}$依次进行如下操作：
$$\{\begin{array}{cc}将超出数额的硬币换成若干个面值次大的硬币& {\lambda }_k>coin[k].numbers\\ continue& {\lambda }_k<=coin[k].numbers\end{array}$$
利用这个方法得到一组除了${\lambda }_1^{{}^{\prime \prime }}$外一定满足贪心性质且合法的替换{${\lambda }_1^{{}^{\prime \prime }},{\lambda }_2^{{}^{\prime \prime }},...,{\lambda }_n^{{}^{\prime \prime }}$}。

因为${\lambda }_k^{{}^{\prime \prime }}$(k>=2)不是coins[k].numbers就是${\lambda }_k^{{}^{\prime }}$，所以一定有$${\lambda }_i<={\lambda }_i^{{}^{\prime \prime }}，i>=2$$
下面利用反证法证明其是一组合法替换。
假设这组替换不合法，则有这组替换$$W_n^{{}^{\prime \prime }}=\sum _{i=2}^n{\lambda }_i^{{}^{\prime \prime }}{B}_i+{\lambda }_1^{{}^{\prime \prime }}{B}_1>=\sum _{i=2}^n{\lambda }_i{B}_i+{\lambda }_1^{{}^{\prime \prime }}{B}_1>\sum _{i=1}^n{\lambda }_i{B}_i=W$$，这与构造过程的W不变相矛盾，故这个替换一定是一组合法替换。

四、贪心选择性质

该Allowence_distribution问题的贪心选择性质为：至少包含一个每次尽可能少给牛但不低于c，且每次选尽可能大但是工资不溢出的构造得到的最优解，

假设$D_i$为第i天的一个分配，其价值为$D_i.value$，其分配为$D_i.\lambda []$。

第一步：

假设有一个解是最优解，但是不满足上述贪心选择性质，由于$D[]$与其天数顺序无关，可以进行排序得到它的一个同构解，并且按照这个顺序进行构造时，每一步的多余工资$$D_i.cost=D_i.value-c>=D_{贪心i}.cost$$
倘若解$D_{贪心}[]$能够发n天工资，则有总价值$$W=nc+\sum _{i=1}^n{D}_{贪心i}.cost+D_{贪心}.last=nc+\sum _{i=1}^n{D}_i.cost+D.last$$
则对于每个n都有：
$$D_{贪心}.last>=D.last$$
故只要$D.last$能够一天工资，$D_{贪心}.last$就至少够发一天工资。

第二步：

则假设有一个最优解$D_i$，由贪心替换引理，可以得到一个最优解的同构$D_i^{{}^{\prime }}$。
证毕。

五、代码实现

可以直接运行的示例：

#include <iostream>
#include <vector>
#include <queue>
#include <list>
#include <math.h>
//#include <assert.h>
using namespace std;
//using Datatype = int;	
typedef long long Datatype;
typedef long long Elemtype;

struct COIN {
	Datatype value;
	Elemtype numbers;
};

Elemtype allowence_distribution_greedy(list<COIN>& coins, Datatype c); // 接受硬币的数据和每天的津贴数，返回够发的天数

int main() {
	list<COIN> coins = { {10,1},{1,100},{5,120} };
	cout << allowence_distribution_greedy(coins, 6);
}


Elemtype allowence_distribution_greedy(list<COIN>& coins, Datatype c) {
	coins.sort([](COIN a, COIN b) { return a.value > b.value; });
	Elemtype days = 0;
	while (1) {
		Datatype today = 0;
		auto it = coins.begin();
		while (it != coins.end())
		{
			if (c == today)break;
			if (c < today + it->value)it++;
			else {
				today += it->value;
				it->numbers--;
				if (0 == it->numbers)coins.erase(it++);
			}
		}
		if (c == today)days++;
		else if (!coins.empty()) {
			// 只要还有硬币，必然再加一块硬币就超过c
			days++;
			auto it = --coins.end();
			it->numbers--;
			//assert(today + it->value > c);
			if (0 == it->numbers)coins.erase(it);
		}
		if (coins.empty())break;
	}
	return days;
}

六、补充
补充部分和本文无任何干系，仅供看个乐子。

1. 本题使用了两次贪心选择性质。
2. 贪心选择性一般而言，都出现在顺序无关的情况下。