The Facets of the Bases Polytope of a Matroid and Two Consequences ()

Brahim Chaourar^{}

Department of Mathematics and Statistics, Al Imam Mohamed Ibn Saud University (IMSIU), Riyadh, Saudi Arabia.

**DOI: **10.4236/ojdm.2018.81002
PDF HTML XML
539
Downloads
846
Views
Citations

Department of Mathematics and Statistics, Al Imam Mohamed Ibn Saud University (IMSIU), Riyadh, Saudi Arabia.

Let *M* be a
matroid defined on a finite set *E* and *L* ⊂ *E* . *L* is
locked in *M* if and are
2-connected, and . In this paper, we prove that the nontrivial facets
of the bases polytope of *M* are
described by the locked subsets. We deduce that finding the maximum-weight
basis of *M* is a polynomial time
problem for matroids with a polynomial number of locked subsets. This class of
matroids is closed under 2-sums and contains the class of uniform matroids, the
Vámos matroid and all the excluded minors of 2-sums of uniform matroids. We
deduce also a matroid oracle for testing uniformity of matroids after one call
of this oracle.

Keywords

Bases Polytope, Facets, Locked Subsets, Maximum-Weight Basis Problem, Polynomially Locked Matroids, Matroid Oracle, Testing Unformity of a Matroid

Share and Cite:

Chaourar, B. (2018) The Facets of the Bases Polytope of a Matroid and Two Consequences. *Open Journal of Discrete Mathematics*, **8**, 14-20. doi: 10.4236/ojdm.2018.81002.

1. Introduction

Sets and their characteristic vectors will not be distinguished. We refer to Oxley [1] and Schrijver [2] about, respectively, matroids and polyhedra terminology and facts.

Let M be a matroid defined on a finite set E. $\mathcal{I}\left(M\right)$ , $\mathcal{B}\left(M\right)$ and the function r are, respectively, the class of independent sets, bases and the rank function of M. ${M}^{\mathrm{*}}$ , ${\mathcal{B}}^{\mathrm{*}}\left(M\right)$ and the function ${r}^{\mathrm{*}}$ are, respectively, the dual matroid, the class of cobases and the dual rank function of M. For any $X\subset E$ , $\mathcal{B}\left(X\right)$ and ${\mathcal{B}}^{\mathrm{*}}\left(X\right)$ are, respectively the class of bases of $M\mathrm{|}X$ and cobases of ${M}^{\mathrm{*}}\mathrm{|}X$ . The polyhedra $Q\left(M\right)$ and $P\left(M\right)$ are, respectively, the convex hulls of the independent sets and the bases of M. Suppose that M (and ${M}^{\mathrm{*}}$ ) is 2-connected. A subset $L\subset E$ is called a locked subset of M if $M\mathrm{|}L$ and ${M}^{*}|\left(E\backslash L\right)$ are 2-connected, and their corresponding ranks are at least 2, i.e., $\mathrm{min}\left\{r\left(L\right),{r}^{*}\left(E\backslash L\right)\right\}\ge 2$ . It is not difficult to see that if L is locked then both L and $E\backslash L$ are closed, respectively, in M and ${M}^{\mathrm{*}}$ (That is why we call it locked). We denote by $\mathcal{L}\left(M\right)$ and $\mathcal{l}\left(M\right)$ , respectively, the class of locked subsets of M and its cardinality, which is called the locked number of M. Given a positive integer k (k does not depend on M or $\left|E\right|$ ), the class of k-locked matroids, denoted by ${\mathcal{L}}_{k}$ , is the class of matroids M such that $\mathcal{l}\left(M\right)\in O\left({\left|E\right|}^{k}\right)$ . ${\mathcal{L}}_{0}$ is the class of matroids M such that $\mathcal{L}\left(M\right)=\varnothing $ , i.e., $\mathcal{l}\left(M\right)=0$ . For a given nonegative integer k, ${\mathcal{L}}_{k}$ is called also a polynomially locked class of matroids and its members are called k-locked or polynomially locked matroids. It is not difficult to see that the class of lockeds subsets of a matroid M is the union of lockeds subsets of the 2-connected components of M. The locked structure of M is the quadruple ( $\mathcal{P}\left(M\right)$ , $\mathcal{S}\left(M\right)$ , $\mathcal{L}\left(M\right)$ , $\rho $ ), where $\mathcal{P}\left(M\right)$ and $\mathcal{S}\left(M\right)$ are, respectively, the class of parallel and coparallel closures, and $\rho $ is the rank function restricted to $\mathcal{P}\left(M\right)\cup \mathcal{S}\left(M\right)\cup \mathcal{L}\left(M\right)\cup \left\{\varnothing \mathrm{,}E\right\}$ .

