To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
Problems of wave interaction with a body with arbitrary shape floating or submerged in water are of immense importance in the literature on the linearized theory of water waves. Wave-free potentials are used to construct solutions to these problems involving bodies with circular geometry, such as a submerged or half-immersed long horizontal circular cylinder (in two dimensions) or sphere (in three dimensions). These are singular solutions of Laplace’s equation satisfying the free surface condition and decaying rapidly away from the point of singularity. Wave-free potentials in two and three dimensions for infinitely deep water as well as water of uniform finite depth with a free surface are known in the literature. The method of constructing wave-free potentials in three dimensions is presented here in a systematic manner, neglecting or taking into account the effect of surface tension at the free surface or for water with an ice cover modelled as a thin elastic plate floating on the water. The forms of the wave motion at the upper surface (free surface or ice-covered surface) related to these wave-free potentials are depicted graphically in a number of figures for all the cases considered.
This work focuses on finding optimal dividend payment and capital injection policies to maximize the present value of the difference between the cumulative dividend payment and the possible capital injections with delays. Starting from the classical Cramér–Lundberg process, using the dynamic programming approach, the value function obeys a quasi-variational inequality. With delays in capital injections, the company will be exposed to the risk of financial ruin during the delay period. In addition, the optimal dividend payment and capital injection strategy should balance the expected cost of the possible capital injections and the time value of the delay period. In this paper, the closed-form solution of the value function and the corresponding optimal policies are obtained. Some limiting cases are also discussed. A numerical example is presented to illustrate properties of the solution. Some economic insights are also given.
It has been known for a long time that the equivariant $2+1$ wave map into the $2$-sphere blows up if the initial data are chosen appropriately. Here, we present numerical evidence for the stability of the blow-up phenomenon under explicit violations of equivariance.
There is much interest within the mathematical biology and statistical physics community in converting stochastic agent-based models for random walkers into a partial differential equation description for the average agent density. Here a collection of noninteracting biased random walkers on a one-dimensional lattice is considered. The usual master equation approach requires that two continuum limits, involving three parameters, namely step length, time step and the random walk bias, approach zero in a specific way. We are interested in the case where the two limits are not consistent. New results are obtained using a Fokker–Planck equation and the results are highly dependent on the simulation update schemes. The theoretical results are confirmed with examples. These findings provide insight into the importance of updating schemes to an accurate macroscopic description of stochastic local movement rules in agent-based models when the lattice spacing represents a physical object such as cell diameter.
A stochastic algorithm for bound-constrained global optimization is described. The method can be applied to objective functions that are nonsmooth or even discontinuous. The algorithm forms a partition on the search region using classification and regression trees (CART), which defines a region where the objective function is relatively low. Further points are drawn directly from the low region before a new partition is formed. Alternating between partition and sampling phases provides an effective method for nonsmooth global optimization. The sequence of iterates generated by the algorithm is shown to converge to an essential global minimizer with probability one under mild conditions. Nonprobabilistic results are also given when random sampling is replaced with points taken from the Halton sequence. Numerical results are presented for both smooth and nonsmooth problems and show that the method is effective and competitive in practice.
A Content Distribution Network (CDN) can be defined as an overlay system that replicatescopies of contents at multiple points of a network, close to the final users, with theobjective of improving data access. CDN technology is widely used for the distribution oflarge-sized contents, like in video streaming. In this paper we address the problem offinding the best server for each customer request in CDNs, in order to minimize theoverall cost. We consider the problem as a transportation problem and a distributedalgorithm is proposed to solve it. The algorithm is composed of two independent phases: adistributed heuristic finds an initial solution that may be later improved by adistributed transportation simplex algorithm. It is compared with the sequential versionof the transportation simplex and with an auction-based distributed algorithm.Computational experiments carried out on a set of instances adapted from the literaturerevealed that our distributed approach has a performance similar or better to itssequential counterpart, in spite of not requiring global information about the contentsrequests. Moreover, the results also showed that the new method outperforms thebased-auction distributed algorithm.
The problem of finding structures with minimum stabbing number has received considerableattention from researchers. Particularly, [10]study the minimum stabbing number of perfect matchings (mspm), spanning trees(msst) and triangulations (mstr) associated to set of points in theplane. The complexity of the mstr remains open whilst the other two are known tobe 𝓝𝓟-hard. This paper presents integer programming(ip) formulations for these three problems, that allowed us to solve them tooptimality through ip branch-and-bound (b&b) or branch-and-cut(b&c) algorithms. Moreover, these models are the basis for the developmentof Lagrangian heuristics. Computational tests were conducted with instances taken from theliterature where the performance of the Lagrangian heuristics were compared with that ofthe exact b&b and b&c algorithms. The results reveal that theLagrangian heuristics yield solutions with minute, and often null, duality gaps forinstances with several hundreds of points in small computation times. To our knowledge,this is the first computational study ever reported in which these three stabbing problemsare considered and where provably optimal solutions are given.
In this study, we consider a scheduling environment with m(m ≥ 1) parallel machines.The set of jobs to schedule is divided into K disjoint subsets. Each subset of jobs isassociated with one agent. The K agents compete to perform their jobs on commonresources. The objective is to find a schedule that minimizes a global objective functionf0, while maintaining the regularobjective function of each agent, fk, at a level nogreater than a fixed value, εk (fk ∈ {fkmax, ∑fk}, k = 0, ..., K). This problem is a multi-agent schedulingproblem with a global objective function. In this study, we consider the casewith preemption and the case without preemption. If preemption is allowed, we propose apolynomial time algorithm based on a network flow approach for the unrelated parallelmachine case. If preemption is not allowed, we propose some general complexity results anddevelop dynamic programming algorithms.
In this article we study the realistic network topology of Synchronous Digital Hierarchy(SDH) networks. We describe how providers fulfill customer connectivity requirements. Weshow that SDH Network design reduces to the Non-Disjoint m-Ring-Star Problem (NDRSP). Wefirst show that there is no two-index integer formulation for this problem. We thenpresent a natural 3-index formulation for the NDRSP together with some classes of validinequalities that are used as cutting planes in a Branch-and-Cut approach. We propose apolyhedral study of a polytope associated with this formulation. Finally, we present ourBranch-and-Cut algorithm and give some experimental results on both random and realinstances.
In this paper we consider the problem of scheduling, on a two-machine flowshop, a set ofunit-time operations subject to time delays with respect to the makespan. This problem isknown to be \hbox{${\cal NP}$}𝒩𝒫-hard in the strong sense. We propose an algorithm basedon a branch and bound enumeration scheme. This algorithm includes the implementation ofnew lower and upper bound procedures, and dominance rules. A computer simulation tomeasure the performance of the algorithm is provided for a wide range of testproblems.
We study an uncapacitated facility location model where customers are served byfacilities of level one, then each level one facility that is opened must be assigned toan opened facility of level two. We identify a polynomially solvable case, and study somevalid inequalities and facets of the associated polytope.
This chapter reduces the inference problem in probabilistic graphical models to an equivalent maximum weight stable set problem on a graph. We discuss methods for recognizing when the latter problem can be solved efficiently by appealing to perfect graph theory. Furthermore, practical solvers based on convex programming and message-passing are presented.
Tractability is the study of computational tasks with the goal of identifying which problem classes are tractable or, in other words, efficiently solvable. The class of tractable problems is traditionally assumed to be solvable in polynomial time by a deterministic Turing machine and is denoted by P.The class contains many natural tasks such as sorting a set of numbers, linear programming (the decision version), determining if a number is prime, and finding a maximum weight matching. Many interesting problems, however, lie in another class that generalizes P and is known as NP: the class of languages decidable in polynomial time on a non-deterministic Turing machine. We trivially have that P isasubset of NP (many researchers also believe that it is a strict subset). It is believed that many problems in the class NP are, in the worst case, intractable and do not admit efficient inference. Problems such as maximum stable set, the traveling salesman problem and graph coloring are known to be NP-hard (at least as hard as the hardest problems in NP). It is, therefore, widely suspected that there are no polynomial-time algorithms for NP-hard problems.
This chapter covers methods for identifying islands of tractability for NP-hard combinatorial problems by exploiting suitable properties of their graphical structure. Acyclic structures are considered, as well as nearly-acyclic ones identified by means of so-called structural decomposition methods. In particular, the chapter focuses on the tree decomposition method, which is the most powerful decomposition method for graphs, and on the hypertree decomposition method, which is its natural counterpart for hypergraphs. These problem-decomposition methods give rise to corresponding notions of width of an instance, namely, treewidth and hypertree width. It turns out that many NP-hard problems can be solved efficiently over classes of instances of bounded treewidth or hypertree width: deciding whether a solution exists, computing a solution, and even computing an optimal solution (if some cost function over solutions is specified) are all polynomial-time tasks. Example applications include problems from artificial intelligence, databases, game theory, and combinatorial auctions.
Many NP-hard problems in different areas such as AI [42], Database Systems [6, 81], Game theory [45, 31, 20], and Network Design [34], are known to be efficiently solvable when restricted to instances whose underlying structures can be modeled via acyclic graphs or acyclic hypergraphs. For such restricted classes of instances, solutions can usually be computed via dynamic programming. However, as a matter of fact, (graphical) structures arising from real applications are in most relevant cases not properly acyclic. Yet, they are often not very intricate and exhibit some rather limited degree of cyclicity, which suffices to retain most of the nice properties of acyclic instances.
Machine learning and data analysis have driven explosive growth in interest in the methods of large-scale optimization. Many commonly used techniques such as stochastic-gradients date back several decades, but owing to their practical success they have gained great importance in machine learning. Before interior point methods totally dominated the field of optimization, first-order methods had already been studied and theoretically analyzed in substantial detail. But interest in these techniques skyrocketed after the prolific rise of applications in machine learning, signal processing, etc. This chapter is a brief introduction to this vast and flourishing area of large-scale optimization.
Introduction
Machine Learning (ML) broadly encompasses a variety of adaptive, autonomous, and intelligent tasks where one must “learn” to predict from observations and feedback. Throughout its evolution, ML has drawn heavily and successfully on optimization algorithms; this relation to optimization is not surprising as “learning” and “adapting” ultimately involve problems where some quality function must be optimized.
But the interaction between ML and optimization is now undergoing rapid change. The increased size, complexity, and variety of ML problems, not only prompts a refinement of existing optimization techniques, but also spurs development of new methods tuned to the specific needs of ML applications.
In particular, ML applications must usually cope with large-scale data, which forces us to prefer “simpler,” perhaps less accurate but more scalable algorithms. Such methods can also crunch through more data, and may actually be better suited for learning – for a more precise characterization see [11]. The use of possibly less accurate methods is also grounded in pragmatic concerns: modeling limitations, observational noise, uncertainty, and computational errors are pervasive in real data.
Optimization problems are often hard to solve precisely. However solutions that are only nearly optimal are often good enough in practical applications. Approximation algorithms can find such solutions efficiently for many interesting problems. Profound theoretical results additionally help us understand what problems are approximable. This chapter gives an overview of existing approximation techniques, along five broad categories: greedy algorithms, linear and semi-definite programming relaxations, metric embeddings and special techniques. It concludes with an overview of the main inapproximability results.
Introduction
NP-hard optimization problems are ubiquitous, and unless P=NP, we cannot expect algorithms that find optimal solutions on all instances in polynomial time. This intractability thus forces us to relax one of the three above mentioned constraints. Approximation algorithms relax the optimality constraint, and aim to do so by as small an amount as possible. We shall concern ourselves with discrete optimization problems, where the goal is to find amongst the set of feasible solutions, the one that minimizes (or maximizes) the value of the objective function. Usually, the space of feasible solutions is defined implicitly, e.g. the set of cuts in a graph on n vertices. The objective function associates with each feasible solution a real value; this usually has a succinct representation as well, e.g. the number of edges in the cut. We measure the performance of an approximation algorithm on a given instance by the ratio of the value of the solution output by the algorithm, to that of the optimal solution.
In this chapter we will introduce submodularity and some of its generalizations, illustrate how it arises in various applications, and discuss algorithms for optimizing submodular functions.
Submodularity is a property of set functions with deep theoretical consequences and far-reaching applications. At first glance it seems very similar to concavity, in other ways it resembles convexity. It appears in a wide variety of applications: in Computer Science it has recently been identified and utilized in domains such as viral marketing [39], information gathering [44], image segmentation [10, 40, 36], document summarization [56], and speeding up satisfiability solvers [73]. Our emphasis in this chapter is on maximization; there are many important results and applications related to minimizing submodular functions that we do not cover.
As a concrete running example, we will consider the problem of deploying sensors in a drinking water distribution network (see Figure 3.1) in order to detect contamination. In this domain, we may have a model of how contaminants, accidentally or maliciously introduced into the network, spread over time. Such a model then allows to quantify the benefit f(A) of deploying sensors at a particular set A of locations (junctions or pipes in the network) in terms of the detection performance (such as average time to detection).
Based on this notion of utility, we then wish to find an optimal subset A ⊆ V of locations maximizing the utility, maxAf(A), subject to some constraints (such as bounded cost). This application requires solving a difficult real-world optimization problem, that can be handled with the techniques discussed in this chapter (Krause et al. [49] show in detail how submodular optimization can be applied in this domain.)