以下是PuLP/GLPK的解决方案。我以前从来没有使用过PuLP,但是它在PyPI上,似乎能完成这项工作。 GLPK非常好,免费。
from collections import defaultdict, namedtuple
from pulp import *
User = namedtuple('User', ('coverage', 'price'))
def solvesetcover(users):
vars = [LpVariable('x{}'.format(i), 0, 1, cat='Binary') for i, user in enumerate(users)]
prob = LpProblem()
totals = defaultdict(int)
for user, var in zip(users, vars):
prob += user.price * var
for elt in user.coverage:
totals[elt] += var
for total in totals.values():
prob += total >= 1
GLPK(msg=0).solve(prob)
return [user for user, var in zip(users, vars) if value(var)]
if __name__ == '__main__':
users = []
users.append(User({1, 2, 3, 4, 5, 6, 7}, 1.16))
users.append(User({8, 9, 10, 11, 12, 13, 14}, 1.08))
users.append(User({1, 8}, 1.04))
users.append(User({2, 3, 9, 10}, 1.02))
users.append(User({4, 5, 6, 7, 11, 12, 13, 14}, 1.01))
print(solvesetcover(users))
看起来像一个难题。我认为暴力是不可能的? – 2014-08-31 15:32:33
http://en.wikipedia.org/wiki/Set_cover_problem – 2014-08-31 16:03:46
它确实设置了封面,因此可能没有多项式时间解决方案。您可以使用指数动态程序或将该问题表述为[ILP](http://en.wikipedia.org/wiki/Integer_programming),并使用通用解算器。 – 2014-08-31 21:18:17