Partition problem

In number theory and computer science, the partition problem, or number partitioning,^[1] is the task of deciding whether a given multiset S of positive integers can be partitioned into two subsets S₁ and S₂ such that the sum of the numbers in S₁ equals the sum of the numbers in S₂. Although the partition problem is NP-complete, there is a pseudo-polynomial time dynamic programming solution, and there are heuristics that solve the problem in many instances, either optimally or approximately. For this reason, it has been called "the easiest hard problem".^[2][3]

There is an optimization version of the partition problem, which is to partition the multiset S into two subsets S₁, S₂ such that the difference between the sum of elements in S₁ and the sum of elements in S₂ is minimized. The optimization version is NP-hard, but can be solved efficiently in practice.^[4]

The partition problem is a special case of two related problems:

• In the subset sum problem, the goal is to find a subset of S whose sum is a certain target number T given as input (the partition problem is the special case in which T is half the sum of S).
• In multiway number partitioning, there is an integer parameter k, and the goal is to decide whether S can be partitioned into k subsets of equal sum (the partition problem is the special case in which k = 2).

However, it is quite different to the 3-partition problem: in that problem, the number of subsets is not fixed in advance – it should be |S|/3, where each subset must have exactly 3 elements. 3-partition is much harder than partition – it has no pseudo-polynomial time algorithm unless P = NP.^[5]

Given S = {3,1,1,2,2,1}, a valid solution to the partition problem is the two sets S₁ = {1,1,1,2} and S₂ = {2,3}. Both sets sum to 5, and they partition S. This solution is not unique. S₁ = {3,1,1} and S₂ = {2,2,1} is another solution.

Not every multiset of positive integers has a partition into two subsets with equal sum. An example of such a set is S = {2,5}.

The partition problem is NP hard. This can be proved by reduction from the subset sum problem.^[6] An instance of SubsetSum consists of a set S of positive integers and a target sum T; the goal is to decide if there is a subset of S with sum exactly T.

Given such an instance, construct an instance of Partition in which the input set contains the original set plus two elements: z₁ and z₂, with z₁ = sum(S) and z₂ = 2T. The sum of this input set is sum(S) + z₁ + z₂ = 2 sum(S) + 2T, so the target sum for Partition is sum(S) + T.

• Suppose there exists a solution S′ to the SubsetSum instance. Then sum(S′) = T, so sum(S′ ∪ z₁) = sum(S) + T, so S′ ∪ z₁ is a solution to the Partition instance.
• Conversely, suppose there exists a solution S′′ to the Partition instance. Then, S′′ must contain either z₁ or z₂, but not both, since their sum is more than sum(S) + T. If S'' contains z₁, then it must contain elements from S with a sum of exactly T, so S'' minus z₁ is a solution to the SubsetSum instance. If S'' contains z₂, then it must contain elements from S with a sum of exactly sum(S) − T, so the other objects in S are a solution to the SubsetSum instance.

As mentioned above, the partition problem is a special case of multiway-partitioning and of subset-sum. Therefore, it can be solved by algorithms developed for each of these problems. Algorithms developed for multiway number partitioning include:

• Greedy number partitioning – loops over the numbers, and puts each number in the set whose current sum is smallest. If the numbers are not sorted, then the runtime is O(n) and the approximation ratio is at most 3/2 ("approximation ratio" means the larger sum in the algorithm output, divided by the larger sum in an optimal partition). Sorting the numbers increases the runtime to O(n log n) and improves the approximation ratio to 7/6. If the numbers are distributed uniformly in [0,1], then the approximation ratio is at most 1 + O ( log ⁡ log ⁡ n n ) {\textstyle 1+O\left({\frac {\log \log n}{n}}\right)} almost surely, and 1 + O ( 1 n ) {\textstyle 1+O\left({\frac {1}{n}}\right)} in expectation.
• Largest Differencing Method (also called the Karmarkar–Karp algorithm) sorts the numbers in descending order and repeatedly replaces numbers by their differences. The runtime complexity is O(n log n). In the worst case, its approximation ratio is similar – at most 7/6. However, in the average case it performs much better than the greedy algorithm: when numbers are distributed uniformly in [0,1], its approximation ratio is at most 1 + 1 / n Θ ( log ⁡ n ) {\textstyle 1+1/n^{\Theta (\log n)}} in expectation. It also performs better in simulation experiments.
• The multifit algorithm uses binary search combined with an algorithm for bin packing. In the worst case, its approximation ratio is 8/7.
• The subset sum problem has an FPTAS which can be used for the partition problem as well, by setting the target sum to sum(S)/2.

