About this Journal Submit a Manuscript Table of Contents
Journal of Applied Mathematics
Volume 2012 (2012), Article ID 487394, 24 pages
doi:10.1155/2012/487394
Research Article

Algorithms for General System of Generalized Resolvent Equations with Corresponding System of Variational Inclusions

1Scientific Computing Key Laboratory of Shanghai Universities and Department of Mathematics, Shanghai Normal University, Shanghai 200234, China
2Center for General Education, Kaohsiung Medical University, Kaohsiung 807, Taiwan

Received 8 January 2012; Accepted 14 January 2012

Academic Editor: Yonghong Yao

Copyright © 2012 Lu-Chuan Ceng and Ching-Feng Wen. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Very recently, Ahmad and Yao (2009) introduced and considered a system of generalized resolvent equations with corresponding system of variational inclusions in uniformly smooth Banach spaces. In this paper we introduce and study a general system of generalized resolvent equations with corresponding general system of variational inclusions in uniformly smooth Banach spaces. We establish an equivalence relation between general system of generalized resolvent equations and general system of variational inclusions. The iterative algorithms for finding the approximate solutions of general system of generalized resolvent equations are proposed. The convergence criteria of approximate solutions of general system of generalized resolvent equations obtained by the proposed iterative algorithm are also presented. Our results represent the generalization, improvement, supplement, and development of Ahmad and Yao corresponding ones.

1. Introduction and Preliminaries

It is well known that the theory of variational inequalities has played an important role in the investigation of a wide class of problems arising in mechanics, physics, optimization and control, nonlinear programming, elasticity, and applied sciences and so on; see, for example, [17] and the references therein. In recent years variational inequalities have been extended and generalized in different directions. A useful and significant generalization of variational inequalities is called mixed variational inequalities involving the nonlinear term [8], which enables us to study free, moving, obstacle, equilibrium problems arising in pure and applied sciences in a unified and general framework. Due to the presence of the nonlinear term, the projection method and its variant forms including the technique of the Wiener-Hopf equations cannot be extended to suggest the iterative methods for solving mixed variational inequalities. To overcome these drawbacks, Hassouni and Moudafi [9] introduced variational inclusions which contain mixed variational inequalities as special cases. They studied the perturbed method for solving variational inclusions. Subsequently, M. A. Noor and K. I. Noor [10] introduced and considered the resolvent equations by virtue of the resolvent operator concept and established the equivalence between the mixed variational inequalities and the resolvent equations. The technique of resolvent equations is being used to develop powerful and efficient numerical techniques for solving mixed (quasi)variational inequalities and related optimization problems. At the same time, some iterative algorithms for approximating a solution of some system of variational inequalities are also introduced and studied in Verma [11]. Pang [12], Cohen and Chaplais [13], Binachi [14], Ansari and Yao [15] considered a system of scalar variational inequalities and Pang showed that the traffic equilibrium problem, the Nash equilibrium, and the general equilibrium programming problem can be modeled as a system of variational inequalities. As generalizations of system of variational inequalities, Agarwal et al. [16] introduced a system of generalized nonlinear mixed quasi-variational inclusions and investigated the sensitivity analysis of solutions for the system of generalized mixed quasi-variational inclusions in Hilbert spaces. In 2007, Peng and Zhu [17] considered and studied a new system of generalized mixed quasi-variational inclusions with ( 𝐻 , 𝜂 ) -monotone operators and Lan et al. [18] studied a new system of nonlinear 𝐴 -monotone multivalued variational inclusions. Furthermore, for more details in the related research work of this field, we invoke the readers to see, for instance, [1930]. Very recently, Ahmad and Yao [31] introduced and considered a new system of variational inclusions in uniformly smooth Banach spaces, which covers the system of variational inclusions in Hilbert spaces considered by [18]. They established an equivalence relation between this system of variational inclusions and a system of generalized resolvent equations, proposed a number of iterative algorithms for this system of variational inclusions, and also gave the convergence criteria.

Let 𝐸 be a real Banach space with its norm , 𝐸 the topological dual of 𝐸 , and 𝑑 the metric induced by the norm . Let 𝐶 𝐵 ( 𝐸 ) (resp., 2 𝐸 ) be the family of all nonempty closed and bounded subsets (resp., all nonempty subsets) of 𝐸 and 𝐷 ( , ) the Hausdorff metric on 𝐶 𝐵 ( 𝐸 ) defined by 𝐷 ( 𝐴 , 𝐵 ) = m a x s u p 𝑥 𝐴 𝑑 ( 𝑥 , 𝐵 ) , s u p 𝑦 𝐵 𝑑 ( 𝐴 , 𝑦 ) , ( 1 . 1 ) where 𝑑 ( 𝑥 , 𝐵 ) = i n f 𝑦 𝐵 𝑑 ( 𝑥 , 𝑦 ) and 𝑑 ( 𝐴 , 𝑦 ) = i n f 𝑥 𝐴 𝑑 ( 𝑥 , 𝑦 ) . We write by 𝐽 𝐸 2 𝐸 the normalized duality mapping defined as 𝐽 ( 𝑥 ) = 𝑓 𝐸 𝑥 , 𝑓 = 𝑥 2 = 𝑓 2 , 𝑥 𝐸 , ( 1 . 2 ) where , denotes the duality pairing between 𝐸 and 𝐸 .

The uniform convexity of a Banach space 𝐸 means that for any 𝜖 > 0 , there exists 𝛿 > 0 , such that for any 𝑥 , 𝑦 𝐸 , 𝑥 1 , 𝑦 1 , 𝑥 𝑦 = 𝜖 ensure the following inequality: 𝑥 + 𝑦 2 ( 1 𝛿 ) . ( 1 . 3 ) The function 𝛿 𝐸 ( 𝜖 ) = i n f 1 𝑥 + 𝑦 2 𝑥 = 1 , 𝑦 = 1 , 𝑥 𝑦 = 𝜖 ( 1 . 4 ) is called the modulus of convexity of 𝐸 .

The uniform smoothness of a Banach space 𝐸 means that for any given 𝜖 > 0 , there exists 𝛿 > 0 such that 𝑥 + 𝑦 + 𝑥 𝑦 2 1 𝜖 𝑦 ( 1 . 5 ) holds. The function 𝜏 𝐸 ( 𝑡 ) = s u p 𝑥 + 𝑦 + 𝑥 𝑦 2 1 𝑥 = 1 , 𝑦 = 𝑡 ( 1 . 6 ) is called the modulus of smoothness of 𝐸 .

It is well known that the Banach space 𝐸 is uniformly convex if and only if 𝛿 𝐸 ( 𝜖 ) > 0 for all 𝜖 > 0 , and it is uniformly smooth if and only if l i m 𝑡 0 𝜏 𝐸 ( 𝑡 ) / 𝑡 = 0 . All Hilbert spaces, 𝐿 𝑝 (or 𝑙 𝑝 ) spaces ( 𝑝 2 ), and the Sobolov spaces 𝑊 𝑝 𝑚 ( 𝑝 2 ) are 2-uniformly smooth, while, for 1 < 𝑝 2 , 𝐿 𝑝 (or 𝑙 𝑝 ) and 𝑊 𝑝 𝑚 spaces are 𝑝 -uniformly smooth.

Proposition 1.1 (see [15]). Let 𝐸 be a uniformly smooth Banach space. Then the normalized duality mapping   𝐽 𝐸 2 𝐸 is single-valued, and for any 𝑥 , 𝑦 𝐸 there holds the following:(i) 𝑥 + 𝑦 2 𝑥 2 + 2 𝑦 , 𝐽 ( 𝑥 + 𝑦 ) ,(ii) 𝑥 𝑦 , 𝐽 ( 𝑥 ) 𝐽 ( 𝑦 ) 2 𝐶 2 𝜏 𝐸 ( 4 𝑥 𝑦 / 𝐶 ) , where 𝐶 = 𝑥 2 + 𝑦 2 / 2 .

Definition 1.2 (see [32]). A mapping 𝑔 𝐸 𝐸 is said to be
(i) 𝑘 -strongly accretive, 𝑘 ( 0 , 1 ) , if for any 𝑥 , 𝑦 𝐸 , there exists 𝑗 ( 𝑥 𝑦 ) 𝐽 ( 𝑥 𝑦 ) such that 𝑔 ( 𝑥 ) 𝑔 ( 𝑦 ) , 𝑗 ( 𝑥 𝑦 ) 𝑘 𝑥 𝑦 2 ; ( 1 . 7 ) (ii)Lipschitz continuous if for any 𝑥 , 𝑦 𝐸 , there exists a constant 𝜆 𝑔 > 0 , such that 𝑔 ( 𝑥 ) 𝑔 ( 𝑦 ) 𝜆 𝑔 𝑥 𝑦 . ( 1 . 8 )

