Stability analysis of solutions to equilibrium problems and applications in economics

Tran Ngoc Tam, Nguyen Minh Hai, Bantaojai Thanatporn

Asian Journal of Economics and Banking

ISSN: 2615-9821

Open Access. Article publication date: 2 December 2020

Issue publication date: 18 December 2020

842

Abstract

Purpose

The purpose of this paper is to study the Hölder calmness of solutions to equilibrium problems and apply it to economics.

Design/methodology/approach

The authors obtain the Hölder calmness by using an effective approach. More precisely, under the key assumption of strong convexity, sufficient conditions for the Hölder continuity of solution maps to equilibrium problems are established.

Findings

A new result in stability analysis for equilibrium problems and applications in economics is archived.

Originality/value

The authors confirm that the paper has not been published previously, is not under consideration for publication elsewhere and is not being simultaneously submitted elsewhere.

Keywords

Citation

Tam, T.N., Hai, N.M. and Thanatporn, B. (2020), "Stability analysis of solutions to equilibrium problems and applications in economics", Asian Journal of Economics and Banking, Vol. 4 No. 3, pp. 121-132. https://doi.org/10.1108/AJEB-06-2020-0023

Publisher

:

Emerald Publishing Limited

Copyright © 2020, Tran Ngoc Tam, Nguyen Minh Hai and Bantaojai Thanatporn.

License