There are exact algorithms, that always find the optimal partition. Since the problem is NP-hard, such algorithms might take exponential time in general, but may be practically usable in certain cases. Algorithms developed for multiway number partitioning include:

• The pseudopolynomial time number partitioning takes O ( n 2 m ) {\textstyle O(n^{2}m)} time and needs O ( n m ) {\textstyle O(nm)} memory, where m is the largest number in the input.
• The Complete Greedy Algorithm (CGA) considers all partitions by constructing a binary tree. Each level in the tree corresponds to an input number, where the root corresponds to the largest number, the level below to the next-largest number, etc. Each branch corresponds to a different set in which the current number can be put. Traversing the tree in depth-first order requires only O ( n ) {\textstyle O(n)} space, but might take O ( 2 n ) {\textstyle O(2^{n})} time. The runtime can be improved by using a greedy heuristic: in each level, develop first the branch in which the current number is put in the set with the smallest sum. This algorithm finds first the solution found by greedy number partitioning, but then proceeds to look for better solutions. Some variations of this idea are fully polynomial-time approximation schemes for the subset-sum problem, and hence for the partition problem as well.^[7][8]
• The Complete Karmarkar–Karp algorithm (CKK) considers all partitions by constructing a binary tree. Each level corresponds to a pair of numbers. The left branch corresponds to putting them in different subsets (i.e., replacing them by their difference), and the right branch corresponds to putting them in the same subset (i.e., replacing them by their sum). This algorithm finds first the solution found by the largest differencing method, but then proceeds to find better solutions. It runs substantially faster than CGA on random instances. Its advantage is much larger when an equal partition exists, and can be of several orders of magnitude. In practice, problems of arbitrary size can be solved by CKK if the numbers have at most 12 significant digits.^[9] CKK can also run as an anytime algorithm: it finds the KK solution first, and then finds progressively better solutions as time allows (possibly requiring exponential time to reach optimality, for the worst instances).^[1] It requires O ( n ) O(n) space, but in the worst case might take O ( 2 n ) O(2^{n}) time.

Algorithms developed for subset sum include:

• Horowitz and Sanhi – runs in time O ( 2 n / 2 ⋅ ( n / 2 ) ) O(2^{n/2}\cdot (n/2)) , but requires O ( 2 n / 2 ) O(2^{n/2}) space.
• Schroeppel and Shamir – runs in time O ( 2 n / 2 ⋅ ( n / 4 ) ) O(2^{n/2}\cdot (n/4)) , and requires much less space – O ( 2 n / 4 ) {\textstyle O(2^{n/4})} .
• Howgrave-Graham and Joux – runs in time O ( 2 n / 3 ) {\textstyle O(2^{n/3})} , but it is a randomized algorithm that only solves the decision problem (not the optimization problem).

Sets with only one, or no partitions tend to be hardest (or most expensive) to solve compared to their input sizes. When the values are small compared to the size of the set, perfect partitions are more likely. The problem is known to undergo a "phase transition"; being likely for some sets and unlikely for others. If m is the number of bits needed to express any number in the set and n is the size of the set then m / n < 1 {\displaystyle m/n<1} tends to have many solutions and m / n > 1 {\displaystyle m/n>1} tends to have few or no solutions. As n and m get larger, the probability of a perfect partition goes to 1 or 0 respectively. This was originally argued based on empirical evidence by Gent and Walsh,^[10] then using methods from statistical physics by Mertens,^[11][12] and later proved by Borgs, Chayes, and Pittel.^[13]