Given a weight function $c\in {R}^{E}$ , the maximum-weight basis problem (MWBP) is the following optimization problem:

Maximize{ $c\left(B\right)$ such that $B\in \mathcal{B}\left(M\right)$ }

The corresponding maximum-weight independent problem is clearly (polynomially time) equivalent to MWBP. MWBP is polynomial on $\left|E\right|$ and $\theta $ , where $\theta $ is the complexity of the used matroid oracle [3] . Even if we use the approach introduced by Mayhew [4] by giving the list of bases (for example) in the input, MWBP is polynomial on the size of the input. However, as Robinson and Welsh [5] note, no matter which of the ways to specify a matroid, the size of the input for a matroid problem on an n-element set is $O\left({2}^{n}\right)$ . It follows that MWBP is not polynomial in its strict sense, that is on $\left|E\right|$ . We prove that MWBP is polynomial on $\left|E\right|$ for polynomially locked classes of matroids, i.e., for any matroid $M\in {\mathcal{L}}_{k}$ (for a fixed k). This class of polynomially locked matroids is closed under 2-sums and contains the class of uniform matroids, the Vámos matroid and all the excluded minors of 2-sums of uniform matroids. These excluded minors are $M\left({K}_{4}\right)$ , ${W}^{3}$ , ${Q}_{6}$ and ${P}_{6}$ [6] . It follows that this class is larger than 2-sums of uniform matroids.

Testing Uniformity of matroids (TUM) is to provide an algorithm in which the matroid is represented by an oracle and which decides whether the given matroid is uniform or not after a number of calls on the oracle which is bounded by a polynomial in the size of the ground set. Jensen and Korte [7] proved that there exists no such algorithm in which the matroid is represented by an independence test oracle (or an oracle polynomially related to an independence test oracle). In this paper, we give a matroid oracle which answers this question.

The remainder of the paper is organized as follows: in Section 2, we give all facets of the bases polytope, then, in Section 3, we deduce two consequences of this characterization. The first one is that MWBP is polynomial (time) for polynomially locked classes of matroids, and the second one is a polynomial time algorithm via a new matroid oracle for testing if a given matroid is uniform or not. In Section 4, we describe some polynomially locked classes of matroids, and we conclude in Section 5.

2. Facets of the Bases Polytope

A description of $Q\left(M\right)$ was given by Edmonds [3] as follows.

Theorem 2.1. $Q\left(M\right)$ is the set of all $x\in {R}^{E}$ such that

$x\left(e\right)\ge 0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}\text{any}\text{\hspace{0.17em}}e\in E$ (1)

$x\left(A\right)\le r\left(A\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}\text{any}\text{\hspace{0.17em}}A\subseteq E$ (2)

Later, a minimal description of $Q\left(M\right)$ was given also by Edmonds [8] as follows.

Theorem 2.2. The inequality (2) is a facet of $Q\left(M\right)$ if and only if A is closed and 2-connected.

It is not difficult to see that $P\left(M\right)$ is the set of all $x\in {R}^{E}$ satisfying the inequalities (1), (2) and

$x\left(E\right)=r\left(E\right)$ (3)

It seems natural to think that the inequality (2) is a facet of $P\left(M\right)$ if and only if A is closed and 2-connected. This is not true because:

Lemma 2.3. If the inequality (2) is a facet of $P\left(M\right)$ then A is a locked subset of M.

