scispace - formally typeset
Search or ask a question

Showing papers in "arXiv: Discrete Mathematics in 2015"


Posted Content
TL;DR: In this paper, the authors survey recent results on temporal graphs and temporal graph problems that have appeared in the Computer Science community, and present a set of modern systems and applications that can be naturally modeled and studied via temporal graphs.
Abstract: A \emph{temporal graph} is, informally speaking, a graph that changes with time. When time is discrete and only the relationships between the participating entities may change and not the entities themselves, a temporal graph may be viewed as a sequence $G_1,G_2\ldots,G_l$ of static graphs over the same (static) set of nodes $V$. Though static graphs have been extensively studied, for their temporal generalization we are still far from having a concrete set of structural and algorithmic principles. Recent research shows that many graph properties and problems become radically different and usually substantially more difficult when an extra time dimension in added to them. Moreover, there is already a rich and rapidly growing set of modern systems and applications that can be naturally modeled and studied via temporal graphs. This, further motivates the need for the development of a temporal extension of graph theory. We survey here recent results on temporal graphs and temporal graph problems that have appeared in the Computer Science community.

149 citations


Posted Content
TL;DR: It is proved that for large (bounded) average degree d, the value of this semidefinite program (SDP) is --with high probability-- 2n*sqrt(d) + n, o(sqrt (d))+o(n), matching a spectral upper bound.
Abstract: Denote by $A$ the adjacency matrix of an Erdos-Renyi graph with bounded average degree. We consider the problem of maximizing $\langle A-E\{A\},X\rangle$ over the set of positive semidefinite matrices $X$ with diagonal entries $X_{ii}=1$. We prove that for large (bounded) average degree $d$, the value of this semidefinite program (SDP) is --with high probability-- $2n\sqrt{d} + n\, o(\sqrt{d})+o(n)$. For a random regular graph of degree $d$, we prove that the SDP value is $2n\sqrt{d-1}+o(n)$, matching a spectral upper bound. Informally, Erdos-Renyi graphs appear to behave similarly to random regular graphs for semidefinite programming. We next consider the sparse, two-groups, symmetric community detection problem (also known as planted partition). We establish that SDP achieves the information-theoretically optimal detection threshold for large (bounded) degree. Namely, under this model, the vertex set is partitioned into subsets of size $n/2$, with edge probability $a/n$ (within group) and $b/n$ (across). We prove that SDP detects the partition with high probability provided $(a-b)^2/(4d)> 1+o_{d}(1)$, with $d= (a+b)/2$. By comparison, the information theoretic threshold for detecting the hidden partition is $(a-b)^2/(4d)> 1$: SDP is nearly optimal for large bounded average degree. Our proof is based on tools from different research areas: $(i)$ A new `higher-rank' Grothendieck inequality for symmetric matrices; $(ii)$ An interpolation method inspired from statistical physics; $(iii)$ An analysis of the eigenvectors of deformed Gaussian random matrices.

109 citations


Posted Content
TL;DR: In this article, the algebraic representation and basic algorithms for MAGs are presented, as well as a set of basic MAG algorithms, constructed from well-known graph algorithms, such as degree computing, Breadth First Search, and Depth First Search (DFS), which can be used as primitives for building other more sophisticated MAG algorithms.
Abstract: We present the algebraic representation and basic algorithms for MultiAspect Graphs (MAGs). A MAG is a structure capable of representing multilayer and time-varying networks, as well as higher-order networks, while also having the property of being isomorphic to a directed graph. In particular, we show that, as a consequence of the properties associated with the MAG structure, a MAG can be represented in matrix form. Moreover, we also show that any possible MAG function (algorithm) can be obtained from this matrix-based representation. This is an important theoretical result since it paves the way for adapting well-known graph algorithms for application in MAGs. We present a set of basic MAG algorithms, constructed from well-known graph algorithms, such as degree computing, Breadth First Search (BFS), and Depth First Search (DFS). These algorithms adapted to the MAG context can be used as primitives for building other more sophisticated MAG algorithms. Therefore, such examples can be seen as guidelines on how to properly derive MAG algorithms from basic algorithms on directed graph. We also make available Python implementations of all the algorithms presented in this paper.

70 citations


Journal ArticleDOI
TL;DR: A new aperiodic tile set containing 11 Wang tiles on 4 colors is presented, which is minimal in the sense that no Wang set with less than 11 tiles is a Periodic.
Abstract: We present a new aperiodic tileset containing 11 Wang tiles on 4 colors, and we show that this tileset is minimal, in the sense that no Wang set with either fewer than 11 tiles or fewer than 4 colors is aperiodic. This gives a definitive answer to the problem raised by Wang in 1961.

67 citations


Posted Content
TL;DR: Improved bounds are obtained for a suite of problems in effective algebraic geometry, including Hilbert nullstellensatz, radical membership and counting rational points in low degree varieties.
Abstract: Let $f$ be a polynomial of degree $d$ in $n$ variables over a finite field $\mathbb{F}$. The polynomial is said to be unbiased if the distribution of $f(x)$ for a uniform input $x \in \mathbb{F}^n$ is close to the uniform distribution over $\mathbb{F}$, and is called biased otherwise. The polynomial is said to have low rank if it can be expressed as a composition of a few lower degree polynomials. Green and Tao [Contrib. Discrete Math 2009] and Kaufman and Lovett [FOCS 2008] showed that bias implies low rank for fixed degree polynomials over fixed prime fields. This lies at the heart of many tools in higher order Fourier analysis. In this work, we extend this result to all prime fields (of size possibly growing with $n$). We also provide a generalization to nonprime fields in the large characteristic case. However, we state all our applications in the prime field setting for the sake of simplicity of presentation. As an immediate application, we obtain improved bounds for a suite of problems in effective algebraic geometry, including Hilbert nullstellensatz, radical membership and counting rational points in low degree varieties. Using the above generalization to large fields as a starting point, we are also able to settle the list decoding radius of fixed degree Reed-Muller codes over growing fields. The case of fixed size fields was solved by Bhowmick and Lovett [STOC 2015], which resolved a conjecture of Gopalan-Klivans-Zuckerman [STOC 2008]. Here, we show that the list decoding radius is equal the minimum distance of the code for all fixed degrees, even when the field size is possibly growing with $n$.

45 citations


Posted Content
TL;DR: It is proved that for large (bounded) average degree $\gamma$, the value of this semidefinite program (SDP) is -with high probability- $2n\sqrt{gamma} + n\, o(\sqrt{\gamma})+ o(n)$.
Abstract: Denote by $A$ the adjacency matrix of an Erd\H{o}s-Renyi graph with bounded average degree. We consider the problem of maximizing $\langle A-{\mathbb E}\{A\},X\rangle$ over the set of positive semidefinite matrices $X$ with diagonal entries $X_{ii}=1$. We prove that for large (bounded) average degree $\gamma$, the value of this semidefinite program (SDP) is -with high probability- $2n\sqrt{\gamma} + n\, o(\sqrt{\gamma})+ o(n)$. Our proof is based on two tools from different research areas. First, we develop a new `higher-rank' Grothendieck inequality for symmetric matrices. In the present case, our inequality implies that the value of the above SDP is arbitrarily well approximated by optimizing over rank-$k$ matrices for $k$ large but bounded. Second, we use the interpolation method from spin glass theory to approximate this problem by a second one concerning Wigner random matrices instead of sparse graphs. As an application of our results, we prove new bounds on community detection via SDP that are substantially more accurate than the state of the art.