Definition 1.3 (see [13]). A set-valued mapping 𝐴 𝐸 2 𝐸 is said to be
(i)accretive, if for any 𝑥 , 𝑦 𝐸 , there exists 𝑗 ( 𝑥 𝑦 ) 𝐽 ( 𝑥 𝑦 ) such that for all 𝑢 𝐴 ( 𝑥 ) and 𝑣 𝐴 ( 𝑦 ) , 𝑢 𝑣 , 𝑗 ( 𝑥 𝑦 ) 0 ; ( 1 . 9 ) (ii) 𝑘 -strongly accretive, 𝑘 ( 0 , 1 ) , if for any 𝑥 , 𝑦 𝐸 , there exists 𝑗 ( 𝑥 𝑦 ) 𝐽 ( 𝑥 𝑦 ) , such that for all 𝑢 𝐴 ( 𝑥 ) and 𝑣 𝐴 ( 𝑦 ) , 𝑢 𝑣 , 𝑗 ( 𝑥 𝑦 ) 𝑘 𝑥 𝑦 2 ; ( 1 . 1 0 ) (iii) 𝑚 -accretive if 𝐴 is accretive and ( 𝐼 + 𝜌 𝐴 ) ( 𝐸 ) = 𝐸 , for every (equivalently, for some) 𝜌 > 0 , where 𝐼 is the identity mapping (equivalently, if 𝐴 is accretive and ( 𝐼 + 𝐴 ) ( 𝐸 ) = 𝐸 ). In particular, it is clear from [9] that if 𝐸 = 𝐻 is a Hilbert space, than 𝐴 𝐸 2 𝐸 is an 𝑚 -accretive mapping if and only if it is a maximal monotone mapping.

Definition 1.4 (see [31]). Let 𝑀 𝐸 2 𝐸 be an 𝑚 -accretive mapping. For any 𝜌 > 0 , the mapping 𝐽 𝜌 𝑀 𝐸 𝐸 associated with 𝑀 defined by 𝐽 𝜌 𝑀 ( 𝑥 ) = ( 𝐼 + 𝜌 𝑀 ) 1 ( 𝑥 ) , 𝑥 𝐸 ( 1 . 1 1 ) is called the resolvent operator.

Definition 1.5 (see [33]). The resolvent operator 𝐽 𝜌 𝑀 𝐸 𝐸 is said to be a retraction if ( 𝐼 + 𝜌 𝑀 ) 1 ( 𝐼 + 𝜌 𝑀 ) 1 ( 𝑥 ) = ( 𝐼 + 𝜌 𝑀 ) 1 ( 𝑥 ) , 𝑥 𝐸 . ( 1 . 1 2 )
It is well known that 𝐽 𝜌 𝑀 is a single-valued and nonexpansive mapping.

Definition 1.6 (see [10]). A set-valued mapping 𝐻 𝐸 𝐶 𝐵 ( 𝐸 ) is said to be 𝐷 -Lipschitz continuous if for any 𝑥 , 𝑦 𝐸 , there exists a constant 𝜆 𝐷 𝐻 > 0 such that 𝐷 ( 𝐻 ( 𝑥 ) , 𝐻 ( 𝑦 ) ) 𝜆 𝐷 𝐻 𝑥 𝑦 . ( 1 . 1 3 )

Let 𝐸 1 and 𝐸 2 be two real Banach spaces, 𝑆 𝐸 1 × 𝐸 2 𝐸 1 and 𝑇 𝐸 1 × 𝐸 2 𝐸 2 single-valued mappings, and 𝐺 𝐸 1 2 𝐸 1 , 𝐹 𝐸 2 2 𝐸 2 , 𝐻 𝐸 1 2 𝐸 1 and 𝑉 𝐸 2 2 𝐸 2 any four multivalued mappings. Let 𝑀 𝐸 1 2 𝐸 1 and 𝑁 𝐸 2 2 𝐸 2 be any nonlinear mappings, 𝑚 𝐸 2 𝐸 1 , 𝑛 𝐸 1 𝐸 2 , 𝑓 𝐸 1 𝐸 1 and 𝑔 𝐸 2 𝐸 2 nonlinear mappings with 𝑓 ( 𝐸 1 ) 𝐷 ( 𝑀 ) and 𝑔 ( 𝐸 2 ) 𝐷 ( 𝑁 ) , respectively. Then we consider the problem of finding ( 𝑥 , 𝑦 ) 𝐸 1 × 𝐸 2 , ( 𝑠 , 𝑣 ) 𝐺 ( 𝑥 ) × 𝐹 ( 𝑦 ) , ( 𝑢 , 𝑡 ) 𝐻 ( 𝑥 ) × 𝑉 ( 𝑦 ) such that 𝑛 𝑚 ( 𝑦 ) 𝑆 ( 𝑠 , 𝑣 ) + 𝑀 ( 𝑓 ( 𝑥 ) ) , ( 𝑥 ) 𝑇 ( 𝑢 , 𝑡 ) + 𝑁 ( 𝑔 ( 𝑦 ) ) , ( 1 . 1 4 ) which is called a general system of variational inclusions. In particular, if 𝑚 ( 𝑦 ) = 0 𝐸 1 , 𝑛 ( 𝑥 ) = 0 𝐸 2 , 𝐺 ( 𝑥 ) = 𝑝 ( 𝑥 ) and 𝑉 ( 𝑦 ) = 𝑞 ( 𝑦 ) , where 𝑝 𝐸 1 𝐸 1 and 𝑞 𝐸 2 𝐸 2 are single-valued mappings, then the general system of variational inclusions (1.14) reduces to the following system of variational inclusions 0 𝑆 ( 𝑝 ( 𝑥 ) , 𝑣 ) + 𝑀 ( 𝑓 ( 𝑥 ) ) , 0 𝑇 ( 𝑢 , 𝑞 ( 𝑦 ) ) + 𝑁 ( 𝑔 ( 𝑦 ) ) , ( 1 . 1 5 ) which was considered by Lan et al. [18] in Hilbert spaces and studied by Ahmad and Yao [31] in Banach spaces, respectively.

Proposition 1.7 (see [31, Lemma  2.1]). ( 𝑥 , 𝑦 ) 𝐸 1 × 𝐸 2 , 𝑢 𝐻 ( 𝑥 ) , 𝑣 𝐹 ( 𝑦 ) is a solution of the system of variational inclusions (1.15) if and only if ( 𝑥 , 𝑦 , 𝑢 , 𝑣 ) satisfies 𝑓 ( 𝑥 ) = 𝐽 𝜌 𝑀 ( 𝑓 ( 𝑥 ) 𝜌 𝑆 ( 𝑝 ( 𝑥 ) , 𝑣 ) ) , 𝜌 > 0 , 𝑔 ( 𝑦 ) = 𝐽 𝛾 𝑁 ( 𝑔 ( 𝑦 ) 𝛾 𝑇 ( 𝑢 , 𝑞 ( 𝑦 ) ) ) , 𝛾 > 0 . ( 1 . 1 6 )

Proposition 1.8 (see [31, Proposition  3.1]). The system of variational inclusions (1.15) has a solution ( 𝑥 , 𝑦 , 𝑢 , 𝑣 ) with ( 𝑥 , 𝑦 ) 𝐸 1 × 𝐸 2 , 𝑢 𝐻 ( 𝑥 ) and 𝑣 𝐹 ( 𝑦 ) if and only if the following system of generalized resolvent equations 𝑆 ( 𝑝 ( 𝑥 ) , 𝑣 ) + 𝜌 1 𝑅 𝜌 𝑀 𝑧 = 0 , 𝑅 𝜌 𝑀 = 𝐼 𝐽 𝜌 𝑀 𝑇 , 𝜌 > 0 , ( 𝑢 , 𝑞 ( 𝑦 ) ) + 𝛾 1 𝑅 𝛾 𝑁 𝑧 = 0 , 𝑅 𝛾 𝑁 = 𝐼 𝐽 𝛾 𝑁 , 𝛾 > 0 , ( 1 . 1 7 ) has a solution ( 𝑧 , 𝑧 , 𝑥 , 𝑦 , 𝑢 , 𝑣 ) with ( 𝑥 , 𝑦 ) 𝐸 1 × 𝐸 2 , 𝑢 𝐻 ( 𝑥 ) , 𝑣 𝐹 ( 𝑦 ) , 𝑧 𝐸 1 and 𝑧 𝐸 2 , where 𝑓 ( 𝑥 ) = 𝐽 𝜌 𝑀 ( 𝑧 ) , 𝑔 ( 𝑦 ) = 𝐽 𝛾 𝑁 ( 𝑧 ) and 𝑧 = 𝑓 ( 𝑥 ) 𝜌 𝑆 ( 𝑝 ( 𝑥 ) , 𝑣 ) , 𝑧 = 𝑔 ( 𝑦 ) 𝛾 𝑇 ( 𝑢 , 𝑞 ( 𝑦 ) ) .

Based on the above Propositions 1.7 and 1.8, Ahmad and Yao [31] presented the following algorithm and established the following strong convergence result for the sequences generated by the algorithm.