Proof. It suffices to prove that if X is closed and 2-connected but $E\backslash L$ is not 2-connected in the dual then the inequality (2) is not a facet. In fact, there exist A and B two disjoint subsets of E such that $E\backslash X=A\cup B$ and ${r}^{\mathrm{*}}\left(E\backslash X\right)={r}^{\mathrm{*}}\left(A\right)+{r}^{\mathrm{*}}\left(B\right)$ , that is, $\left|E\backslash X\right|-r\left(E\right)+r\left(X\right)=\left|A\right|-r\left(E\right)+r\left(E\backslash A\right)$ $+\left|B\right|-r\left(E\right)+r\left(E\backslash B\right)$ . It follows that: $r\left(E\right)+r\left(X\right)=r\left(E\backslash A\right)+r\left(E\backslash B\right)\ge $ $x\left(E\backslash A\right)+x\left(E\backslash B\right)=x\left(E\right)+x\left(X\right)$ , which implies the inequality (2). So the inequality (2) is redundant and cannot be a facet. ,

We give now a minimal description of $P\left(M\right)$ . A part of the proof is inspired from a proof given by Pulleyblank [9] to describe the nontrivial facets of $Q\left(M\right)$ . Independently, Fujishige [10] , and Feichtner and Sturmfels [11] , gave a characterization of nontrivial facets of $P\left(M\right)$ . We give here below a new and complete formulation with a new proof.

Theorem 2.4. A minimal description of $P\left(M\right)$ is the set of all $x\in {R}^{E}$ satisfying the equality (3) and the following inequalities:

$x\left(P\right)\le 1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}\text{any}\text{\hspace{0.17em}}\text{parallel}\text{\hspace{0.17em}}\text{closure}\text{\hspace{0.17em}}\text{\hspace{0.17em}}P\subseteq E$ (4)

$x\left(S\right)\ge \left|S\right|-1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}\text{any}\text{\hspace{0.17em}}\text{coparallel}\text{\hspace{0.17em}}\text{closure}\text{\hspace{0.17em}}S\subseteq E$ (5)

$x\left(L\right)\le r\left(L\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}\text{any}\text{\hspace{0.17em}}\text{locked}\text{\hspace{0.17em}}\text{subset}\text{\hspace{0.17em}}L\subseteq E$ (6)

Proof. Without loss of generality, we can suppose that M is without parallel or coparallel elements so the inequalities (5) become as (1) and the inequalities (4) become the following:

$x\mathrm{(}e\mathrm{)}\le 1\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}\text{any}\text{\hspace{0.17em}}e\in E$ (7)

Let $C\left(M\right)$ be the cone generated by the incidence vectors of the bases of M. It suffices to prove that the minimal description of $C\left(M\right)$ is given by (1) and the following inequalities:

$x\left(e\right)\le x\left(E\right)/r\left(E\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}\text{any}\text{\hspace{0.17em}}e\in E$ (8)

$x\left(L\right)/r\left(L\right)\le x\left(E\right)/r\left(E\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}\text{any}\text{\hspace{0.17em}}\text{locked}\text{\hspace{0.17em}}\text{subset}\text{\hspace{0.17em}}L\subseteq E$ (9)

It is not difficult to see via induction and operations of deletion and contraction that the inequalities (1) and (8) are facets of $C\left(M\right)$ . It remains to prove that the inequality (9) is a facet of $C\left(M\right)$ if and only if L is a locked subset of M. According to Lemma 2.3, it suffices to prove the inverse way. Note that (9) is equivalent to the following inequality:

$\left(r\left(L\right)-r\left(E\right)\right)x\left(L\right)+r\left(L\right)x\left(E\backslash L\right)\ge 0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}\text{any}\text{\hspace{0.17em}}\text{locked}\text{\hspace{0.17em}}\text{subset}\text{\hspace{0.17em}}L\subseteq E$ (10)

Let $ax\ge 0$ be a valid inequality for $C\left(M\right)$ which is tight for all $B\in \mathcal{B}\left(L\right)$ .

Claim 1: ${a}_{j}={a}_{k}$ for all j and k of L.

