A polynomial-time approximation algorithm for the number of k-matchings in bipartite graphs