Algorithm 1.9 (see [31, Algorithm  3.1]). For given ( 𝑥 0 , 𝑦 0 ) 𝐸 1 × 𝐸 2 , 𝑢 0 𝐻 ( 𝑥 0 ) , 𝑣 0 𝐹 ( 𝑦 0 ) , 𝑧 0 𝐸 1 and 𝑧 0 𝐸 2 , compute { 𝑧 𝑘 } , { 𝑧 𝑘 } , { 𝑥 𝑘 } , { 𝑦 𝑘 } , { 𝑢 𝑘 } , and { 𝑣 𝑘 } by the iterative scheme: 𝑓 𝑥 𝑘 = 𝐽 𝜌 𝑀 𝑧 𝑘 , 𝑔 𝑦 𝑘 = 𝐽 𝛾 𝑁 𝑧 𝑘 , 𝑢 𝑘 𝑥 𝐻 𝑘 𝑢 𝑘 + 1 𝑢 𝑘 𝐻 𝑥 𝐷 𝑘 + 1 𝑥 , 𝐻 𝑘 , 𝑣 𝑘 𝑦 𝐹 𝑘 𝑣 𝑘 + 1 𝑣 𝑘 𝐹 𝑦 𝐷 𝑘 + 1 𝑦 , 𝐹 𝑘 , 𝑧 𝑘 + 1 𝑥 = 𝑓 𝑘 𝑝 𝑥 𝜌 𝑆 𝑘 , 𝑣 𝑘 , 𝑧 𝑘 + 1 𝑦 = 𝑔 𝑘 𝑢 𝛾 𝑇 𝑘 𝑦 , 𝑞 𝑘 , 𝑘 = 0 , 1 , 2 , . ( 1 . 1 8 )

Theorem 1.10 (see [31, Theorem  3.1]). Let 𝐸 1 and 𝐸 2 be two real uniformly smooth Banach spaces with modulus of smoothness 𝜏 𝐸 1 ( 𝑡 ) 𝐶 1 𝑡 2 and 𝜏 𝐸 2 ( 𝑡 ) 𝐶 2 𝑡 2 for 𝐶 1 , 𝐶 2 > 0 , respectively. Let 𝐻 𝐸 1 𝐶 𝐵 ( 𝐸 1 ) , 𝐹 𝐸 2 𝐶 𝐵 ( 𝐸 2 ) be 𝐷 -Lipschitz continuous mappings with constants 𝜆 𝐷 𝐻 and 𝜆 𝐷 𝐹 , respectively, and let 𝑀 𝐸 1 2 𝐸 1 , 𝑁 𝐸 2 2 𝐸 2 be 𝑚 -accretive mappings such that the resolvent operators associated with 𝑀 and 𝑁 are retractions. Let 𝑓 𝐸 1 𝐸 1 , 𝑔 𝐸 2 𝐸 2 be both strong accretive with constants 𝛼 and 𝛽 , respectively, and Lipschitz continuous with constants 𝛿 1 and 𝛿 2 , respectively. Let 𝑝 𝐸 1 𝐸 1 , 𝑞 𝐸 2 𝐸 2 be Lipschitz continuous with constants 𝜆 𝑝 and 𝜆 𝑞 , respectively, and let 𝑆 𝐸 1 × 𝐸 2 𝐸 1 , 𝑇 𝐸 1 × 𝐸 2 𝐸 2 be Lipschitz continuous in the first and second arguments with constants 𝜆 𝑆 1 , 𝜆 𝑆 2 and 𝜆 𝑇 1 , 𝜆 𝑇 2 , respectively.
If there exist constants 𝜌 > 0 and 𝛾 > 0 , such that 0 < 𝐵 + 𝜃 1 + 𝜃 4 𝐵 1 𝐵 < 1 , 0 < + 𝜃 2 + 𝜃 3 1 𝐵 < 1 , ( 1 . 1 9 ) where 𝐵 = 1 2 𝛼 + 6 4 𝐶 1 𝛿 2 1 , 𝐵 = 1 2 𝛽 + 6 4 𝐶 2 𝛿 2 2 and 𝜃 1 = ( 1 + 𝜌 𝜆 𝑆 1 𝜆 𝑝 ) / ( 1 𝜌 ( 𝜆 𝑆 1 𝜆 𝑝 + 𝜆 𝑆 2 𝜆 𝐷 𝐹 ) ) , 𝜃 2 = ( 1 + 𝜌 𝜆 𝑆 2 𝜆 𝐷 𝐹 ) / ( 1 𝜌 ( 𝜆 𝑆 1 𝜆 𝑝 + 𝜆 𝑆 2 𝜆 𝐷 𝐹 ) ) , 𝜃 3 = ( 1 + 𝛾 𝜆 𝑇 2 𝜆 𝑞 ) / ( 1 𝛾 ( 𝜆 𝑇 1 𝜆 𝐷 𝐻 + 𝜆 𝑇 2 𝜆 𝑞 ) ) , 𝜃 4 = 𝛾 𝜆 𝑇 1 𝜆 𝐷 𝐻 / ( 1 𝛾 ( 𝜆 𝑇 1 𝜆 𝐷 𝐻 + 𝜆 𝑇 2 𝜆 𝑞 ) ) , then there exist ( 𝑥 , 𝑦 ) 𝐸 1 × 𝐸 2 , 𝑢 𝐻 ( 𝑥 ) , 𝑣 𝐹 ( 𝑦 ) and ( 𝑧 , 𝑧 ) 𝐸 1 × 𝐸 2 satisfying the system of generalized resolvent equations (1.17) (in this case, ( 𝑥 , 𝑦 , 𝑢 , 𝑣 ) is a solution of system of variational inclusions (1.15)), and the iterative sequences { 𝑧 𝑘 } , { 𝑧 𝑘 } , { 𝑥 𝑘 } , { 𝑦 𝑘 } , { 𝑢 𝑘 } , and { 𝑣 𝑘 } generated by Algorithm 1.9 converge strongly to 𝑧 , 𝑧 , 𝑥 , 𝑦 , 𝑢 , and 𝑣 , respectively.

In this paper we introduce and study a general system of generalized resolvent equations with corresponding general system of variational inclusions in uniformly smooth Banach spaces. Motivated and inspired by the above Proposition 1.8, we establish an equivalence relation between general system of generalized resolvent equations and general system of variational inclusions. By using Nadler [34] we propose some new iterative algorithms for finding the approximate solutions of general system of generalized resolvent equations, which include Ahmad and Yao’s corresponding algorithms as special cases to a great extent. Furthermore, the convergence criteria of approximate solutions of general system of generalized resolvent equations obtained by the proposed iterative algorithm are also presented. There is no doubt that our results represent the generalization, improvement, supplement, and development of Ahmad and Yao corresponding ones [31].

2. Main Results

Let 𝐸 1 and 𝐸 2 be two real Banach spaces, let 𝑆 𝐸 1 × 𝐸 2 𝐸 1 and 𝑇 𝐸 1 × 𝐸 2 𝐸 2 be single-valued mappings, and let 𝐺 𝐸 1 2 𝐸 1 , 𝐹 𝐸 2 2 𝐸 2 , 𝐻 𝐸 1 2 𝐸 1 and 𝑉 𝐸 2 2 𝐸 2 be any four multivalued mappings. Let 𝑀 𝐸 1 2 𝐸 1 and 𝑁 𝐸 2 2 𝐸 2 be any nonlinear mappings, 𝑚 𝐸 2 𝐸 1 , 𝑛 𝐸 1 𝐸 2 , 𝑓 𝐸 1 𝐸 1 and 𝑔 𝐸 2 𝐸 2 nonlinear mappings with 𝑓 ( 𝐸 1 ) 𝐷 ( 𝑀 ) and 𝑔 ( 𝐸 2 ) 𝐷 ( 𝑁 ) , respectively. Then we consider the problem of finding ( 𝑥 , 𝑦 ) 𝐸 1 × 𝐸 2 , ( 𝑠 , 𝑣 ) 𝐺 ( 𝑥 ) × 𝐹 ( 𝑦 ) , ( 𝑢 , 𝑡 ) 𝐻 ( 𝑥 ) × 𝑉 ( 𝑦 ) , 𝑧 𝐸 1 , 𝑧 𝐸 2 such that 𝑆 ( 𝑠 , 𝑣 ) + 𝜌 1 𝑅 𝜌 𝑀 𝑧 𝑇 = 𝑚 ( 𝑦 ) , 𝜌 > 0 , ( 𝑢 , 𝑡 ) + 𝛾 1 𝑅 𝛾 𝑁 𝑧 = 𝑛 ( 𝑥 ) , 𝛾 > 0 , ( 2 . 1 ) where 𝑅 𝜌 𝑀 = 𝐼 𝐽 𝜌 𝑀 , 𝑅 𝛾 𝑁 = 𝐼 𝐽 𝛾 𝑁 and 𝐽 𝜌 𝑀 , 𝐽 𝛾 𝑁 are the resolvent operators associated with 𝑀 and 𝑁 , respectively.