A related problem, somewhat similar to the Birthday paradox, is that of determining the size of the input set so that we have a probability of one half that there is a solution, under the assumption that each element in the set is randomly selected with uniform distribution between 1 and some given value. The solution to this problem can be counter-intuitive, like the birthday paradox.

Equal-cardinality partition is a variant in which both parts should have an equal number of items, in addition to having an equal sum. This variant is NP-hard too.^[5]: SP12 Proof. Given a standard Partition instance with some n numbers, construct an Equal-Cardinality-Partition instance by adding n zeros. Clearly, the new instance has an equal-cardinality equal-sum partition iff the original instance has an equal-sum partition. See also Balanced number partitioning.

Product partition is the problem of partitioning a set of integers into two sets with the same product (rather than the same sum). This problem is strongly NP-hard.^[14]

Kovalyov and Pesch^[15] discuss a generic approach to proving NP-hardness of partition-type problems.

One application of the partition problem is for manipulation of elections. Suppose there are three candidates (A, B and C). A single candidate should be elected using a voting rule based on scoring, e.g. the veto rule (each voter vetoes a single candidate and the candidate with the fewest vetoes wins). If a coalition wants to ensure that C is elected, they should partition their votes among A and B so as to maximize the smallest number of vetoes each of them gets. If the votes are weighted, then the problem can be reduced to the partition problem, and thus it can be solved efficiently using CKK. The same is true for any other voting rule that is based on scoring.^[16]

• Borgs, Christian; Chayes, Jennifer; Pittel, Boris (2001), "Phase transition and finite-size scaling for the integer partitioning problem", Random Structures and Algorithms, 19 (3–4): 247–288, CiteSeerX 10.1.1.89.9577, doi:10.1002/rsa.10004, S2CID 6819493
• Gent, Ian; Walsh, Toby (August 1996). "Phase Transitions and Annealed Theories: Number Partitioning as a Case Study". In Wolfgang Wahlster (ed.). Proceedings of 12th European Conference on Artificial Intelligence. ECAI-96. John Wiley and Sons. pp. 170–174. CiteSeerX 10.1.1.2.4475.
• Gent, Ian; Walsh, Toby (1998), "Analysis of Heuristics for Number Partitioning", Computational Intelligence, 14 (3): 430–451, CiteSeerX 10.1.1.149.4980, doi:10.1111/0824-7935.00069, S2CID 15344203
• Korf, Richard E. (1998), "A complete anytime algorithm for number partitioning", Artificial Intelligence, 106 (2): 181–203, CiteSeerX 10.1.1.90.993, doi:10.1016/S0004-3702(98)00086-1, ISSN 0004-3702
• Mertens, Stephan (November 1998), "Phase Transition in the Number Partitioning Problem", Physical Review Letters, 81 (20): 4281–4284, arXiv:cond-mat/9807077, Bibcode:1998PhRvL..81.4281M, doi:10.1103/PhysRevLett.81.4281, S2CID 119541289
• Mertens, Stephan (2001), "A physicist's approach to number partitioning", Theoretical Computer Science, 265 (1–2): 79–108, arXiv:cond-mat/0009230, doi:10.1016/S0304-3975(01)00153-0, S2CID 16534837
• Mertens, Stephan (2006). "The Easiest Hard Problem: Number Partitioning". In Allon Percus; Gabriel Istrate; Cristopher Moore (eds.). Computational complexity and statistical physics. USA: Oxford University Press. pp. 125–140. arXiv:cond-mat/0310317. Bibcode:2003cond.mat.10317M. ISBN 9780195177374.
• Mertens, Stephan (1999), "A complete anytime algorithm for balanced number partitioning", arXiv:cs/9903011