Matroid intersection

In combinatorial optimization, the matroid intersection problem is to find a largest common independent set in two matroids over the same ground set. If the elements of the matroid are assigned real weights, the weighted matroid intersection problem is to find a common independent set with the maximum possible weight. These problems generalize many problems in combinatorial optimization including finding maximum matchings and maximum weight matchings in bipartite graphs and finding arborescences in directed graphs.

The matroid intersection theorem, due to Jack Edmonds,^[1] says that there is always a simple upper bound certificate, consisting of a partitioning of the ground set amongst the two matroids, whose value (sum of respective ranks) equals the size of a maximum common independent set. Based on this theorem, the matroid intersection problem for two matroids can be solved in polynomial time using matroid partitioning algorithms.

Examples

Let G = (U,V,E) be a bipartite graph. One may define a partition matroid M_U on the ground set E, in which a set of edges is independent if no two of the edges have the same endpoint in U. Similarly one may define a matroid M_V in which a set of edges is independent if no two of the edges have the same endpoint in V. Any set of edges that is independent in both M_U and M_V has the property that no two of its edges share an endpoint; that is, it is a matching. Thus, the largest common independent set of M_U and M_V is a maximum matching in G.

Similarly, if each edge has a weight, then the maximum-weight independent set of M_U and M_V is a Maximum weight matching in G.

Algorithms

There are several polynomial-time algorithms for weighted matroid intersection, with different run-times. The run-times are given in terms of $n$ - the number of elements in the common base-set, $r$ - the maximum between the ranks of the two matroids, $T$ - the number of operations required for a circuit-finding oracle, and $k$ - the number of elements in the intersection (in case we want to find an intersection of a specific size $k$ ).

Edmonds' algorithm uses linear programming and polyhedra.^[1]
Lawler's algorithm.^[2]
Iri and Tomizawa's algorithm^[3]
Andras Frank's algorithm^[4] uses $O(n^{3}T)$ arithmetic operations.
Orlin and Vande-Vate's algorithm.^[5]
Cunningham's algorithm^[6] requires $O(r^{1.5}nT)$ operations on general matroids, and $O(nr^{2}\log {r})$ operations on linear matroid, for two r-by-n matrices.
Brezovec, Cornuejos and Glover^[7] present two algorithms for weighted matroid intersection.
- The first algorithm requires that all weights be integers, and finds an intersection of cardinality $k$ in time $O((n^{3}k^{2}+nkT)\cdot (1+\log {W}))$ .
- The second algorithm runs in time $O(nr\cdot (\log {n}+r+T))$ .
Huang, Kakimura and Kamiyama^[8] show that the weighted matroid intersection problem can be solved by solving W instances of the unweighted matroid intersection problem, where W is the largest given weight, assuming that all given weights are integral. This algorithm is faster than previous algorithms when W is small. They also present an approximation algorithm that finds an e-approximate solution by solving $O(\epsilon ^{-1}\log {r})$ instances of the unweighted matroid intersection problem, where r is the smaller rank of the two input matroids.
Ghosh, Gurjar and Raj^[9] study the run-time complexity of matroid intersection in the parallel computing model.
Bérczi, Király, Yamaguchi and Yokoi^[10] present strongly polynomial-time algorithms for weighted matroid intersection using more restricted oracles.

Extensions

Maximizing weight subject to cardinality

In a variant of weighted matroid intersection, called "(P_k)", the goal is to find a common independent set with the maximum possible weight among all such sets with cardinality k, if such a set exists. This variant, too, can be solved in polynomial time.^[7]

Three matroids

The matroid intersection problem becomes NP-hard when three matroids are involved, instead of only two.

One proof of this hardness result uses a reduction from the Hamiltonian path problem in directed graphs. Given a directed graph G with n vertices, and specified nodes s and t, the Hamiltonian path problem is the problem of determining whether there exists a simple path of length n − 1 that starts at s and ends at t. It may be assumed without loss of generality that s has no incoming edges and t has no outgoing edges. Then, a Hamiltonian path exists if and only if there is a set of n − 1 elements in the intersection of three matroids on the edge set of the graph: two partition matroids ensuring that the in-degree and out-degree of the selected edge set are both at most one, and the graphic matroid of the undirected graph formed by forgetting the edge orientations in G, ensuring that the selected edge set has no cycles.^[11]

Matroid parity

Another computational problem on matroids, the matroid parity problem, was formulated by Lawler^[12] as a common generalization of matroid intersection and non-bipartite graph matching. However, although it can be solved in polynomial time for linear matroids, it is NP-hard for other matroids, and requires exponential time in the matroid oracle model.^[13]

Valuated matroids

A valuated matroid is a matroid equipped with a value function v on the set of its bases, with the following exchange property: for any two distinct bases $A$ and $B$ , if $a\in A\setminus B$ , then there exists an element $b\in B\setminus A$ such that both $(A\setminus \{a\})\cup \{b\}$ and $(B\setminus \{b\})\cup \{a\}$ are bases, and : $v(A)+v(B)\leq v(B\setminus \{b\}\cup \{a\})+v(A\setminus \{a\}\cup \{b\})$ .

