Arising from structural graph theory, treewidth has become a focus of study in fixed-parameter tractable algorithms. Many NP-hard problems are known to be solvable in (O) over tilde (n . 2(O(tau))) time, where tau is ...
详细信息
ISBN:
(纸本)9781450380539
Arising from structural graph theory, treewidth has become a focus of study in fixed-parameter tractable algorithms. Many NP-hard problems are known to be solvable in (O) over tilde (n . 2(O(tau))) time, where tau is the treewidth of the input graph. Analogously, many problems in P should be solvable in (O) over tilde (n . tau(O(1))) time;however, due to the lack of appropriate tools, only a few such results are currently known. In our paper, we show this holds for linear programs: Given a linear program of the form min(Ax=b,l <= x <= u) c(T)x whose dual graph G(A) has treewidth tau, and a corresponding width-tau tree decomposition, we show how to solve it in time (O) over tilde (n . tau(2) log(1/epsilon)), where n is the number of variables and epsilon is the relative accuracy. When a tree decomposition is not given, we use existing techniques in vertex separators to obtain algorithms with (O) over tilde (n . tau(4)log(1/epsilon)) and (O) over tilde (n . tau(2)log(1/epsilon) + n(1.5)) run-times. Besides being the first of its kind, our algorithm has run-time nearly matching the fastest run-time for solving the sub-problem Ax = b (under the assumption that no fast matrix multiplication is used). We obtain these results by combining recent techniques in interior-point methods (IPMs), sketching, and a novel representation of the solution under a multiscale basis similar to the wavelet basis. This representation further yields the first IPM with o (rank(A)) time per iteration when the treewidth is small.
In this paper, we present a new algorithm for solving linear programs that requires only (O) over tilde (root rank(A)L) iterations where A is the constraint matrix of a linear program with m constraints, n variables, ...
详细信息
ISBN:
(纸本)9781479965175
In this paper, we present a new algorithm for solving linear programs that requires only (O) over tilde (root rank(A)L) iterations where A is the constraint matrix of a linear program with m constraints, n variables, and bit complexity L. Each iteration of our method consists of solving (O) over tilde (1) linear systems and additional nearly linear time computation. Our method improves upon the previous best iteration bounds by factor of (Omega) over tilde((m/rank (A)))(1/4)) for methods with polynomial time computable iterations and by (Omega) over tilde (m/rank (A)) 1/2) for methods which solve at most (O) over tilde (1) linear systems in each iteration each achieved over 20 years ago. Applying our techniques to the linear program formulation of maximum flow yields an (O) over tilde (vertical bar E vertical bar root vertical bar V vertical bar log(2) U) time algorithm for solving the maximum flow problem on directed graphs with vertical bar E vertical bar edges, vertical bar V vertical bar vertices, and capacity ratio U. This improves upon the previous fastest running time of O(vertical bar E vertical bar min{vertical bar E vertical bar(1/2), vertical bar V vertical bar(2/3)} log(vertical bar V vertical bar(2)/vertical bar E vertical bar) log(U)) achieved over 15 years ago by Goldberg and Rao and improves upon the previous best running times for solving dense directed unit capacity graphs of O(vertical bar E vertical bar min{vertical bar E vertical bar(1/2), vertical bar V vertical bar(2/3)}) achieved by Even and Tarjan over 35 years ago and a running time of O(vertical bar E vertical bar(10/7)) achieved recently by Madry.
We recently reported a criterion for blind separation of non-negative sources, using a new concept called convex analysis for mixtures of non-negative sources (CAMNS). Under some assumptions that are considered realis...
详细信息
ISBN:
(纸本)9781424414833
We recently reported a criterion for blind separation of non-negative sources, using a new concept called convex analysis for mixtures of non-negative sources (CAMNS). Under some assumptions that are considered realistic for sparse or high-contrast signals, the criterion is that the true source signals can be perfectly recovered by finding the extreme points of some observation-constructed convex set. In our last work we also developed methods for fulfilling the CAMNS criterion, but only for two to three sources. In this paper we propose a systematic linear programming (LP) based method that is applicable to any number of sources. The proposed method has two advantages. First, its dependence on LP means that the method does not suffer from local minima. Second, the maturity of LP solvers enables efficient implementation of the proposed method in practice. Simulation results are provided to demonstrate the efficacy of the proposed method.
This paper proposes a simple row-action type iterative algorithm which is appropriate to reconstruct sparse objects from a limited number of projections. The main idea is to use the L-1 norm to pick up a sparse soluti...
详细信息
ISBN:
(纸本)0780376366
This paper proposes a simple row-action type iterative algorithm which is appropriate to reconstruct sparse objects from a limited number of projections. The main idea is to use the L-1 norm to pick up a sparse solution from a set of feasible solutions to the measurement equation. By perturbing the linear program to a quadratic program, we use the duality of the nonlinear programming to construct a row-action type iterative algorithm to find a solution, we also prove that the algorithm converges for any initial image. We show that this method works well in the 3D blood-vessel reconstruction and its computation time is shorter compared to our previous method.
This paper shows how to solve linear programs of the form min(Ax=b, x >= 0) c (inverted perpendicular) x with n variables in time O*((n(omega) + n(2.5-alpha/2) + n(2+1/6)) log(n/delta)) where omega is the exponent ...
详细信息
ISBN:
(纸本)9781450367059
This paper shows how to solve linear programs of the form min(Ax=b, x >= 0) c (inverted perpendicular) x with n variables in time O*((n(omega) + n(2.5-alpha/2) + n(2+1/6)) log(n/delta)) where omega is the exponent of matrix multiplication, alpha is the dual exponent of matrix multiplication, and delta is the relative accuracy. For the current value of omega similar to 2.37 and alpha similar to 0.31, our algorithm takes O*(n(omega) log(n/delta)) time. When omega = 2, our algorithm takes O*(n(2+1/6) log(n/delta)) time. Our algorithm utilizes several new concepts that we believe may be of independent interest: (1) We define a stochastic central path method. (2) We show how to maintain a projection matrix root WA(inverted perpendicular)(AWA(inverted perpendicular))(-1)A root W in sub-quadratic time under l(2) multiplicative changes in the diagonal matrix W.
A computer program based on Lemke's complementary pivot algorithm is presented. This can be used to solve linear and quadratic programming problems. The program has been extensively tested on a wide range of probl...
详细信息
In chemical processes, we often take advantage of customizing daily problem-solving when a final solution will be obtained after solving a family of optimization problems repeatedly. In aid of the neural network, in t...
详细信息
In chemical processes, we often take advantage of customizing daily problem-solving when a final solution will be obtained after solving a family of optimization problems repeatedly. In aid of the neural network, in this note we have engaged in improving and customizing the algorithm of progressive linear programs (PROLP) after introducing manifold criteria to estimate active constraints. Through a few numerical experiments, the effect of the approach has been examined.
linear programs (LP) have been used popularly as practical optimization methods in many fields for many years. Due to the recent diversification of problem-solving, however, the problem must be formulated in a difficu...
详细信息
linear programs (LP) have been used popularly as practical optimization methods in many fields for many years. Due to the recent diversification of problem-solving, however, the problem must be formulated in a difficult and complicated manner. It is thus necessary to develop special methods that can solve large complex LP effectively as well as flexibly. A new approach for the simplex method, termed PAPA (pivot and probe algorithm), is considered promising for such situations. Since its effectiveness is confined to a special form of problem, we generalized it in this paper by applying a dual algorithm in the framework of the two-phase method for LP solution. In comparison with the revised simplex method, we examined numerically its efficiency by solving a variety of randomly generated test problems. Then, we paid special attention to application of the method in chemical processes. Thereupon, we pointed out the proposed method is especially effective where certain reference information is available to solve a family of problems repeatedly before arriving at a final result. This is why we propose calling the method ''progressive linear programs'' (PROLP).
Static traffic assignment aims to disclose the spatial distribution of vehicular flow over a transportation network subject to given traffic demands, and plays an essential role in transportation engineering. User-opt...
详细信息
Static traffic assignment aims to disclose the spatial distribution of vehicular flow over a transportation network subject to given traffic demands, and plays an essential role in transportation engineering. User-optimal pattern adheres the individual rationality of motorists, in which everyone chooses a route that minimizes his own travel cost, while considering congestion effects influenced by the aggregated movement of vehicles. User optimal traffic assignment, which is also known as the user equilibrium, entails solving an optimization problem with a strictly convex objective function and linear constraints. However, the performances of general-purpose solvers are quite disappointing. This paper proposes two highly efficient computation models for the user equilibrium problem. The first one exploits the second-order cone reformulation of a convex power function, resulting in a second-order cone program, and no approximation is incurred. The second one approximates the convex objective function using a piece-wise linear function, and comes down to a linear program. An adaptive path generation oracle is devised in order to circumvent path enumeration in problem setup. Case studies demonstrate that the proposed method can deal with large-scale transportation systems, and outperforms the most popular iterative algorithm in the literature.
Both the combinatorial and the circuit diameter of polyhedra are of interest to the theory of linear programming for their intimate connection to a best-case performance of linear programming algorithms. We study the ...
详细信息
Both the combinatorial and the circuit diameter of polyhedra are of interest to the theory of linear programming for their intimate connection to a best-case performance of linear programming algorithms. We study the diameters of dual network flow polyhedra associated to b-flows on directed graphs and prove quadratic upper bounds for both of them: the minimum of and for the combinatorial diameter, and for the circuit diameter. Previously, bounds on these diameters have only been known for bipartite graphs. The situation is much more involved for general graphs. In particular, we construct a family of dual network flow polyhedra with members that violate the circuit diameter bound for bipartite graphs by an arbitrary additive constant.
暂无评论