Abstract

We study an important case of integer linear programs (ILPs) of the form $\max\{c^Tx \ \vert\ \mathcal Ax = b, l \leq x \leq u,\, x \in \mathbb{Z}^{n t} \} $ with $n t$ variables and lower and upper bounds $\ell, u\in\mathbb Z^{nt}$. In $n$-fold ILPs nonzero entries only appear in the first $r$ rows of the matrix $\mathcal A$ and in small blocks of size $s\times t$ along the diagonal underneath. Despite this restriction, many optimization problems can be expressed in this form. It is known that $n$-fold ILPs are fixed-parameter tractable (FPT) regarding the parameters $s, r,$ and $\Delta$, where $\Delta$ is the greatest absolute value of any entry in $\mathcal A$. The state-of-the-art technique is a local search algorithm that subsequently moves in an improving direction where the number of iterations and the search for such an improving direction each take time $\Omega(n)$. This leads to a running time quadratic in $n$. We introduce a technique based on color coding which allows us to compute these improving directions in logarithmic time after a single initialization step. This yields an algorithm for $n$-fold ILPs with a running time that is near-linear in $nt$, the number of variables. More precisely, our algorithm runs in time $(rs\Delta)^{\mathcal{O}(r^2s + s^2)} L^2 nt \log^{\mathcal{O}(1)}(nt)$, where $L$ is the encoding length of the largest integer in the input. Further, in contrast to the algorithms in recent literature, we do not need to solve the LP relaxation in order to handle unbounded variables. Instead we give a structural lemma to introduce appropriate bounds. On the other hand, if we are given such an LP solution, the running time can be decreased by a factor of $L$.

Highlights

  • Solving integer linear programs of the form max {cT x | Ax = b, x ∈ Z≥0} is one of the most fundamental tasks in optimization

  • Handling unbounded variables in an n-fold ILP is a non-trivial issue in the previous algorithms from literature

  • We present an algorithm, which solves n-fold ILPs in timeO(r2s+s2)L · nt log5(nt) + LP, where LP is the time to solve the LP relaxation of the n-fold ILP

Read more

Summary

Introduction

LP is the running time required for solving the corresponding LP relaxation This augmentation algorithm is the last one in a line of research, where local improvement/augmenting steps are used to converge to an optimal solution. For example when analyzing big data, large real world graphs as in telecommunication networks or DNA strings in biology, the duration of the computation would go far beyond the scope of an acceptable running time [3, 6, 13] For this reason even problems which have an algorithm of quadratic running time are still studied from the viewpoint of approximation algorithms with the objective to obtain results in subquadratic time, even at the cost of a worse quality [3, 6, 13]. Handling unbounded variables in an n-fold ILP is a non-trivial issue in the previous algorithms from literature They had to solve the corresponding LP relaxation and use proximity results.

Summary of Results
Related Work
Preliminaries
Efficient Computation of Improving Directions
The Augmenting Step Algorithm
Bounds on 1-norm
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call