Computational Complexity of Normalizing Constants for the Product of Determinantal Point Processes
Tatsuya Matsuoka,Naoto Ohsaka
DOI: https://doi.org/10.1016/j.tcs.2024.114513
IF: 1.002
2024-03-26
Theoretical Computer Science
Abstract:We consider the product of determinantal point processes (DPPs), a point process whose probability mass is proportional to the product of principal minors of multiple matrices, as a natural, promising generalization of DPPs. We study the computational complexity of computing its normalizing constant , which is among the most essential probabilistic inference tasks. Our complexity-theoretic results (almost) rule out the existence of efficient algorithms for this task unless the input matrices are forced to have favorable structures. In particular, we prove the following: • Computing ∑Sdet(AS,S)p exactly for every (fixed) positive even integer p is UP -hard and Mod 3 P -hard, which gives a negative answer to an open question posed by Kulesza and Taskar [51] . • ∑Sdet(AS,S)det(BS,S)det(CS,S) is NP -hard to approximate within a factor of 2O(|I|1−ε) or 2O(n1/ε) for any ε>0 , where |I| is the input size and n is the order of the input matrix. This result is stronger than the # P -hardness for the case of two matrices derived by Gillenwater [36] . • There exists a kO(k)nO(1) -time algorithm for computing ∑Sdet(AS,S)det(BS,S) , where k is the maximum rank of A and B or the treewidth of the graph formed by nonzero entries of A and B . Such parameterized algorithms are said to be fixed-parameter tractable . These results can be extended to the fixed-size case. Further, we present two applications of fixed-parameter tractable algorithms given a matrix A of treewidth w : • We can compute a 2n2p−1 -approximation to ∑Sdet(AS,S)p for any fractional number p>1 in wO(wp)nO(1) time. • We can find a 2n -approximation to unconstrained maximum a posteriori inference in wO(wn)nO(1) time.
computer science, theory & methods