The corresponding general system of variational inclusions of (2.1) is the problem (1.14), that is, find ( 𝑥 , 𝑦 ) 𝐸 1 × 𝐸 2 , ( 𝑠 , 𝑣 ) 𝐺 ( 𝑥 ) × 𝐹 ( 𝑦 ) , ( 𝑢 , 𝑡 ) 𝐻 ( 𝑥 ) × 𝑉 ( 𝑦 ) such that 𝑛 𝑚 ( 𝑦 ) 𝑆 ( 𝑠 , 𝑣 ) + 𝑀 ( 𝑓 ( 𝑥 ) ) , ( 𝑥 ) 𝑇 ( 𝑢 , 𝑡 ) + 𝑁 ( 𝑔 ( 𝑦 ) ) . ( 2 . 2 )

Proposition 2.1. ( 𝑥 , 𝑦 ) 𝐸 1 × 𝐸 2 , ( 𝑠 , 𝑣 ) 𝐺 ( 𝑥 ) × 𝐹 ( 𝑦 ) , ( 𝑢 , 𝑡 ) 𝐻 ( 𝑥 ) × 𝑉 ( 𝑦 ) are solutions of general system of variational inclusions (1.14) if and only if ( 𝑥 , 𝑦 , 𝑢 , 𝑣 , 𝑠 , 𝑡 ) satisfies 𝑓 ( 𝑥 ) = 𝐽 𝜌 𝑀 [ ] 𝑓 ( 𝑥 ) 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) , 𝜌 > 0 , 𝑔 ( 𝑦 ) = 𝐽 𝛾 𝑁 [ ] 𝑔 ( 𝑦 ) 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) , 𝛾 > 0 . ( 2 . 3 )

Proof. The proof of Proposition 2.1 is a direct consequence of the definition of resolvent operator, and hence, is omitted.

Next we first establish an equivalence relation between general system of generalized resolvent equations (2.1) and general system of variational inclusions (1.14) and then prove the existence of a solution of (2.1) and convergence of sequences generated by the proposed algorithms.

Proposition 2.2. The general system of variational inclusions (1.14) has a solution ( 𝑥 , 𝑦 , 𝑢 , 𝑣 , 𝑠 , 𝑡 ) with ( 𝑥 , 𝑦 ) 𝐸 1 × 𝐸 2 , ( 𝑠 , 𝑣 ) 𝐺 ( 𝑥 ) × 𝐹 ( 𝑦 ) and ( 𝑢 , 𝑡 ) 𝐻 ( 𝑥 ) × 𝑉 ( 𝑦 ) if and only if general system of generalized resolvent equations (2.1) has a solution ( 𝑧 , 𝑧 , 𝑥 , 𝑦 , 𝑢 , 𝑣 , 𝑠 , 𝑡 ) with ( 𝑥 , 𝑦 ) 𝐸 1 × 𝐸 2 , ( 𝑠 , 𝑣 ) 𝐺 ( 𝑥 ) × 𝐹 ( 𝑦 ) , ( 𝑢 , 𝑡 ) 𝐻 ( 𝑥 ) × 𝑉 ( 𝑦 ) , ( 𝑧 , 𝑧 ) 𝐸 1 × 𝐸 2 , where 𝑓 ( 𝑥 ) = 𝐽 𝜌 𝑀 𝑧 , 𝑔 ( 𝑦 ) = 𝐽 𝛾 𝑁 𝑧 , ( 2 . 4 ) and 𝑧 = 𝑓 ( 𝑥 ) 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) and 𝑧 = 𝑔 ( 𝑦 ) 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) .

Proof. Let ( 𝑥 , 𝑦 ) 𝐸 1 × 𝐸 2 , ( 𝑠 , 𝑣 ) 𝐺 ( 𝑥 ) × 𝐹 ( 𝑦 ) , ( 𝑢 , 𝑡 ) 𝐻 ( 𝑥 ) × 𝑉 ( 𝑦 ) be a solution of general system of variational inclusions (1.14). Then, by Proposition 2.1, it satisfies the following system of equations 𝑓 ( 𝑥 ) = 𝐽 𝜌 𝑀 [ ] , 𝑓 ( 𝑥 ) 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) 𝑔 ( 𝑦 ) = 𝐽 𝛾 𝑁 [ ] . 𝑔 ( 𝑦 ) 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) ( 2 . 5 )
Let 𝑧 = 𝑓 ( 𝑥 ) 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) and 𝑧 = 𝑔 ( 𝑦 ) 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) . Then we have 𝑓 ( 𝑥 ) = 𝐽 𝜌 𝑀 𝑧 , 𝑔 ( 𝑦 ) = 𝐽 𝛾 𝑁 𝑧 , ( 2 . 6 ) and hence 𝑧 = 𝐽 𝜌 𝑀 ( 𝑧 ) 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) and 𝑧 = 𝐽 𝛾 𝑁 ( 𝑧 ) 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) . Thus it follows that 𝐼 𝐽 𝜌 𝑀 𝑧 = 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) , 𝐼 𝐽 𝛾 𝑁 𝑧 = 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) , ( 2 . 7 ) that is, 𝑆 ( 𝑠 , 𝑣 ) + 𝜌 1 𝑅 𝜌 𝑀 𝑧 𝑇 = 𝑚 ( 𝑦 ) , ( 𝑢 , 𝑡 ) + 𝛾 1 𝑅 𝛾 𝑁 𝑧 = 𝑛 ( 𝑥 ) . ( 2 . 8 ) Therefore, ( 𝑧 , 𝑧 , 𝑥 , 𝑦 , 𝑢 , 𝑣 , 𝑠 , 𝑡 ) is a solution of general system of generalized resolvent equations (2.1).
Conversely, let ( 𝑧 , 𝑧 , 𝑥 , 𝑦 , 𝑢 , 𝑣 , 𝑠 , 𝑡 ) be a solution of general system of generalized resolvent equations (2.1). Then 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) = 𝑅 𝜌 𝑀 𝑧 , 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) = 𝑅 𝛾 𝑁 𝑧 . ( 2 . 9 )
Now observe that 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) = 𝑅 𝜌 𝑀 𝑧 = 𝐼 𝐽 𝜌 𝑀 𝑧 = 𝐽 𝜌 𝑀 𝑧 𝑧 = 𝐽 𝜌 𝑀 [ ] [ ] , 𝑓 ( 𝑥 ) 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) 𝑓 ( 𝑥 ) 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) ( 2 . 1 0 ) which leads to 𝑓 ( 𝑥 ) = 𝐽 𝜌 𝑀 [ ] , 𝑓 ( 𝑥 ) 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) ( 2 . 1 1 ) and also that 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) = 𝑅 𝛾 𝑁 𝑧 = 𝐼 𝐽 𝛾 𝑁 𝑧 = 𝐽 𝛾 𝑁 𝑧 𝑧 = 𝐽 𝛾 𝑁 [ ] [ ] , 𝑔 ( 𝑦 ) 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) 𝑔 ( 𝑦 ) 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) ( 2 . 1 2 ) which leads to 𝑔 ( 𝑦 ) = 𝐽 𝛾 𝑁 [ ] . 𝑔 ( 𝑦 ) 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) ( 2 . 1 3 ) Consequently, we have 𝑓 ( 𝑥 ) = 𝐽 𝜌 𝑀 [ ] , 𝑓 ( 𝑥 ) 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) 𝑔 ( 𝑦 ) = 𝐽 𝛾 𝑁 [ ] . 𝑔 ( 𝑦 ) 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) ( 2 . 1 4 ) Therefore, by Proposition 2.1, ( 𝑥 , 𝑦 , 𝑢 , 𝑣 , 𝑠 , 𝑡 ) is a solution of general system of variational inclusions (1.14).

Proof (Alternative). Let 𝑧 = 𝑓 ( 𝑥 ) 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) , 𝑧 = 𝑔 ( 𝑦 ) 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) . ( 2 . 1 5 ) Then, utilizing (2.4), we can write 𝑧 = 𝐽 𝜌 𝑀 𝑧 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) , 𝑧 = 𝐽 𝛾 𝑁 𝑧 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) ( 2 . 1 6 ) which yield that 𝑆 ( 𝑠 , 𝑣 ) + 𝜌 1 𝑅 𝜌 𝑀 𝑧 𝑇 = 𝑚 ( 𝑦 ) , ( 𝑢 , 𝑡 ) + 𝛾 1 𝑅 𝛾 𝑁 𝑧 = 𝑛 ( 𝑥 ) , ( 2 . 1 7 ) the required general system of generalized resolvent equations.

