Approximation Algorithm Research Articles

Exact and approximation algorithms for routing a convoy through a graph

Exact and approximation algorithms for routing a convoy through a graph

Probabilistic Answer Set Programming with Discrete and Continuous Random Variables

Abstract Probabilistic Answer Set Programming under the credal semantics extends Answer Set Programming with probabilistic facts that represent uncertain information. The probabilistic facts are discrete with Bernoulli distributions. However, several real-world scenarios require a combination of both discrete and continuous random variables. In this paper, we extend the PASP framework to support continuous random variables and propose Hybrid Probabilistic Answer Set Programming. Moreover, we discuss, implement, and assess the performance of two exact algorithms based on projected answer set enumeration and knowledge compilation and two approximate algorithms based on sampling. Empirical results, also in line with known theoretical results, show that exact inference is feasible only for small instances, but knowledge compilation has a huge positive impact on performance. Sampling allows handling larger instances but sometimes requires an increasing amount of memory.

Algorithmic Persuasion with Evidence

In a game of persuasion with evidence, a sender has private information. By presenting evidence on the information, the sender wishes to persuade a receiver to take a single action (e.g., hire a job candidate, or convict a defendant). The sender’s utility depends solely on whether the receiver takes the action. The receiver’s utility depends on both the action and the sender’s private information. We study three natural variations. First, we consider the problem of computing an equilibrium of the game without commitment power. Second, we consider a persuasion variant, where the sender commits to a signaling scheme and the receiver, after seeing the evidence, takes the action or not. Third, we study a delegation variant, where the receiver first commits to taking the action if being presented certain evidence, and the sender presents evidence to maximize the probability the action is taken. We study these variants through the computational lens, and give hardness results, optimal approximation algorithms, and polynomial-time algorithms for special cases. Among our results is an approximation algorithm that rounds a semidefinite program that might be of independent interest, since, to the best of our knowledge, it is the first such approximation algorithm in algorithmic economics.

Machine Learning-enhanced Signature of Metastasis-related T Cell Marker Genes for Predicting Overall Survival in Malignant Melanoma.

In this study, we aimed to investigate disparities in the tumor immune microenvironment (TME) between primary and metastatic malignant melanoma (MM) using single-cell RNA sequencing (scRNA-seq) and to identify metastasis-related T cell marker genes (MRTMGs) for predicting patient survival using machine learning techniques. We identified 6 distinct T cell clusters in 10×scRNA-seq data utilizing the Uniform Manifold Approximation and Projection (UMAP) algorithm. Four machine learning algorithms highlighted SRGN, PMEL, GPR143, EIF4A2, and DSP as pivotal MRTMGs, forming the foundation of the MRTMGs signature. A high MRTMGs signature was found to be correlated with poorer overall survival (OS) and suppression of antitumor immunity in MM patients. We developed a nomogram that combines the MRTMGs signature with the T stage and N stage, which accurately predicts 1-year, 3-year, and 5-year OS probabilities. Furthermore, in an immunotherapy cohort, a high MRTMG signature was associated with an unfavorable response to anti-programmed death 1 (PD-1) therapy. In conclusion, primary and metastatic MM display distinct TME landscapes with different T cell subsets playing crucial roles in metastasis. The MRTMGs signature, established through machine learning, holds potential as a valuable biomarker for predicting the survival of MM patients and their response to anti-PD-1 therapy.

Application of Bivariate Reproducing Kernel-Based Best Interpolation Method in Electrical Tomography

Electrical Tomography (ET) technology is widely used in multiphase flow detection due to its advantages of low cost, visualization, fast response, non-radiation, and non-invasiveness. However, ill-posed solutions lead to low image reconstruction resolution, which limits its practical engineering applications. Although existing interpolation approximation algorithms can alleviate the effects of the ill-posed solutions to some extent, the imaging results remain suboptimal due to the limited approximation capability of these methods. This paper proposes a Bivariate Reproducing Kernel-Based Best Interpolation (BRKBI) method, which offers smaller approximation errors and clearer image reconstruction quality compared to existing methods. The effectiveness of the BRKBI method is validated through theoretical analysis and experimental comparisons.