Suppose this is not true. Let $X=\left\{j\in L\text{\hspace{0.17em}}\text{such}\text{\hspace{0.17em}}\text{that}\text{\hspace{0.17em}}{a}_{j}\text{\hspace{0.17em}}\text{takes}\text{\hspace{0.17em}}\text{minimum}\text{\hspace{0.17em}}\text{over}\text{\hspace{0.17em}}L\right\}$ , $Y=L\backslash X$ and $B\in \mathcal{B}\left(L\right)\cap \mathcal{B}\left(Y\right)$ . Since L is 2-connected in M, and since, by assumption, X is a strict subset of L, then $r\left(X\right)>\left|B\cap X\right|$ . Thus there exists $e\in X\backslash B$ such that $\left(B\cap X\right)\cup e$ is an independent set of M. It follows that there exists $f\in B\cap Y$ such that $\stackrel{\u02dc}{B}=\left(B\backslash f\right)\cup e\in \mathcal{B}\left(L\right)$ . But: $a\left(\stackrel{\u02dc}{B}\right)=a\left(B\right)$ $-a\left(f\right)+a\left(e\right)<a\left(B\right)$ , a contradiction.

Claim 2: For any $X\subseteq E$ , $B\in \mathcal{B}\left(X\right)$ if and only if $E\backslash B\in {\mathcal{B}}^{\mathrm{*}}\left(E\backslash X\right)$ .

It suffices to prove one way and use duality for the other way.

Let $B\in \mathcal{B}\left(X\right)$ then

$\begin{array}{c}\left|B\cap X\right|=r\left(X\right)=\left|X\right|-{r}^{*}\left(E\right)+{r}^{*}\left(E\backslash X\right)\\ =\left|E\right|-\left|E\backslash X\right|-{r}^{*}\left(E\right)+{r}^{*}\left(E\backslash X\right)\\ =r\left(E\right)-\left|E\backslash X\right|+{r}^{*}\left(E\backslash X\right).\end{array}$

Thus,

$\begin{array}{c}\left|\left(E\backslash B\right)\cap \left(E\backslash X\right)\right|=\left|E\backslash X\right|-\left|B\cap \left(E\backslash X\right)\right|\\ =\left|E\backslash X\right|-\left|B\right|+\left|B\cap X\right|\\ =\left|E\backslash X\right|-\left|B\right|+r\left(E\right)-\left|E\backslash X\right|+{r}^{*}\left(E\backslash X\right)\\ ={r}^{*}\left(E\backslash X\right).\end{array}$

Since $E\backslash B$ is a basis in the dual, then $E\backslash B\in {\mathcal{B}}^{*}\left(E\backslash X\right)$ .

Claim 3: ${a}_{j}={a}_{k}$ for all j and k of $E\backslash L$ .

Using claim 2, $E\backslash L$ being 2-connected and a similar argument on $E\backslash B$ as in claim 1, we conclude.

Claim 4: $ax\ge 0$ is a multiple of inequality (10).

By claims 1 and 3, $ax\ge 0$ becomes: ${a}_{L}x\left(L\right)+{a}_{E\backslash L}x\left(E\backslash L\right)\ge 0$ . Thus, for $B\in \mathcal{B}\left(L\right)$ , we have:

${a}_{L}\left|B\cap L\right|+{a}_{E\backslash L}\left|B\cap \left(E\backslash L\right)\right|=\mathrm{0,}$

that is,

${a}_{L}r\left(L\right)+{a}_{E\backslash L}\left(r\left(E\right)r\left(L\right)\right)=0.$

But ${a}_{L}=r\left(L\right)r\left(E\right)$ and ${a}_{E\backslash L}=r\left(L\right)$ is a solution of this equation, so we conclude. ,

3. MWBP and TUM

Since the bases polytope is completly described by the locked structure of the matroid, so a natural matroid oracle follows.

The k-locked oracle

Input: a nonegative integer k and a matroid M defined on E.

Output: 1) No if $M\notin {\mathcal{L}}_{k}$ , and

2) ( $\mathcal{P}\left(M\right)$ , $\mathcal{S}\left(M\right)$ , $\mathcal{L}\left(M\right)$ , $\rho $ ) if $M\in {\mathcal{L}}_{k}$ .

