Semidefinite Programming and Ramsey Numbers
Bernard Lidick\'y, Florian Pfender

TL;DR
This paper develops a novel technique using flag algebra to determine exact small Ramsey numbers, successfully establishing new bounds and exact values for several specific cases, challenging the notion that the method is only for large graphs.
Contribution
The paper introduces a new approach with flag algebras to find exact small Ramsey numbers, expanding the method's applicability beyond large graphs.
Findings
Established exact values for several small Ramsey numbers.
Developed a technique to adapt flag algebra for small graph problems.
Provided new upper bounds for multiple small graph and hypergraph Ramsey numbers.
Abstract
Finding exact Ramsey numbers is a problem typically restricted to relatively small graphs. The flag algebra method was developed to find asymptotic results for very large graphs, so it seems that the method is not suitable for finding small Ramsey numbers. But this intuition is wrong, and we will develop a technique to do just that in this paper. We find new upper bounds for many small graph and hypergraph Ramsey numbers. As a result, we prove the exact values , , , , , , and . We hope that this technique will be adapted to address other questions for smaller graphs with the flag algebra method.
| Previous bounds | Order | Graphs | Our upper bound |
| 8 | 1418 | ||
| 10 | 37133 | ||
| 10 | 38322 | ||
| 10 | 38440 | ||
| 10 | 38450 | ||
| 9 | 134037 | ||
| 8 | 11667 | ||
| 8 | 11765 | ||
| 8 | 11773 | ||
| 8 | 8722 | ||
| 8 | 18503 | ||
| 8 | 18601 | ||
| 8 | 9795 | ||
| 9 | 23398 | ||
| 9 | 23427 | ||
| 9 | 150078 | ||
| 8 | 11747 | ||
| 8 | 14889 | ||
| 8 | 15286 | ||
| 8 | 15311 | ||
| 8 | 14169 | ||
| 8 | 18186 | ||
| 8 | 18583 | ||
| 8 | 9478 | ||
| 8 | 19339 | ||
| 8 | 15206 | ||
| 8 | 19259 | ||
| 8 | 19656 | ||
| 9 | 23315 | ||
| 8 | 15304 | ||
| 8 | 19357 | ||
| 9 | 23419 | ||
| 9 | 23428 | ||
| 8 | 16649 | ||
| 8 | 14529 | ||
| 8 | 8836 | ||
| 8 | 14113 | (tight) | |
| 8 | 12327 | ||
| 8 | 17591 | ||
| 8 | 15471 | ||
| 8 | 18600 | ||
| 8 | 9778 | ||
| 8 | 9837 | ||
| 9 | 15067 | (tight) | |
| 9 | 15076 | ||
| 9 | 25482 | (tight) | |
| 9 | 49758 | ||
| 8 | 10114 | (tight) | |
| 8 | 10361 | ||
| 8 | 13780 | (tight) | |
| 8 | 8048 | (tight) | |
| 8 | 1398 | ||
| 8 | 11391 | ||
| 8 | 11748 | ||
| 8 | 15217 | ||
| 8 | 17547 | ||
| 8 | 9519 | ||
| 8 | 14456 | ||
| 8 | 9568 | ||
| 8 | 18543 | ||
| 8 | 12024 | ||
| 8 | 12122 | ||
| 8 | 15439 | ||
| 8 | 15591 | ||
| 9 | 116054 | (tight)111We provide the tight bound 13 in this paper, but it was not obtained by direct FA computation. | |
| 9 | 147411 | ||
| 7 | 120737 | ||
| 7 | 141516 | ||
| 6 | 15625 | ||
| 6 | 16272 | ||
| 6 | 18571 | ||
| 7 | 102305 | (tight) | |
| 7 | 7283 | ||
| 6 | 11193 | ||
| 6 | 120853 | ||
| 6 | 155664 | ||
| 6 | 88612 | ||
| 7 | 192287 | ||
| 6 | 9928 | ||
| 6 | 9386 | ||
| 6 | 170041 | ||
| 5 | 4418 | ||
| 6 | 56381 | ||
| 6 | 2589 | (tight) | |
| 7 | 145774 | ||
| 6 | 9476 | ||
| 6 | 11410 | ||
| 6 | 15170 | ||
| 6 | 12554 | ||
| 7 | 16169 | ||
| 7 | 5802 | ||
| 6 | 1345 | ||
| 8 | 1432 | (tight) | |
| 9 | 126456 | ||
| 8 | 5848 |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Semidefinite Programming and Ramsey Numbers
Bernard Lidický Department of Mathematics, Iowa State University, Ames, IA, E-mail: [email protected].
Florian Pfender Department of Mathematical and Statistical Sciences, University of Colorado Denver, E-mail: [email protected].
Abstract
Finding exact Ramsey numbers is a problem typically restricted to relatively small graphs. The flag algebra method was developed to find asymptotic results for very large graphs, so it seems that the method is not suitable for finding small Ramsey numbers. But this intuition is wrong, and we will develop a technique to do just that in this paper.
We find new upper bounds for many small graph and hypergraph Ramsey numbers. As a result, we prove the exact values , , , , , , and .
We hope that this technique will be adapted to address other questions for smaller graphs with the flag algebra method.
1 Introduction
Let be graphs. Ramsey’s celebrated Theorem [37] implies that for every edge coloring of a large enough complete graph with colors from , there exists some such that the contains a copy of with all edges colored . The Ramsey number is the smallest for which we are guaranteed to find such a monochromatic copy. A Ramsey graph is an extremal example for this number, i.e. a -edge-coloring of which does not contain a copy of in color for any .
The theory of flag algebras was developed by Razborov [38]. It has been used to find new results on graphs [3, 12, 13, 39], hypergraphs [2, 19, 22, 25], graphons [21], permutations [4], discrete geometry [20, 23], and even phylogenetic trees [1], to name a few.
The easiest and most popular usage of the theory is the plain flag algebra method. Formally, the method works with homomorphisms from linear combinations of combinatorial structures (graphs) to real numbers. The homomorphisms can be viewed as subgraph densities of (small) graphs in a very large graph, or more precisely, a graph limit.
The core of the plain method is to use the Cauchy-Schwarz Inequality to generate valid inequalities which hold for the subgraph densities of a large number of small graphs in the extremal graph (limit). Combinations of these inequalities are then used to produce the desired bounds. The right combination of the inequalities is usually found via semidefinite programming.
Finding exact Ramsey numbers is a problem typically restricted to relatively small graphs. The flag algebra method is designed to find asymptotic results for very large graphs, so it seems that the method is not suitable for finding small Ramsey numbers. But this intuition is wrong, and we will develop a technique to do just that in this paper. This technique may be adapted to address other questions for smaller graphs with the flag algebra method. So far, we have used variants of these ideas in [30] and [31].
We give a summary of new results in Section 2. We provide a very brief introduction to the theory of flag algebras in Section 3. We describe how to use the theory to obtain bounds on Ramsey numbers in Section 4. For better exposition, we describe the technique on a toy example proving that in Section 5, and expand on this example to a more general situation in Section 6. In the Appendix, we summarize the results of all the computations we tried.
The proofs involve extensive computations, and it is impractical to provide the actual solutions here. Even the certificates are impractically large to be provided as ancillary files. Instead, we provide the computer programs we used to obtain the results. This gives the interested reader the opportunity to recreate our results, and to try the methods on related questions. The programs and brief descriptions can be found in electronic form at http://lidicky.name/pub/ramsey and on arXiv https://arxiv.org/abs/1704.03592 as ancillary files.
2 Results
Here, we only present the new upper bounds we achieved together with the previously best known bounds referenced in a dynamical survey by Radziszowski [36]. Results presented here will be included in 2020 version of [36]. We use standard notation for all graphs and hypergraphs appearing here. In particular, stands for a complete (hyper)graph on vertices, minus one edge.
2.1 Graphs
We establish the following graph Ramsey numbers.
Theorem 1**.**
.
A Ramsey graph is the balanced complete -partite graph on vertices. Previously, the best upper bound was from [28].
Theorem 2**.**
.
A Ramsey graph is the balanced complete -partite graph on vertices. We are not aware of a previous non-trivial upper bound.
Theorem 3**.**
.
Here, stands for the graph of a -dimensional cube. Our flag algebra computations give an upper bound of , the previous lower bound was from [24]. In this case, the problem is small enough for a complete enumeration, and we found the exact number and all Ramsey graphs this way.
Theorem 4**.**
.
The flag algebra computation gives an upper bound for the order of a Ramsey graph barely above 16. Assuming this to be the correct bound, we examine the solution more closely. The flag algebra computation gives a list of graphs on 8 vertices that are unlikely to appear in a Ramsey graph on vertices, so we further assume that this graph does not contain any such subgraphs. This provides a significant restriction on the possible graphs on 9 or more vertices and we can enumerate all such graphs on up to 16 vertices. We find one Ramsey graph on 16 vertices this way, the Clebsch graph.
Theorem 5**.**
.
Previously, the best upper bound was by Piwakowski [35]. A Ramsey graph (which was not known to be Ramsey at the time) was constructed by Exoo [17].
Theorem 6**.**
.
Here, we improve the upper bound from to . The lower bound is by Tse [45].
We are able to improve the following bounds. Bounds without citations come from general theorems about Ramsey numbers. We denote the wheel on vertices by and a book on vertices by . That is, and .
Theorem 7**.**
New upper bounds on graph Ramsey numbers.
[TABLE]
Here we note that some of the upper bounds for wheels turned out to be tight later shown in [34].
Theorem 8**.**
New upper bounds on multi-color graph Ramsey numbers.
[TABLE]
2.2 3-uniform hypergraphs
In a couple cases, we are able to improve bounds on Ramsey numbers for -uniform hypergraphs.
Theorem 9**.**
* and .*
Both lower bounds are from [16], and we are not aware of a previous upper bound for the first quantity. The second quantity was previously bounded by .
We establish one new hypergraph Ramsey number.
Theorem 10**.**
.
To the best of our knowledge, this number has not been studied before. Using the computations for the upper bound similarly to the proof of Theorem 4, we construct the unique Ramsey -graph on vertices. This Ramsey -graph is highly symmetric, and we describe it here.
The -graph has edges, it is vertex and vertex-pair transitive with degree and co-degree . In fact, every vertex link (the -graph spanned by the edges incident to a vertex after deleting that vertex) is isomorphic to a -vertex Möbius ladder, i.e., with the antipodal chords added.
With vertex set , the edge set is
[TABLE]
In this case, the flag algebra computations result in a sharp bound. From this, we can use standard arguments to show that a large set of subgraphs (other than and the complement of ) can not occur in an -vertex Ramsey graph. The computer is then used to enumerate all such -graphs up to vertices, and finds that there is only one allowed -graph on vertices. Thus, in any Ramsey graph on vertices, all -vertex subgraphs must be isomorphic to this -graph. With this information, constructing is easy, either by hand or by computer.
2.3 Tournaments, directed graphs and further directions
Erdős and Moser [15] noted that Ramsey’s Theorem implies that for every , there exists a minimum number , such that every tournament on vertices contains a transitive tournament on vertices as a subtournament. The number is known for . Our method is applicable for this problem as well. We were able to improve the upper bound for .
Theorem 11**.**
.
The previous best upper bound was from [41]. As an immediate consequence, we can improve the upper bound for all larger , which was previously .
Corollary 12**.**
* for .*
It is also possible to use our method for Ramsey numbers of directed graphs in tournaments but we have not explored this direction. See the appendix for all bounds we have tried to improve, and for more information on the size of the computations.
3 Flag algebra terminology
Let us now introduce the terminology related to flag algebras needed in this paper. For more details about the method, see [38]. This section is included in order to make the paper self-contained. A reader familiar with the theory may wish to skip to the next section.
For a list , an edge colored graph is -free if it does not contain a copy of as a subgraph in color for any . Since we deal mostly with blow-ups of edge colored -free complete graphs, we restrict our attention to this particular case of Razborov’s much more general theory. We say that a graph is a blow-up of an edge colored complete graph if it can be obtained from an edge-colored complete graph by a blow-up of the vertices, i.e., vertices are replaced by independent sets, and edges are replaced by complete bipartite graphs between the sets, and all edges inherit the given color. For brevity, we will just write blow-up graph for these objects. The central notions we are going to introduce are an algebra and algebras , where is a fixed blow-up graph.
In order to precisely describe the algebras and , we first need to introduce some additional notation. Let be the set of all finite blow-up graphs. Next, for every , let be the set of blow-up graphs on exactly vertices. For and , we denote by the probability that a randomly chosen subset of vertices in induces a subgraph isomorphic to . Note that if . Let be the set of all formal linear combinations of elements of with real coefficients. Furthermore, let be the linear subspace of generated by all linear combinations of the form
[TABLE]
Finally, we define to be the space factorized by .
The space has naturally defined linear operations of addition and scalar multiplication by real numbers. To introduce a multiplication inside , we first define it on the elements of in the following way. For , and , we define to be the probability that a randomly chosen subset of of size and its complement induce in subgraphs isomorphic to and , respectively. We set
[TABLE]
The multiplication on has a unique linear extension to , which yields a well-defined multiplication also in the factor algebra . A formal proof of this can be found in [38, Lemma 2.4].
Let us now move to the definition of an algebra , where is an arbitrary blow-up graph with a fixed labelling of its vertex set. The labelled graph is usually called a type within the flag algebra framework. Without loss of generality, we will assume that the vertices of are labelled by . Now we follow almost the same lines as in the definition of . We define to be the set of all finite blow-up graphs with a fixed embedding of , i.e., an injective mapping from to such that induces in a subgraph isomorphic to . The elements of are usually called -flags and the subgraph induced by is called the root of a -flag.
Again, for every , we define to be the set of the -flags from that have size (i.e., the -flags with the underlying blow-up graph having vertices). Analogously to the case for , for two blow-up graphs with the embeddings of given by , we set to be the probability that a randomly chosen subset of vertices in together with induces a subgraph that is isomorphic to through an isomorphism that preserves the embedding of . In other words, the isomorphism has to satisfy . Let be the set of all formal linear combinations of elements of with real coefficients, and let be the linear subspace of generated by all the linear combinations of the form
[TABLE]
We define to be factorized by .
We now describe the multiplication of two elements from . Let , , and be the fixed embedding of in . As in the definition of multiplication for , we define to be the probability that a randomly chosen subset of of order and its complement in of order , extend in to subgraphs isomorphic to and , respectively. This definition naturally extends to .
Now consider an infinite sequence of blow-up graphs of increasing orders. We say that the sequence is convergent if the probability has a limit for every . A standard compactness argument (e.g., using Tychonoff’s theorem) yields that every infinite sequence of blow-up graphs has a convergent subsequence. All the following results can be found in [38]. Fix a convergent increasing sequence of blow-up graphs. For every , we set and linearly extend to . We usually refer to the mapping as the limit of the sequence. The obtained mapping is a homomorphism from to . Moreover, for every , we obtain . Let be the set of all such homomorphisms, i.e., the set of all homomorphisms from the algebra to such that for every . It is an interesting result that this set is exactly the set of all limits of convergent sequences of blow-up graphs [38, Theorem 3.3].
Let be a convergent sequence of blow-up graphs and be its limit. For and an embedding of in , we define to be the blow-up graph rooted on the copy of that corresponds to . For every and , we define . Picking at random gives rise to a probability distribution on mappings from to , for every . Since converges (as tends to infinity) for every , the sequence of these probability distributions on mappings from to also converges [38, Theorems 3.12 and 3.13]. We denote the limit probability distribution by . In fact, for any such that , the homomorphism itself fully determines the random distribution [38, Theorem 3.5]. Furthermore, any mapping from the support of the distribution is in fact a homomorphism from to such that for all [38, Proof of Theorem 3.5].
The last notion we introduce is the averaging (or downward) operator . It is a linear operator defined on the elements of by , where is the (unlabeled) blow-up graph from corresponding to , and is the probability that a random injective mapping from to is an embedding of in yielding a -flag isomorphic to . The key relation between and is the following:
[TABLE]
where the integration is over the probability space given by the random distribution on . Therefore, if almost surely for some , then \phi\left(\left\llbracket{A^{\sigma}}\right\rrbracket_{\sigma}\right)\geq 0. In particular,
[TABLE]
The plain method is a tool from the flag algebra framework that, for a given density problem of the form
[TABLE]
where , systematically searches for ‘best possible’ inequalities of the form (2). If we fix in advance an upper bound on the size of graphs in the terms of inequalities we will be using, we can find the best inequalities of the form (2) using semidefinite programming.
To reduce the size of and with it the size of all required computations, it is often benefitial to use a partially color-blind setting. In this setting, the colors are partitioned into classes, and two blow-up graphs are considered to be the same if they differ only by a permutation of colors inside the classes. All of the theory described in this chapter naturally works for this setting as well.
4 Using flag algebra to bound Ramsey numbers
For some , start with a -free -edge-coloring of a . Now replace every vertex by a large independent set of size , say. If this blow-up graph contains a copy of in color , then two of the vertices in this copy are in the same independent set. Making larger and larger, this graph sequence becomes an object that can be analyzed by the plain flag algebra method.
Formally, we consider the model of blow-ups of -edge-colored complete graphs, for which every copy of in color contains at least two vertices in the same independent set. This model can easily be described in the theory of flag algebras. For readers familiar with the language of graph limits, we look at the -colored graphon of , i.e. a step function , where every contains exactly one for off-diagonal steps, and all [math]s for the diagonal steps.
In this model, we find a lower bound for the density of non-edges via the plain flag algebra method. The minimum is achieved exactly by a balanced blow-up of any Ramsey graph. Therefore, if is a lower bound for the density of non-edges, then
[TABLE]
More generally, we can look at lower bounds for the density of independent sets of size . Again, the minimum is achieved exactly by a balanced blow-up of any Ramsey graph on vertices, and it follows that
[TABLE]
Notice that we can make use of the integrality of . If we want to show that for some , all we need to show is that for some . In most cases, we found the same bounds by using different , but in some cases, the bounds were different.
The application of the plain flag algebra method requires the enumeration of all small graphs in the model. A computer is used to enumerate the small graphs of the prescribed order, set up the inequalities and then solve the resulting semidefinite program. The process of setting up the semidefinite program and processing the solution of the program is by now standard in the community and it is briefly described in the next two sections. It has been automated for graphs, 3-graphs and oriented graphs by the software package Flagmatic [46]. While we can not use Flagmatic in our specific application to blow-up graphs, our computations follow the same lines.
The semidefinite program can be solved by state of the art solvers CSDP [5] and SDPA [48]. These solvers use floating point arithmetic, and in most applications of the plain flag algebra method the following rounding step requires some thought, sometimes ingenuity, to turn the results into a proof. In our application, though, we are usually not interested in sharp bounds as we can use the integrality of , and the rounding is easy. Round the result to a desired level of precision, while keeping the resulting matrix positive semidefinite. Due to continuity, the resulting bounds are almost unchanged. We end up with a certificate consisting of several (sometimes very large) rational positive semidefinite matrices.
5 Illustration of the method:
In this section we illustrate our method on the smallest non trivial Ramsey number . This may be the most complicated proof of this fact ever published. In fact, at an early point in the proof we determine all -colorings of without monochromatic triangles, from which it is easy to find the unique Ramsey graph on five vertices. For larger Ramsey numbers a similar complete enumeration is not feasible, and our method, which only uses relatively small graphs, can find new upper bounds.
Recall that there is a 2-edge-coloring of without monochromatic triangles, see Figure 2, so all we need to show is that .
Proof of .
Let , and suppose that is a 2-edge-colored with no monochromatic triangle. Let be a blow-up of on vertices where every vertex of is replaced by an independent set of size for . Clearly, . The number of non-edges in is . This is minimized if for all . Hence, the number of non-edges is at least , which gives an asymptotic density of non-edges of at least .
Denote by the minimum asymptotic density of non-edges over all 2-colored blow-up graphs with no monochromatic triangles. Therefore, and hence . In order to prove that the largest graph with no monochromatic triangles has at most vertices, it is enough to show that . If there was a complete graph on vertices with no monochromatic triangle, then there would be a blow-up graph with .
We work in : the class of -colored blow-up graphs with no monochromatic triangles. In figures, we will use solid and dotted lines to distinguish the two colors. We use the color-blind setting, so for example,
is considered to be the same graph as
.
Forbidden subgraphs in are monochromatic triangles
(and
, but this already follows from color-blindness). Since all graphs in are blow-up graphs, triples inducing exactly one edge
and triples inducing exactly two edges with different colors
are also forbidden subgraphs.
This leaves exactly seven graphs on vertices in , taking color-blindness into account:
[TABLE]
With a slight abuse of notation, we use the drawing of a graph also for the asymptotic density , making our equalities and inequalities much more intuitive. As a first equality, we have in :
[TABLE]
We use one type of size two, the edge . We use flags of size three. In the figures, we use a gray square and a white square to distinguish the two labeled vertices. We have three flags for in a vector
[TABLE]
Using (1), we have
[TABLE]
and we want to show that \vbox{\hbox{\includegraphics[page=17]{fig-R33}}}>\frac{1}{6}.
Let be the following matrix, which is the rounded solution of a suitable semidefinite program which we will describe later.
[TABLE]
The matrix is positive semidefinite, the smallest eigenvalue is greater than , and thus
[TABLE]
We explicitely compute the right side of (5). Here is an example for the required computations.
[TABLE]
Performing similar computations for all required multiplications, we get the following table, in which we omitted all zeros and multiplied all entries by to avoid fractions.
[TABLE]
This gives
[TABLE]
We subtract the result from (4) and obtain
[TABLE]
∎
Now we give an explanation on how to formulate a semidefinite program to find . First we expand (5).
[TABLE]
We combine this with (4) and use (3) to obtain the following.
[TABLE]
To find the best lower bound t\leq~{}\vbox{\hbox{\includegraphics[page=17]{fig-R33}}}~{}, we formulate this as a semidefinite program.
[TABLE]
This program can be solved with CSDP, which provides a numerical approximate solution
[TABLE]
Now we need to round to a rational matrix we can use in the proof. The numerical eigenvalues of are
[TABLE]
If we round entry-wise, this can easily result in a matrix that is no longer positive semidefinite since has an eigenvalue that is very close to zero. In the proof above, we rounded the entries by hand and verified that we get positive eigenvalues. In other instances, in particular for larger , we use an automated approach. We numerically diagonalize , and round the eigenvectors and eigenvalues to compute the rounded matrix . Numerically, we get , where is a diagonal matrix with the eigenvalues of and
[TABLE]
If we round both and to two decimal places we obtain
[TABLE]
which is positive definite by its construction, and the calculation can be performed with rational numbers without introducing any numerical errors. From this, we recompute the value of . With this approximation, we get . We can get closer to the numerical solution if we round to more decimal places. But since the last set of inequalities in the proof of is far from sharp, this is not necessary.
6 More details on using the method
In this section, we will run through a more general example to show how to adapt the illustration from the previous section. Our goal is to show that for two graphs and a bound .
Again, is a -colored , and contains neither a red or a blue . If we replace the vertices of by large independent sets, then the density of non-edges in the resulting blow-up graph is at least . We are computing in the model of red/blue-colored graphs with the following forbidden subgraphs:
any complete graph with a red , 2. 2.
any complete graph with a blue , 3. 3.
,
,
.
Here, 1 and 2 are from the particular Ramsey problem, and 3 ensures that we are considering blow-up graphs only. The manual input for the next step are the adjacency matrices of these graphs. As standard in the plain flag algebra method, the computer then generates all graphs in this model up to a given order , where is part of the input. Typically, we want to pick as large as possible while keeping the total number of generated graphs under 200,000. Larger numbers of graphs greatly slow down the subsequent computations and increase the memory demands. A personal computer can handle about 10,000 graphs. The main bottleneck in computation and memory is later in the process when solving the semidefinite program, so we are not very worried about implementation here. Notice that the number of constraints in the semidefinite program is about the number of graphs. Any reasonable implementation of generating graphs of order from graphs of order while testing for duplicates is sufficient. We could likely achieve a small speed-up here by using nauty or a similar software. Next, the computer generates all possible types of flags so that products of any two flags of the same type have order . By the definition of the product, we restrict ourselves to types whose order has the same parity as . This allows us to gather these flags in vectors of flag densities such that for every positive semidefinite matrix .
We then formulate a semidefinite program which optimizes the matrices used to give us the best possible lower bound of the density of non-edges, which then in turn gives an upper bound for . The formulation of the semidefinite program is performed by the computer, along the lines in the illustration in the previous section, mainly computing the coefficients from . Up to this point, the computer code was developed by ourselves.
For solving the semidefinite program, we use CSDP [5], a numerical solver using floating point arithmetic. As mentioned above, this is the main bottleneck in computation and memory. We then use yet another self developed computer program to change the numerical solution to an exact one. We take the numerical positive semidefinite matrix, find numerical eigenvectors and eigenvalues, round them to rational numbers, and then reconstruct a rational positive semidefinite matrix. This process puts out a positive semidefinite matrix which can be used as a certificate for the correctness of the bound. The rounding may slightly worsen the resulting bound on . But in the end, due to the discrete nature of the problem, this small difference typically does not change the result.
More colors or hypergraphs can easily be encoded in a similar model, at the cost of many more graphs on the same number of vertices. As mentioned above, and exhibited in the previous section, if for some , we can reduce the number of graphs through a color-blind process, in which the two colors and may be permuted.
Acknowledgements
The research of the first author was supported in part by NSF grants DMS-1600390 and DMS-1855653 and the research of the second author was supported in part by NSF grants DMS-1600483 and DMS-1855622.
The research reported in this paper is partially supported by the HPC@ISU equipment at Iowa State University, some of which has been purchased through funding provided by NSF under MRI grant number CNS 1229081 and CRI grant number 1205413.
This work used computer clusters at the Center for Computational Mathematics, University of Colorado Denver, purchased thanks to NSF Grants AGS 0835579 (Gross) and CNS 0958354 (Colibri).
This work utilized the Janus supercomputer, which is supported by the NSF (award number CNS-0821794) and the University of Colorado Boulder. The Janus supercomputer is a joint effort of the University of Colorado Boulder, the University of Colorado Denver and the National Center for Atmospheric Research.
We are also grateful to the Charles University in Prague for providing access to the computers Kaminka, and Kamenozrout.
The first author gratefully acknowledges the support and hospitality and access to computing resources provided by the IMA during his visit which took place from Sep, 2014 to Dec, 2014.
The authors would like to thank Stanisław Radziszowski for providing valuable comments about the paper.
Appendix A All attempted bounds
The following table provides a summary of computations we performed, where shaded rows correspond to improved upper bounds. The purpose of the table is to illustrate the size of the computations, and to also show our attempts where the method provided an upper bound that did not improve on the best known one. The basic parameter of computations is the order of graphs in . A bigger value of typically gives a better result. On the other hand, the number of graphs of order grows quickly and becomes unmanageable soon. For every computation, we list both and the number of graphs in . We still have some calculations in progress. We plan to update the arXiv preprint at https://arxiv.org/abs/1704.03592 when they finish.
One of the main issues is the memory needed by CSDP when solving the semidefinite program. The memory demands grow quickly with the number of graphs in . If the number of graphs is around 10,000, the instance is solvable on a desktop. Numbers under 100,000 will fit in 128G of memory, which requires a supercomputer. Numbers above 100,000 require high memory super computers. All instances we have tried fit in about 300G of memory. Even larger instances than we tried could be solvable as very high memory nodes may have even terabytes of memory, but one would have to be very patient.
The running time also depends heavily on the number of graphs. The CSDP solver runs in iterations and it took 30 to 60 iterations to solve most of the problems in this class. The larger instances compute a few iterations per day to a few days per iteration on the supercomputers we use. Let us mention that we obtained a significant speedup (10) of the CSDP solver by compiling it with Intel Math Kernel Library.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] N. Alon, H. Naves, and B. Sudakov. On the maximum quartet distance between phylogenetic trees. SIAM J. Discrete Math. , 30(2):718–735, 2016.
- 2[2] R. Baber and J. Talbot. Hypergraphs do jump. Combin. Probab. Comput. , 20(2):161–171, 2011.
- 3[3] J. Balogh, P. Hu, B. Lidický, and F. Pfender. Maximum density of induced 5-cycle is achieved by an iterated blow-up of 5-cycle. European J. Combin. , 52(part A):47–58, 2016.
- 4[4] J. Balogh, P. Hu, B. Lidický, O. Pikhurko, B. Udvari, and J. Volec. Minimum number of monotone subsequences of length 4 in permutations. Combin. Probab. Comput. , 24(4):658–679, 2015.
- 5[5] B. Borchers. CSDP, A C library for semidefinite programming. Optimization Methods and Software , 11(1-4):613–623, 1999.
- 6[6] L. Boza. Personal communication. 2012.
- 7[7] L. Boza. Sobre el. número de Ramsey R ( K 6 − e , K 4 ) 𝑅 subscript 𝐾 6 𝑒 subscript 𝐾 4 R(K_{6}-e,K_{4}) , 2013. VIII Encuentro Andaluz de Matemática Discreta, Sevilla, Spain.
- 8[8] L. Boza. Sobre los números de Ramsey R ( K 5 − e , K 5 ) 𝑅 subscript 𝐾 5 𝑒 subscript 𝐾 5 R(K_{5}-e,K_{5}) y R ( K 6 − e , K 4 ) 𝑅 subscript 𝐾 6 𝑒 subscript 𝐾 4 R(K_{6}-e,K_{4}) , 2014. IX Jornada de Matemática Discreta y Algoritmica, JMDA 2014, Tarragona, Spain.