Algorithm 2.3. For given ( 𝑥 0 , 𝑦 0 ) 𝐸 1 × 𝐸 2 , ( 𝑠 0 , 𝑣 0 ) 𝐺 ( 𝑥 0 ) × 𝐹 ( 𝑦 0 ) , ( 𝑢 0 , 𝑡 0 ) 𝐻 ( 𝑥 0 ) × 𝑉 ( 𝑦 0 ) , ( 𝑧 0 , 𝑧 0 ) 𝐸 1 × 𝐸 2 , compute 𝑧 1 𝑥 = 𝑓 0 𝑆 𝑠 𝜌 0 , 𝑣 0 𝑦 𝑚 0 , 𝑧 1 𝑦 = 𝑔 0 𝑇 𝑢 𝛾 0 , 𝑡 0 𝑥 𝑛 0 . ( 2 . 1 8 ) For ( 𝑧 1 , 𝑧 1 ) 𝐸 1 × 𝐸 2 , we take ( 𝑥 1 , 𝑦 1 ) 𝐸 1 × 𝐸 2 such that 𝑓 ( 𝑥 1 ) = 𝐽 𝜌 𝑀 ( 𝑧 1 ) and 𝑔 ( 𝑦 1 ) = 𝐽 𝛾 𝑁 ( 𝑧 1 ) . Then, by Nadler [34], there exist ( 𝑠 1 , 𝑣 1 ) 𝐺 ( 𝑥 1 ) × 𝐹 ( 𝑦 1 ) , ( 𝑢 1 , 𝑡 1 ) 𝐻 ( 𝑥 1 ) × 𝑉 ( 𝑦 1 ) such that 𝑢 1 𝑢 0 𝐻 𝑥 ( 1 + 1 ) 𝐷 1 𝑥 , 𝐻 0 , 𝑣 1 𝑣 0 𝐹 𝑦 ( 1 + 1 ) 𝐷 1 𝑦 , 𝐹 0 , 𝑠 1 𝑠 0 𝐺 𝑥 ( 1 + 1 ) 𝐷 1 𝑥 , 𝐺 0 , 𝑡 1 𝑡 0 𝑉 𝑦 ( 1 + 1 ) 𝐷 1 𝑦 , 𝑉 0 , ( 2 . 1 9 ) where 𝐷 ( , ) is the Hausdorff metric on 𝐶 𝐵 ( 𝐸 1 ) (for the sake of convenience, we also denote by 𝐷 ( , ) the Hausdorff metric on 𝐶 𝐵 ( 𝐸 2 ) ). Compute 𝑧 2 𝑥 = 𝑓 1 𝑆 𝑠 𝜌 1 , 𝑣 1 𝑦 𝑚 1 , 𝑧 2 𝑦 = 𝑔 1 𝑇 𝑢 𝛾 1 , 𝑡 1 𝑥 𝑛 1 . ( 2 . 2 0 ) By induction, we can obtain sequences ( 𝑥 𝑘 , 𝑦 𝑘 ) 𝐸 1 × 𝐸 2 , ( 𝑠 𝑘 , 𝑣 𝑘 ) 𝐺 ( 𝑥 𝑘 ) × 𝐹 ( 𝑦 𝑘 ) , ( 𝑢 𝑘 , 𝑡 𝑘 ) 𝐻 ( 𝑥 𝑘 ) × 𝑉 ( 𝑦 𝑘 ) , ( 𝑧 𝑘 , 𝑧 𝑘 ) 𝐸 1 × 𝐸 2 by the iterative scheme: 𝑓 𝑥 𝑘 = 𝐽 𝜌 𝑀 𝑧 𝑘 , 𝑔 𝑦 𝑘 = 𝐽 𝛾 𝑁 𝑧 𝑘 , 𝑢 ( 2 . 2 1 ) 𝑘 𝑥 𝐻 𝑘 𝑢 𝑘 + 1 𝑢 𝑘 1 1 + 𝐷 𝐻 𝑥 𝑘 + 1 𝑘 + 1 𝑥 , 𝐻 𝑘 , 𝑣 𝑘 𝑦 𝐹 𝑘 𝑣 𝑘 + 1 𝑣 𝑘 1 1 + 𝐷 𝐹 𝑦 𝑘 + 1 𝑘 + 1 𝑦 , 𝐹 𝑘 , 𝑠 𝑘 𝑥 𝐺 𝑘 𝑠 𝑘 + 1 𝑠 𝑘 1 1 + 𝐷 𝐺 𝑥 𝑘 + 1 𝑘 + 1 𝑥 , 𝐺 𝑘 , 𝑡 𝑘 𝑦 𝑉 𝑘 𝑡 𝑘 + 1 𝑡 𝑘 1 1 + 𝐷 𝑉 𝑦 𝑘 + 1 𝑘 + 1 𝑦 , 𝑉 𝑘 , 𝑧 ( 2 . 2 2 ) 𝑘 + 1 𝑥 = 𝑓 𝑘 𝑆 𝑠 𝜌 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 , 𝑧 𝑘 + 1 𝑦 = 𝑔 𝑘 𝑇 𝑢 𝛾 𝑘 , 𝑡 𝑘 𝑥 𝑛 𝑘 , ( 2 . 2 3 ) for 𝑘 = 0 , 1 , 2 , .
The general system of generalized resolvent equations (2.1) can also be rewritten as 𝑧 = 𝑓 ( 𝑥 ) 𝑆 ( 𝑠 , 𝑣 ) + 𝑚 ( 𝑦 ) + 𝐼 𝜌 1 𝑅 𝜌 𝑀 𝑧 , 𝑧 = 𝑔 ( 𝑦 ) 𝑇 ( 𝑢 , 𝑡 ) + 𝑛 ( 𝑥 ) + 𝐼 𝛾 1 𝑅 𝛾 𝑁 𝑧 . ( 2 . 2 4 )
Utilizing this fixed-point formulation, we suggest the following iterative algorithm.

Algorithm 2.4. For given ( 𝑥 0 , 𝑦 0 ) 𝐸 1 × 𝐸 2 , ( 𝑠 0 , 𝑣 0 ) 𝐺 ( 𝑥 0 ) × 𝐹 ( 𝑦 0 ) , ( 𝑢 0 , 𝑡 0 ) 𝐻 ( 𝑥 0 ) × 𝑉 ( 𝑦 0 ) , ( 𝑧 0 , 𝑧 0 ) 𝐸 1 × 𝐸 2 , compute 𝑧 1 𝑥 = 𝑓 0 𝑠 𝑆 0 , 𝑣 0 𝑦 + 𝑚 0 + 𝐼 𝜌 1 𝑅 𝜌 𝑀 𝑧 0 , 𝑧 1 𝑦 = 𝑔 0 𝑢 𝑇 0 , 𝑡 0 𝑥 + 𝑛 0 + 𝐼 𝛾 1 𝑅 𝛾 𝑁 𝑧 0 . ( 2 . 2 5 ) For ( 𝑧 1 , 𝑧 1 ) 𝐸 1 × 𝐸 2 , we take ( 𝑥 1 , 𝑦 1 ) 𝐸 1 × 𝐸 2 such that 𝑓 ( 𝑥 1 ) = 𝐽 𝜌 𝑀 ( 𝑧 1 ) and 𝑔 ( 𝑦 1 ) = 𝐽 𝛾 𝑁 ( 𝑧 1 ) . Then, by Nadler [34], there exist ( 𝑠 1 , 𝑣 1 ) 𝐺 ( 𝑥 1 ) × 𝐹 ( 𝑦 1 ) , ( 𝑢 1 , 𝑡 1 ) 𝐻 ( 𝑥 1 ) × 𝑉 ( 𝑦 1 ) such that 𝑢 1 𝑢 0 𝐻 𝑥 ( 1 + 1 ) 𝐷 1 𝑥 , 𝐻 0 , 𝑣 1 𝑣 0 𝐹 𝑦 ( 1 + 1 ) 𝐷 1 𝑦 , 𝐹 0 , 𝑠 1 𝑠 0 𝐺 𝑥 ( 1 + 1 ) 𝐷 1 𝑥 , 𝐺 0 , 𝑡 1 𝑡 0 𝑉 𝑦 ( 1 + 1 ) 𝐷 1 𝑦 , 𝑉 0 , ( 2 . 2 6 ) where 𝐷 ( , ) is the Hausdorff metric on 𝐶 𝐵 ( 𝐸 1 ) (for the sake of convenience, we also denote by 𝐷 ( , ) the Hausdorff metric on 𝐶 𝐵 ( 𝐸 2 ) ). Compute 𝑧 2 𝑥 = 𝑓 1 𝑠 𝑆 1 , 𝑣 1 𝑦 + 𝑚 1 + 𝐼 𝜌 1 𝑅 𝜌 𝑀 𝑧 1 , 𝑧 2 𝑦 = 𝑔 1 𝑢 𝑇 1 , 𝑡 1 𝑥 + 𝑛 1 + 𝐼 𝛾 1 𝑅 𝛾 𝑁 𝑧 1 . ( 2 . 2 7 ) By induction, we can obtain sequences ( 𝑥 𝑘 , 𝑦 𝑘 ) 𝐸 1 × 𝐸 2 , ( 𝑠 𝑘 , 𝑣 𝑘 ) 𝐺 ( 𝑥 𝑘 ) × 𝐹 ( 𝑦 𝑘 ) , ( 𝑢 𝑘 , 𝑡 𝑘 ) 𝐻 ( 𝑥 𝑘 ) × 𝑉 ( 𝑦 𝑘 ) , ( 𝑧 𝑘 , 𝑧 𝑘 ) 𝐸 1 × 𝐸 2 by the iterative scheme: 𝑓 𝑥 𝑘 = 𝐽 𝜌 𝑀 𝑧 𝑘 , 𝑔 𝑦 𝑘 = 𝐽 𝛾 𝑁 𝑧 𝑘 , 𝑢 𝑘 𝑥 𝐻 𝑘 𝑢 𝑘 + 1 𝑢 𝑘 1 1 + 𝐷 𝐻 𝑥 𝑘 + 1 𝑘 + 1 𝑥 , 𝐻 𝑘 , 𝑣 𝑘 𝑦 𝐹 𝑘 𝑣 𝑘 + 1 𝑣 𝑘 1 1 + 𝐷 𝐹 𝑦 𝑘 + 1 𝑘 + 1 𝑦 , 𝐹 𝑘 , 𝑠 𝑘 𝑥 𝐺 𝑘 𝑠 𝑘 + 1 𝑠 𝑘 1 1 + 𝐷 𝐺 𝑥 𝑘 + 1 𝑘 + 1 𝑥 , 𝐺 𝑘 , 𝑡 𝑘 𝑦 𝑉 𝑘 𝑡 𝑘 + 1 𝑡 𝑘 1 1 + 𝐷 𝑉 𝑦 𝑘 + 1 𝑘 + 1 𝑦 , 𝑉 𝑘 , 𝑧 𝑘 + 1 𝑥 = 𝑓 𝑘 𝑠 𝑆 𝑘 , 𝑣 𝑘 𝑦 + 𝑚 𝑘 + 𝐼 𝜌 1 𝑅 𝜌 𝑀 𝑧 𝑘 , 𝑧 𝑘 + 1 𝑦 = 𝑔 𝑘 𝑢 𝑇 𝑘 , 𝑡 𝑘 𝑥 + 𝑛 𝑘 + 𝐼 𝛾 1 𝑅 𝛾 𝑁 𝑧 𝑘 , ( 2 . 2 8 ) for 𝑘 = 0 , 1 , 2 , .
For positive stepsize 𝛿 , 𝛿 , the general system of generalized resolvent equations (2.1) can also be rewritten as 𝑓 𝑥 , 𝑧 = 𝑓 𝑥 , 𝑧 𝛿 𝑧 𝐽 𝜌 𝑀 𝑧 + 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) = 𝑓 𝑥 , 𝑧 𝛿 𝑓 ( 𝑥 ) 𝐽 𝜌 𝑀 ( , 𝑔 𝑓 ( 𝑥 ) ) + 𝜌 ( 𝑆 ( 𝑠 , 𝑣 ) 𝑚 ( 𝑦 ) ) 𝑦 , 𝑧 = 𝑔 𝑦 , 𝑧 𝛿 𝑧 𝐽 𝛾 𝑁 𝑧 + 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) = 𝑔 𝑦 , 𝑧 𝛿 𝑔 ( 𝑦 ) 𝐽 𝛾 𝑁 ( . 𝑔 ( 𝑦 ) ) + 𝛾 ( 𝑇 ( 𝑢 , 𝑡 ) 𝑛 ( 𝑥 ) ) ( 2 . 2 9 )
This fixed point formulation enables us to propose the following iterative algorithm.

