我知道现在回答已经晚了,但我想我已经找到了解决这个问题的简单方法。我们使用回溯按字典顺序枚举S 的子集,并检查到目前为止生成的子集sum。
当sum 超过k 时,有趣的部分来了:我们需要检查生成的subset 是否是先前报告的项目的适当子集。
一种解决方案是保留所有报告的子集并检查是否包含,但这很浪费。
相反,我们计算k 和sum 之间的差异。如果S 中有一个元素e 使得e not in subset 和e <= (k - sum),那么我们生成的集合是先前报告的子集的适当子集,我们可以安全地跳过它。
下面是用普通的旧 C++ 编写的完整工作程序,展示了这个想法:
#include <iostream>
#include <vector>
#include <set>
#include <algorithm>
typedef std::set<int> Set;
typedef std::vector<int> SubSet;
bool seen_before(const Set &universe, const SubSet &subset, int diff) {
Set::const_iterator i = std::mismatch(universe.begin(), universe.end(),
subset.begin()).first;
return i != universe.end() && *i <= diff;
}
void process(const SubSet &subset) {
if (subset.empty()) {
std::cout << "{}\n";
return;
}
std::cout << "{" << subset.front();
for (SubSet::const_iterator i = subset.begin() + 1, e = subset.end();
i != e; ++i) {
std::cout << ", " << *i;
}
std::cout << "}\n";
}
void generate_max_subsets_rec(const Set &universe, SubSet &subset,
long sum, long k) {
Set::const_iterator i = subset.empty()
? universe.begin()
: universe.upper_bound(subset.back()),
e = universe.end();
if (i == e) {
if (!seen_before(universe, subset, k - sum))
process(subset);
return;
}
for (; i != e; ++i) {
long new_sum = sum + *i;
if (new_sum > k) {
if (!seen_before(universe, subset, int(k - sum)))
process(subset);
return;
} else {
subset.push_back(*i);
if (new_sum == k)
process(subset);
else
generate_max_subsets_rec(universe, subset, new_sum, k);
subset.pop_back();
}
}
}
void generate_max_subsets(const Set &universe, long k) {
SubSet subset;
subset.reserve(universe.size());
generate_max_subsets_rec(universe, subset, 0, k);
}
int main() {
int items[] = {1, 2, 3, 4, 5};
Set u(items, items + (sizeof items / sizeof items[0]));
generate_max_subsets(u, 7);
return 0;
}
输出是按字典顺序排列的所有最大子集,每行一个:
{1, 2, 3}
{1, 2, 4}
{1, 5}
{2, 5}
{3, 4}