39 citations


Posted Content
TL;DR: This work gives two efficient algorithms for computing shortest time-respecting paths on a temporal network, and proves that there is a natural analogue of Menger’s theorem holding for arbitrary temporal networks.
Abstract: In this work we consider \emph{temporal networks}, i.e. networks defined by a \emph{labeling} $\lambda$ assigning to each edge of an \emph{underlying graph} $G$ a set of \emph{discrete} time-labels. The labels of an edge, which are natural numbers, indicate the discrete time moments at which the edge is available. We focus on \emph{path problems} of temporal networks. In particular, we consider \emph{time-respecting} paths, i.e. paths whose edges are assigned by $\lambda$ a strictly increasing sequence of labels. We begin by giving two efficient algorithms for computing shortest time-respecting paths on a temporal network. We then prove that there is a \emph{natural analogue of Menger's theorem} holding for arbitrary temporal networks. Finally, we propose two \emph{cost minimization parameters} for temporal network design. One is the \emph{temporality} of $G$, in which the goal is to minimize the maximum number of labels of an edge, and the other is the \emph{temporal cost} of $G$, in which the goal is to minimize the total number of labels used. Optimization of these parameters is performed subject to some \emph{connectivity constraint}. We prove several lower and upper bounds for the temporality and the temporal cost of some very basic graph families such as rings, directed acyclic graphs, and trees.

39 citations


Posted Content
TL;DR: In this paper, the authors define a combinatorial class of native clique minors in Chimera graphs with vertex images of uniform, near minimal size, and provide a polynomial-time algorithm that finds a maximum clique minor in a given induced subgraph of a Chimera graph.
Abstract: The current generation of D-Wave quantum annealing processor is designed to minimize the energy of an Ising spin configuration whose pairwise interactions lie on the edges of a {\em Chimera} graph $\mathcal C_{M,N,L}$. In order to solve an Ising spin problem with arbitrary pairwise interaction structure, the corresponding graph must be minor-embedded into a Chimera graph. We define a combinatorial class of {\em native clique minors} in Chimera graphs with vertex images of uniform, near minimal size, and provide a polynomial-time algorithm that finds a maximum native clique minor in a given induced subgraph of a Chimera graph. These minors allow improvement over recent work and have immediate practical applications in the field of quantum annealing.

33 citations


Posted Content
TL;DR: Extensions and modifications of the classical degree sequence graphic realization problem studied by Erd\H{o}s-Gallai and Havel-Hakimi are introduced, as well as of the corresponding connected graphic realization version.
Abstract: In this paper we introduce extensions and modifications of the classical degree sequence graphic realization problem studied by Erd\H{o}s-Gallai and Havel-Hakimi, as well as of the corresponding connected graphic realization version. We define the joint-degree matrix graphic (resp. connected graphic) realization problem, where in addition to the degree sequence, the exact number of desired edges between vertices of different degree classes is also specified. We give necessary and sufficient conditions, and polynomial time decision and construction algorithms for the graphic and connected graphic realization problems. These problems arise naturally in the current topic of graph modeling for complex networks. From the technical point of view, the joint-degree matrix realization algorithm is straightforward. However, the connected joint-degree matrix realization algorithm involves a novel recursive search of suitable local graph modifications. Also, we outline directions for further work of both theoretical and practical interest. In particular, we give a Markov chain which converges to the uniform distribution over all realizations. We show that the underline state space is connected, and we leave the question of the mixing rate open.

30 citations


Journal ArticleDOI
TL;DR: A meta-heuristic solver is developed that utilizes D-Wave Systems’ quantum annealer (or any other QUBO problem optimizer) to solve larger or denser problems, by iteratively solving subproblems, while keeping the rest of the variables fixed.
Abstract: A quantum annealer heuristically minimizes quadratic unconstrained binary optimization (QUBO) problems, but is limited by the physical hardware in the size and density of the problems it can handle. We have developed a meta-heuristic solver that utilizes D-Wave Systems' quantum annealer (or any other QUBO problem optimizer) to solve larger or denser problems, by iteratively solving subproblems, while keeping the rest of the variables fixed. We present our algorithm, several variants, and the results for the optimization of standard QUBO problem instances from OR-Library of sizes 500 and 2500 as well as the Palubeckis instances of sizes 3000 to 7000. For practical use of the solver, we show the dependence of the time to best solution on the desired gap to the best known solution. In addition, we study the dependence of the gap and the time to best solution on the size of the problems solved by the underlying optimizer.

30 citations


Posted Content
TL;DR: In this paper, the authors examined the trees extremal with respect to the $F$-index and showed that the second sum has an exceptional applicative potential and is known as the "em forgotten index".
Abstract: In a study on the structure--dependency of the total $\pi$-electron energy from 1972, Trinajstic and one of the present authors have shown that it depends on the sums $\sum_{v\in V}d(v)^2$ and $\sum_{v\in V}d(v)^3$, where $d(v)$ is the degree of a vertex $v$ of the underling molecular graph $G$. The first sum was later named {\it first Zagreb index} and over the years became one of the most investigated graph--based molecular structure descriptors. On the other hand, the second sum, except in very few works on the general first Zagreb index and the zeroth--order general Randic index, has been almost completely neglected. Recently, this second sum was named {\em forgotten index}, or shortly the $F$-{\em index}, and shown to have an exceptional applicative potential. In this paper we examine the trees extremal with respect to the $F$-index.

Posted Content
TL;DR: This manuscript provides a more complete picture on the relationship between submodular with convexity and concavity, by extending many of the results connecting submodularity with conveXity to the concave aspects of sub Modular maximization.
Abstract: Seminal work by Edmonds and Lovasz shows the strong connection between submodularity and convexity. Submodular functions have tight modular lower bounds, and subdifferentials in a manner akin to convex functions. They also admit poly-time algorithms for minimization and satisfy the Fenchel duality theorem and the Discrete Seperation Theorem, both of which are fundamental characteristics of convex functions. Submodular functions also show signs similar to concavity. Submodular maximization, though NP hard, admits constant factor approximation guarantees. Concave functions composed with modular functions are submodular, and they also satisfy diminishing returns property. This manuscript provides a more complete picture on the relationship between submodularity with convexity and concavity, by extending many of the results connecting submodularity with convexity to the concave aspects of submodularity. We first show the existence of superdifferentials, and efficiently computable tight modular upper bounds of a submodular function. While we show that it is hard to characterize this polyhedron, we obtain inner and outer bounds on the superdifferential along with certain specific and useful supergradients. We then investigate forms of concave extensions of submodular functions and show interesting relationships to submodular maximization. We next show connections between optimality conditions over the superdifferentials and submodular maximization, and show how forms of approximate optimality conditions translate into approximation factors for maximization. We end this paper by studying versions of the discrete seperation theorem and the Fenchel duality theorem when seen from the concave point of view. In every case, we relate our results to the existing results from the convex point of view, thereby improving the analysis of the relationship between submodularity, convexity, and concavity.

Posted Content
TL;DR: In this article, the authors extend the graph uncertainty principle to weighted graphs and show that a naive extension of their definitions leads to inconsistent results such as discontinuity of graph spread when regarded as a function of the graph structure.
Abstract: The uncertainty principle states that a signal cannot be localized both in time and frequency. With the aim of extending this result to signals on graphs, Agaskar&Lu introduce notions of graph and spectral spreads. They show that a graph uncertainty principle holds for some families of unweighted graphs. This principle states that a signal cannot be simultaneously localized both in graph and spectral domains. In this paper, we aim to extend their work to weighted graphs. We show that a naive extension of their definitions leads to inconsistent results such as discontinuity of the graph spread when regarded as a function of the graph structure. To circumvent this problem, we propose another definition of graph spread that relies on an inverse similarity matrix. We also discuss the choice of the distance function that appears in this definition. Finally, we compute and plot uncertainty curves for families of weighted graphs.

Journal ArticleDOI
TL;DR: In this article, the authors presented a constant approximation algorithm for the minimum weight (k,m)$-CDS problem in unit disk graphs under the assumption that $k and $m$ are two fixed constants with $m\geq k.
Abstract: In a wireless sensor network, the virtual backbone plays an important role. Due to accidental damage or energy depletion, it is desirable that the virtual backbone is fault-tolerant. A fault-tolerant virtual backbone can be modeled as a $k$-connected $m$-fold dominating set ($(k,m)$-CDS for short). In this paper, we present a constant approximation algorithm for the minimum weight $(k,m)$-CDS problem in unit disk graphs under the assumption that $k$ and $m$ are two fixed constants with $m\geq k$. Prior to this work, constant approximation algorithms are known for $k=1$ with weight and $2\leq k\leq 3$ without weight. Our result is the first constant approximation algorithm for the $(k,m)$-CDS problem with general $k,m$ and with weight. The performance ratio is $(\alpha+2.5k\rho)$ for $k\geq 3$ and $(\alpha+2.5\rho)$ for $k=2$, where $\alpha$ is the performance ratio for the minimum weight $m$-fold dominating set problem and $\rho$ is the performance ratio for the subset $k$-connected subgraph problem (both problems are known to have constant performance ratios.)

Posted Content
TL;DR: In this article, the authors studied the problem of low-pattern complexity in finite-word configurations and subshifts of low pattern complexity using tools of algebraic geometry, and showed that the configuration must be a sum of periodic configurations over integers, possibly with unbounded values.
Abstract: We study multidimensional configurations (infinite words) and subshifts of low pattern complexity using tools of algebraic geometry. We express the configuration as a multivariate formal power series over integers and investigate the setup when there is a non-trivial annihilating polynomial: a non-zero polynomial whose formal product with the power series is zero. Such annihilator exists, for example, if the number of distinct patterns of some finite shape D in the configuration is at most the size |D| of the shape. This is our low pattern complexity assumption. We prove that the configuration must be a sum of periodic configurations over integers, possibly with unbounded values. As a specific application of the method we obtain an asymptotic version of the well-known Nivat's conjecture: we prove that any two-dimensional, non-periodic configuration can satisfy the low pattern complexity assumption with respect to only finitely many distinct rectangular shapes D.

Posted Content
TL;DR: It is shown that determining whether a graph has a cograph 2-decomposition, resp.
Abstract: Symbolic ultrametrics define edge-colored complete graphs K_n and yield a simple tree representation of K_n. We discuss, under which conditions this idea can be generalized to find a symbolic ultrametric that, in addition, distinguishes between edges and non-edges of arbitrary graphs G=(V,E) and thus, yielding a simple tree representation of G. We prove that such a symbolic ultrametric can only be defined for G if and only if G is a so-called cograph. A cograph is uniquely determined by a so-called cotree. As not all graphs are cographs, we ask, furthermore, what is the minimum number of cotrees needed to represent the topology of G. The latter problem is equivalent to find an optimal cograph edge k-decomposition {E_1,...,E_k} of E so that each subgraph (V,E_i) of G is a cograph. An upper bound for the integer k is derived and it is shown that determining whether a graph has a cograph 2-decomposition, resp., 2-partition is NP-complete.

Posted Content
TL;DR: The main technical result is that computing the maximum number of redundant labels is APX-hard, i.e., there is no PTAS unless P = NP, and the existence of minimal temporal designs with at least nlogn labels is shown.
Abstract: We study the design of small cost temporally connected graphs, under various constraints. We mainly consider undirected graphs of $n$ vertices, where each edge has an associated set of discrete availability instances (labels). A journey from vertex $u$ to vertex $v$ is a path from $u$ to $v$ where successive path edges have strictly increasing labels. A graph is temporally connected iff there is a $(u,v)$-journey for any pair of vertices $u,v,~u ot= v$. We first give a simple polynomial-time algorithm to check whether a given temporal graph is temporally connected. We then consider the case in which a designer of temporal graphs can \emph{freely choose} availability instances for all edges and aims for temporal connectivity with very small \emph{cost}; the cost is the total number of availability instances used. We achieve this via a simple polynomial-time procedure which derives designs of cost linear in $n$. We also show that the above procedure is (almost) optimal when the underlying graph is a tree, by proving a lower bound on the cost for any tree. However, there are pragmatic cases where one is not free to design a temporally connected graph anew, but is instead \emph{given} a temporal graph design with the claim that it is temporally connected, and wishes to make it more cost-efficient by removing labels without destroying temporal connectivity (redundant labels). Our main technical result is that computing the maximum number of redundant labels is APX-hard, i.e., there is no PTAS unless $P=NP$. On the positive side, we show that in dense graphs with random edge availabilities, there is asymptotically almost surely a very large number of redundant labels. A temporal design may, however, be \emph{minimal}, i.e., no redundant labels exist. We show the existence of minimal temporal designs with at least $n \log{n}$ labels.

Posted Content
TL;DR: In this paper, it was shown that the conjecture holds for any constant k ≥ 4 and that for any k ≥ 3, all k-colourings of a cubic graph that is neither complete nor the triangular prism are Kempe equivalent.
Abstract: Let $G$ be a graph with a vertex colouring $\alpha$. Let $a$ and $b$ be two colours. Then a connected component of the subgraph induced by those vertices coloured either $a$ or $b$ is known as a Kempe chain. A colouring of $G$ obtained from $\alpha$ by swapping the colours on the vertices of a Kempe chain is said to have been obtained by a Kempe change. Two colourings of $G$ are Kempe equivalent if one can be obtained from the other by a sequence of Kempe changes. A conjecture of Mohar (2007) asserts that, for $k \geq 3$, all $k$-colourings of a $k$-regular graph that is not complete are Kempe equivalent. It was later shown that all $3$-colourings of a cubic graph that is neither $K_4$ nor the triangular prism are Kempe equivalent. In this paper, we prove that the conjecture holds for each $k\geq 4$. We also report the implications of this result on the validity of the Wang-Swendsen-Kotecký algorithm for the antiferromagnetic Potts model at zero-temperature.

Posted Content
TL;DR: In this article, a method based on a topological entropy for addressing the comparison between discrete piecewise linear functions is presented, which is certified by the stability theorem for persistent entropy.
Abstract: In this paper we present a novel methodology based on a topological entropy, the so-called persistent entropy, for addressing the comparison between discrete piecewise linear functions. The comparison is certified by the stability theorem for persistent entropy. The theorem is used in the implementation of a new algorithm. The algorithm transforms a discrete piecewise linear function into a filtered simplicial complex that is analyzed with persistent homology and persistent entropy. Persistent entropy is used as discriminant feature for solving the supervised classification problem of real long length noisy signals of DC electrical motors. The quality of classification is stated in terms of the area under receiver operating characteristic curve (AUC=94.52%).

Journal ArticleDOI
TL;DR: A gossip protocol is a procedure for spreading secrets among a group of agents, using a connection graph, and it is defined in terms of the topology of the graphs on which they are successful, wherein strong success (the protocol always terminates, possibly assuming fair scheduling) from weak success ( the protocol sometimes terminates).
Abstract: A gossip protocol is a procedure for spreading secrets among a group of agents, using a connection graph. The goal is for all agents to get to know all secrets, in which case we call the execution of the protocol successful. We consider distributed and dynamic gossip protocols. In distributed gossip the agents themselves instead of a global scheduler determine whom to call. In dynamic gossip not only secrets are exchanged but also telephone numbers (agent identities). This results in increased graph connectivity. We define six such distributed dynamic gossip protocols, and we characterize them in terms of the topology of the graphs on which they are successful, wherein we distinguish strong success (the protocol always terminates, possibly assuming fair scheduling) from weak success (the protocol sometimes terminates). For five of these protocols strong (fair) and weak success are characterized by weakly connected graphs. This result is surprising because the protocols are fairly different. In the sixth protocol an agent may only call another agent if it does not know the other agent's secret. Strong success for this protocol is characterized by graphs for which the set of non-terminal nodes is strongly connected. Weak success for this protocol is characterized by weakly connected graphs satisfying further topological constraints that we define in the paper. One direction of this characterization is surprisingly harder to prove than the other results in this contribution.

Posted Content
TL;DR: This work argues that an additional step can help: reassembling the spanning trees before sampling, and develops a deterministic polynomial-time algorithm that improves on Sebo's previously best approximation ratio of $\frac{8}{5}$.
Abstract: Many recent approximation algorithms for different variants of the traveling salesman problem (asymmetric TSP, graph TSP, s-t-path TSP) exploit the well-known fact that a solution of the natural linear programming relaxation can be written as convex combination of spanning trees. The main argument then is that randomly sampling a tree from such a distribution and then completing the tree to a tour at minimum cost yields a better approximation guarantee than simply taking a minimum cost spanning tree (as in Christofides' algorithm). We argue that an additional step can help: reassembling the spanning trees before sampling. Exchanging two edges in a pair of spanning trees can improve their properties under certain conditions. We demonstrate the usefulness for the metric s-t-path TSP by devising a deterministic polynomial-time algorithm that improves on Sebő's previously best approximation ratio of 8/5.

Posted Content
TL;DR: Bundala et al. as discussed by the authors improved upon the known upper bounds for minimal depth sorting networks on $17, 19,$ and $20$ channels and showed that no sorting network using less layers exists.
Abstract: We present new parallel sorting networks for $17$ to $20$ inputs. For $17, 19,$ and $20$ inputs these new networks are faster (i.e., they require less computation steps) than the previously known best networks. Therefore, we improve upon the known upper bounds for minimal depth sorting networks on $17, 19,$ and $20$ channels. Furthermore, we show that our sorting network for $17$ inputs is optimal in the sense that no sorting network using less layers exists. This solves the main open problem of [D. Bundala & J. Zavodný. Optimal sorting networks, Proc. LATA 2014].

Posted Content
TL;DR: Practical methods to estimate Laplacian spectra of direct and strong product graphs from spectral properties of their factor graphs using a few heuristic assumptions are developed.
Abstract: Calculating a product of multiple graphs has been studied in mathematics, engineering, computer science, and more recently in network science, particularly in the context of multilayer networks. One of the important questions to be addressed in this area is how to characterize spectral properties of a product graph using those of its factor graphs. While several such characterizations have already been obtained analytically (mostly for adjacency spectra), characterization of Laplacian spectra of direct product and strong product graphs has remained an open problem. Here we develop practical methods to estimate Laplacian spectra of direct and strong product graphs from spectral properties of their factor graphs using a few heuristic assumptions. Numerical experiments showed that the proposed methods produced reasonable estimation with percentage errors confined within a +/-10% range for most eigenvalues.

Posted Content
TL;DR: The characterization of symbolic ultrametrics that additionally distinguishes between edges and non-edges of arbitrary edge-colored graphs G and thus, yielding a tree representation of G, by means of so-called cographs is characterized, resulting in several new open problems, and NP-hardness results.
Abstract: Tree representations of (sets of) symmetric binary relations, or equivalently edge-colored undirected graphs, are of central interest, e.g.\ in phylogenomics. In this context symbolic ultrametrics play a crucial role. Symbolic ultrametrics define an edge-colored complete graph that allows to represent the topology of this graph as a vertex-colored tree. Here, we are interested in the structure and the complexity of certain combinatorial problems resulting from considerations based on symbolic ultrametrics, and on algorithms to solve them. This includes, the characterization of symbolic ultrametrics that additionally distinguishes between edges and non-edges of \emph{arbitrary} edge-colored graphs $G$ and thus, yielding a tree representation of $G$, by means of so-called cographs. Moreover, we address the problem of finding "closest" symbolic ultrametrics and show the NP-completeness of the three problems: symbolic ultrametric editing, completion and deletion. Finally, as not all graphs are cographs, and hence, don't have a tree representation, we ask, furthermore, what is the minimum number of cotrees needed to represent the topology of an arbitrary non-cograph $G$. This is equivalent to find an optimal cograph edge $k$-decomposition $\{E_1,\dots,E_k\}$ of $E$ so that each subgraph $(V,E_i)$ of $G$ is a cograph. We investigate this problem in full detail, resulting in several new open problems, and NP-hardness results. For all optimization problems proven to be NP-hard we will provide integer linear program (ILP) formulations to efficiently solve them.

Posted Content
TL;DR: It is proved that for any finite simple graph G, the spectral functions F(G(m) of successive refinements converge for m to infinity uniformly on compact subsets of (0,1) and exponentially fast to a universal limiting eigenvalue distribution function F which only depends on the clique number respectively the dimension d of the largest complete subgraph of G and not on the starting graph G.
Abstract: Given a finite simple graph G, let G' be its barycentric refinement: it is the graph in which the vertices are the complete subgraphs of G and in which two such subgraphs are connected, if one is contained into the other. If L(0)=0 = dim(G). Let G(m) be the sequence of barycentric refinements of G=G(0). We prove that for any finite simple graph G, the spectral functions F(G(m)) of successive refinements converge for m to infinity uniformly on compact subsets of (0,1) and exponentially fast to a universal limiting eigenvalue distribution function F which only depends on the clique number respectively the dimension d of the largest complete subgraph of G and not on the starting graph G. In the case d=1, where we deal with graphs without triangles, the limiting distribution is the smooth function F(x) = 4 sin^2(pi x/2). This is related to the Julia set of the quadratic map T(z) = 4z-z^2 which has the one dimensional Julia set [0,4] and F satisfies T(F(k/n))=F(2k/n) as the Laplacians satisfy such a renormalization recursion. The spectral density in the d=1 case is then the arc-sin distribution which is the equilibrium measure on the Julia set. In higher dimensions, where the limiting function F still remains unidentified, F' appears to have a discrete or singular component.

Book ChapterDOI
TL;DR: This paper examines finite field trigonometry as a tool to construct trigonometric digital transforms by using properties of the k-cosine function over GF(p), and the Finite Field Discrete Cosine Transform (FFDCT) is introduced.
Abstract: This paper examines finite field trigonometry as a tool to construct trigonometric digital transforms. In particular, by using properties of the k-cosine function over GF(p), the Finite Field Discrete Cosine Transform (FFDCT) is introduced. The FFDCT pair in GF(p) is defined, having blocklengths that are divisors of (p+1)/2. A special case is the Mersenne FFDCT, defined when p is a Mersenne prime. In this instance blocklengths that are powers of two are possible and radix-2 fast algorithms can be used to compute the transform.

Posted Content
TL;DR: In this article, the Lovasz Local Lemma is used for the analysis of focused stochastic local search algorithms, which search a state space probabilistically by repeatedly selecting a constraint that is violated in the current state and moving to a random nearby state which, hopefully, addresses the violation without introducing new ones.
Abstract: We develop tools for analyzing focused stochastic local search algorithms. These are algorithms which search a state space probabilistically by repeatedly selecting a constraint that is violated in the current state and moving to a random nearby state which, hopefully, addresses the violation without introducing many new ones. A large class of such algorithms arise from the algorithmization of the Lovasz Local Lemma, a non-constructive tool for proving the existence of satisfying states. Here we give tools that provide a unified analysis of such algorithms and of many more, expressing them as instances of a general framework.

Posted Content
TL;DR: The main theorem says that the spanning trees in the convex combination can be chosen such that many of them are such "Gao trees" simultaneously at all sufficiently narrow cuts.
Abstract: We consider the $s$-$t$-path TSP: given a finite metric space with two elements $s$ and $t$, we look for a path from $s$ to $t$ that contains all the elements and has minimum total distance. We improve the approximation ratio for this problem from 1.599 to 1.566. Like previous algorithms, we solve the natural LP relaxation and represent an optimum solution $x^*$ as a convex combination of spanning trees. Gao showed that there exists a spanning tree in the support of $x^*$ that has only one edge in each narrow cut (i.e., each cut $C$ with $x^*(C)<2$). Our main theorem says that the spanning trees in the convex combination can be chosen such that many of them are such "Gao trees''.

Posted Content
TL;DR: A polynomial time algorithm is given which determines if a given graph containing a triangle and no induced seven-vertex path is 3-colorable, and gives an explicit coloring if one exists.
Abstract: In this paper, we give a polynomial time algorithm which determines if a given graph containing a triangle and no induced seven-vertex path is 3-colorable, and gives an explicit coloring if one exists. In previous work, we gave a polynomial time algorithm for three-coloring triangle-free graphs with no induced seven-vertex path. Combined, our work shows that three-coloring a graph with no induced seven-vertex path can be done in polynomial time.

Posted Content
TL;DR: A new approach to the computation of the hypervolume indicator, based on partitioning the dominated region into a set of axis-parallel hyperrectangles or boxes, is proposed and the nonincremental algorithm competes with the currently most practically efficient algorithms.
Abstract: We propose a new approach to the computation of the hypervolume indicator, based on partitioning the dominated region into a set of axis-parallel hyperrectangles or boxes. We present a nonincremental algorithm and an incremental algorithm, which allows insertions of points, whose time complexities are $O(n^{\lfloor \frac{p-1}{2} \rfloor+1})$ and $O(n^{\lfloor \frac{p}{2} \rfloor+1})$, respectively. While the theoretical complexity of such a method is lower bounded by the complexity of the partition, which is, in the worst-case, larger than the best upper bound on the complexity of the hypervolume computation, we show that it is practically efficient. In particular, the nonincremental algorithm competes with the currently most practically efficient algorithms. Finally, we prove an enhanced upper bound of $O(n^{p-1})$ and a lower bound of $\Omega (n^{\lfloor \frac{p}{2}\rfloor} \log n )$ for $p \geq 4$ on the worst-case complexity of the WFG algorithm.