Given a weighted bipartite graph G = (X+Y, E) and two valuated matroids, one on X with bases set B_X and valuation v_X, and one on Y with bases B_Y and valuation v_Y, the valuated independent assignment problem is the problem of finding a matching M in G, such that M_X (the subset of X matched by M) is a base in B_X , M_Y is a base in B_Y , and subject to this, the sum $w(M)+v_{X}(M_{X})+v_{Y}(M_{Y})$ is maximized. The weighted matroid intersection problem is a special case in which the matroid valuations are constant, so we only seek to maximize $w(M)$ subject to M_X is a base in B_X and M_Y is a base in B_Y .^[14] Murota presents a polynomial-time algorithm for this problem.^[15]

References

^ ^a ^b Edmonds, Jack (1970), "Submodular functions, matroids, and certain polyhedra", in R. Guy; H. Hanam; N. Sauer; J. Schonheim (eds.), Combinatorial structures and their applications (Proc. 1969 Calgary Conference), Gordon and Breach, New York, pp. 69–87. Reprinted in M. Jünger et al. (Eds.): Combinatorial Optimization (Edmonds Festschrift), LNCS 2570, pp. 1126, Springer-Verlag, 2003.
^ Lawler, Eugene L. (1975), "Matroid intersection algorithms", Mathematical Programming, 9 (1): 31–56, doi:10.1007/BF01681329, S2CID 206801650
^ Iri, Masao; Tomizawa, Nobuaki (1976). "An Algorithm for Finding an Optimal "Independent Assignment"". 日本オペレーションズ・リサーチ学会論文誌. 19 (1): 32–57. doi:10.15807/jorsj.19.32.
^ Frank, András (1981), "A weighted matroid intersection algorithm", Journal of Algorithms, 2 (4): 328–336, doi:10.1016/0196-6774(81)90032-8
^ Orlin, James B.; VandeVate, John (1983). On a "primal" matroid intersection algorithm (Sloan School of Management Working Paper #1446-83). hdl:1721.1/2050.
^ Cunningham, William H. (1986-11-01). "Improved Bounds for Matroid Partition and Intersection Algorithms". SIAM Journal on Computing. 15 (4): 948–957. doi:10.1137/0215066. ISSN 0097-5397.
^ ^a ^b Brezovec, Carl; Cornuéjols, Gérard; Glover, Fred (1986), "Two algorithms for weighted matroid intersection", Mathematical Programming, 36 (1): 39–53, doi:10.1007/BF02591988, S2CID 34567631
^ Huang, Chien-Chung; Kakimura, Naonori; Kamiyama, Naoyuki (2019-09-01). "Exact and approximation algorithms for weighted matroid intersection". Mathematical Programming. 177 (1): 85–112. doi:10.1007/s10107-018-1260-x. hdl:2324/1474903. ISSN 1436-4646. S2CID 254138118.
^ Ghosh, Sumanta; Gurjar, Rohit; Raj, Roshan (2022-01-01), "A Deterministic Parallel Reduction from Weighted Matroid Intersection Search to Decision", Proceedings of the 2022 Annual ACM-SIAM Symposium on Discrete Algorithms (SODA), Proceedings, Society for Industrial and Applied Mathematics, pp. 1013–1035, doi:10.1137/1.9781611977073.44, ISBN 978-1-61197-707-3, S2CID 245799113, retrieved 2022-11-28
^ Bérczi, Kristóf; Király, Tamás; Yamaguchi, Yutaro; Yokoi, Yu (2022-09-28). "Matroid Intersection under Restricted Oracles". arXiv:2209.14516 [cs.DS].
^ Welsh, D. J. A. (2010) [1976], Matroid Theory, Courier Dover Publications, p. 131, ISBN 9780486474397.
^ Lawler, Eugene L. (1976), "Chapter 9: The Matroid Parity Problem", Combinatorial Optimization: Networks and Matroids, New York: Holt, Rinehart and Winston, pp. 356–367, MR 0439106.
^ Jensen, Per M.; Korte, Bernhard (1982), "Complexity of matroid property algorithms", SIAM Journal on Computing, 11 (1): 184–190, doi:10.1137/0211014, MR 0646772.
^ Murota, Kazuo (1996-11-01). "Valuated Matroid Intersection I: Optimality Criteria". SIAM Journal on Discrete Mathematics. 9 (4): 545–561. doi:10.1137/S0895480195279994. ISSN 0895-4801.
^ Murota, Kazuo (November 1996). "Valuated Matroid Intersection II: Algorithms". SIAM Journal on Discrete Mathematics. 9 (4): 562–576. doi:10.1137/S0895480195280009. ISSN 0895-4801.