Note that this oracle has time complexity $O\left({\left|E\right|}^{k+1}\right)$ because we need to count at most ${\left|E\right|}^{k+1}$ members of $\mathcal{L}\left(M\right)$ in order to know that M is not k-locked, even if the memory complexity can be $O\left(\left|E\right|+\mathcal{l}\left(M\right)\right)$ . Actually this matroid oracle permits to recognize if a given matroid is k-locked or not for a given nonegative integer k (which does not depend on M or $\left|E\right|$ ).

The first consequence of Theorem 2.4 then follows.

Corollary 3.1. Given a nonegative integer k, a matroid $M\in {\mathcal{L}}_{k}$ , the k-locked oracle to acess M and a weight function $c\in {R}^{E}$ . Then there exists a polynomial time algorithm on the size of $E\left(M\right)$ for solving MWBP in M.

Proof. Let M be a such matroid. Since $M\in {\mathcal{L}}_{k}$ then it can be described by its locked structure in the input of MWBP by using the k-locked oracle. MWBP is equivalent for optimizing on $P\left(M\right)$ , which is also equivalent to separating on $P\left(M\right)$ . Since the number of facets of $P\left(M\right)$ is $2\left|E\right|+\mathcal{l}\left(M\right)$ then separating can be done on $O\left(\left|E\right|+\mathcal{l}\left(M\right)\right)$ . But M is polynomially locked, hence $\mathcal{l}\left(M\right)\in $ $O\left({\left|E\right|}^{k}\right)$ and separating on $P\left(M\right)$ can be done on $O\left({\left|E\right|}^{k}\right)$ . ,

The k-locked oracle is stronger than the rank and the independence oracles for polynomially locked matroids:

We can get the rank of any subset $X\subseteq E$ by choosing the weight function c equal to the characteristic vector of X and optimizing on $P\left(M\right)$ , which can be done in polynomial time. The obtained optimum value of c is the requested rank.

For the independence oracle, we can decide if a subset $X\subseteq E$ is independent or not by choosing the same previous weight function and decide that X is independent if it is included in the optimum basis, and not independent otherwise.

For testing uniformity of matroids, we need the following result [6] .

Theorem 3.2. A 3-connected matroid M is uniform if and only if $\mathcal{l}\left(M\right)=0$ .

We can see through the proof of this theorem that $\mathcal{l}\left(M\right)=0$ if M is uniform whatever its connectivity. For disconnected matroids, we have the following result [1] .

Theorem 3.3. A disconnected matroid M is uniform if and only if $r\left(M\right)=\left|E\right|$ or $r\left(M\right)=0$ .

For 2-connected matroids, we can write:

Proposition 3.4. A 2-connected matroid M is uniform if and only if one of the following properties holds:

i) $\mathcal{l}\left(M\right)=0$ and $\left|\mathcal{P}\left(M\right)\right|=\left|E\right|=\left|\mathcal{S}\left(M\right)\right|$ ;

ii) $\left|\mathcal{P}\left(M\right)\right|=1$ ;

iii) $\left|\mathcal{S}\left(M\right)\right|=1$ .

So we can now characterize uniform matroids as follows.

Corollary 3.5. M is uniform if and only if one of the following properties holds:

i) $\mathcal{l}\left(M\right)=0$ and $\left|\mathcal{P}\left(M\right)\right|=\left|E\right|=\left|\mathcal{S}\left(M\right)\right|$ ;

ii) $\left|\mathcal{P}\left(M\right)\right|=1$ ;

iii) $\left|\mathcal{S}\left(M\right)\right|=1$ ;

iv) $r\left(M\right)=\left|E\right|$ ;

v) $r\left(M\right)=0$ .

A natural matroid oracle follows.

The locked number oracle

Input: a matroid M defined on E.

Output: $\mathcal{l}\left(M\right)$ , $r\left(M\right)$ , $\left|\mathcal{P}\left(M\right)\right|$ , $\left|\mathcal{S}\left(M\right)\right|$ .

Note that, except for $\mathcal{l}\left(M\right)$ , all other outputs of this oracle can be computed in a polynomial time given a locked structure of M. We can now give an algorithm which tests if a given matroid is uniform after one call of the locked number oracle.

Testing Uniformity of Matroids