Published in Asian Journal of Economics and Banking. Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) licence. Anyone may reproduce, distribute, translate and create derivative works of this article (for both commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this licence may be seen at http://creativecommons.org/licences/by/4.0/legalcode


1. Introduction

Many important problems such as optimization problems, variational inequality problems, complementarity problems, Nash equilibrium problems, minimax problems, fixed-point and coincidence-point problems and traffic network problems are considered as special cases of an equilibrium problem (). This implies possibility of a wide application of results in the equilibrium problem theory to several important fields, including economics, physics (especially, mechanics), engineering, transportation, sociology, chemistry, biology, etc. ().

Existence conditions of solutions for the equilibrium problem and its generalizations are the first and most developed topic. There are a lot of works devoted to this topic in the literature, see e.g. , , ,, , and, and the references therein. The second one is stability analysis of solutions. Stability conditions, especially the lower semicontinuity and Hölder continuity of solution maps, for such problems have been extensively studied and received increasing attention from many researchers so far (; Anh and Khanh, ; Anh and Khanh, ; Kimura and Yao, ; Kimura and Yao, 2013; Anh et al., 2018; and references therein).

The paper aims at investigating the stability analysis in the sense of Hölder calmness of the solution maps to equilibrium problems. To be more precise, sufficient conditions for the solution maps to be Hölder calm are established. At the end of the paper, we present applications of the main results in economics. Namely, the Hölder calmness of solution maps to mean-variance portfolio and Nash equilibrium problems is derived.

The rest of the paper is organized as follows. Section 2 states the equilibrium problem setting and recalls some definitions and their properties needed in what follows. Sufficient conditions for the Hölder calmness of the solution maps are established in Section 3. In Section 4, we discuss the Hölder calmness of solution maps to mean-variance portfolio and Nash equilibrium problems.

2. Preliminaries

In this paper, we use d(·,·) for the metric in metric spaces. For two subsets A,B⊆X, we use the following notations:

d(a,B):=infbBd(a,b),
H*(A,B):=supaAd(a,B),
H(A,B):=max{H*(A,B),H*(B,A)},
ρ(A,B):=supaA,bBd(a,b).

Denote ℝ+ is the set of the nonnegative real numbers and B(x,r) is the closed ball of radius r ≥ 0 and is centered at x. intA, conv(A) and diam A=sup⁡x,z∈Ad(x,z) stand for the interior, the convex hull and the diameter, respectively, of a subset A. For a set-valued map G:XY, gph G:={(x,y)X×Y:yG(x)} is the graph of G. Recall that X is called a metric linear space if and only if it is both a metric space and a linear space and the metric d of X is translation invariant (i.e. d(x+z,y+z)=d(x,y),∀x,y,zX) and, for any convergent sequences (λm) in ℝ and (xm) in X, we have limm(λmxm)=(limmλm)(limmxm).

From now on, unless otherwise stated, let X be a metric linear space, Λ,M be metric spaces and AX be nonempty. Let K:Λ⇉A be nonempty-convex-valued and φ:A×A×M→ℝ(λ,μ)∈Λ×M, we consider the following parametric equilibrium problem, i.e. the family of the corresponding problems when (λ,μ) varies in Λ×M,

(EP) find x¯∈K(λ) such that φ(x¯,y,μ)≥0 for all yK(λ).

For each (λ,μ)Λ×M, denote the solution set of (EP) by S(λ,μ), i.e.:

S(λ,μ):={xK(λ):φ(x,y,μ)0,yK(λ)}.

As existence of solutions has been studied much in the literature, we do not include existence investigations and always assume that S(λ,μ) is nonempty in the neighborhood of the considered point.

We first recall some notions needed in the sequel.

Definition 2.1. Let nγ>0 and θ≥0. It is said that:

  • a function g:X→ℝ is nγ-Hölder continuous at x¯∈X if there is a neighborhood U of x¯ such that, for all x1,x2∈U:

    d(g(x1),g(x2))ndγ(x1,x2);

  • a function g:X→ℝ is n.γ-Hölder calm at x¯∈X if there is a neighborhood U of x¯ such that, for all xU:

    d(g(x),g(x¯))ndγ(x,x¯);

  • a function g:X×X×M→ℝ is n.γ-Hölder calm at μ¯∈M, θ-uniformly over a subset BX if there is a neighborhood U ofμ¯ such that, for all μ∈U and x,yB:xy:

    d(g(x,y,μ),g(x,y,μ¯))ndγ(μ,μ¯)dθ(x,y);

if θ = 0, we say that g(x,y,·) is n.γ-Hölder calm at μ¯, uniformly over B; and

  • a multifunction KΛ⇉X is n.γ-Hölder calm at λ¯∈Λ if there exists a neighborhood N of λ¯ such that, for all λ∈N:

    H(K(λ),K(λ¯))dγ(λ,λ¯).

If γ = 1, then Hölder continuity is called Lipschitz continuity.

We say that a certain property is satisfied in a subset BX if and only if it is satisfied at every point of B.

Definition 2.2. Let g:X→ℝ, BX, and h,β be positive.

  • g is called h.β-strongly convex on a convex subset B if and only if, for all x1,x2∈B and t∈ (0,1):

    g((1t)x1+tx2)(1t)g(x1)+tg(x2)ht(1t)dβ(x1,x2).

  • g is called h.β-strongly convex-like in B (B not necessarily convex) if and only if, for all x1,x2∈B and t∈ (0,1), there is zB such that:

    g(z)(1t)g(x1)+tg(x2)ht(1t)dβ(x1,x2).

Recall that a function g:X×X→ℝ is called monotone on BX if:

g(x,y)+g(y,x)0,x,yB.

In what follows, we use the following assumptions that play an important role in investigating conditions for the Hölder calmness of solution maps to the equilibrium problems:

H1.

K is 1.α-Hölder calm on Λ.

H2.

The map y↦φ(x,y,μ) is h.β-strongly convex as well as m.1-Hölder continuous in conv(K(Λ)).

H3.

The map (x,Y)↦φ(x,Y,μ) is monotone on K(Λ)×K(Λ) for all μ∈M.

H4.

The map μ↦φ(x,Y,μ) is n.γ-Hölder calm on M, θ-uniformly over K(Λ).

3. Hölder calmness of solutions

In this section, we state the main results of the paper. Namely, we establish sufficient conditions for the Hölder calmness of the solution maps to equilibrium problems.

Theorem 3.1. Assume that all assumptions (H1–H4) are satisfied with θ<β. Then, the solution map S is single-valued and Hölder calm on Λ × M.

Proof.

Let (λ¯,μ¯)∈Λ×M be arbitrary, we need to prove that S is Hölder calm at (λ¯,μ¯). We divide the proof into three steps.

Step 1. For all x11∈S(λ¯,μ¯) and x21∈S(λ,μ¯), we claim that:

(1) d1:=d(x11,x21)(4mlh)1βdαβ(λ¯,λ).

By the definition of the solution set, we have, for all yK(λ¯) and zK(λ):

(2) min{φ(x11,y,μ¯),φ(x21,z,μ¯)}0.

The Hölder calmness of K leads to the existence of x1∈K(λ¯) and x2∈K(λ) such that:

(3) max{d(x11,x2),d(x21,x1)}ldα(λ¯,λ).

Letting x^=12(x11+x21), it follows from the strong convexity in H2 that:

(4) φ(x11,x^,μ¯)12φ(x11,x11,μ¯)+12φ(x11,x21,μ¯)14hd1β.

By virtue of the monotonicity of φ and x11∈S(λ¯,μ¯), one gets:

f(x11,x11,μ¯)=0
f(x11,x21,μ¯)f(x21,x11,μ¯).

Hence, implies that:

(5) 14hd1β12φ(x21,x11,μ¯)φ(x11,x11+x212,μ¯).

Letting z = x2 and y=12(x11+x1) in , one has:

(6) min{12φ(x21,x2,μ¯),φ(x11,x11+x12,μ¯)}0

From and , we arrive at:

14hd1β12[φ(x21,x2,μ¯)φ(x21,x11,μ¯)]+[φ(x11,x11+x12,μ¯)φ(x11,x11+x212,μ¯)].

Hence, using the m.1-Hölder calmness in H2 and , we obtain:

14hd1β12md(x11,x2)+12md(x21,x1)
    12mldα(λ¯,λ)+12mldα(λ¯,λ),
that is:
14hd1βmldα(λ¯,λ).

Hence, we obtain (1).

Step 2. We will show that, for all x21∈S(λ,μ¯) and x22∈S(λ,μ):

(7) d2:=d(x21,x22)(nh)1βθdγβθ(μ¯,μ).

By the definition of the solution, one has, for all y,zK(λ),

(8) min{φ(x21,y,μ¯),φ(x22,z,μ)}0.

Putting y=12(x22+x21) in , we have:

φ(x21,x22+x212,μ¯)0.

By the strong convexity given in H2, we have:

(9) f(x21,x22,μ¯)+f(x21,x21,μ¯)12hd2β0.

The monotonicity of φ assumed in H3 and x21∈S(λ,μ¯) imply that:

φ(x21,x21,μ¯)=0
φ(x21,x22,μ¯)φ(x22,x21,μ¯).

Hence, combining with , we have:

(10) 12hd2βφ(x22,x21,μ¯).

Now, letting z=12(x22+x21) in derives:

φ(x22,x22+x212,μ)0.

Also by the strong convexity given in H2, we have:

φ(x22,x22,μ)+φ(x22,x21,μ)12hd2β0.

As φ(x22,x22,μ)=0, this yields that:

12hd2βφ(x22,x21,μ).

Adding and this inequality, and using assumption H4, we get:

hd2βf(x22,x21,μ)f(x22,x21,μ¯)ndγ(μ¯,μ)d2θ,
i.e. we obtain .

Step 3. We are now ready to complete the proof. For all x11∈S(λ¯,μ¯) and x22∈S(λ,μ), we have:

d(x11,x22)d1+d2.

Hence, from and , we get, with k1=(4mlh)1β and k2=(nh)1βθ:

ρ(S(λ¯,μ¯),S(λ,μ))k1dαβ(λ¯,λ)+k2dγβθ(μ¯,μ).

Letting λ¯=λ and μ¯=μ in this inequality, we see that the diameter of S(λ¯,μ¯) is 0 [for arbitrary (λ¯,μ¯)], i.e. the solution map of (EP) is single-valued in Λ × M. The proof is complete.

In the special case where K(λ)≡ (K is a nonempty set), we have the following result.

Theorem 3.2. For (EP) with K(λ)≡ (K, assume that assumptions H3 and H4 are satisfied and H2 is replaced by the following condition:

H′2: The map yφ(x,y,μ) is h.β-strongly convex-like on K.

Then, S is Hölder calm on Λ × M.

Proof. For any x1S(μ¯),x2S(μ) and y,zK, one has:

min{φ(x1,y,μ¯),φ(x2,z,μ)}0.

By the strong convex-likeness of f in K, there isz¯∈K such that:

(11) φ(x1,z¯,μ¯)(1t)φ(x1,x2,μ¯)+tφ(x1,x1,μ¯)ht(1t)dβ(x1,x2).

This and the monotonicity of φ imply that:

ht(1t)dβ(x1,x2)(1t)φ(x2,x1,μ¯)φ(x1,z¯,μ¯).

As φ(x1,z¯,μ¯)≥0, this inequality leads to:

(12) htdβ(x1,x2)φ(x2,x1,μ¯).

With the same arguments but with x1 replaced by x2 and μ¯ by μ in , we have:

(13) h(1t)dβ(x1,x2)φ(x2,x1,μ).

Adding and , and using assumption H3, one gets:

hdβ(x1,x2)φ(x2,x1,μ)φ(x2,x1,μ¯)
ndγ(μ¯,μ)dθ(x1,x2).

This means the required Hölder condition. From this, we obtain the uniqueness of the solution.

4. Applications

4.1 Mean-variance portfolio

In this subsection, we present an application of our results to a practical situation. Namely, we study Hölder properties of solution maps to parametric quadratic programing, a special case of the equilibrium problem, which has well-known applications in the view of practice. A good sample is the model of mean-variance portfolio based on investor’s utility maximization. This problem can be formulated as parametric quadratic optimization problems. Consider a universe of n assets with the known data is c=(c1,…,cn)T and Q=[σij], where ci is the expected return for asset i and σij is the covariance of returns for assets i and j. So, c is the vector of expected returns and Q is the n × n variance–covariance matrix of asset returns. Herein, Q is a positive definite matrix which follows from the properties of variance–covariance matrices. Denote the vector of asset holdings by x=(x1,…,xn)T. Then, the expected return of the portfolio x is cTx and its variance is σ2=xTQx. A portfolio is said to be efficient if for some fixed level of expected return no other portfolio gives smaller variance (). In other words, an efficient portfolio can be considered as the one for which at some fixed level of variance no other portfolio gives larger expected return. The determination of the efficient portfolio frontier in this mean-variance model is equivalent to solving the following parametric problem because of ():

(MP)minf(x,μ):=12xTQxμcTxs.t.Ax=b,x+n,
where μ ≥ 0 is an investor’s variance aversion parameter. The linear constraints Ax = b can represent budget constraints, bounds on asset holdings, etc. Nonnegativity constraints x ≥ 0 are short-sale constraints (non-negative asset holdings).

Lemma 4.1. If u(x) is h.β-strongly convex and v(x) is convex on a convex set A⊂ℝn, then g(x)=u(x)+v(x) is h.β-strongly convex on A.

Proof. For all x1,x2∈A and t∈[0,1], we have:

g(tx1+(1t)x2)=u(tx1+(1t)x2)+v(tx1+(1t)x2)
tu(x1)+(1t)u(x2)ht(1t)dβ(x1,x2)+tv(x1)
+(1t)v(x2)
t(u(x1)+v(x1))+(1t)(u(x2)+v(x2))
ht(1t)dβ(x1,x2)
tg(x1)+(1t)g(x2)ht(1t)dβ(x1,x2).

Hence, g is a strongly convex function.

Lemma 4.2. The function g(x)=xT*Q is strongly convex with respect to the Euclidean norm in ℝn, where Q is a positive definite matrix.

g(tx1+(1t)x2)=(tx1+(1t)x2)TQ(tx1+(1t)x2)
=t2x1TQx1+t(1t)x1TQx2+t(1t)x2TQx1
+(1t)2x2TQx2
=t2x1TQx1+(1t)2x2TQx2+2t(1t)x1TQx2
=tx1TQx1+(1t)x2TQx2
t(1t)(x1TQx1+x2TQx22x1TQx2)
=tx1TQx1+(1t)x2TQx2t(1t)||x1x2||2
=tg(x1)+(1t)g(x2)t(1t)||x1x2||2.

We now show that g(x)=x*TQ*x is strongly convex with respect to the Euclidean norm ||·||E. As any two norms on a finite dimensional normed space are equivalent, there are a, b > 0 such that a||·||E||·||b||·||E. Therefore, for all x1,x2∈ℝn and t∈[0,1], we get:

g(tx1+(1t)x2)tg(x1)+(1t)g(x2)a2t(1t)||x1x2||E2.

The proof is complete.

If we set φ(x,y,μ):=f(y,μ)f(x,μ), then (MP) becomes the parametric equilibrium problem. For μ≥0, we denote the solution set of (MP) by S1(μ).

Lemma 4.3. If the function y ↦ g(y) is strongly convex, then the function yφ(x,y):=f(y)f(x) is also strongly convex.

Proof. We omit the proof as it is trivial.

We have the result for the Hölder calmness of the solution map S1(μ) as follows:

Corollary 4.1. If Q is a positive definite matrix then the solution map S1 is Hölder calm on (0, +∞).

Proof. We prove this corollary by checking all assumptions of Theorem 3.1. We see that assumptions H1, H3 and H4 are obviously fulfilled. The strong convexity and Hölder calmness in H2 is derived from Lemmas 4.1, 4.2 and 4.3.

4.2 Nash equilibrium problem

In this subsection, we discuss an application of the main results to a Nash equilibrium problem (; Facchinei and Kanzow, ). Let us state the model of the problem as follows. Let I⊂ℕ be an index set, and Λ,M,Xi (iI) be linear metric vector spaces. We use the notations X=iIXi and Xi=jI,jiXj. For each xX, we denote its ith coordinate and its projection on X-i by xi and x-i, respectively. We also denote an element of X by (xi,x−i). For iI, let fi be extended real-valued functions defined on X × M and Ki be set-valued mappings from Λ into Xi. For any fixed (λ,μ)∈Λ×M, an abstract economy, introduced by , is the set of data:

Γ(λ,μ)={Xi,Ki(λ),fi(·,μ):iI}.

A profile of strategies x*∈X is said to be a social Nash equilibrium of Γ(λ,μ) if for any iI,x*∈K(λ) and:

(14) fi(x*,μ)=maxyiK(λ)fi((yi,xi*),μ).

For each (λ,μ)∈Λ×M, the set of all social Nash equilibria of the parametric abstract economy Γ(λ,μ) is denoted by N(λ,μ). The set-valued mapping NΛ×MX is called a social Nash equilibria set-valued mapping. The following lemma gives us the equivalence between the social Nash equilibrium problem and (EP).

Lemma 4.4. For given (λ¯,μ¯)Λ×M,x* is a Nash equilibrium of abstract economy Γ(λ¯,μ¯) if and only if x*∈K(λ¯) and:

(15) φ(x*,y,μ¯)0,yK(λ¯),
i.e. x* is a solution to (EP), where:
K(λ¯):=iIKi(λ¯)andφ(x,y,μ¯)=iI[fi((xi,xi),μ¯)fi((yi,xi),μ¯)].

Proof. If holds for λ¯∈Λ, then it is obvious that x* is a solution to (EP). Conversely, let x*∈X be such that x*∈K(λ¯) and holds. For each iI, we choose y^∈X in such a way satisfying jI{i}. Then, y^∈K(λ¯) and:

φ(x*,y^,μ¯)=iI[fi((xi*,xi*),μ¯)fi((y^i,xi*),μ¯)]
=jIji[fj((xj*,xj*),μ¯)fj((y^j,xj*),λ¯)]+fi((xi*,xi*),μ¯)fi((y^i,xi*),μ¯)
=fi((xi*,xi*),μ¯)fi((y^i,xi*),μ¯).

Combining this with , we derive fi((xi*,xi*),μ¯)fi((yi,xi*),μ¯) for all yiKi(λ¯ and iI. Consequently, we conclude that for each iI, x*∈K(λ¯) and fi(x*,μ¯)=maxyiKi(λ¯)fi((yi,xi*),μ¯), i.e. x* is a social Nash equilibrium of abstract economy Γ(λ¯,μ¯).

Lemma 4.4 allows us to apply Theorem 3.1 to the mapping N, obtaining the following result.

Theorem 4.1. Assume that the set-valued map K:Λ⇉X,  and the function φ:X×X×Λ{} defined by:

K(λ):=iIKi(λ),
and:
φ(x,y,λ):=iI[fi((xi,xi),μ)fi((yi,xi),μ)],
respectively, satisfy conditions (i)-(iv) of Theorem 3.1 on Λ × M. Then, the set-valued mapping N is Hölder calm on Λ × M.

5. Conclusions

In this paper, we are successful in establishing the Hölder calmness of solution maps to equilibrium problems where both objective functions and constraints are perturbed. Our results are new. These main results are applied to some special cases including mean-variance portfolio and Nash equilibrium problem. Tools used in this paper can be used for investigating more generalized settings.

References

Alleche, B. and Rădulescu, V. (2016), “Solutions and approximate solutions of quasi-equilibrium problems in Banach spaces”, Journal of Optimization Theory and Applications, Vol. 170 No. 2, pp. 629-649.

Anh, L.Q. and Khanh, P.Q. (2008), “Various kinds of semicontinuity and solution sets of parametric multivalued symmetric vector quasi equilibrium problems”, Journal of Global Optimization, Vol. 41 No. 4, pp. 539-558.

Anh, L.Q. and Khanh, P.Q. (2010), “Continuity of solution maps of parametric quasi equilibrium problems”, Journal of Global Optimization, Vol. 46 No. 2, pp. 247-259.

Anh, L.Q., Duoc, P.T. and Tam, T.N. (2018), “On Hölder continuity of solution maps to parametric vector primal and dual equilibrium problems”, Optimization, Vol. 67 No. 8, pp. 1169-1182.

Ansari, Q.H., Konnov, I.V. and Yao, J.C. (2001), “Existence of a solution and variational principles for vector equilibrium problems”, Journal of Optimization Theory and Applications, Vol. 110 No. 3, pp. 481-492.

Bianchi, M. and Rita, P. (2006), “Sensitivity for parametric vector equilibria”, Optimization, Vol. 55 No. 3, pp. 221-230.

Blum, E. and Oettli, W. (1994), “From optimization and variational inequalities to equilibrium problems”, Math Student, Vol. 63, pp. 123-145.

Castellani, M., Pappalardo, M. and Passacantando, M. (2010), “Existence results for nonconvex equilibrium problems”, Optim Meth Software, Vol. 25 No. 1, pp. 49-58.

Debreu, G. (1952), “A social equilibrium existence theorem”, Proceedings of the National Academy of Sciences, Vol. 38 No. 10, pp. 886-893.

Facchinei, F. and Kanzow, C. (2007), “Generalized Nash equilibrium problems”, 4OR, Vol. 5 No. 3, pp. 173-210.

Farrar, D.E. (1965), The Investment Decision under Uncertainty, Prentice-Hall, New York, NY.

Hai, N.X., Khanh, P.Q. and Quan, N.H. (2009), “On the existence of solutions to quasi variational inclusion problems”, Journal of Global Optimization, Vol. 45 No. 4, pp. 565-581.

Hai, N.X. and Khanh, P.Q. (2007), “Existence of solutions to general quasi equilibrium problems and applications”, Journal of Optimization Theory and Applications, Vol. 133 No. 3, pp. 317-327.

Jafari, S., Farajzadeh, A.P., Moradi, S. and Khanh, P.Q. (2017), “Existence results for φ-quasi monotone equilibrium problems in convex metric spaces”, Optimization, Vol. 66 No. 3, pp. 293-310.

Kassay, G. and Radulescu, V.D. (2018), Equilibrium Problems and Applications, Academic Press, New York, NY.

Kimura, K. and Yao, J.C. (2008), “Sensitivity analysis of solution mappings of parametric vector quasi-equilibrium problems”, Journal of Global Optimization, Vol. 41 No. 2, pp. 187-202.

Kimura, K. and Yao, J.C. (2008), “Semicontinuity of solution mappings of parametric generalized vector equilibrium problems”, Journal of Optimization Theory and Applications, Vol. 138 No. 3, pp. 429-443.

Li, S.J., Liu, H.M., Zhang, Y. and Fang, Z.M. (2013), “Continuity of solution mappings to parametric generalized strong vector equilibrium problems”, Journal of Global Optimization, Vol. 55 No. 3, pp. 597-610.

Markowitz, H.M. (1956), “The optimization of a quadratic function subject to linear constraints”, Naval Research Logistics Quarterly, Vol. 3 Nos 1/2, pp. 111-133.

Sadeqi, I. and Alizadeh, C.G. (2011), “Existence of solutions of generalized vector equilibrium problems in reflexive Banach spaces”, Nonlinear Anal, Vol. 74 No. 6, pp. 2226-2234.

Acknowledgements

The authors are very grateful to the referee for his valuable remarks and suggestions, which helped us improve the paper.

Corresponding author

Tran Ngoc Tam can be contacted at: tntam@ctu.edu.vn

Related articles