Algorithm 2.5. For given ( 𝑥 0 , 𝑦 0 ) 𝐸 1 × 𝐸 2 , ( 𝑠 0 , 𝑣 0 ) 𝐺 ( 𝑥 0 ) × 𝐹 ( 𝑦 0 ) , ( 𝑢 0 , 𝑡 0 ) 𝐻 ( 𝑥 0 ) × 𝑉 ( 𝑦 0 ) , ( 𝑧 0 , 𝑧 0 ) 𝐸 1 × 𝐸 2 , compute ( 𝑥 1 , 𝑦 1 ) 𝐸 1 × 𝐸 2 and ( 𝑧 1 , 𝑧 1 ) 𝐸 1 × 𝐸 2 such that 𝑓 𝑥 1 , 𝑧 1 𝑥 = 𝑓 0 , 𝑧 0 𝛿 𝑓 𝑥 0 𝐽 𝜌 𝑀 𝑓 𝑥 0 𝑆 𝑠 + 𝜌 0 , 𝑣 0 𝑦 𝑚 0 , 𝑔 𝑦 1 , 𝑧 1 𝑦 = 𝑔 0 , 𝑧 0 𝛿 𝑔 𝑦 0 𝐽 𝛾 𝑁 𝑔 𝑦 0 𝑇 𝑢 + 𝛾 0 , 𝑡 0 𝑥 𝑛 0 . ( 2 . 3 0 ) Then, by Nadler [34], there exist ( 𝑠 1 , 𝑣 1 ) 𝐺 ( 𝑥 1 ) × 𝐹 ( 𝑦 1 ) , ( 𝑢 1 , 𝑡 1 ) 𝐻 ( 𝑥 1 ) × 𝑉 ( 𝑦 1 ) such that 𝑢 1 𝑢 0 𝐻 𝑥 ( 1 + 1 ) 𝐷 1 𝑥 , 𝐻 0 , 𝑣 1 𝑣 0 𝐹 𝑦 ( 1 + 1 ) 𝐷 1 𝑦 , 𝐹 0 , 𝑠 1 𝑠 0 𝐺 𝑥 ( 1 + 1 ) 𝐷 1 𝑥 , 𝐺 0 , 𝑡 1 𝑡 0 𝑉 𝑦 ( 1 + 1 ) 𝐷 1 𝑦 , 𝑉 0 , ( 2 . 3 1 ) where 𝐷 ( , ) is the Hausdorff metric on 𝐶 𝐵 ( 𝐸 1 ) (for the sake of convenience, we also denote by 𝐷 ( , ) the Hausdorff metric on 𝐶 𝐵 ( 𝐸 2 ) ). Compute ( 𝑥 2 , 𝑦 2 ) 𝐸 1 × 𝐸 2 and ( 𝑧 2 , 𝑧 2 ) 𝐸 1 × 𝐸 2 such that 𝑓 𝑥 2 , 𝑧 2 𝑥 = 𝑓 1 , 𝑧 1 𝛿 𝑓 𝑥 1 𝐽 𝜌 𝑀 𝑓 𝑥 1 𝑆 𝑠 + 𝜌 1 , 𝑣 1 𝑦 𝑚 1 , 𝑔 𝑦 2 , 𝑧 2 𝑦 = 𝑔 1 , 𝑧 1 𝛿 𝑔 𝑦 1 𝐽 𝛾 𝑁 𝑔 𝑦 1 𝑇 𝑢 + 𝛾 1 , 𝑡 1 𝑥 𝑛 1 . ( 2 . 3 2 ) By induction, we can obtain sequences ( 𝑥 𝑘 , 𝑦 𝑘 ) 𝐸 1 × 𝐸 2 , ( 𝑠 𝑘 , 𝑣 𝑘 ) 𝐺 ( 𝑥 𝑘 ) × 𝐹 ( 𝑦 𝑘 ) , ( 𝑢 𝑘 , 𝑡 𝑘 ) 𝐻 ( 𝑥 𝑘 ) × 𝑉 ( 𝑦 𝑘 ) , ( 𝑧 𝑘 , 𝑧 𝑘 ) 𝐸 1 × 𝐸 2 by the iterative scheme: 𝑢 𝑘 𝑥 𝐻 𝑘 𝑢 𝑘 + 1 𝑢 𝑘 1 1 + 𝐷 𝐻 𝑥 𝑘 + 1 𝑘 + 1 𝑥 , 𝐻 𝑘 , 𝑣 𝑘 𝑦 𝐹 𝑘 𝑣 𝑘 + 1 𝑣 𝑘 1 1 + 𝐷 𝐹 𝑦 𝑘 + 1 𝑘 + 1 𝑦 , 𝐹 𝑘 , 𝑠 𝑘 𝑥 𝐺 𝑘 𝑠 𝑘 + 1 𝑠 𝑘 1 1 + 𝐷 𝐺 𝑥 𝑘 + 1 𝑘 + 1 𝑥 , 𝐺 𝑘 , 𝑡 𝑘 𝑦 𝑉 𝑘 𝑡 𝑘 + 1 𝑡 𝑘 1 1 + 𝐷 𝑉 𝑦 𝑘 + 1 𝑘 + 1 𝑦 , 𝑉 𝑘 , 𝑓 𝑥 𝑘 + 1 , 𝑧 𝑘 + 1 𝑥 = 𝑓 𝑘 , 𝑧 𝑘 𝑓 𝑥 𝛿 𝑘 𝐽 𝜌 𝑀 𝑓 𝑥 𝑘 𝑆 𝑠 + 𝜌 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 , 𝑔 𝑦 𝑘 + 1 , 𝑧 𝑘 + 1 𝑦 = 𝑔 𝑘 , 𝑧 𝑘 𝛿 𝑔 𝑦 𝑘 𝐽 𝛾 𝑁 𝑔 𝑦 𝑘 𝑇 𝑢 + 𝛾 𝑘 , 𝑡 𝑘 𝑥 𝑛 𝑘 , ( 2 . 3 3 ) for 𝑘 = 0 , 1 , 2 , .
Note that for 𝛿 = 𝛿 = 1 , 𝑓 ( 𝑥 𝑘 , 𝑧 𝑘 ) = 𝑓 ( 𝑥 𝑘 ) , 𝑔 ( 𝑦 𝑘 , 𝑧 𝑘 ) = 𝑔 ( 𝑦 𝑘 ) , Algorithm 2.5 reduces to the following algorithm which solves the general system of variational inclusions (1.14).