Input: a matroid M defined on E.

Output: a) M is uniform if one of the following properties holds:

i) $\mathcal{l}\left(M\right)=0$ and $\left|\mathcal{P}\left(M\right)\right|=\left|E\right|=\left|\mathcal{S}\left(M\right)\right|$ ;

ii) $\left|\mathcal{P}\left(M\right)\right|=1$ ;

iii) $\left|\mathcal{S}\left(M\right)\right|=1$ ;

iv) $r\left(M\right)=\left|E\right|$ ;

v) $r\left(M\right)=0$ .

b) Else, M is not uniform.

4. Some Polynomially Locked Matroids

Since 2-sums preserve k-lockdness for $k\ge 1$ [12] , $\mathcal{l}\left(M\left({K}_{4}\right)\right)=4$ , $\mathcal{l}\left({W}^{3}\right)=3$ , $\mathcal{l}\left({Q}_{6}\right)=2$ , $\mathcal{l}\left({P}_{6}\right)=1$ , $\mathcal{l}\left({V}_{8}\right)=5$ , then we can say:

Theorem 4.1. If $k\ge 1$ then ${\mathcal{L}}_{k}$ is closed under 2-sums, contains all the excluded minors of 2-sums of uniform matroids and the Vámos matroid.

In particular, ${\mathcal{L}}_{1}$ is closed under 2-sums and contains all the above matroids.

It follows that ${\mathcal{L}}_{1}$ contains strictly 2-sums of uniform matroids.

5. Conclusion

We have given a complete description of all facets of the bases polytope of a matroid and deduce two consequences. One about MWBP and the second about TUM. Future investigations can be characterizing some or all polynomially locked classes of matroids.

Acknowledgements

The author is grateful to anonymous referees in a previous version of this paper.

Conflicts of Interest

The authors declare no conflicts of interest.

[1] | Oxley, J.G. (1992) Matroid Theory. Oxford University Press, Oxford. |

[2] | Schrijver, A. (1986) Theory of Linear and Integer Programming. John Wiley and Sons, Chichester. |

[3] |
Edmonds, J. (1971) Matroids and the Greedy Algorithm. Mathematical Programming, 1, 127-136. https://doi.org/10.1007/BF01584082 |

[4] |
Mayhew, D. (2008) Matroid Complexity and Nonsuccinct Descriptions. SIAM Journal on Discrete Mathematics, 22, 455466. https://doi.org/10.1137/050640576 |

[5] |
Robinson, G.C. and Welsh, D.J.A. (1980) The Computational Complexity of Matroid Properties. Mathematical Proceedings of the Cambridge Philosophical Society, 87, 29-45. https://doi.org/10.1017/S0305004100056498 |

[6] |
Chaourar, B. (2011) A Characterization of Uniform Matroids. ISRN Algebra, 2011, Article ID: 208478. https://doi.org/10.5402/2011/208478 |

[7] |
Jensen, P.M. and Korte, B. (1982) Complexity of Matroid Property Algorithms. SIAM Journal on Computing, 11, 184-190. https://doi.org/10.1137/0211014 |

[8] | Giles, R. (1975) Submodular Functions, Graphs and Integer Polyhedra. PhD Thesis, University of Waterloo, Waterloo. |

[9] | Pulleyblank, W.R. (1989) Polyhedral Combinatorics. In: Nemhauser, G.L., et al., Eds., Handbooks in Operations Research and Management Science, Elsevier, Amsterdam, 371-446. |

[10] |
Fujishige, S. (1984) A Characterization of Faces of the Base Polyhedron Associated with a Sub Modular System. Journal of the Operations Research Society of Japan, 27, 112-128. https://doi.org/10.15807/jorsj.27.112 |

[11] | Feichtner, E.M. and Sturmfels, B. (2005) Matroid Polytopes, Nested Sets and Bergman Fans. Portugaliae Mathematica, 62, 437-468. |

[12] |
Chaourar, B. (2008) On the Kth Best Basis of a Matroid. Operations Research Letters, 36, 239-242. https://doi.org/10.1016/j.orl.2007.05.007 |

Journals Menu

Copyright © 2021 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.