Edge Manipulations for the Maximum Vertex-Weighted Bipartite \(b\) -matching

In this paper, we explore the Mechanism Design aspects of the Maximum Vertex-weighted \(b\) -Matching (MVbM) problem on bipartite graphs \((A\cup T,E)\) . The set \(A\) comprises agents, while \(T\) represents tasks. The set \(E\) , which connects \(A\) and \(T\) , is the private information of either agents or tasks. In this framework, we investigate three mechanisms – \(\mathbb{M}_{BFS}\) , \(\mathbb{M}_{DFS}\) , and \(\mathbb{M}_{G}\) . We examine scenarios in which either agents or tasks are strategic and report their adjacent edges to one of the three mechanisms. In both cases, we assume that the strategic entities are bounded by their statements: they can hide edges, but they cannot report edges that do not exist. First, we consider the case in which agents can manipulate. In this framework, \(\mathbb{M}_{BFS}\) and \(\mathbb{M}_{DFS}\) are optimal but not truthful. By characterizing the Nash Equilibria induced by \(\mathbb{M}_{BFS}\) and \(\mathbb{M}_{DFS}\) , we reveal that both mechanisms have a Price of Anarchy ( \(PoA\) ) and Price of Stability ( \(PoS\) ) of \(2\) . These efficiency guarantees are tight; no deterministic mechanism can achieve a lower \(PoA\) or \(PoS\) . In contrast, the third mechanism, \(\mathbb{M}_{G}\) , is not optimal, but truthful and its approximation ratio is \(2\) . We demonstrate that this ratio is optimal; no deterministic and truthful mechanism can outperform it. We then shift our focus to scenarios where tasks can exhibit strategic behaviour. In this case, \(\mathbb{M}_{BFS}\) , \(\mathbb{M}_{DFS}\) , and \(\mathbb{M}_{G}\) all maintain truthfulness, making \(\mathbb{M}_{BFS}\) and \(\mathbb{M}_{DFS}\) truthful and optimal mechanisms. In conclusion, we investigate the manipulability of \(\mathbb{M}_{BFS}\) and \(\mathbb{M}_{DFS}\) through experiments on randomly generated graphs. We observe that (i) \(\mathbb{M}_{BFS}\) is less prone to be manipulated by the first agent than \(\mathbb{M}_{DFS}\) , and (ii) \(\mathbb{M}_{BFS}\) is more manipulable on instances in which the total capacity of the agents is equal to the number of tasks. 1

Scaling whole-chip QAOA for higher-order ising spin glass models on heavy-hex graphs

We show that the quantum approximate optimization algorithm (QAOA) for higher-order, random coefficient, heavy-hex compatible spin glass Ising models has strong parameter concentration across problem sizes from 16 up to 127 qubits for p = 1 up to p = 5, which allows for computationally efficient parameter transfer of QAOA angles. Matrix product state (MPS) simulation is used to compute noise-free QAOA performance. Hardware-compatible short-depth QAOA circuits are executed on ensembles of 100 higher-order Ising models on noisy IBM quantum superconducting processors with 16, 27, and 127 qubits using QAOA angles learned from a single 16-qubit instance using the JuliQAOA tool. We show that the best quantum processors find lower energy solutions up to p = 2 or p = 3, and find mean energies that are about a factor of two off from the noise-free distribution. We show that p = 1 QAOA energy landscapes remain very similar as the problem size increases using NISQ hardware gridsearches with up to a 414 qubit processor.

A Novel Noise-Aware Classical Optimizer for Variational Quantum Algorithms

