Constructions of optimal orthogonal arrays with repeated rows
Charles J. Colbourn, Douglas R. Stinson, Shannon Veitch

TL;DR
This paper presents new methods for constructing optimal orthogonal arrays with repeated rows, maximizing the ratio of repeated rows to the array's index, and explores the existence of basic OAs linked to Hadamard matrices.
Contribution
It introduces constructions for optimal orthogonal arrays with large repeated row ratios and establishes a link between basic OAs with n=2 and Hadamard matrices.
Findings
Constructed optimal OAs for any k ≥ n+1 with large λ
Developed basic OAs with n=2 and k=4t+1 based on Hadamard matrices
Solved the problem of constructing basic OAs with n=2 modulo Hadamard conjecture
Abstract
We construct orthogonal arrays OA (of strength two) having a row that is repeated times, where is as large as possible. In particular, we consider OAs where the ratio is as large as possible; these OAs are termed optimal. We provide constructions of optimal OAs for any , albeit with large . We also study basic OAs; these are optimal OAs in which . We construct a basic OA with and , provided that a Hadamard matrix of order exists. This completely solves the problem of constructing basic OAs wth , modulo the Hadamard matrix conjecture.
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.
Taxonomy
Topicsgraph theory and CDMA systems · Optimal Experimental Design Methods · Optimization and Packing Problems
Constructions of optimal orthogonal arrays with repeated rows
Charles J. Colbourn
School of Computing, Informatics and Decision Systems Engineering, Arizona State University, Tempe, Arizona 85287, U.S.A.
Douglas R. Stinson D.R. Stinson’s research is supported by NSERC discovery grant RGPIN-03882. David R. Cheriton School of Computer Science, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
Shannon Veitch
Dept. of Combinatorics and Optimization, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
Abstract
We construct orthogonal arrays (of strength two) having a row that is repeated times, where is as large as possible. In particular, we consider OAs where the ratio is as large as possible; these OAs are termed optimal. We provide constructions of optimal OAs for any , albeit with large . We also study basic OAs; these are optimal OAs in which . We construct a basic OA with and , provided that a Hadamard matrix of order exists. This completely solves the problem of constructing basic OAs wth , modulo the Hadamard matrix conjecture.
1 Introduction
Let , and be integers. An orthogonal array is a by array, , with entries from a set of cardinality such that, within any two columns of , every ordered pair of symbols from occurs in exactly rows of . For much information on orthogonal arrays, see [2].
In this note, we are interested in that contain a row that is repeated times, where is as large as possible. We observe that, by relabelling the symbols in the orthogonal array, these rows can all be assumed to be of the form for some fixed symbol . We will denote this particular symbol either by [math] or by .
The following theorem, along with an elementary combinatorial proof, can be found in [6, Theorem 2.2]. However, we should note that this result is also an immediate consequence of a much more general result from [3].
Theorem 1.1**.**
Let , and be integers. If there is an containing a row that is repeated times, then
[TABLE]
An , say , containing a row that is repeated
[TABLE]
times will be termed optimal. Another way to view the optimality property is to observe that the ratio is as large as possible in an optimal OA.
We note that, in a recent paper, Culus and Toulouse [1] discuss an application where it is beneficial to construct optimal orthogonal arrays. They also construct several small examples of optimal OAs using linear programs.
The rest of this paper is organized as follows. Section 2 establishes some basic results about optimal OAs. In Section 2.1, we present some small examples of basic OAs, which are optimal OAs in which . In Section 3, we give a complete solution (modulo the Hadamard matrix conjecture) to the the problem of constructing basic OAs with . Section 4 gives constructions of optimal OAs for arbitrary values of and . Section 5 examines the effect of deleting a small number of columns from an optimal OA. Finally, Section 6 is a brief summary.
2 Preliminary results
From the proof of [6, Theorem 2.2], the following result can be derived immediately.
Lemma 2.1**.**
[6, Corollary 2.4]** Let , and be integers. Suppose there is an optimal , say , containing a row that is repeated times. Then every other row of contains exactly
[TABLE]
occurrences of the symbol [math].
We can simplify the formula (2) for , as follows.
Corollary 2.2**.**
Let , and be integers. Suppose there is an optimal , say , containing a row that is repeated times. Then every other row of contains exactly occurrences of the symbol [math] and thus .
Proof.
In an optimal , equation (1) holds. Suppose we substitute this expression for into equation (2). We obtain
[TABLE]
∎
We will define a quadruple of positive integers to be feasible if the following conditions are satisfied:
- •
and ,
- •
, and
- •
is a positive integer.
The above conditions are all necessary for the existence of an optimal .
Lemma 2.3**.**
Suppose is a feasible quadruple, and , where . Then is a feasible quadruple.
Proof.
Define and . It suffices to observe that
[TABLE]
∎
If there exists an optimal and we take copies of every row, then we obtain an optimal , where . The most interesting parameter cases are those where we cannot just take multiple copies of a smaller OA. Therefore, we define a feasible quadruple to be basic if . Similarly, an optimal is basic if .
As we already mentioned, an optimal OA is any that has an -times repeated row, where the ratio is as large as possible. A basic OA is an optimal OA where the value of (or equivalently, ) is as small as possible. Note that basic OAs include the case , which of course do not contain repeated rows.
The possible basic quadruples are quite constrained if is prime and .
Theorem 2.4**.**
If is a basic quadruple, and is prime, then . Further, and for some integer such that .
Proof.
We have that from (1). Further, for some integer , by Corollary 2.2. Therefore
[TABLE]
and hence
[TABLE]
Then and , so . Since is prime and , we have .
Referring again to (1), it follows that
[TABLE]
Finally, a basic quadruple requires that . Thus, it is necessary that , which simplifies to . ∎
2.1 Some examples of basic OAs
An optimal has a total of rows. If we delete rows of the form , then the number of remaining rows is , which is divisible by . This suggests that we might attempt to construct the by cyclically rotating “starting rows” times.
We illustrate the technique in a small example.
Example 2.1**.**
We construct a basic from the following two starting rows:
[TABLE]
We cyclically rotate these starting rows five times, and then adjoin rows of [math]’s. This yields the desired orthogonal array, which is presented in Figure 1.
It is possible to verify that this process will yield an without actually constructing the whole array. It suffices to look all the ordered pairs that are (cyclically) a fixed distance apart in the starting rows. Each row has five entries, so we only have to consider pairs at distance one and two, because .
- •
For distance one, we have from the first starting row and from the second starting row. We see we have three occurrences of each of and , and one occurrence of .
- •
For distance two, we have from the first starting row and from the second starting row. Again, we have three occurrences of each of and , and one occurrence of .
This means that, when we rotate the starting rows and adjoin two rows of [math]’s, we are guaranteed to get the desired orthogonal array.
Finally, we note that existence of a basic is also reported in [1, Table 4].
We now give starting rows for a few other small examples.
Example 2.2**.**
A basic with can be constructed from the following two starting rows:
[TABLE]
Existence of a basic is also reported in [1, Table 4].
Example 2.3**.**
A basic with can be constructed from the following two starting rows:
[TABLE]
Example 2.4**.**
A basic with can be constructed from the following two starting rows:
[TABLE]
Example 2.5**.**
We use a slightly different technique to obtain a basic with . We have three starting rows, consisting of symbols from the set :
[TABLE]
First, cyclically rotate each starting row seven times. Then develop each row modulo (the point is fixed). Finally, adjoin three rows of ’s. The resulting rows form a basic .
Existence of a basic is also reported in [1, Table 4].
Example 2.6**.**
We construct a basic with from four starting rows, consisting of symbols from the set :
[TABLE]
First, cyclically rotate each starting row nine times. Then develop each row modulo (the point is fixed). Finally, adjoin four rows of ’s. The resulting rows form a basic .
3 Basic OAs for
In this section, we determine all the basic quadruples with . For these quadruples, we can construct a basic provided that a suitable Hadamard matrix exists.
Lemma 3.1**.**
If is a basic quadruple with , then . Further, and for some positive integer .
Proof.
Take in Theorem 2.4. Then and we have and , where . Hence is even. Writing , we obtain and . ∎
Theorem 3.2**.**
There exists a basic if and only if there is a -BIBD.
Proof.
First suppose that a -BIBD exists. This BIBD has blocks and replication number . Let be the by incidence matrix of this BIBD. Construct the matrix
[TABLE]
We claim that is a basic . Clearly the first two rows are identical, so and we just need to verify that is an OA with the stated parameters.
Choose any two distinct columns of . These columns correspond to two points in the BIBD, say and , where . The number of occurrences of in these two columns is . The number of occurrences of is , as is the number of occurrences of . Finally, the number of occurrences of is .
Conversely, suppose is a basic . We can assume that the symbols in the OA are [math] and . Without loss of generality, suppose that there are rows consisting entirely of zeroes, and then delete them, creating a by matrix . We will show that is the incidence matrix of a -BIBD.
Clearly is the incidence matrix of a set system on points. Given any two points and , the number of blocks containing these two points is the same as the number of occurrences of in the two associated columns of the OA, which is . To complete the proof, we show that every block in the set system has size . This can be seen easily by recalling from Lemma 2.1 that every row of contains exactly zeroes. Therefore, the number of ones in a row of is . This completes the proof. ∎
Now we show how to construct a basic from a certain Hadamard matrix. We use a few standard results, all of which can be found in [5], for example.
Theorem 3.3**.**
If there exists a Hadamard matrix of order , then there exists a basic .
Proof.
It is well-known that a Hadamard matrix of order is equivalent to a symmetric -BIBD. The derived BIBD is a -BIBD. If we then complement every block in this BIBD, we obtain a -BIBD. Finally, apply Theorem 3.2. ∎
It is known that Hadamard matrices exist for all orders , , and it is conjectured that Hadamard matrices exist for all orders , .
4 General constructions for optimal OAs
Suppose we fix and , where . Denote ; then . Our goal is to find an optimal for some value of . Note that in such an OA. Also, .
One possible approach would be to take all possible -tuples that contain precisely occurrences of [math], and then adjoin an appropriate number of rows consisting entirely of [math]’s.
We illustrate the idea using a small example.
Example 4.1**.**
Suppose we take and . Here we have and . There are
[TABLE]
-tuples on the symbol set that contain precisely two zeroes. If we then adjoin rows of [math]’s, it is not hard to check that we obtain an optimal . The ratio , as required.
Of course, we know from Example 2.5 that a basic exists. The value of in the above-constructed optimal is much larger.
As in Example 4.1, we take all possible -tuples that contain precisely occurrences of 0. In order to have an orthogonal array, any two columns must contain every ordered pair of symbols exactly times, for some . Since every possible -tuple containing [math]’s is used, it suffices to consider the first two columns. If these two columns contain every ordered pair the same number of times, then so will every other pair of columns in the array.
There are four cases to consider for the first two elements in a row: , , , and , where . We consider these cases in turn.
For each , there are
[TABLE]
rows beginning with . This result is the same for rows beginning with . For each , there are
[TABLE]
rows beginning with . Since , we find that
[TABLE]
Using (5), it is now easy to show that (3) and (4) are equal:
[TABLE]
Therefore, each ordered pair other than appears the same number of times (given by (3) or (4)), which we denote by . Thus, adjoining the appropriate number of rows consisting entirely of [math]’s will result in an orthogonal array.
The number of rows beginning with is
[TABLE]
It then follows that we need to adjoin
[TABLE]
rows of [math]’s. Substituting the value of obtained from (3), we have
[TABLE]
Therefore, substituting , we have
[TABLE]
as desired.
Thus we have proven the following result.
Theorem 4.1**.**
Suppose and suppose is an integer. Then there is an optimal , where
[TABLE]
4.1 An improvement
We next show that the constructed in Theorem 4.1 can be partitioned into optimal . In order to describe how this is done, it is useful to change the set of symbols on which the OAs are defined. Suppose we begin with the above-mentioned , constructed on symbols . Delete the rows of [math]’s. Then we replace all occurrences of [math] in the remaining rows by , and the symbols are replaced by , respectively. We consider the new symbol set as . Denote the resulting array by .
For any row of , let denote the sum of the non-infinite elements in row , reduced modulo . Then, for any , let consist of all the rows of such that . Clearly every row of is in precisely one of .
It is obvious that every is fixed by any permutation of the columns. Therefore, the number of occurrences of a particular pair of symbols in two given columns does not depend on the two columns that are chosen. Hence, we can restrict our attention to the first two columns of these arrays.
For two (not necessarily distinct) symbols and for , let denote the number of occurrences of the ordered pair in the first two columns of . Also, let denote the number of occurrences of in the first two columns of . Therefore
[TABLE]
where
[TABLE]
and
[TABLE]
We will show, for any , that is independent of . First however, we state and prove a small technical lemma.
Lemma 4.2**.**
Suppose that , is an integer and . Then unless , and .
Proof.
First, suppose , so . The following inequalities are equivalent:
[TABLE]
However,
[TABLE]
because .
Now, suppose and , We have , so it follows that . ∎
Remark 4.1**.**
The exception to Lemma 4.2 is , and . But this is a trivial case, as a basic exists, and this OA can trivially be “partitioned” into OAs.
In the following discussion, we assume that , is an integer, , and .
We next define a mapping on the rows of . Let be any row of , where . Let (there are occurrences of in , so Lemma 4.2 ensures that and hence exists).
Let . Define , where
[TABLE]
for .
The process above can also be described as follows. Find the first entry in row , past the second column, that is not equal to . Then add modulo to that entry.
It is clear that the mapping gives a bijection from the rows in to the rows in , for . Also, leaves the points in the first two columns of any row of unaltered. Since
[TABLE]
we have for all . Therefore, if we adjoin rows of ’s to any , we obtain an optimal , where
[TABLE]
The above discussion, along with Remark 4.1, proves the following.
Theorem 4.3**.**
Suppose and suppose is an integer. Then there is an optimal , where
[TABLE]
that can be partitioned into optimal .
Corollary 4.4**.**
Suppose and suppose is an integer. Then there is an optimal , where
[TABLE]
Remark 4.2**.**
We have noted that when . A basic OA with these parameters is in fact an , which is equivalent to a projective plane of order . However, in cases where a projective plane of order is known not to exist, Corollary 4.4 provides examples of for certain large values of .
4.2 A further improvement
We now prove an extension of Theorem 4.3 where we can sometimes reduce the value of by additional factors of , depending on the parameters and .
Suppose we can write , where are integers such that for . Evidently, we can take
[TABLE]
As before, our starting point is the optimal obtained from Theorem 4.1 in which the symbol set is . The value of is given by (6) and there are rows consisting only of ’s. Delete these rows and call the resulting array .
Now, we consider the columns of to be partitioned into classes, where the th class consists of columns. Denote these classes as .
For any row of , let , where, for , is the modulo sum of the non-infinite elements in row in the columns in . We will say that the -tuple is the type of row .
Then, for any possible type , let consist of all the rows of such that . (That is, comprises all the rows of having type .) This yields a partition of the rows of into subsets. We will show that each , when augmented with an appropriate number of rows of ’s, is an optimal , where
[TABLE]
The number of rows of ’s to be added to each to construct an orthogonal array is .
We use a “bijection” proof similar to Theorem 4.3. First, we note that any permutation of the columns within any column class is an automorphism of every . Therefore, to prove that yields an orthogonal array (as described above), we just need to consider the ordered pairs of symbols occurring in pairs of columns of the following types:
(a)
the first two columns of any , and
(b)
the first column of and the the first column of , where .
We can use a bijection similar to Theorem 4.3, but we apply the bijection independently for every column class . Let be any row of , where . For each column class , let be the first column within , past the second column, that is not equal to (note that exists because ).
Fix any -tuple . Define the mapping as follows: , where
[TABLE]
for .
The mapping is a bijection from the rows in to the rows in , for any . Also, for any row , leaves the points in the first two columns of every unaltered. Therefore, the ordered pairs occurring in pairs of columns of types (a) and (b) in any is independent of . Hence, we obtain the following theorem and corollary.
Theorem 4.5**.**
Suppose and suppose is an integer. Suppose
[TABLE]
Then there is an optimal , where
[TABLE]
that can be partitioned into optimal .
Corollary 4.6**.**
Suppose and suppose is an integer. Suppose
[TABLE]
Then there is an optimal , where
[TABLE]
Example 4.2**.**
Suppose we take and . Then and . From Corollary 4.6, we obtain an optimal , where
[TABLE]
5 Deleting columns from optimal OAs
An optimal can exist only when . However, it is certainly of interest to determine the maximum possible ratio for values of and where . We might reasonably expect that deleting a small number of columns, say columns, from an optimal could yield an where the ratio is as large as possible.
5.1 Deleting a single column
We first prove that this approach works for by proving a modification of Theorem 1.1. The modified bound uses a similar proof technique to [4, Theorem 3.1] (see also [5, Theorem 8.7]).
Theorem 5.1**.**
Let , , be integers. Suppose is an containing a row that is repeated times, and let be a positive integer. Then
[TABLE]
Further, equality occurs if and only if every row has either or occurrences of the symbol 0.
Proof.
Suppose the last rows of are . Let denote the number of occurrences of the symbol [math] in row of . Define . It is clear from elementary counting that
[TABLE]
Let be a positive integer. Then
[TABLE]
Solving for , we obtain (7).
In the case of equality, every term in the sum
[TABLE]
must equal [math]. Therefore, for every . ∎
We examine a special case of the bound proven in Theorem 5.1.
Corollary 5.2**.**
Let , , be integers. If and there is an containing a row that is repeated times, then
[TABLE]
Proof.
Write for some integer and take in Theorem 5.1. We find that
[TABLE]
∎
Theorem 5.3**.**
Suppose there is an optimal containing a row that is repeated times. Then and there is an containing a row that is repeated times, where (8) is met with equality.
Proof.
Let be an optimal containing a row that is repeated times. If we delete any one column from , then the resulting array is an containing a row that is repeated times. Since is optimal, we have
[TABLE]
from (1), so (8) is met with equality. Further, must be a positive integer from Corollary 2.2. ∎
5.2 Deleting multiple columns
We now consider a different approach. Suppose we return to our original bound, Theorem 1.1. Since must be an integer, the following variation is immediate.
Theorem 5.4**.**
Let , and be integers. If there is an containing a row that is repeated times, then
[TABLE]
An orthogonal array for which the bound in Theorem 5.4 is met with equality will be termed -optimal.
Assume we start with an optimal , so
[TABLE]
We determine how many columns we can remove, denoted by , such that
[TABLE]
The resulting will be -optimal.
The following numerical lemma will be useful.
Lemma 5.5**.**
Suppose that
[TABLE]
is an integer and
[TABLE]
is a positive integer. Then
[TABLE]
Proof.
From (9), we find that
[TABLE]
Rearranging this inequality, we see that
[TABLE]
and therefore
[TABLE]
Since is an integer, this proves that
[TABLE]
∎
The following theorem is now an immediate consequence of Lemma 5.5.
Theorem 5.6**.**
If there is an optimal and is a positive integer such that (9) holds, then there is an -optimal .
Proof.
Let be an optimal containing a row that is repeated times. If we delete any columns from , then the resulting array is an containing a row that is repeated times.
From Lemma 5.5, we have that
[TABLE]
so the resulting is -optimal. ∎
To illustrate the application of Theorem 5.6, we consider the case .
Theorem 5.7**.**
Suppose there is a basic and suppose
[TABLE]
is a positive integer. Then there is an -optimal .
Proof.
From Lemma 3.1, a basic has and , where is a positive integer, and . If we take , , and in (9), the inequality becomes
[TABLE]
The stated result now follows from Theorem 5.6. ∎
Corollary 5.8**.**
For , there is an -optimal .
Proof.
Begin with a basic (Example 2.3) and take in Theorem 5.7. ∎
We now discuss some examples of -optimal OAs from Culus and Toulouse [1].
Example 5.1**.**
Some OAs reported in [1, Table 4] are in fact -optimal:
- •
an with
- •
an with
- •
an with .
All of these OAs can be obtained by deleting a column from a basic OA.
Finally, the with that is depicted in [1, Table 2] is also -optimal.
We make a few observations about the inequality in Theorem 5.4.
In Section 2, we noted that taking multiple copies of an optimal OA would yield another optimal OA. However, a similar result does not hold for -optimal OAs. To illustrate, we consider OAs with and . For an , Theorem 5.4 asserts that
[TABLE]
The with mentioned in Example 5.1 is -optimal because . If we take two copies of this OA, we obtain an with . This OA is -optimal because . However, if we take three copies of the , we obtain an with , which is not -optimal because . (We do not know if an with exists.)
Finally, consider a hypothetical . Here, Theorem 5.4 yields . An with would be optimal, because (1) holds. However, an optimal cannot exist by Corollary 2.2, because . Therefore there is no -optimal .
6 Summary
Clearly there is much work to be done on the problem of constructing optimal and basic orthogonal arrays. At present, we only have a couple of examples of basic OAs with . Thus, it is of particular interest to construct additional examples, or better yet, infinite classes of these arrays.
Acknowledgements
This work benefitted from the use of the CrySP RIPPLE Facility at the University of Waterloo.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] J.-F. Culus and S. Toulouse. How far from a worst solution a random solution of a k 𝑘 k CSP instance can be? Lecture Notes in Computer Science 10979 (2018), 374–386 (IWOCA 2018).
- 2[2] A.S. Hedayat, N.J.A. Sloane and John Stufken. Orthogonal Arrays, Theory and Applications. Springer, 1999.
- 3[3] R. Mukerjee, P. Qian and J. Wu. On the existence of nested orthogonal arrays. Discrete Math. 308 (2008), 4635–4642.
- 4[4] D.R. Stinson. Applications and generalizations of the variance method in combinatorial designs. Utilitas Math. 22 (1982), 323–333.
- 5[5] D.R. Stinson. Combinatorial Designs, Constructions and Analysis , Springer, 2004.
- 6[6] D.R. Stinson. Bounds for orthogonal arrays with repeated rows. Bulletin of the ICA , to appear.