Algorithm 2.6. For given ( 𝑥 0 , 𝑦 0 ) 𝐸 1 × 𝐸 2 , ( 𝑠 0 , 𝑣 0 ) 𝐺 ( 𝑥 0 ) × 𝐹 ( 𝑦 0 ) , ( 𝑢 0 , 𝑡 0 ) 𝐻 ( 𝑥 0 ) × 𝑉 ( 𝑦 0 ) , compute ( 𝑥 1 , 𝑦 1 ) 𝐸 1 × 𝐸 2 such that 𝑓 𝑥 1 = 𝐽 𝜌 𝑀 𝑓 𝑥 0 𝑆 𝑠 𝜌 0 , 𝑣 0 𝑦 𝑚 0 , 𝑔 𝑦 1 = 𝐽 𝛾 𝑁 𝑔 𝑦 0 𝑇 𝑢 𝛾 0 , 𝑡 0 𝑥 𝑛 0 . ( 2 . 3 4 ) Then, by Nadler [34], there exist ( 𝑠 1 , 𝑣 1 ) 𝐺 ( 𝑥 1 ) × 𝐹 ( 𝑦 1 ) , ( 𝑢 1 , 𝑡 1 ) 𝐻 ( 𝑥 1 ) × 𝑉 ( 𝑦 1 ) such that 𝑢 1 𝑢 0 𝐻 𝑥 ( 1 + 1 ) 𝐷 1 𝑥 , 𝐻 0 , 𝑣 1 𝑣 0 𝐹 𝑦 ( 1 + 1 ) 𝐷 1 𝑦 , 𝐹 0 , 𝑠 1 𝑠 0 𝐺 𝑥 ( 1 + 1 ) 𝐷 1 𝑥 , 𝐺 0 , 𝑡 1 𝑡 0 𝑉 𝑦 ( 1 + 1 ) 𝐷 1 𝑦 , 𝑉 0 , ( 2 . 3 5 ) where 𝐷 ( , ) is the Hausdorff metric on 𝐶 𝐵 ( 𝐸 1 ) (for the sake of convenience, we also denote by 𝐷 ( , ) the Hausdorff metric on 𝐶 𝐵 ( 𝐸 2 ) ). Compute ( 𝑥 2 , 𝑦 2 ) 𝐸 1 × 𝐸 2 such that 𝑓 𝑥 2 = 𝐽 𝜌 𝑀 𝑓 𝑥 1 𝑆 𝑠 𝜌 1 , 𝑣 1 𝑦 𝑚 1 , 𝑔 𝑦 2 = 𝐽 𝛾 𝑁 𝑔 𝑦 1 𝑇 𝑢 𝛾 1 , 𝑡 1 𝑥 𝑛 1 . ( 2 . 3 6 ) By induction, we can obtain sequences ( 𝑥 𝑘 , 𝑦 𝑘 ) 𝐸 1 × 𝐸 2 , ( 𝑠 𝑘 , 𝑣 𝑘 ) 𝐺 ( 𝑥 𝑘 ) × 𝐹 ( 𝑦 𝑘 ) , ( 𝑢 𝑘 , 𝑡 𝑘 ) 𝐻 ( 𝑥 𝑘 ) × 𝑉 ( 𝑦 𝑘 ) by the iterative scheme: 𝑓 𝑥 𝑘 + 1 = 𝐽 𝜌 𝑀 𝑓 𝑥 𝑘 𝑆 𝑠 𝜌 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 , 𝑔 𝑦 𝑘 + 1 = 𝐽 𝛾 𝑁 𝑔 𝑦 𝑘 𝑇 𝑢 𝛾 𝑘 , 𝑡 𝑘 𝑥 𝑛 𝑘 , 𝑢 𝑘 𝑥 𝐻 𝑘 𝑢 𝑘 + 1 𝑢 𝑘 1 1 + 𝐷 𝐻 𝑥 𝑘 + 1 𝑘 + 1 𝑥 , 𝐻 𝑘 , 𝑣 𝑘 𝑦 𝐹 𝑘 𝑣 𝑘 + 1 𝑣 𝑘 1 1 + 𝐷 𝐹 𝑦 𝑘 + 1 𝑘 + 1 𝑦 , 𝐹 𝑘 , 𝑠 𝑘 𝑥 𝐺 𝑘 𝑠 𝑘 + 1 𝑠 𝑘 1 1 + 𝐷 𝐺 𝑥 𝑘 + 1 𝑘 + 1 𝑥 , 𝐺 𝑘 , 𝑡 𝑘 𝑦 𝑉 𝑘 𝑡 𝑘 + 1 𝑡 𝑘 1 1 + 𝐷 𝑉 𝑦 𝑘 + 1 𝑘 + 1 𝑦 , 𝑉 𝑘 , ( 2 . 3 7 ) for 𝑘 = 0 , 1 , 2 , .

We now study the convergence analysis of Algorithm 2.3. In a similar way, one can study the convergence of other algorithms.

Theorem 2.7. Let 𝐸 1 and 𝐸 2 be two real uniformly smooth Banach spaces with modulus of smoothness 𝜏 𝐸 1 ( 𝑡 ) 𝐶 1 𝑡 2 and 𝜏 𝐸 2 ( 𝑡 ) 𝐶 2 𝑡 2 for 𝐶 1 , 𝐶 2 > 0 , respectively. Let 𝐺 𝐸 1 𝐶 𝐵 ( 𝐸 1 ) , 𝐹 𝐸 2 𝐶 𝐵 ( 𝐸 2 ) , 𝐻 𝐸 1 𝐶 𝐵 ( 𝐸 1 ) , 𝑉 𝐸 2 𝐶 𝐵 ( 𝐸 2 ) be 𝐷 -Lipschitz continuous mappings with constants 𝜆 𝐷 𝐺 , 𝜆 𝐷 𝐹 , 𝜆 𝐷 𝐻 , and 𝜆 𝐷 𝑉 , respectively, and let 𝑀 𝐸 1 2 𝐸 1 , 𝑁 𝐸 2 2 𝐸 2 be 𝑚 -accretive mappings such that the resolvent operators associated with 𝑀 and 𝑁 are retractions. Let 𝑓 𝐸 1 𝐸 1 , 𝑔 𝐸 2 𝐸 2 be both strong accretive with constants 𝛼 and 𝛽 , respectively, and Lipschitz continuous with constants 𝛿 1 and 𝛿 2 , respectively. Let 𝑚 𝐸 2 𝐸 1 , 𝑛 𝐸 1 𝐸 2 be Lipschitz continuous with constants 𝜆 𝑚 and 𝜆 𝑛 , respectively, and 𝑆 𝐸 1 × 𝐸 2 𝐸 1 , 𝑇 𝐸 1 × 𝐸 2 𝐸 2 Lipschitz continuous in the first and second arguments with constants 𝜆 𝑆 1 , 𝜆 𝑆 2 and 𝜆 𝑇 1 , 𝜆 𝑇 2 , respectively.
If there exist constants 𝜌 > 0 and 𝛾 > 0 , such that 0 < 𝐵 + 𝜃 1 + 𝜃 4 𝐵 1 𝐵 < 1 , 0 < + 𝜃 2 + 𝜃 3 1 𝐵 < 1 , ( 2 . 3 8 ) where 𝐵 = 1 2 𝛼 + 6 4 𝐶 1 𝛿 2 1 , 𝐵 = 1 2 𝛽 + 6 4 𝐶 2 𝛿 2 2 , and 𝜃 1 = ( 1 + 𝜌 𝜆 𝑆 1 𝜆 𝐷 𝐺 ) / ( 1 𝜌 ( 𝜆 𝑆 1 𝜆 𝐷 𝐺 + 𝜆 𝑆 2 𝜆 𝐷 𝐹 + 𝜆 𝑚 ) ) , 𝜃 2 = 𝜌 ( 𝜆 𝑆 2 𝜆 𝐷 𝐹 + 𝜆 𝑚 ) / ( 1 𝜌 ( 𝜆 𝑆 1 𝜆 𝐷 𝐺 + 𝜆 𝑆 2 𝜆 𝐷 𝐹 + 𝜆 𝑚 ) ) , 𝜃 3 = ( 1 + 𝛾 𝜆 𝑇 2 𝜆 𝐷 𝑉 ) / ( 1 𝛾 ( 𝜆 𝑇 1 𝜆 𝐷 𝐻 + 𝜆 𝑇 2 𝜆 𝐷 𝑉 + 𝜆 𝑛 ) ) , 𝜃 4 = 𝛾 ( 𝜆 𝑇 1 𝜆 𝐷 𝐻 + 𝜆 𝑛 ) / ( 1 𝛾 ( 𝜆 𝑇 1 𝜆 𝐷 𝐻 + 𝜆 𝑇 2 𝜆 𝐷 𝑉 + 𝜆 𝑛 ) ) , then there exist ( 𝑥 , 𝑦 ) 𝐸 1 × 𝐸 2 , ( 𝑠 , 𝑣 ) 𝐺 ( 𝑥 ) × 𝐹 ( 𝑦 ) , ( 𝑢 , 𝑡 ) 𝐻 ( 𝑥 ) × 𝑉 ( 𝑦 ) and ( 𝑧 , 𝑧 ) 𝐸 1 × 𝐸 2 satisfying the general system of generalized resolvent equations (2.1) (in this case, ( 𝑥 , 𝑦 , 𝑢 , 𝑣 , 𝑠 , 𝑡 ) is a solution of general system of variational inclusions (1.14)), and the iterative sequences { 𝑧 𝑘 } , { 𝑧 𝑘 } , { 𝑥 𝑘 } , { 𝑦 𝑘 } , { 𝑢 𝑘 } , { 𝑣 𝑘 } , { 𝑠 𝑘 } , and { 𝑡 𝑘 } generated by Algorithm 2.3 converge strongly to 𝑧 , 𝑧 , 𝑥 , 𝑦 , 𝑢 , 𝑣 , 𝑠 , and 𝑡 , respectively.

