This paper introduces a method to extend bounds on the minimum distance of cyclic codes to abelian codes using the concept of $ ext{B}$-apparent distance, improving bounds and identifying conditions for exact distances.
Contribution
It develops a novel technique to generalize ds-bounds from cyclic to abelian codes and provides criteria for when the apparent distance equals the true minimum distance.
Findings
01
Enhanced bounds for abelian code distances
02
Conditions for $ ext{B}$-apparent distance to match true minimum distance
03
Constructed tables of codes with verified minimum distances
Abstract
In this paper we develop a technique to extend any bound for the minimum distance of cyclic codes constructed from its defining sets (ds-bounds) to abelian (or multivariate) codes through the notion of B-apparent distance. We use this technique to improve the searching for new bounds for the minimum distance of abelian codes. We also study conditions for an abelian code to verify that its B-apparent distance reaches its (true) minimum distance. Then we construct some tables of such codes as an application
ε1(M)=Δ(M(g1,k))=Z(g1,k) and ω1(M)=D1,k.
ε1(M)=Δ(M(g1,k))=Z(g1,k) and ω1(M)=D1,k.
Δ(M(mk))
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
TopicsCoding theory and cryptography · graph theory and CDMA systems · Finite Group Theory Research
Full text
From ds-bounds for cyclic codes to true distance for abelian codes.
J.J. Bernal, M. Guerreiro111M. Guerreiro is with Departamento de Matemática, Universidade Federal de Viçosa, 36570-900 Viçosa-MG, Brazil. Supported by CNPq-Brazil, Processo PDE 233497/2014-5.
E-mail: [email protected]
and J. J. Simón
222J. J. Bernal and J. J. Simón are with the Departamento de Matemáticas, Universidad de Murcia, 30100 Murcia, Spain. Partially supported by MINECO project MTM2016-77445-P and Fundación Séneca of Murcia.
E-mail: {josejoaquin.bernal, jsimon}@um.es
333This work was partially presented at the ISITA 2016, Monterey, CA, USA.
Abstract
In this paper we develop a technique to extend any bound for the minimum distance of cyclic codes constructed from its defining sets (ds-bounds) to abelian (or multivariate) codes through the notion of B-apparent distance. We use this technique to improve the searching for new bounds for the minimum distance of abelian codes. We also study conditions for an abelian code to verify that its B-apparent distance reaches its (true) minimum distance. Then we construct some tables of such codes as an application.
The study of abelian codes is an important topic in Coding Theory, having an extensive literature, because they have good algebraic properties that allow one to construct good codes with efficient encoding and decoding algorithms. More precisely, regarding decoding, the two most known general techniques are permutation decoding [2] and the so-called locator decoding [6] that uses the Berlekamp-Massey algorithm [21] (see also [16]).
Even though the mentioned decoding methods require to know the minimum distance, or a bound for it, there are not much literature or studies on its computation and properties, or it does exist only for specific families of abelian codes (see [6]). Concerning BCH bound, in [10], P. Camion introduces an extension from cyclic to abelian codes which is computed through the apparent distance of such codes. Since then there have appeared some papers improving the original computation and giving a notion of multivariate BCH bound and codes (see, for example, [3, 20]).
These advances lead us to some natural questions about the extension to the multivariate case of all generalizations and improvements of the BCH bound known for cyclic codes; specifically, those bounds on the minimum distance for cyclic codes which are defined from defining sets.
There are dozens of papers on this topic regarding approaches from matrix methods ([5, 25]) through split codes techniques ([12, 15]) until arriving at the most classical generalizations based on computations over the
defining set, as the Hartmann-Tzeng (HT) bound [13], the Ross (R) bound [19] and the improvements by Van Lint and Wilson, as the shifting bound (SB) [22].
Having so many references on the subject, it seems very necessary to find a general method that allows one to extend any bound for the minimum distance of cyclic codes based on the defining set to the multivariate case. This is our first goal. We shall show a method to extend to the multivariate case any bound of the mentioned type via associating an apparent distance to such bound.
The second target of our study is to improve the searching for new bounds for abelian codes. At this point, we must honestly say that these searches may only have interest for codes whose minimum distances are not known (in fact, if one knows the minimum distance of a code one does not need a bound for it), so that, our examples consider codes of lengths necessarily large.
In our opinion, long abelian codes are not so bad (see [1]) in terms of performance.
As we work with long codes, certainly it seems impossible at the moment to compute their minimum distance, and so, it is natural to ask for conditions on them to ensure that a founded bound is in fact the minimum distance. This is the last goal of this paper. We found conditions for a bivariate abelian code to reach the mentioned equality and we write these conditions in terms of its defining set from the notion of composed polynomial matrices (CP-matrices, for short). We comment the extension of this results to several variables. We illustrate with some examples (of large codes) how this technique works.
In order to achieve our goals, we give in Section 3 a notion of defining set bound (ds-bound) for the minimum distance of cyclic codes. In Section 4, we revisit the relation between the weight of codewords of abelian codes and the apparent distance of their discrete Fourier transforms. In Section 5, we use this technique
to define the apparent distance of an abelian code with respect to a set of ds-bounds. In Section 6, we adapt a known algorithm given in [3] (of linear complexity by Remark 23) to compute the B-apparent distance of an abelian code. Finally, we study the abelian codes which verify the equality between its BCH bound and its minimum distance. For two variables, we find some sufficient conditions that are easy to extend to several variables.
2 Preliminaries
Let Fq be a finite field with q elements, with q a power of a prime p, ri be positive integers, for all i∈{1,…,s}, and n=r1⋯rs. We denote by Zri the ring of integers modulo ri and we shall always write its elements as canonical representatives.
An abelian code of length n is an ideal in the algebra Fq(r1,…,rs)=Fq[X1,…,Xs]/⟨X1r1−1,…,Xsrs−1⟩ and throughout this work we assume that this algebra is semisimple; that is, gcd(ri,q)=1, for all i∈{1,…,s}. Abelian codes are also called multidimensional cyclic codes (see, for example, [14]).
The codewords are identified with polynomials f(X1,…,Xs) in which, for each monomial, the degree of the indeterminate Xk belongs to Zrk. We denote by I the set Zr1×⋯×Zrs and we write the elements f∈Fq(r1,…,rs) as f=f(X1,…,Xs)=∑aiXi, where i=(i1,…,is)∈I and Xi=X1i1⋯Xsis. Given a polynomial f∈Fq[X1,…,Xs] we denote by f its image under the canonical projection onto Fq(r1,…,rs).
For each i∈{1,…,s}, we denote by Rri (resp. Uri) the set of all ri-th roots of unity (resp. all ri-th primitive roots of unity) and define R=∏i=1sRri (U=∏i=1sUri).
For f=f(X1,…,Xs)∈Fq[X1,…,Xs] and αˉ∈R, we write f(αˉ)=f(α1,…,αs). For i=(i1,…,is)∈I, we write αˉi=(α1i1,…,αsis).
It is a known fact that every abelian code C in Fq(r1,…,rs) is totally determined by its root set or set of zeros, namely
[TABLE]
The set of non zeros is denoted by Z(C)=R∖Z(C). For a fixed αˉ∈U, the code C is determined by its defining set, with respect to αˉ, which is defined as
[TABLE]
Given an element a=(a1,…,as)∈I, we shall define its q-orbit modulo (r1,…,rs) as Q(a)={(a1⋅qi,…,as⋅qi)∈I∣i∈N}. In the case of a semisimple algebra, it is known that any defining set Dαˉ(C) is a disjoint union of q-orbits modulo (r1,…,rs). Conversely, every union of q-orbits modulo (r1,…,rs) determines an abelian code (an ideal) in Fq(r1,…,rs) (see, for example, [3] for details). We note that the notions of root set and defining set also apply to polynomials. Moreover, if C is the ideal generated by the polynomial f in Fq(r1,…,rn), then Dαˉ(C)=Dαˉ(f).
We recall that the notion of defining set also applies to cyclic codes. For s=1 and r1=n, a q-orbit is called a ** q-cyclotomic coset** of a positive integer b modulo n and it is the set Cq(b)={b⋅qi∈Zn∣i∈N}.
Throughout this paper, we fix the notation L∣Fq for an extension field containing Uri, for all i∈{1,…,s}. The discrete Fourier transform of a polynomial f∈Fq(r1,…,rs) with respect to αˉ∈U (also called Mattson-Solomon polynomial in [20]) is the polynomial
φαˉ,f(X)=∑j∈If(αˉj)Xj∈L(r1,…,rs).
It is known that the discrete Fourier transform may be viewed as an isomorphism of algebras
φαˉ:L(r1,…,rs)⟶(L∣I∣,⋆),
where the multiplication “⋆” in L∣I∣ is defined coordinatewise. Thus, we may see φαˉ,f as a vector in L∣I∣ or as a polynomial in L(r1,…,rs) (see [10, Section 2.2]).
The inverse of the discrete Fourier transform is φαˉ,g−1(X)=r1r2⋯rs1∑j∈Ig(αˉ−j)Xj.
3 Defining set bounds for cyclic codes
In this section we deal with cyclic codes; that is, r1=n. By P(Zn) we denote the power set of Zn. We take an arbitrary α∈Un.
Definition 1
A defining set bound (or ds-bound, for short) for the minimum distance of cyclic codes is a family of relations δ={δn}n∈N such that, for each n∈N, δn⊆P(Zn)×N and it satisfies the following conditions:
If C is a cyclic code in F(n) such that N⊆Dα(C), then 1≤a≤d(C), for all (N,a)∈δn.
2. 2.
If N⊆M are subsets of Zn then
(N,a)∈δn implies (M,a)∈δn.
3. 3.
For all N∈P(Zn), (N,1)∈δn.
From now on, sometimes we write simply δ to denote a ds-bound or any of its elements independently on the length n of the code. It will be clear from the context which one is being used.
Remarks 2
(1) For example, the BCH bound states that, for any cyclic code in Fq(n) which in its set of zeros has a string of t−1 consecutive powers of some α∈Un, the minimum distance of the code is at least t [17, Theorem 7.8].
Now, define δ⊂P(Zn)×N as follows: for any a≥2, (N,a)∈δ if and only if there exist i0,i1,…,ia−2 in N which are consecutive integers modulo n. Then the BCH bound says that δ is a ds-bound for any cyclic code (we only have to state Condition 3 as a convention; so that (∅,1)∈δBCH).
(2) It is easy to check that all extensions of the BCH bound, all new bounds from the defining set of a cyclic code as in [5, 13, 18, 19, 25] and the new bounds and improvements arising from Corollary 1, Theorem 5 and results in Section 4 and Section 5 in [22], also verify Definition 1.
In general, for any bound for the minimum distance of a cyclic code, say b, we denote the corresponding ds-bound by δb. In order to relate the idea of ds-bound with the Camion’s apparent distance, which will be defined later, we consider the following family of maps.
Definition 3
Let δ be a ds-bound for the minimum distance of cyclic codes. The optimal ds-bound associated to δ is the family δ={δn}n∈N of maps δn:P(Zn)⟶N defined as
δn(N)=max{b∈N∣(N,b)∈δn}.
The following result is immediate.
Lemma 4
Let δ be a ds-bound for the minimum distance of cyclic codes. Then, for each n∈Z:
If C is a cyclic code in F(n) such that N⊆Dα(C), then 1≤δn(N)≤d(C).
2. 2.
If N⊆M⊆Zn, then
δn(N)≤δn(M).
As we noted above, we may omit the index of the map δˉn, because it will be clear from the context
for which value it is being taken.
4 Apparent distance of matrices
We begin this section recalling the notion and notation of a hypermatrix that will be used hereby, as it is described in [3]. For any i∈I=Zr1×⋯×Zrs, we write its k-th coordinate as i(k). A hypermatrix with entries in a field S indexed by I (or an I-hypermatrix over S) is an s-dimensional I-array, denoted by M=(ai)i∈I, with ai∈S [23]. The set of indices, the dimension and the ground field will be omitted if they are clear from the context. For s=2, M is a matrix and when s=1, M is a vector. We write M=0 when all its entries are [math] and M=0, otherwise. As usual, a hypercolumn is defined as HM(j,k)={ai∈M∣i(j)=k}, with 1≤j≤s and 0≤k<rj, where ai∈M means that ai is an entry of M. A hypercolumn can be seen as an (s−1)-dimensional hypermatrix. In the case s=2, we refer to hypercolumns as rows or columns and, when s=1, we say entries.
For any I-hypermatrix M with entries in a field, we define the support of M as the set supp(M)={i∈I∣ai=0}. Its complement with respect to I will be denoted by D(M). When D(M) (or supp(M)) is an union of q-orbits we say that M is a q-orbits hypermatrix.
Let D⊆I. The hypermatrix afforded by D is defined as M=(ai)i∈I, where ai=1 if i∈D and ai=0 otherwise; it will be denoted by M=M(D). Note that if D is union of q-orbits then M(D) is a q-orbits hypermatrix. To define and compute the apparent distance of an abelian code we will use the hypermatrix afforded by its defining set, with respect to αˉ∈U.
We define a partial ordering on the set \left\{M(D)\,\mid\,D\text{ is union of q−orbitsinI}\right\} as follows:
[TABLE]
Clearly, this condition is equivalent to D′⊆D.
We begin with the apparent distance of a vector in Ln.
Definition 5
Let δ be a ds-bound for the minimum distance of cyclic codes and v∈Ln a vector. The apparent distance of v with respect to δ (or δ-apparent distance of v, for short), denoted by δ∗(v), is defined as
If v=0, then δ∗(v)=0.
2. 2.
If v=0, then δ∗(v)=δ(Zn∖supp(v)).
From now on we denote by B a set of ds-bounds which are used to proceed a computation of the apparent distances of matrices, hypermatrices or abelian codes.
Definition 6
Let v∈Ln. The apparent distance of v with respect to B denoted by ΔB(v), is:
If v=0, then ΔB(v)=0.
2. 2.
If v=0, then ΔB(v)=max{δ∗(v)∣δ∈B}.
Remarks 7
The following properties arise straightforward from the definition above, for any v∈Ln.
If v=0 then ΔB(v)≥1.
2. 2.
*If supp(v)⊆supp(w) then ΔB(v)≥ΔB(w). *
Proposition 8
Let f∈L(n) and v be the vector of its coefficients. Fix any α∈Un. Then ΔB(v)≤ω(φα,f−1)=∣Z(f)∣.
Proof. Set N=Zn∖supp(v) and let C be the abelian code generated by φα,f−1 in L(n). Then d(C)≤ω(φα,f−1). By properties of the discrete Fourier transform, we have N=Dα(φα,f−1)=Dα(C) hence, by Lemma 4 and the definition of apparent distance, ΔB(v)≤d(C). This gives the desired inequality. The last equality is obviuos.
The notion of apparent distance appeared for the first time in [10] and originally it was defined for polynomials. Its computation reflects a bound of the nonzeros (in the sense given in the preliminaries) of a multivariate polynomial. The aim of the apparent distance was to extend the notion of BCH bound, from cyclic to abelian codes as we will comment in the following paragraphs. The first algorithm for its computation was made in terms of coefficients of polynomials. Later, in [20], R. E. Sabin gave an algorithm in terms of matrices. The notion of strong apparent distance, that appeared in [3], is a slight but powerful modification of the original one, defined for multivariate polynomials and hypermatrices, and it is the predecessor of the current apparent distance defined with respect to a list of ds-bounds.
Remark 9
To identify notations from previous works with the ones used here, given a polynomial f∈Fq(r1,…,rs), if we denote by M(f) its hypermatrix of coefficients (in the obvious sense), then the strong apparent distance of f in [3] is sd∗(f)=ΔδBCH(M(f)); that is, B={δBCH}, together with the convention that δBCH(∅)=1.
Now let us show how the notion of apparent distance for abelian codes works as the BCH bound for cyclic codes. All results in the following corollary are proved in [10] and [20].
Corollary 10
Let C be a cyclic code in Fq(n) and α∈Un. Then
If g,e∈C are the generating polynomial and the idempotent generator of C, respectively, then ΔB[M(φα,g)]=ΔB[M(φα,e)]≤ΔB[M(φα,c)], for all c∈C.
2. 2.
If c∈C is a codeword with φα,c=f∈L(n), then ω(c)≥ΔB(M(f)) and consequently
3. 3.
The number on the left of the last inequality is known as the apparent distance of the cyclic code C with respect to the set B and α∈Un or the B-apparent distance of C with respect to α∈Un.
Proof.(1) comes from the fact that, for all c∈C, we have supp(M(φα,c))⊆supp(M(φα,g))=supp(M(φα,e)), together with Remark 7. (2) is Proposition 8. (3) is immediate from (1) and (2).
Now we shall define the apparent distance of matrices and hypermatrices with respect to a set B of ds-bounds.
Definition 11
Let M be an s-dimensional I-hypermatrix over a field L. The apparent distance of M with respect to B, denoted by ΔB(M), is defined as follows:
ΔB(0)=0* and, for s=1, Definition 6 applies.*
2. 2.
For s=2 and a nonzero matrix M, note that HM(1,i) is the i-th row and HM(2,j) is the j-th column of M. Define the row support of M as supp1(M)={i∈{0,…r1−1}∣HM(1,i)=0} and the column support of M as supp2(M)={k∈{0,…r2−1}∣HM(2,k)=0}.
Then put
[TABLE]
and set Δ1(M)=ω1(M)⋅ϵ1(M).
Analogously, we compute the apparent distance Δ2(M) for the other variable and finally we define the apparent distance of M with respect to B by
[TABLE]
3. 3.
For s>2, proceed as follows: suppose that one knows how to compute the apparent distance
ΔB(N), for all non zero hypermatrices N of dimension s−1. Then first compute the “hypermatrix support” of M=0 with respect to the j-th hypercolumn, that is,
[TABLE]
Now put
[TABLE]
and set Δj(M)=ωj(M)⋅ϵj(M).
Finally, define the apparent distance of M with respect to B (or the B-apparent distance) as:
[TABLE]
As we have already commented in Remark 9, by taking B={δBCH}, ΔB(M) is the strong apparent distance in [3].
Now, as in the case of cyclic codes, we relate the apparent distance to the weight of codewords. For each multivariate polynomial f=∑i∈IaiXi, consider the hypermatrix of the coefficients of f, denoted by M(f)=(ai)i∈I. For any j∈{0,…,s}, if we write f=∑k=0rj−1fj,kXjk, where fj,k=fj,k(Xj) and Xj=X1⋯Xj−1⋅Xj+1⋯Xs, then M(fj,k)=HM(j,k). This means that “fixed” the variable Xj in f, for each power k of Xj, the coefficient fj,k is a polynomial in Xj, and HM(j,k) is the hypermatrix obtained from its coefficients. Now we extend Proposition 8 to several variables.
Theorem 12
Let f∈L(r1,…,rs) and M=M(f) be the hypermatrix of its coefficients. Fix any αˉ∈U. Then ΔB(M(f))≤ω(φα,f−1)=∣Z(f)∣.
Proof. For f=0, the result is obvious. Consider f=0. The case n=1 is Proposition 8. We prove the theorem for matrices; that is, for s=2. The general case follows directly by induction.
Set M=M(f)=0, αˉ=(α1,α2) and write f=∑k=0r2−1f2,kX2k. Then HM(2,k) is the vector of coefficients of f2,k. Clearly, supp2(M)={k∈{0,…,r2−1}∣f2,k=0} and, for any k∈supp2(M), we have ΔB(HM(2,k))≥1. Then ω(φα1,f2,k−1)≥1, by Proposition 8.
Now, for each fixed k∈Zr2, by the definition of discrete Fourier transform,
∣Z(f2,k)∣=ω(φα1,f2,k−1), hence, if k∈supp2(M) there exists t∈Zr1 (at least one) such that α1t is a non zero of f2,k .
Set g(X2)=f(α1t,X2)=∑k=0r2−1f2,k(α1t)X2k and note that it is a non zero polynomial. Let v(g) be the vector of coefficients of g(X2). Then supp(v(g))⊆supp2(M) and so ΔB(v(g))≥max{δˉ(Zr2∖supp2(M))∣δ∈B}=ω2(M).
As ∣Z(g)∣=ω(φα2,g−1)≥ΔB(v(g)) then, for any k∈Zr2,
[TABLE]
Finally, as in the univariate case, it is clear that ∣Z(f)∣=ω(φαˉ,f−1). The extension to more variables is clear and this completes the proof.
Example 13
Set n=96=4×24 and q=5. Fix α1∈U4 and α2∈U24 and consider the 5-orbits matrix M afforded by D=Q(0,0)∪Q(0,1)∪Q(0,2)∪Q(0,3)∪Q(0,6)∪Q(0,7)∪Q(0,9)∪Q(1,1)∪Q(1,2)∪Q(1,3)∪Q(2,1)∪Q(2,2)∪Q(3,6). Choose B={δBS,δBCH}, where
δBS is the Betti-Sala bound in [5].
One may check that supp1(M)={0,1,2,3} and then ω1(M)=1. On the other hand, supp2(M)=Z24 so ω2(M)=1. Now ΔB(HM(1,0))=8, by using δBS (see [5, Example 4.2]) which is the maximum of the values of the two bounds considered, hence ϵ1(M)=8. It is clear that ϵ2(M)=4, so that Δ1(M)=8 and Δ2(M)=4. Hence Δ{δBS,δBCH}(M)=8.
The computation of the apparent distance in several variables is a natural extension of that of one variable, and, moreover, the relationship between apparent distance and weight of a codeword is essentially the same in any case. However, condition (2) of Remark 7(2) does not necessarily hold in two or more variables and so we cannot extend directly the results of Corollary 10. Let us show the situation in the following example.
Example 14
Let M be the 2-orbits matrix of order 5×7 such that supp(M)=Q(0,0)∪Q(1,0)∪Q(1,3) and N be the 2-orbits matrix such that supp(N)=Q(1,0)∪Q(1,3). Then N<M, however, one may check that ΔδBCH(N)=6 and ΔδBCH(M)=7. So Remark 7(2) does not hold in this case.
As we will see in the next section, the notion of B-apparent distance of an abelian code with respect to some roots of unity will be a natural extension of that of cyclic codes, while its computation will be an interesting problem to solve.
5 The B-apparent distance of an abelian code
The following definition changes a little the usual way to present the notions of apparent distance from [10] and the strong apparent distance from [3] (see also [20]). We recall that B denotes a set of ds-bounds, which are used to proceed a concrete computation of the apparent distances.
Definition 15
Let C be an abelian code in Fq(r1,…,rs).
1) The apparent distance of C with respect to αˉ∈U and B (or the (B,αˉ)-apparent distance) is
[TABLE]
2) The apparent distance of C with respect to B is
[TABLE]
The following result is a consequence of Theorem 12.
Corollary 16
For any abelian code C in Fq(r1,…,rs) and any B as above, ΔB(C)≤d(C).
Proof. Let g∈C such that ω(g)=d(C). By Theorem 12, ΔB(M(φαˉ,g))≤ω(g), for any αˉ∈U. From this, the result follows directly.
It is certain that to compute the apparent distance for each element of a code in order to obtain its apparent distance can be as hard work as to compute the minimum distance of such a code. Thus, to improve the efficiency of the computation the following result tells us that we may restrict our attention to the idempotents of the code. It also allows us to reformulate Definition 15 as it is presented in the previously mentioned papers.
Proposition 17
Let C be an abelian code in Fq(r1,…,rs). The apparent distance of C with respect to αˉ∈U and B verifies the equality
[TABLE]
Proof. Consider any c∈C. Since Fq(r1,…,rs) is semisimple, then there exists an idempotent e∈C such that the ideals generated by c and e in Fq(r1,…,rs) coincide; that is, ⟨c⟩=⟨e⟩, and so Dαˉ(c)=Dαˉ(e). This means that supp(M(φαˉ,c))=supp(M(φαˉ,e)). Note that the computation of the apparent distance is based on the fact that the entries (of the matrices) are zero or not; that is, once an entry is non zero, its specific value is irrelevant. From this fact, it is easy to see that ΔB(M(φαˉ,c))=ΔB(M(φαˉ,e)) and so we get the desired equality.
Let e∈Fq(r1,…,rs) be an idempotent and E be the ideal generated by e. Then φαˉ,e⋆φαˉ,e=φαˉ,e, for any αˉ∈U and thus, if φαˉ,e=∑i∈IaiXi, we have ai∈{1,0}⊆Fq and ai=0 if and only if i∈Dαˉ(E). Hence M(φαˉ,e)=M(Dαˉ(E)). Conversely, let M be a hypermatrix afforded by a set D which is a union of q-orbits. We know that D determines a unique ideal C in Fq(r1,…,rs) such that Dαˉ(C)=D. Let e∈C be its generating idempotent. Clearly, M(φαˉ,e)=M(D).
Now let C be an abelian code, αˉ∈U and let M be the hypermatrix afforded by Dαˉ(C). For any q-orbits hypermatrix P≤M [see the ordering (1)] there exists a unique idempotent e′∈C such that P=M(φαˉ,e′) and, for any codeword f∈C, there is a unique idempotent e(f) such that ΔB(M(φαˉ,f))=ΔB(M(φαˉ,e(f))). Therefore,
[TABLE]
This fact drives us to give the following definition.
Definition 18
For a q-orbits hypermatrix M, its minimum B-apparent distance is
[TABLE]
Finally, in the next theorem we set the relationship between the apparent distance of an abelian code and the minimum apparent distance of hypermatrices.
Theorem 19
Let C be an abelian code in Fq(r1,…,rs) and let e be its generating idempotent. For any αˉ∈U, we have ΔB,αˉ(C)=B−mad(M(φαˉ,e)). Therefore, ΔB(C)=max{B−mad(M(φαˉ,e))∣αˉ∈U}.
Proof. It follows directly from the preceding paragraphs.
6 Computing minimum apparent distance
In [3] it is presented an algorithm to find, for any abelian code, a list of matrices (or hypermatrices in case of more than 2 variables) representing some of its idempotents whose apparent distances based on the BCH bound (called the strong apparent distance) go decreasing until the minimum value is reached. It is a kind of “suitable idempotents chase through hypermatrices” [3, p. 2]. This algorithm is based on certain manipulations of the (q-orbits) hypermatrix afforded by the defining set of the abelian code. It is not so hard to see that it is possible to obtain an analogous algorithm in our case.
We reproduce here the result and the algorithm in the case of two variables under our notation. Then we will use the mentioned algorithm to improve the searching for new bounds for abelian codes.
Definition 20
With the notation of the previous sections, let D be a union of q-orbits and M=M(D) the hypermatrix afforded by D. We say that HM(j,k) is an involved hypercolumn (row or column for two variables) in the computation of ΔB(M), if ΔB(HM(j,k))=ϵj(M) and Δj(M)=ΔB(M).
We denote the set of indices of involved hypercolumns by Ip(M).
Note that the involved hypercolumns are those which contribute in the computation of the B-apparent distance.
The next result shows a sufficient condition to get at once the minimum B-apparent distance of a hypermatrix.
Proposition 21
With the notation as above, let D be a union of q-orbits and M=M(D) the hypermatrix afforded by D. If ΔB(HM(j,k))=1, for some (j,k)∈Ip(M), then B−mad(M)=ΔB(M).
Proof. It is a modification of that in [3, Proposition 23] having in mind the use of different ds-bounds.
Theorem 22
Let Q be the set of all q-orbits modulo (r1,r2), μ∈{1,…,∣Q∣−1} and {Qj}j=1μ a subset of Q. Set D=∪j=1μQj and M=M(D). Then there exist two sequences: the first one is formed by nonzero q-orbits matrices, M=M0>⋯>Ml=0 and the second one is formed by positive integers
m0≥⋯≥ml, with l≤μ and mi≤ΔB(Mi), verifying the following property:
If P is a q-orbits matrix such that 0=P≤M, then ΔB(P)≥ml and if ΔB(P)<mi−1 then P≤Mi, where 0<i≤l.
Moreover, if l′∈{0,…,l} is the first element satisfying ml′=ml then
ΔB(Ml′)=B−mad(M).
Proof. It follows the same lines of that in [3, Proposition 25] having in mind the use of different ds-bounds.
Algorithm for matrices.
Set I=Zr1×Zr2. Consider the q-orbits matrix M=(mij)(i,j)∈I and a set B of ds-bounds.
Step 1. Compute the apparent distance of M with respect to B and set m0=ΔB(M).
Step 2.
a)
If there exists (j,k)∈Ip(M) (see notation below Definition 20) such that ΔB(HM(j,k))=1, then we finish giving the sequences M=M0 and m0=ΔB(M) (because of Proposition 21).
b)
If ΔB(HM(j,k))=1, for all (k,b)∈Ip(M), we set
[TABLE]
and construct the q-orbits matrix M1=(aij)(i,j)∈I such that
[TABLE]
In other words, M1<M is the (q-orbits) matrix with maximum support such that the involved rows and columns of M
are replaced by zero. One may prove that if 0=P<M and ΔB(P)<m0 then P≤M1.
Step 3.
a)
If M1=0, then we finish giving the sequences M=M0 and m0=ΔB(M).
b)
If M1=0, we set m1=min{m0,ΔB(M1)}, and we get the sequences M=M0>M1 and m0≥m1. Then, we go back to Step 1 with M1 in the place of M and m1 in the place of m0. ■
Remark 23
If the q-orbits matrix has μq-orbits, the algorithm has at most μ steps. ■
Example 24
We take the setting of Example 13 and consider the abelian code C with D(C)=D. In this case, the matrix M=M(D(C)) is the same as that in the mentioned example. Choose again B={δBCH,δBS}. This code has dimF5(C)=73 and ΔB(C)=8.
As Ip(M)={(1,0)}, the matrix M1 has the first row all zero and the others equal to the ones of M; that is D(M1)=(∪i∈Z24Q(0,i))∪Q(1,1)∪Q(1,2)∪Q(1,3)∪Q(2,1)∪Q(2,2)∪Q(3,6).
Now supp1(M1)={1,2,3}, ω1(M1)=2 and ϵ1(M1)=4 and we also have supp2(M1)={0,…,23}, ω2(M1)=1 and ϵ2(M1)=4. Hence ΔB(M1)=8. Here Ip(M1)={(1,1)} and we get the matrix M2 having its first and second rows all zero and the others equal to the ones of M and M1; that is, D(M2)=(∪i∈Z24Q(0,i))∪(∪i∈Z24Q(1,i))∪Q(2,1)∪Q(2,2)∪Q(3,6).
Here supp1(M2)={2,3}, ω1(M2)=3 and ϵ1(M2)=3 and we also have supp2(M2)={0,…,23}, ω2(M2)=1 and ϵ2(M2)=4. Hence ΔB(M2)=9 and Ip(M2)={(1,2)} and we get the matrix M3 having its first, second and third rows all zero and the others equal to the ones of M, M1 and M2; that is, D(M2)=(∪i∈Z24Q(0,i))∪(∪i∈Z24Q(1,i))∪(∪i∈Z24Q(2,i))∪Q(3,6).
Finally, supp1(M3)={3}, ω1(M2)=4 and ϵ1(M2)=2; so we also have supp2(M2)=Z24∖{6}, ω2(M3)=2 and ϵ2(M3)=4. Hence ΔB(M3)=8 and Ip(M3)={(1,3)} and we get M4=0. Therefore, Δ{δBCH,δBS}(M)=8.
The closest code to C we know is a (105,51,7) binary cyclic code in [12, Table II]. The known bounds for linear codes with the same length and dimension are between 10 and 15.
The reader may find some tables with examples of this kind in [4].
7 True minimum distance in abelian codes
In this section we study the problem of find abelian codes such that its apparent distance or its multivariate BCH bound reaches its minimum distance. We keep all the notation from the preceding sections. In [8, 9] it is presented a characterization of cyclic and BCH codes whose apparent distance reaches their minimum distance. Our aim is to extend those results for multivariate codes.
Theorem 25
Let C be an abelian code in Fq(r1,…,rs). The following conditions are equivalent:
ΔB(C)=d(C).
2. 2.
There exist an element α∈U and a codeword c∈C such that its image under the discrete Fourier transform, g=φα,c, verifies:
(a)
ΔB(M(g))=ΔB,α(C)=min{ΔB(M(φα,v))∣v∈C}.
2. (b)
ΔB(M(g))=Z(g).
Proof.[1⇒2] Let c∈C a be codeword with ω(c)=d(C) and α∈U such that
ΔB,α(C)=ΔB(C). Set g=φα,c. Then
[TABLE]
where the first equality is given by hypothesis. Thus, the inequality becomes equality.
[2⇒1] Suppose that there is a codeword c∈C satisfying the hypotheses. Then
[TABLE]
and, again, equalities hold.
Remarks 26
The conditions in statement (2) of Theorem 25 refers only to a single element in U. This is an important reduction that will be very useful later. On the other hand, we recall that if M is the hypermatrix afforded by Dα(C) then ΔB,α(C)=B−mad(M).
For a given abelian code, the problem of finding, if any, a codeword verifying Condition (2) of Theorem 25 is in general difficult to solve. In the case that the codeword is an idempotent we will be able to find it through the computation of the minimum apparent distance. So far we only know an actual way to find the desired idempotent codeword of the mentioned theorem and it is given in the following result.
Proposition 27
Let C be a code in Fq(r1,r2), α∈U and M the matrix afforded by its defining set
Dα(C). Let P≤M be a q-orbits matrix and e∈L(r1,r2) be the idempotent such that that P=M(e). If P verifies
Δ(P)B=B−mad(M)* and*
2. 2.
Δ(P)B=Z(e).
Then d(C)=B−mad(M)=ΔB(C).
Proof. Since P≤M, then φα,e−1∈C, hence ω(φα,e−1)≥d(C)≥ΔB(C). On the other hand, by the hypothesis (2),
[TABLE]
Therefore, d(C)=B−mad(M)=ΔB(C).
So, for a given code C with afforded matrix M, if we want to know whether d(C)=ΔB(C) by using
Proposition 27, in the case that the requested codeword is an idempotent, we have to analyze all (q-orbits) matrices P≤M. If Dα(C)=t and such idempotent exists, we have to do at most 2t steps. This is a search on the set of idempotents of C. The reader may note that the original computation of the apparent distance in [10] and [20] requires to compute the apparent distance of exactly the same set of q-orbits matrices. This might be an important reduction in some cases.
We wonder if it is possible to simplify the procedure to find a q-orbits matrix P, as in Proposition 27, by analyzing the sequence of matrices in the algorithm for the computation of the minimum apparent distance of the matrix M afforded by Dα(C); i.e. the computation of B−mad(M). The algorithm gives us an interesting reduction.
In the algorithm for the computation of the strong apparent distance as in Theorem 22, we consider the sequence of matrices
[TABLE]
and let j0 be the first index such that ΔB(Mj0)=mℓ=B−mad(M). If m0=ΔB(M) equals mℓ, then P=M and we do not have any reduction. However, if m0>mℓ, then ΔB(P)=mℓ<mj0−1 which implies P≤Mj0, hence we can start our search from Mj0;
that is, we have to check only at most 2t−j0 matrices in order to find the hypothetical matrix of Theorem 25
We wonder if the existence of a matrix P≤M satisfying the conditions of Proposition 27 implies the existence of a matrix in the sequence (2) also satisfying those conditions. The answer is negative, as the following very simple example shows.
Example 28
Set Δ=ΔB, with B={δBCH}. There exists an abelian code C, with matrix M afforded by Dα(C) with respect to α∈U such that:
For every q-orbits matrix in the sequence M=M0>⋯>0 we have Δ(Mj)=Z(ej), where ej∈L(r1,r2) is the idempotent that verifies Mj=M(ej).
2. 2.
d(C)=Δ(C)**
Proof. Set q=2 and r1=r2=7. Let C be the code such that Dα(C)=Q(0,3)∪Q(1,3)∪Q(1,5)∪Q(1,6)∪Q(3,0)∪(3,2)∪Q(3,3)∪Q(3,4)∪Q(3,5)∪Q(3,6) with respect to α∈U. The matrix afforded by Dα(C) is
[TABLE]
Let a(X1)=(1+X1)(1+X12+X13), b(X2)=(1+X2)(1+X22+X23). If e∈C is the idempotent generator then φα,e(X1,X2)=a(X1)b(X2)+X13X2+X16X22+X13X24. One may compute
Z(φα,e)=25, by using GAP.
On the other hand, computing B−mad(M), we obtain the chain M0>0 and Δ(M)=B−mad(M)=9. Now consider
the q-orbits matrix
[TABLE]
Note that P does not belong to the sequence M0>0; however, as a∣X17−1 and b∣X27−1, then g1=ab satisfies the hypothesis of Proposition 34 which means that C satisfies the condition (2b) of Theorem 25. Now, as Δ(P)=9 then condition (2a) of the same theorem is satisfied and thus d(C)=Δ(C).
As we have seen, although Proposition 27 gives us a sufficient condition, it does not guarantee that we can find the desired codeword if it is not an idempotent, not even by using the algorithm. Now, in order to construct codes C satisfying d(C)=ΔB(C), we try to move forward into a different direction: we firstly characterize those polynomials that verify (2b); that is, ΔB(M(g))=∣Z(g)∣. In the univariate case, those polynomials were characterized in [9]. Before to extend the results to multivariate polynomials, we need to put some restrictions on the election of ds-bounds that we may use. In the case of polynomials in one variable, one may see that the condition (2b) in Theorem 25 forces us to use exclusively the BCH bound as, to the best of our knowledge, the computation of ∣Z(g)∣ is only known to be obtained in terms of the degree of Xhg (viewed as a polynomial); that is, a list of consecutive exponents of the monomials with coefficient zero of the highest degrees.
In the reminder of this section, consider B={δBCH} and denote Δ=ΔδBCH, for the sake of simplicity. Let us recall some facts from univariate polynomials that will be used herein. For 0=g∈L(n), let mg=gcd(Xn−1,Xhg), which does not depend on h∈N. As we pointed out in Remark 9 about notation, sd∗=Δ{δBCH}=Δ. The proof of the following result is essentially the same as [9, Proposition 1]; so we ommit it.
Proposition 29
Let g∈L(n). Then
Δ(M(g))=Z(g)* if and only if Xhg∣Xn−1, for some h∈N.*
Consider g=g(X1,X2)∈L(r1,r2) and write M=M(g). In general, as we have seen in Theorem 12, Δ(M(g))≤Z(g). We want to describe the polynomials g in L(r1,r2) such that Δ(M(g))=Z(g) so we assume that the equality holds; moreover, we impose the following condition
[TABLE]
where, as in Definition 11, Δi(M)=ωi(M)⋅ϵi(M), for i∈{1,2}.
We also write, following the notation of paragraph prior to Theorem 12,
[TABLE]
For all u∈L, the polynomials of the form g(u,X2) and g(X1,u) have the obvious meaning. For j=1,2,
denote by Zj=πj(Z(g)) the projection of Z(g)
onto the j-th coordinate and we also write Zj=πj(Z(g)).
For each j∈{1,2}, we set Mj={k∈{0,…,ri−1}∣(j,k)∈Ip(M)}; that is, Δ(M(gj,k))=εj, for any k∈Mj. Note that Ip(M)=(1,M1)∪(2,M2), (see Definition 20).
Now, for each k∈M1, define
[TABLE]
and analogously, for each k∈M2, define
[TABLE]
So, if we set, for j∈{1,2}, Dj,k=min{Z(g(u,Xj))∣u∈Z(gj,k)} then Z(gj,k)⋅Dj,k≤∣Dj,k∣.
We know that, for each k∈M1, it happens ε1(M)=Δ(M(g1,k))≤Z(g1,k) and for each v∈Z(g1,k), we have ω1(M)≤Δ(M(g(X1,v))) and so ω1(M)≤D1,k, hence ε1(M)ω1(M)≤∣D1,k∣. However, by the condition (3), ε1(M)ω1(M)=Δ1(M)=Z(g) and, by definition, ∣D1,k∣≤Z(g). Therefore, for all k∈M1,
[TABLE]
Analogously, D2,k=Z(g), for k∈M2.
In fact, for j∈{1,2} and any k∈Mj, we have Z(gj,k)⋅Dj,k=∣Dj,k∣, hence
[TABLE]
Keeping in mind the equalities obtained in the previous paragraphs, we get the following two results.
Lemma 30
Let g=g(X1,X2)∈L(r1,r2) be a polynomial such that M=M(g) satisfies the
condition (3). Then:
For each k∈M1, Δ(M(g1,k))=Z(g1,k)=Z2 and Δ(M(g(X1,v)))=Z(g(X1,v)), for any
v∈Z(g1,k).
2. 2.
For each k∈M2, Δ(M(g2,k))=Z(g2,k)=Z1 and Δ(M(g(u,X2)))=Z(g(u,X2)), for any
u∈Z(g2,k).
Proof. We prove (1.) as the the proof of (2.) is entirely analogous. As we have already seen, if condition (3) is satisfied then (5) and (6) also hold; so that, ε1(M)=Δ(M(g1,k))=Z(g1,k).
Once we have the first equality, if ω1(M)=Δ(M(g(X1,v)))<Z(g(X1,v)), for some v∈Z(g1,k), then it must happen ε1(M)ω1(M)<∣D1,j∣, a contradiction. Finally, if v∈Z2, then there exists u∈Rr1 such that (u,v)∈Z(g)=D1,k, hence u∈Z(g1,k). This proves the other equalities of this lemma.
Proposition 31
Let g=g(X1,X2)∈L(r1,r2) be a polynomial such that M=M(g) satisfies the condition (3). Then there exist a=a(X1)∈L(r1),b=b(X2)∈L(r2) and
F=F(X1,X2)∈L(r1,r2) such that
g=abF and
X1h1a∣(X1r1−1), for some h1∈Z1, with Δ(M(a))=ε2(M).
2. 2.
X2h2b∣(X2r2−1), for some h2∈Z2, with Δ(M(b))=ε1(M).
Proof. By Lemma 30.1 and by Proposition 29, for each k∈M2, if we denote mk=gcd(g2,k,X1r1−1), then
[TABLE]
By definition, Δ(M(g2,k))=r1−deg(X1k′g2,k), for some k′∈N. As Δ(M(mk))=Δ(M(g2,k)) and, by [9, Lemma 2],
Δ(M(mk))=r1−degmk, then X1k′g2,k and mk are associated.
Now we claim that mk∣g2,j, for all j∈{0,…,r2−1}. Indeed, for a fixed
k∈M2, by (5), we have D2,k=Z(g) which implies Z(g2,j)⊆Z(g2,k) or rather Z(g2,k)⊆Z(g2,j). Hence, mk∣g2,j, for all j∈{0,…,r2−1}.
Denoting by g2,j′=mkg2,j, for all j∈{0,…,r2−1} and a(X1)=mk, we may write
[TABLE]
with Δ(M(a(X1)))=ε2(M).
Analogously, we get
[TABLE]
with Δ(M(b(X2)))=ε1(M) and b(X2)g1,i′=g1,i, for any i∈{0,…,r1−1}.
It is important to note that 1=gcd(X1r1−1,g2,0′,…,g2,r2−1′) and
1=gcd(X2r2−1,g1,0′,…,g1,r1−1′).
Now by writing
[TABLE]
we get g(X1,X2)=a(X1)f(X1,X2)=b(X2)h(X1,X2).
Recall that Z1=π1(Z(g)) and Z2=π2(Z(g)).
First note that if v∈Z2, then there exists u∈Z1 such that
(u,v)∈Z(g). This implies (u,v)∈D1,k, for k∈M1,
hence v∈Z(g2,k)=Z(b).
Therefore, Z2⊆Z(b) and, by Lemma 30.1,
Z2=Z(b).
Consider g(X1,v)=b(v)h(X1,v)=b(v)∑i=0r1−1g2,i′(v)X1i, for any v∈Rr2. If b(v)=0, then g(X1,v)=0, otherwise all g2,i′ would have a common zero, which is not possible. Conversely,
if g(X1,v)=0, then b(v)=0. This proves that v∈Z(b)=Z2 if and only if g(X1,v)=0.
Obviously g(X1,v)=0 also implies b(v)=0, hence v∈Z2, as Z(b)⊆Z2. Now let us write
[TABLE]
If v∈Z(b), then g(X1,v)=0 and we have f(X1,v)=0. Since a(X1)=0 we must have f2,i(v)=0, for all i∈{0,…,r1−1}. Hence Z(b)⊂Z(fi) and b(X2)∣f2,i, for all i∈{0,…,r1−1}. Now if f2,i(v)=0, for all i∈{0,…,r1−1}, then f(X1,v)=0 and g(X1,v)=0, which implies v∈Z(b), as we have seen before. Hence b(X2)=gcd(X2r2−1,f2,i), for all i∈{0,…,r1−1}. Therefore, f(X1,X2)=b(X2)f′(X1,X2) and
[TABLE]
Analogously, one may prove that a(X1)∣h1,j, for all j∈{0,…,r2−1}, and get h(X1,X2)=a(X1)h′(X1,X2), hence
[TABLE]
Finally, note that the decompositions g=abf′ and g=abh′ from (9) and (10), has been done in L[X1,X2], which is a domain, and so, we have f′(X1,X2)=h′(X1,X2). By writing F(X1,X2)=f′(X1,X2)=h′(X1,X2), we get
[TABLE]
It is clear that the condition (3) plays an important role in all the previous proofs. Recall that a polynomial g∈L(r1,r2), with coefficient matrix M=M(g), satisfies such condition if
[TABLE]
For those polynomials, we have obtained a factorization g=abF, which describes them, where Δ(M(a))=ε2(M) and Δ(M(b))=ε1(M).
At this point, two questions arise for an abelian code satisfying Theorem 25 with a codeword image g=φα,c as in such theorem.
Is it always true that g satisfies also condition (3)?
2. 2.
Suppose a polynomial g∈L(r1,r2) already satisfies condition (3) and so we have a decomposition g=abF. What can we say about F? More specifically, is it true that M(F) is a q-orbits matrix? And is it true that Δ(M(a))=Δ(M(b))=Δ(M(g))?
We shall answer all these questions in the following examples.
Example 32
There exists an abelian code C generated by an idempotent e∈Fq(r1,r2), with image g=φα,f, satisfying the following properties.
Δ1(M(g))<Δ2(M(g))* (so the condition (3) is not fully satisfied)*
2. 2.
Δ(M(g))=Z(g).
3. 3.
d(C)=Δ(C).
Proof. Set q=2, r1=5, r2=9 and C be the code with D(C)=Q(1,3), a minimal code with generator idempotent e(X1,X2)=X14X27+X13X28+X14X26+X12X28+X13X26+X14X24+X13X25+X12X26+X1X27+X14X23+X12X25+X1X26+X13X23+X14X2+X13X22+X12X23+X1X24+X14+X12X22+X1X23+X13+X12+X1X2+X1. Using the program GAP, we computed
d(C)=24. One may check that
φα,e=g(X1,X2)=X1X23+X14X23+X12X26+X13X26. Some direct computations yield Δ1(M(g))=18 and Δ(M(g))=24, so assertion (1) of this example is satisfied. As ω(e)=24 we also get assertion (2). Since C is minimal, we have, 24=Δ(M(g))=B−mad(M(g)). On the other hand Δ(M(g))≤Δ(C)≤d(C)=24. Thus d(C)=Δ(C) and we get assertion (3).
The previous example gives a negative answer to question 1. The following example answers question 2.
Example 33
Under the same notation of Proposition 31, there exists an abelian code C generated by an
idempotent e∈Fq(r1,r2), with image g=φα,e,
such that the following properties hold.
M(g)* satisfies the condition (3) and then g=abF as in the mentioned proposition (see paragraph prior to Example 32).*
2. 2.
d(C)=Δ(C).
3. 3.
F(X1,X2)* has at least two nonzero monomials and M(F) is not a q-orbits matrix.*
4. 4.
Δ(a)=Δ(b), but Δ(a)Δ(b)=Δ(M(g))
5. 5.
Z(g)=Z1×Z2**
Proof. Let q=2, r1=r2=5 and C be the code with D(C)=Q(1,1)∪Q(1,3).
In this case, g(X1,X2)=X14X24+X13X24+X14X22+X13X23+X12X22+X1X23+X12X2+X1X2 and φα,g−1=X13X24+X14X22+X14X2+X13X22+X12X23+X1X24+X1X23+X12X2, so that ∣Z(g)∣=8.
By using GAP, we computed d(C)=8 and one may check that Δ1(M(g))=Δ2(M(g))=Δ(M(g))=8; so that (1.) and (2.) hold.
Let us factorize g. In the case
[TABLE]
M2={1,4} and a=a(X1)=1+X1. Note that (1+X1)X1 is a common factor in all summands of g(X1)(X2).
On the other hand,
[TABLE]
M1={2,3} and b=b(X2)=1+X2. Here (1+X2)X2 is a common factor in all summands of g(X2)(X1).
Moreover,
[TABLE]
This gives us (3.)
Now, one may easily check that Δ(M(a))=Δ(M(b))=4, hence Δ(M(a))⋅Δ(M(b))=Δ(M(g)). This gives us (4).
Finally, by using GAP we compute Z(g)=Q(1,3)∪Q(1,4) and clearly Z(g)=Z1×Z2.
To finish our argumentation from Theorem 25 we prove that a polynomial that satisfies condition (3), and so factorizes g=abF with F(X1,X2) a monomial in L(r1,r2), verifies that its image under the discrete Fourier transform satisfy condition (2) of the mentioned theorem.
Proposition 34
Suppose g∈L(r1,r2) is such that g(X1,X2)=a(X1)b(X2), where a and b satisfy Proposition 29. Set M=M(g). Then
Z(g)=Z1×Z2**
2. 2.
Δ(M)=Δ(M(a))⋅Δ(M(b))=Z(g).
3. 3.
Δ1(M)=Δ2(M)=Δ(M)=Z(g)* (the condition (3)).*
4. 4.
Δ(M(a))=ε2(M)=ω1(M)* and*
5. 5.
Δ(M(b))=ε1(M)=ω2(M).
Proof. Assertion (1) and the equality Δ(M(a))⋅Δ(M(b))=Z(g) come directly from the decomposition of g together with the hypothesis that a and b satisfy Proposition 29.
Now set M=M(g). Since g=ab then Hm(2,j)=M(a), for all j∈supp2(M) and supp2(M)=supp(M(b)); so that ε2(M)=Δ(M(a)) and ω2(M)=Δ(M(b)). Analogously, ε1(M)=Δ(M(a)) and ω1(M)=Δ(M(b)). From this, we get all assertions.
To sum up, from Proposition 31 we have obtained what kind of polynomials we have to use to reach condition (3). This will be the main idea in order to construct abelian codes C, with d(C)=Δ(C). We address this problem in the following section.
7.1 Application 1: construction of abelian codes for which its multivariate BCH bound, apparent distance and minimum distance coincide.
In this section, we continue considering B={δBCH} and denoting Δ=ΔδBCH, for the sake of simplicity by the same reasons given in the paragraphs prior Proposition 29. Bearing in mind Proposition 34 and Proposition 27, we introduce the following definition.
Definition 35
A matrix P of order r1×r2, with entries in L is called a composed polynomial matrix (CP-matrix, for short) if there exist polynomials a=a(X1)∈L(r1) and b=b(X2)∈L(r2) such that P=M(ab), where ab∈L(r1,r2).
Note that, for a CP-matrix P, its support is a direct product supp(P)=supp(a(X1))×supp(b(X2)). The polynomials a and b are called the polynomial factors of P. The reader may see that to check if a matrix is a CP-matrix is a trivial task, because it must happen π1(supp(P))=supp(a)=supp(M(a)) and π2(supp(P))=supp(b)=supp(M(b)). The following result is an immediate consequence of Proposition 34.
Corollary 36
Let P=M(g) be a CP-matrix of order r1×r2 with polynomial factors a and b; that is, g=ab.
If X1h1a∣X1r1−1 and X2h2b∣X2r2−1, for some h1,h2∈N, then
Z(ab)=Z(a)×Z(b).**
2. 2.
Δ(P)=Δ(M(a))⋅Δ(M(b))=Z(g).
3. 3.
Δ1(P)=Δ2(P)=Δ(P)=Z(g)* (the condition (3)).*
4. 4.
Δ(M(a))=ε2(P)=ω1(P)* and*
5. 5.
Δ(M(b))=ε1(P)=ω2(P).
Example 37
Set q=2, r1=3 and r2=7; so that n=21. Let P be the CP-matrix with polynomial factors a=X1+X12 and b=X2+X22+X24, and g=ab. In this case, X12a∣X13−1 and X26b∣X27−1. Now Δ(M(a))=Δ(M(X12a))=2,
Δ(M(b))=Δ(M(X26b))=4; hence Δ(P)=8=∣Z(g)∣.
The next example shows that the hypothesis on the polynomials a and b of Corollary 36 are not superfluous.
Example 38
Set q=2, r1=5 and r2=7; so that n=35. Let P be the CP-matrix with factors a=X1+X12+X13+X14 and b=X2+X22+X24. In this case, X1h1a∤X13−1, for all h1∈Z5. On the other hand X26b∣X27−1. Now Δ(M(a))=2, Δ(M(b))=4. Although Δ1(P)=Δ2(P)=Δ(P)=8, one may check that Z(ab)=16.
Now we give a method for constructing the desired abelian codes. First, a technical lemma.
Lemma 39
Let D⊂Zr1×Zr2 be union of q-orbits and M=M(D), the matrix afforded by D. If supp(M)=π1(supp(M))×π2(supp(M)) then B−mad=ΔB(M), where B is any set of ds-bounds.
In the case B={δBCH}, the above equality coincides with the multivariate BCH bound in [3, Theorem 30].
Proof. Clearly, in this case all rows (columns) have the same support and so if one row or column is involved then all of them are too. The last assertion comes directly from the computation of the multivariate BCH bound.
Remark 40
We have already mentioned that for any CP-matrix, M, one has supp(M)=π1(supp(M))×π2(supp(M)). The converse is true for those matrices satisfying hypothesis of lemma above; that is, if D⊂Zr1×Zr2 is union of q-orbits and M=M(D) is the matrix afforded by D with supp(M)=π1(supp(M))×π2(supp(M)) then M is a CP-matrix.
Theorem 41
Let K be an intermediate field Fq⊆K⊆L, a=a(X1)∈K(r1) and b=b(X2)∈K(r2) be such that a∣X1r1−1 and b∣X2r2−1. If there exist (α1,α2)∈U, h1∈Zr1 and h2∈Zr2 for which φα1,X1h1a−1∈Fq(r1) and φα2,X2h2b−1∈Fq(r2), then the abelian code
C=(φα1,X1h1a−1⋅φα2,X2h2b−1)
in Fq(r1,r2) verifies Δ(M(ab))=Δ(C)=d(C).
Moreover, in this case, for any β1∈Ur1 and β2∈Ur2 the abelian code C(β1,β2)=(φβ1,X1h1a−1⋅φβ2,X2h2b−1) is an ideal of Fq(r1,r2) and verifies Δ(M(ab))=Δ(C(β1,β2))=d(C(β1,β2))=d(C).
Proof. Set g(X1,X2)=X1h1a(X1)⋅X2h2b(X2) and α=(α1,α2). By definition of the discrete Fourier transform, it is easy to see that the particular factorization of g implies that φ(α1,α2),g−1=φα1,X1h1a−1⋅φα2,X2h2b−1. On the other hand, it is clear that M(g) is a CP-matrix satisfying the hypothesis of Corollary 36. This, in turn, implies that statement 2(b) of Theorem 25 is satisfied.
Let M be the matrix afforded by Dα(C). Since C=(φα,g−1) then supp(M)=supp(M(g)), hence M is also a CP-matrix and Δ(M)=Δ(M(g)). By Lemma 39, B−mad(M)=Δ(M) and so statement 2(a) of Theorem 25 is also satisfied. Thus Δ(C)=d(C). The final assertion is a direct consequence of [9, Remark 2] together with the fact that, under these hypothesis, all afforded matrices are CP-matrices.
Now, we may apply all known criteria for univariate polynomials to have inverse of the discrete Fourier transform in an specific quotient ring. The following corollary concretizes the proposed construction.
It comes from [9, Remark 2] and the theorem above.
Corollary 42
Let K be an intermediate field Fq⊆K⊆L, a=a(X1)∈K(r1) and b=b(X2)∈K(r2) be such that a∣X1r1−1 and b∣X2r2−1. If there exist (α1,α2)∈U, h1∈Zr1 and h2∈Zr2 for which
[(X1h1a)(α1i)]q=(X1h1a)(α1i),
for all i∈{0,…,r1−1} and [(X2h2b)(α2j)]q=(X2h2b)(α2j), for all j∈{0,…,r2−1}, then the
family of abelian codes
[TABLE]
in Fq(r1,r2) verifies Δ(M(ab))=Δ(C(β1,β2))=d(C(β1,β2)).
The following example shows how to use Corollary 42.
Example 43
Set q=2, r1=3 and r2=45 (so n=135). Fix α1∈U3 and α2∈U45. Consider the polynomials a=X+1 and b=Y40+Y39+Y38+Y36+Y35+Y32+Y30+Y25+Y24+Y23+Y21+Y20+Y17+Y15+Y10+Y9+Y8+Y6+Y5+Y2+1. Then a∣X3−1. Note that supp(X⋅a(X))={1,2}=C2(1) modulo 3. By [9, Lemma 1] we have
that h1=1 works. Now, the polynomial b appears in [9, Example 5] where it was mentioned that b∣x45−1 in F2[x] (so that K=F2). In that example, it is shown that, for α2∈U45 (for instance, the one with minimal polynomial Y12+Y3+1), since b(1)=1 and b(α23)=α230, then (Y5b)(1)=1, (Y5b)(α23)=(α23)5α230=α245=1. So h2=5 will work because (Y5b)(α26)=(Y5b)(α212)=(Y5b)(α224)=1; note that C2(3)={3,6,12,24} modulo 45. Now set C=(φα1,Xa−1⋅φα2,Y5b−1)⊂F2(r1,r2)=F2(3,45). Then D(α1,α2)(C)=C2(1)×(C2(1)∪C2(3)∪C2(9)∪C2(21)).
One may check that 10=Δ(M(ab)); so that d(C)=10 and dimF2(C)=87.
The next example shows that from a code satisfying the conditions of Theorem 25, we can obtain a code with better parameters by making slight modifications on the defining set in such a way that the new code verifies the same conditions, but it has higher dimension, for example.
Example 44
Set q=2, r1=3 and r2=45. Fix α1∈U3 and α2∈U45. Consider the code C in
Example 43; that is D(α1,α2)(C)=C2(1)×(C2(1)∪C2(3)∪C2(9)∪C2(21)) and set g=Xa⋅Y5b. As one may check there are three subsets determining Δ(M(g)); to witt
[TABLE]
If one computes Δ(M(g)) by considering S1 then clearly C2(1)×C2(21) will have no influence in the computation. Hence one may construct the new code C′ for which D(α1,α2)(C′)=C2(1)×(C2(1)∪C2(3)∪C2(9)) such that Δ(C′)=Δ(C)=Δ(M(g)).
Note that the matrix afforded by D=D(α1,α2)(C′) is also a CP-matrix and so B−mad(M(D))=Δ(M(D))=10. Since C is a subcode of C′ then c=φα1,g−1∈C′ and, clearly, g satisfies conditions (2a) and (2b) of Theorem 25 for C′; hence Δ(C′)=Δ(M(g))=d(C′)=10=d(C)=Δ(C). Since dimF2(C′)=95 and dimF2(C)=87, C′ is a code with better parameters than those of C.
Let K be an intermediate field Fq⊆K⊆L and a=a(X1)∈K(r1)
be such that a∣X1r1−1 with φα1,X1h1a−1∈Fq(r1), for some α1∈Ur1 and h1∈Zr1.
Let g be an irreducible factor of X2r2−1 in K[X2] with defining set Dα2(g),
for some α2∈Ur2. Set b=(X2n−1)/g. If there are positive integers j,t such that b(α2j)=α2t and gcd(j,gcd(q−1,r2)r2)∣t, then there exists h2∈Zr2 such that the abelian code C=(φα1,X1h1a−1⋅φα2,X2h2b−1) in Fq(r1,r2) verifies Δ(M(ab))=Δ(C)=d(C).
Our last application of this section is the following result that comes from [9, Corollary 7].
Corollary 46
Let a=a(X1)∈L(r1) be such that a∣X1r1−1 with φα1,X1h1a−1∈F2(r1), for some α1∈Ur1 and h1∈Zr1, and suppose r2=2m−1, for some m∈N. Then there exist at least mϕ(r2) binary codes C of length n=r1r2 such that Δ(M(ab))=Δ(C)=d(C).
In the following examples we take the advantage of the information for cyclic codes from the tables appearing in [9]
Example 47
We first show in Table 1 some abelian codes of lenght 7×15=105 constructed from a list of divisors ai of X17−1 in F2[X1] and divisors bj of X215−1 in F2[X2] as in Corollary 45. The divisors are a1=1+X1, a2=1+X1+X13, a3=1+X12+X13, b1=1+X2+X22X215−1, b2=1+X2+X24X215−1 and b3=1+X23+X24X215−1.
In Table 2 we also have abelian codes of lenght 105, but in this case we consider r1=5 and r2=21.
Here we choose only one divisor of X15−1 in F2[X1]; the 5-th cyclotomic polynomial, Φ5. The other divisors bj′ of X221−1 in F2[X2] from which we may construct abelian codes as in Corollary 45 are b1′=1+X2+X22X221−1, b2′=1+X2+X23X221−1, b3′=1+X22+X23X221−1, b4′=1+X2+X22+X24+X26X221−1 and b5′=1+X22+X24+X25+X26X221−1.
7.2 Application 2: True distance in BCH multivariate codes
In [3, Definition 33], the notion of BCH multivariate code appears. Let us recall this definition focused on the bivariate case.
Definition 48
Let γˉ⊆{1,2} and δˉ={(δk)k∈γˉ∣2≤δk≤rk}. An abelian code C in Fq(r1,r2) is a bivariate BCH code of designed distance δˉ if there exists a list of positive integers bˉ=(bk)k∈γˉ such that
[TABLE]
for some α∈U, where {bk,…,bk+δk−2} is a list of consecutive integers modulo rk and I(k,u)={i∈I∣i(k)=u}.
The BCH multivariate codes are denoted Bq(α,γˉ,δˉ,bˉ).
Let C be an abelian code in Fq(r1,r2) with M=M(Dα(C)) the matrix afforded by its defining set with respect to some αˉ=(α1,α2)∈U. If M satisfies supp(M)=π1(supp(M))×π2(supp(M)) then Dα(C)=π1(supp(M))×π2(supp(M)). We set S1=π1(supp(M)) and S2=π2(supp(M)). Then, one may consider the cyclic codes C1 and C2 with defining sets D1=Zr1∖S1 and D2=Zr2∖S2 with respect to α1 and α2, respectively (note that it may happen Dα(C)=D1×D2).
Now suppose that the code C is an abelian code as described in Theorem 41 keeping the notation for the polynomials a and b and having in mind Remark 40. By the proof of this theorem one also may deduce that viewing φα1,X1h1a−1 in Fq(r1) and φα2,X2h2b−1 in Fq(r2) it happens that C1=(φα1,X1h1a−1)⊆Fq(r1) and C2=(φα2,X2h2b−1)⊆Fq(r2). It is also clear that the cyclic codes C1 and C2 verify that their minimum distances equal their respective maximum BCH bounds, as a and b satisfy the conditions in [9, Corollary 5].
In some sense we may consider C1 and C2 as “projected codes”, but clearly C is not a product of them under any classical algebraic operation; however, for their defining sets the equality under the product occurs, and all properties of ds-bounds are related, as the following results show.
Lemma 49
Under the same notation from previous paragraphs, let C be an abelian code in Fq(r1,r2), with
M=M(Dα(C)) and suppose supp(M)=π1(supp(M))×π2(supp(M)). Consider D1=Zr1∖π1(supp(M)), D2=Zr2∖π2(supp(M)) and let Ci be the
cyclic code with Dαi(Ci)=Di, for i∈{1,2}. Then
For any set B of ds-bounds, ΔB,α(C)=ΔB,α1(C1)⋅ΔB,α2(C2).
2. 2.
C* is a nonzero BCH multivariate code if and only if C1 and C2 are BCH cyclic codes in the classical sense (see [17]).*
Moreover, if case (2) holds, with Ci=(αi,δi,bi), for i∈{1,2}, then C=Bq((α1,α2),{1,2},{δ1,δ2},{b1,b2}).
Proof. Assertion (1) comes from Corollary 36, having in mind Remark 40.
Now we prove assertion (2). First, suppose that C is a multivariate BCH code. Assume that 1∈γˉ, and let B={b1,…,b1+δ1−2} be its list of consecutive integers
modulo r1. Consider a q-cyclotomic coset T⊆D1 and take t∈T. Since t∈D1
then the set (t,Zr2)⊆Dα(C) (with the obvious meaning).
If (t,Zr2)∩Q(b1+l,Zr2)=∅, for some l∈{0,…,δ1−2}, then we are done. Otherwise, it must happen 2∈γˉ and (t,Zr2)⊂⋃l=0δ2−2Q(Zr1,b2+l). Since C=0 then we may take an element u∈Zr2∖⋃l=0δ2−2Cq(b2+l) and clearly (t,u)∈⋃l=0δ2−2Q(Zr1,b2+l), a contradiction.
The final assertion comes immediately from the fact that M is a CP-matrix.
Theorem 50
Let K be an intermediate field Fq⊆K⊆L, a=a(X1)∈K(r1) and
b=b(X2)∈K(r2) be such that a∣X1r1−1 and b∣X2r2−1. If there exist (α1,α2)∈U,
h1∈Zr1 and h2∈Zr2 for which φα1,X1h1a−1∈Fq(r1) and
φα2,X2h2b−1∈Fq(r2), with at least one of the inverses different from 1,
then there exists in Fq(r1,r2) a family of permutation equivalent BCH multivariate codes {Cβ=Bq(β,γˉ,δˉ,bˉ)∣β∈U} such that
γˉ⊆{1,2}* and*
(a)
If supp(a)=Zr1 then 1∈γˉ.
2. (b)
If supp(b)=Zr2 then 2∈γˉ.
2. 2.
δˉ={δk∣k∈γˉ}* with δ1=Δ(M(a)) and δ2=Δ(M(b)).*
3. 3.
k∈γˉ∏δk=Δ(Cβ)=d(Cβ), for each
βˉ∈U.
4. 4.
φβ1,X1h1a−1⋅φβ2,X2h2b−1=φβ,X1h1aX2h2b−1∈Cβ, where β=(β1,β2).
Proof. Set bˉ={bk∣k∈γˉ}. First, suppose supp(M(a))=supp(a)=Zr1 and supp(M(b))=supp(b)=Zr2. Then by [9, Theorem 2], for each
β∈U, there exist BCH cyclic codes B(a)=Bq(β1,δ1,b1) and B(b)=Bq(β2,δ2,b2) such that δ1=Δ(M(a)), δ2=Δ(M(b)), φβ1,X1h1a−1∈B(a) and φβ2,X2h2b−1∈B(b). Now let C=Cβ be the abelian code with Dβ(C)=Dβ1(B(a))×Dβ2(B(b)) and set M=M(Dβ(C)). It is clear that M is a CP-matrix and, moreover, following the notation of Lemma 49, D1=Dβ1(B(a)) and D2=Dβ2(B(b)). Then C is a bivariate BCH code with γˉ={1,2}, δˉ={δ1,δ2} and b={b1,b2}. Now by statement (1) of Lemma 49, we have Δβ(C)=δ1δ2 and clearly statement * (5)* of this theorem holds.
It remains to prove the equality Δ(Cβ)=d(Cβ). On the one hand, we have Δβ(C)=Δ(M(a))Δ(M(b))=Δ(M(ab)) and, on the other hand, by hypothesis, Δ(M(ab))=Z(ab). Hence, by applying Theorem 25, we are done.
We may take, again, advantage from cyclic codes to transform a given abelian code C=(g), with d(C)=Δ(C) into another abelian code with higher dimension, as in Example 44, until to get a new BCH code.
Example 51
We continue with the code C from Example 43. Recall that q=2, r1=3, r2=45 and we have fixed α1∈U3 and α2∈U45. We have polynomials a=X+1 and b=Y40+Y39+Y38+Y36+Y35+Y32+Y30+Y25+Y24+Y23+Y21+Y20+Y17+Y15+Y10+Y9+Y8+Y6+Y5+Y2+1 such that a∣X13−1 and h1=1 works, and b∣X245−1 and h2=5 works in the sense of Theorem 50. Hence the hypothesis of this theorem are satisfied.
Now we have to follow the proof to construct our multivariate BCH code. The proof of Theorem 50 uses a construction from [9, Theorem 2]. Clearly B(a) is the BCH code in F2(3) with defining set Dα1(B(a))=C2(1). On the other hand, by [9, Example 8], the code B(b) has defining set Dα2(B(b))=C2(1)∪C2(3) so that it is a BCH code. In fact, B(b)=B2(α2,5,1), following the usual notation for BCH codes.
Thus C(α1,α2)=B2((α1,α2),{1,2},{2,5},{1,1}), d(C(α1,α2))=10 and dimF2(C(α1,α2))=58. This code has better parameters than the code C from Example 43 and C′ from Example 44.
We finish by extending Corollary 45 to bivariate BCH codes.
Corollary 52
Let K be an intermediate field Fq⊆K⊆L and a=a(X1)∈L(r1) be
such that a∣X1r1−1, with φα1,X1h1a−1∈Fq(r1), for some α1∈Ur1
and h1∈Zr1.
Let g be an irreducible factor of X2r2−1 in K[X2] with defining set Dα2(g),
for some α2∈Ur2. Set b=(X2n−1)/h. If there are positive integers j,t such that
b(α2j)=α2t and gcd(j,gcd(q−1,r2)r2)∣t then
there exists a bivariate BCH code C=Bq(α,γˉ,δˉ,bˉ)
in Fq(r1,r2) verifying Δ(M(ab))=Δ(C)=d(C), for certain α,γˉ,δˉ,bˉ.
Proof. Comes immediately from Corollary 45 together with Theorem 50.
Example 53
We shall extend to bivariate BCH codes those abelian codes on Table 1 from
Example 47. Recall that we had a list of divisors ai of X17−1 in
F2[X1] and divisors bj of X215−1 in F2[X2], namely, a1=1+X1, a2=1+X1+X13,
a2=1+X12+X13, b1=1+X2+X22X215−1, b2=1+X2+X24X215−1 and
b3=1+X23+X24X215−1 from which we constructed the mentioned table.
Now, one may check easily that the codes determined by a2, a3 and a1a3 are all BCH. Specifically,
B(a2)=B2(α1,4,5), B(a3)=B2(α1,4,0) and B(a1a3)=B2(α1,3,5), while the code determined by a2a3 is all F2(7).
On the other hand, as it is shown in [9, Example 9] one may construct from b1 the code B(b1)=B2(α1,2,0), of dimension 14, from b2 the code B2(α1,4,13) of dimension 10, and from b3 the code B2(α1,4,0) which also has dimension 10.
Thus, Table 3 is the new table of bivariate BCH codes in F2(7,15):
In the case of codes in F2(5,21), the code determined by a=Φ5 is F2(5), so we construct B(b1′)=B2(α2,2,0), B(b2′)=B2(α2,3,19), B(b3′)=B2(α2,3,1), B(b4′)=B2(α2,6,17) and B(b5′)=B2(α2,6,0). Table 4 is the new table of bivariate BCH codes in F2(5,21).
8 Conclusion
We have developed a technique to extend any bound for the minimum distance of cyclic codes constructed from its defining sets (ds-bounds) to abelian (or multivariate) codes through the notion of B-apparent distance. We used this technique to improve the searching for new bounds for abelian codes having unknown minimum distance. We have also studied conditions for an abelian code to verify that its B-apparent distance reaches its (true) minimum distance and we have constructed some tables of such codes as an application.
Bibliography25
The reference list from the paper itself. Each links out to its DOI / PubMed record.
1[1] E.R. Berlekamp and J. Justesen, Some long cyclic linear binary codes are not so bad . IEEE Transactions on Information Theory 20 (1974), 351-356.
2[2] J.J. Bernal, J.J. Simón, Partial permutation decoding for abelian codes . IEEE Transactions on Information Theory 59 (8) (2013), 5152-5170.
3[3] J.J. Bernal, D.H. Bueno-Carreño, J.J. Simón, Apparent distance and a notion of BCH multivariate codes . IEEE Transactions on Information Theory, 62 (2) (2016), 655-668.
4[4] J.J. Bernal, M. Guerreiro, J.J. Simón, Ds-bounds for cyclic codes: new bounds for abelian codes . Proceedings of ISITA 2016, Monterey, CA, USA, 712-716.
5[5] E. Betti, M. Sala, A new bound for the minimum distance of a cyclic code from its defining set . IEEE Transactions on Information Theory 52 (8) (2006) 3700-3706.
6[6] R.E. Blahut, Decoding of cyclic codes and codes on curves . In W.C. Huffman and V. Pless (Eds.), Handbook of Coding Theory . Vol. II, 1569-1633, 1998.
7[7] A. E. Brouwer, T. Verhoeff, An updated table of the minimum-distance bounds for binary linear codes . IEEE Transactions on Information Theory 39 (2) (1993) 662–677.
8[8] D. H. Bueno-Carreño, J.J. Bernal and J.J. Simón, A characterization of cyclic codes whose minimum distance equals their maximum BCH bound . Proceedings ACA 2013, Málaga, 109-113.