A key component of variational quantum algorithms (VQAs) is the choice of classical optimizer employed to update the parameterization of an ansatz. It is well recognized that quantum algorithms will, for the foreseeable future, necessarily be run on noisy devices with limited fidelities. Thus, the evaluation of an objective function (e.g., the guiding function in the quantum approximate optimization algorithm (QAOA) or the expectation of the electronic Hamiltonian in variational quantum eigensolver (VQE)) required by a classical optimizer is subject not only to stochastic error from estimating an expected value but also to error resulting from intermittent hardware noise. Model-based derivative-free optimization methods have emerged as popular choices of a classical optimizer in the noisy VQA setting, based on empirical studies. However, these optimization methods were not explicitly designed with the consideration of noise. In this work we adapt recent developments from the “noise-aware numerical optimization” literature to these commonly used derivative-free model-based methods. We introduce the key defining characteristics of these novel noise-aware derivative-free model-based methods that separate them from standard model-based methods. We study an implementation of such noise-aware derivative-free model-based methods and compare its performance on demonstrative VQA simulations to classical solvers packaged in scikit-quant. History: Accepted by Giacomo Nannicini, Area Editor for Quantum Computing and Operations Research. Accepted for Special Issue. Funding: This material is based upon work supported by the U.S. Department of Energy, Office of Science, National Quantum Information Science Research Centers and the Office of Advanced Scientific Computing Research, Accelerated Research for Quantum Computing program under contract number DE-AC02-06CH11357. Supplemental Material: The software that supports the findings of this study is available within the paper and its Supplemental Information ( https://pubsonline.informs.org/doi/suppl/10.1287/ijoc.2023.0177 ) as well as from the IJOC GitHub software repository ( https://github.com/INFORMSJoC/2023.0177 ). The complete IJOC Software and Data Repository is available at https://informsjoc.github.io/ .

Fixed-Ratio Approximation Algorithm for the Minimum Cost Cover of a Digraph by Bounded Number of Cycles

We consider polynomial time approximation for the minimum cost cycle cover problem of an edge-weighted digraph, where feasible covers are restricted to have at most k disjoint cycles. In the literature this problem is referred to as Min-k-SCCP. The problem is closely related to classic Traveling Salesman Problem (TSP) and Vehicle Routing Problem (VRP) and has many important applications in algorithms design and operations research. Unlike its unconstrained variant, the Min-k-SCCP is strongly NP-hard even on undirected graphs and remains intractable in very specific settings. For any metric, the problem can be approximated in polynomial time within ratio 2, while in fixed-dimensional Euclidean spaces it admits Polynomial Time Approximation Schemes (PTAS). In the same time, approximation of the more general asymmetric Min-k-SCCP still remains weakly studied. In this paper, we propose the first fixed-ratio approximation algorithm for this problem, which extends the recent breakthrough Svensson-Tarnawski-Vegh and Traub-Vygen results for the Asymmetric Traveling Salesman Problem.

Efficiently Constructing Convex Approximation Sets in Multiobjective Optimization Problems

Convex approximation sets for multiobjective optimization problems are a well-studied relaxation of the common notion of approximation sets. Instead of approximating each image of a feasible solution by the image of some solution in the approximation set up to a multiplicative factor in each component, a convex approximation set only requires this multiplicative approximation to be achieved by some convex combination of finitely many images of solutions in the set. This makes convex approximation sets efficiently computable for a wide range of multiobjective problems: even for many problems for which (classic) approximations sets are hard to compute. In this article, we propose a polynomial-time algorithm to compute convex approximation sets that builds on an exact or approximate algorithm for the weighted sum scalarization and is therefore applicable to a large variety of multiobjective optimization problems. The provided convex approximation quality is arbitrarily close to the approximation quality of the underlying algorithm for the weighted sum scalarization. In essence, our algorithm can be interpreted as an approximate version of the dual variant of Benson’s outer approximation algorithm. Thus, in contrast to existing convex approximation algorithms from the literature, information on solutions obtained during the approximation process is utilized to significantly reduce both the practical running time and the cardinality of the returned solution sets while still guaranteeing the same worst-case approximation quality. We underpin these advantages by the first comparison of all existing convex approximation algorithms on several instances of the triobjective knapsack problem and the triobjective symmetric metric traveling salesman problem. History: Accepted by Andrea Lodi, Area Editor for Design & Analysis of Algorithms–Discrete. Funding: This research was supported by the German Research Foundation [Project 398572517]. Supplemental Material: The software that supports the findings of this study is available within the paper and its Supplemental Information ( https://pubsonline.informs.org/doi/suppl/10.1287/ijoc.2023.0220 ) as well as from the IJOC GitHub software repository ( https://github.com/INFORMSJoC/2023.0220 ). The complete IJOC Software and Data Repository is available at https://informsjoc.github.io/ .

Efficient quantum circuits based on the quantum natural gradient

Efficient preparation of arbitrary entangled quantum states is crucial for quantum computation. This is particularly important for noisy intermediate-scale quantum simulators relying on variational hybrid quantum-classical algorithms. To that end, we propose symmetry-conserving modified quantum approximate optimization algorithm (SCom-QAOA) circuits. The depths of these circuits depend not only on the desired fidelity to the target state but also on the amount of entanglement the state contains. The parameters of the SCom-QAOA circuits are optimized using the quantum natural gradient method based on the Fubini-Study metric. The SCom-QAOA circuit transforms an unentangled state into a ground state of a gapped one-dimensional Hamiltonian with a circuit depth that depends not on the system size but rather on the finite correlation length. In contrast, the circuit depth grows proportionally to the system size for preparing low-lying states of critical one-dimensional systems. Even in the latter case, SCom-QAOA circuits with depth less than the system size were sufficient to generate states with fidelity in excess of 99%, which is relevant for near-term applications. The proposed scheme enlarges the set of the initial states accessible for variational quantum algorithms and widens the scope of investigation of nonequilibrium phenomena in quantum simulators. Published by the American Physical Society 2024

Improved Approximation Algorithms for Relational Clustering

Improved Approximation Algorithms for Relational Clustering

Estimating Asset Pricing Models in the Presence of Cross-Sectionally Correlated Pricing Errors

In this study, we propose an adversarial learning approach to the asset pricing model estimation problem which aims to find estimates of factors and loadings that capture time-series covariations while minimizing the worst-case cross-sectional pricing errors. The proposed estimator is defined by a novel min-max optimization problem in which finding a solution is known to be difficult. This contrasts with other related estimators that admit a well-defined analytic solution but do not effectively account for correlations among the pricing errors. To this end, we propose an approximate algorithm based on the alternating optimization procedure and empirically demonstrate that our proposed adversarial estimation framework outperforms other existing factor models, especially when the explanatory power of the pricing model is limited.

An approximation algorithm for the prize-collecting connected dominating set problem

An approximation algorithm for the prize-collecting connected dominating set problem

Multi-objective Route Planning of an Unmanned Air Vehicle in Continuous Terrain: An Exact and an Approximation Algorithm

Multi-objective Route Planning of an Unmanned Air Vehicle in Continuous Terrain: An Exact and an Approximation Algorithm

Algorithms for solving the inverse scattering problem for the Manakov model

The paper considers algorithms for solving inverse scattering problems based on the discretization of the Gelfand–Levitan–Marchenko integral equations, associated with the system of nonlinear Schrödinger equations of the Manakov model. The numerical algorithm of the first order approximation for solving the scattering problem is reduced to the inversion of a series of nested block Toeplitz matrices using the Levinson-type bordering method. Increasing the approximation accuracy violates the Toeplitz structure of block matrices. Two algorithms are described that solve this problem for second order accuracy. One algorithm uses a block version of the Levinson bordering algorithm, which recovers the Toeplitz structure of the matrix by moving some terms of the systems of equations to the right-hand side. Another algorithm is based on the Toeplitz decomposition of an almost block-Toeplitz matrix and the Tyrtyshnikov bordering algorithm. The speed and accuracy of calculations using the presented algorithms are compared on an exact solution (the Manakov vector soliton).

Continuous Approximate Dynamic Programming Algorithm to Promote Multiple Battery Energy Storage Lifespan Benefit in Real-Time Scheduling

Continuous Approximate Dynamic Programming Algorithm to Promote Multiple Battery Energy Storage Lifespan Benefit in Real-Time Scheduling

Parameter Optimization of Tuned Inerter Damper for Vibration Suppression in Structures with Damping

This study investigates the parameter optimization of tuned inerter damper (TID) in a single degree of freedom (SDOF) system, with a focus on vibration mitigation. As a type of inerter based dynamic vibration absorber (IDVA), TID achieves vibration suppression of the primary structure by replacing the mass block in traditional dynamic vibration absorber (DVA) with an inerter, thereby minimizing the increase in physical mass. With the displacement response of the primary structure as the objective function, this study first employs the classical fixed point theory (FPT) to derive the analytical formulations for the optimal parameters of TID following the H∞ optimization approach and neglecting the inherent damping of primary structure. The influence of optimal parameter deviations on the vibration mitigation effect of TID is also analyzed, revealing that the deviation of the optimal natural frequency ratio has a significant impact on the vibration mitigation performance. By flattening the amplitude curve of the displacement transfer function of primary structure between the fixed points, this study derives the analytical solution for the optimal damping ratio of TID using the extended fixed point theory (EFPT). When the inherent damping of primary structure is considered, this study employs the approximate extended fixed point theory (AEFPT) to derive the approximate optimal parameters of TID. A comparative study of the optimal natural frequency ratios obtained using FPT, AEFPT and numerical searches reveals that the discrepancies among the three methods are minimal for structures with low inherent damping. However, as the inherent damping ratio of the primary structure increases, the optimal natural frequency ratio obtained using FPT deviates significantly from the exact value, whereas the approximate optimal natural frequency ratio derived using AEFPT can significantly reduce this deviation. Combining the conclusion that the vibration mitigation effect of TID is significantly influenced by the natural frequency ratio, this study suggests that for highly damped primary structures, the use of AEFPT can yield more optimal TID parameters, thereby enhancing the robustness of vibration mitigation performance.

Refining uniform approximation algorithm for low-rank Chebyshev embeddings

Abstract Nowadays, low-rank approximations are a critical component of many numerical procedures. Traditionally the problem of low-rank approximation of matrices is solved in unitary invariant norms such as Frobenius or spectral norm due to the existence of efficient methods for constructing approximations. However, recent results discover the potential of low-rank approximations in the Chebyshev norm, which naturally arises in many applications. In this paper, we investigate the problem of uniform approximation of vectors, which is the main component in the low-rank approximations of matrices in the Chebyshev norm. The principal novelty of this paper is the accelerated algorithm for solving uniform approximation problems. We also analyze the iterative procedure of the proposed algorithm and demonstrate that it has a geometric convergence rate. Finally, we provide an extensive numerical evaluation, which demonstrates the effectiveness of the proposed procedures.

Optimized Noise Suppression for Quantum Circuits

Quantum computation promises to advance a wide range of computational tasks. However, current quantum hardware suffers from noise and is too small for error correction. Thus, accurately utilizing noisy quantum computers strongly relies on noise characterization, mitigation, and suppression. Crucially, these methods must also be efficient in terms of their classical and quantum overhead. Here, we efficiently characterize and mitigate crosstalk noise, which is a severe error source in, for example, cross-resonance based superconducting quantum processors. For crosstalk characterization, we develop a simplified measurement experiment. Furthermore, we analyze the problem of optimal experiment scheduling and solve it for common hardware architectures. After characterization, we mitigate noise in quantum circuits by a noise-aware qubit routing algorithm. Our integer programming algorithm extends previous work on optimized qubit routing by swap insertion. We incorporate the measured crosstalk errors in addition to other, more easily accessible noise data in the objective function. Furthermore, we strengthen the underlying integer linear model by proving a convex hull result about an associated class of polytopes, which has applications beyond this work. We evaluate the proposed method by characterizing crosstalk noise for two chips with up to 127 qubits and leverage the resulting data to improve the approximation ratio of the Quantum Approximate Optimization Algorithm by up to 10% compared with other established noise-aware routing methods. Our work clearly demonstrates the gains of including noise data when mapping abstract quantum circuits to hardware native ones. History: Accepted by Giacomo Nannicini, Area Editor for Quantum Computing and Operations Research. Accepted for Special Issue on Quantum Computing. Funding: This work was supported by Bavarian state government; Bayerische Staatsministerium für Wirtschaft, Landesentwicklung und Energie. Supplemental Material: The software that supports the findings of this study is available within the paper and its Supplemental Information ( https://pubsonline.informs.org/doi/suppl/10.1287/ijoc.2024.0551 ) as well as from the IJOC GitHub software repository ( https://github.com/INFORMSJoC/2024.0551 ). The complete IJOC Software and Data Repository is available at https://informsjoc.github.io/ .