Proof. From Algorithm 2.3 we have 𝑧 𝑘 + 1 𝑧 𝑘 = 𝑓 𝑥 𝑘 𝑆 𝑠 𝜌 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 𝑓 𝑥 𝑘 1 𝑆 𝑠 𝜌 𝑘 1 , 𝑣 𝑘 1 𝑦 𝑚 𝑘 1 𝑥 𝑘 𝑥 𝑘 1 𝑓 𝑥 𝑘 𝑥 𝑓 𝑘 1 + 𝑥 𝑘 𝑥 𝑘 1 𝑆 𝑠 𝜌 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 𝑆 𝑠 𝑘 1 , 𝑣 𝑘 1 𝑦 𝑚 𝑘 1 . ( 2 . 3 9 ) By Proposition 1.1, we have (see, e.g., the proof of [32, Theorem  3]) 𝑥 𝑘 𝑥 𝑘 1 𝑓 𝑥 𝑘 𝑥 𝑓 𝑘 1 2 1 2 𝛼 + 6 4 𝐶 𝛿 2 1 𝑥 𝑘 𝑥 𝑘 1 2 . ( 2 . 4 0 ) Since 𝑆 is Lipschitz continuous in both arguments, 𝐺 , 𝐹 are 𝐷 -Lipschitz continuous, and 𝑚 is Lipschitz continuous, we have 𝑆 𝑠 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 𝑆 𝑠 𝑘 1 , 𝑣 𝑘 1 𝑦 𝑚 𝑘 1 = 𝑆 𝑠 𝑘 , 𝑣 𝑘 𝑠 𝑆 𝑘 1 , 𝑣 𝑘 𝑠 + 𝑆 𝑘 1 , 𝑣 𝑘 𝑠 𝑆 𝑘 1 , 𝑣 𝑘 1 𝑚 𝑦 𝑘 𝑦 𝑚 𝑘 1 𝑆 𝑠 𝑘 , 𝑣 𝑘 𝑠 𝑆 𝑘 1 , 𝑣 𝑘 + 𝑆 𝑠 𝑘 1 , 𝑣 𝑘 𝑠 𝑆 𝑘 1 , 𝑣 𝑘 1 + 𝑚 𝑦 𝑘 𝑦 𝑚 𝑘 1 𝜆 𝑆 1 𝑠 𝑘 𝑠 𝑘 1 + 𝜆 𝑆 2 𝑣 𝑘 𝑣 𝑘 1 + 𝜆 𝑚 𝑦 𝑘 𝑦 𝑘 1 𝜆 𝑆 1 1 1 + 𝑘 𝐷 𝐺 𝑥 𝑘 𝑥 , 𝐺 𝑘 1 + 𝜆 𝑆 2 1 1 + 𝑘 𝐷 𝐹 𝑦 𝑘 𝑦 , 𝐹 𝑘 1 + 𝜆 𝑚 𝑦 𝑘 𝑦 𝑘 1 1 1 + 𝑘 𝜆 𝑆 1 𝜆 𝐷 𝐺 𝑥 𝑘 𝑥 𝑘 1 + 1 1 + 𝑘 𝜆 𝑆 2 𝜆 𝐷 𝐹 𝑦 𝑘 𝑦 𝑘 1 + 𝜆 𝑚 𝑦 𝑘 𝑦 𝑘 1 = 1 1 + 𝑘 𝜆 𝑆 1 𝜆 𝐷 𝐺 𝑥 𝑘 𝑥 𝑘 1 + 1 1 + 𝑘 𝜆 𝑆 2 𝜆 𝐷 𝐹 + 𝜆 𝑚 𝑦 𝑘 𝑦 𝑘 1 . ( 2 . 4 1 ) Utilizing (2.41) and Proposition 1.1, we have 𝑥 𝑘 𝑥 𝑘 1 𝑆 𝑠 𝜌 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 𝑆 𝑠 𝑘 1 , 𝑣 𝑘 1 𝑦 𝑚 𝑘 1 2 𝑥 𝑘 𝑥 𝑘 1 2 𝑆 𝑠 2 𝜌 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 𝑆 𝑠 𝑘 1 , 𝑣 𝑘 1 𝑦 𝑚 𝑘 1 , 𝐽 𝑥 𝑘 𝑥 𝑘 1 𝑆 𝑠 𝜌 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 𝑆 𝑠 𝑘 1 , 𝑣 𝑘 1 𝑦 𝑚 𝑘 1 𝑥 𝑘 𝑥 𝑘 1 2 𝑆 𝑠 + 2 𝜌 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 𝑆 𝑠 𝑘 1 , 𝑣 𝑘 1 𝑦 𝑚 𝑘 1 × 𝑥 𝑘 𝑥 𝑘 1 𝑆 𝑠 𝜌 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 𝑆 𝑠 𝑘 1 , 𝑣 𝑘 1 𝑦 𝑚 𝑘 1 𝑥 𝑘 𝑥 𝑘 1 2 1 + 2 𝜌 1 + 𝑘 𝜆 𝑆 1 𝜆 𝐷 𝐺 𝑥 𝑘 𝑥 𝑘 1 + 1 1 + 𝑘 𝜆 𝑆 2 𝜆 𝐷 𝐹 + 𝜆 𝑚 𝑦 𝑘 𝑦 𝑘 1 × 𝑥 𝑘 𝑥 𝑘 1 𝑆 𝑠 𝜌 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 𝑆 𝑠 𝑘 1 , 𝑣 𝑘 1 𝑦 𝑚 𝑘 1 𝑥 𝑘 𝑥 𝑘 1 2 1 + 𝜌 1 + 𝑘 𝜆 𝑆 1 𝜆 𝐷 𝐺 × 𝑥 𝑘 𝑥 𝑘 1 2 + 𝑥 𝑘 𝑥 𝑘 1 𝑆 𝑠 𝜌 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 𝑆 𝑠 𝑘 1 , 𝑣 𝑘 1 𝑦 𝑚 𝑘 1 2 1 + 𝜌 1 + 𝑘 𝜆 𝑆 2 𝜆 𝐷 𝐹 + 𝜆 𝑚 × 𝑦 𝑘 𝑦 𝑘 1 2 + 𝑥 𝑘 𝑥 𝑘 1 𝑆 𝑠 𝜌 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 𝑆 𝑠 𝑘 1 , 𝑣 𝑘 1 𝑦 𝑚 𝑘 1 2 = 1 1 + 𝜌 1 + 𝑘 𝜆 𝑆 1 𝜆 𝐷 𝐺 𝑥 𝑘 𝑥 𝑘 1 2 1 + 𝜌 1 + 𝑘 𝜆 𝑆 1 𝜆 𝐷 𝐺 + 𝜆 𝑆 2 𝜆 𝐷 𝐹 + 𝜆 𝑚 × 𝑥 𝑘 𝑥 𝑘 1 𝑆 𝑠 𝜌 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 𝑆 𝑠 𝑘 1 , 𝑣 𝑘 1 𝑦 𝑚 𝑘 1 2 1 + 𝜌 1 + 𝑘 𝜆 𝑆 2 𝜆 𝐷 𝐹 + 𝜆 𝑚 𝑦 𝑘 𝑦 𝑘 1 2 , ( 2 . 4 2 ) which implies that 𝑥 𝑘 𝑥 𝑘 1 𝑆 𝑠 𝜌 𝑘 , 𝑣 𝑘 𝑦 𝑚 𝑘 𝑆 𝑠 𝑘 1 , 𝑣 𝑘 1 𝑦 𝑚 𝑘 1 2 1 + 𝜌 ( 1 + 1 / 𝑘 ) 𝜆 𝑆 1 𝜆 𝐷 𝐺