Higher Order Strongly Biconvex Functions and Biequilibrium Problems

Show more

1. Introduction

Variational inequality theory, which was introduced and studied by Stampacchia [1] can be viewed as an important and significant extension of the variational principles, the origin of which can be traced back to Euler, Lagrange, Newton, and Bernoulli brothers. This theory has had a great impact and influence in the development of several branches of pure, applied, and engineering sciences. It has been shown that variational inequalities include complementarity problems, fixed-point, optimization problems, and game theory as special cases. Variational inequalities have drawn much attention from researchers because of their importance as a core model for studying many mathematical problems which include convex programming, equilibrium problem, inclusion problem, split feasibility problem, complementarity problem, minimization problem, see [2] [3].

It is well known that the minimum of differentiable convex function on a convex set can be characterized by an inequality, which is called the variational inequality. This simple fact has influenced several branches of pure and applied sciences. For the formulation, applications, numerical methods, and other aspects of the variational inequalities, see [4] - [17] and the references therein. Variational inequalities can be viewed as a novel and significant extension of the variational principles. Blum and Oettlie [18] and Noor and Oettlie [19] introduced the concept of equilibrium, which has appeared as an interesting generalization of the variational inequalities. Equilibrium problems include variational inequalities and fixed point problems as special cases. Several equilibrium problems provide us a unified framework to study unrelated arising in transportation, economics, and management sciences in a unified framework.

In recent years, several extensions and generalizations of the convex sets and convex functions have been considered and investigated. Noor *et al*. [20] [21] introduced the concept of biconvex sets and biconvex functions. It is known that the differentiable biconvex functions are biconvex functions. The converse also holds under certain conditions, see [20]. Noor *et al*. [20] [21] proved that the minimum of the differentiable biconvex functions on the biconvex set can be characterized by a class of variational inequalities, which is known as the bivariational inequality. Noor [20] [21] studied the biconvex equilibrium problems, which can be viewed as an interesting and significant extension of the variational inequalities. For the recent developments in bivariational inequalities and biconvex equilibrium problems, see [20] [21] [22] and the references therein. These results have inspired a great deal of subsequent work, which has expanded the role and applications of convexity in nonlinear optimization and engineering sciences.

In many problems, a set may not be a convex set. To overcome this drawback, the underlying set can be made a *k*-convex set with respect to an arbitrary function. Micherda *et al*. [23], Hazy [24], and Crestescu *et al*. [25] defined the so- called
$\left(h\mathrm{,}k\right)$ convex function which is a natural generalization of the usual convexity, the *s*-convexity in the first and second sense. It is worth mentioning that
$k\left(t\right)=t{\text{e}}^{i\phi}$, the *φ*-convex functions were introduced and studied by Noor [26].

We would like to point out that biconvex sets, biconvex functions, *k*-convex sets, and *k*-convex functions are different generalizations of convex sets and convex functions in various directions. These types of functions have played a leading role in the developments of various branches of pure and applied sciences. It is natural to unify these classes of functions and to investigate their characterizations. Motivated and inspired by the recent activities in these areas, we introduce some new classes of biconvex sets and biconvex functions involving an arbitrary function *k*, which are called modified *k*-biconvex sets and *k*-biconvex functions. These new classes of *k*-biconvex sets and *k*-biconvex functions include the *φ*-biconvex sets, *φ*-biconvex and Toader type *k*-convex sets and *k*-convex functions. These *k*-biconvex functions can be viewed as modified refinement of the
$\left(h\mathrm{,}k\right)$ convex functions of Micherda *et al*. [23] and Hazy [24]. Several new concepts are defined and their properties have been studied. We prove that the minimum of the differential *k*-biconvex functions on the *k*- biconvex sets can be characterized by a class of equilibrium problems. This result inspired us to consider the biequilibrium problems.

As a result of interaction between different branches of mathematical and engineering sciences, we now have a variety of techniques to suggest and analyze various iterative algorithms for solving equilibrium problems, variational inequalities and related optimization problems. Analysis of these problems requires a blend of techniques and ideas from convex analysis, functional analysis, numerical analysis and nonsmooth analysis. There are several methods for solving variational inequalities and equilibrium problems. Due to the nature of the biequilibrium problems, projection and resolvent methods can not be applied for solving directional equilibrium-like problems. To overcome these difficulties, one usually uses the auxiliary principle technique, which is mainly due to Glowinski *et al*. [4]. The main idea involving this technique is to first consider an auxiliary problem and then to show that the solution of the auxiliary problem is the solution of the original problem by using the fixed-point approach. Noor [9] [11] [12] [19] and Noor *et al*. [13] [14] [15] [16] have used this approach to suggest and analyze some iterative methods for solving various classes of variational inequalities and equilibrium problems. It is worth mentioning that one can use several principles to suggest auxiliary problem associated with the original problem. Bregman [27] introduced the convex functions involving the distance functions and discussed its applications in convex programming and optimization problems. Noor *et al*. [13] [14] [15] [16] and Zhu and Marcote [17] used the auxiliary principle involving the Bregman functions. They have suggested and analyzed a wide class of iterative method for solving variational inequalities and equilibrium problems. They have shown [2] [3] that the Bregman functions have some practical important types of functions such as Burg entropy, which is very important in information theory and Shannon entropy, having applications in several areas of applied mathematics such as machine learning. This shows that the entropy methods have been used in the development of methods for solving equilibrium problems in various ways. Thus one concluded the entropy technique has been used effectively in studying the variational inequalities and equilibrium problems. We again show that the auxiliary principle technique involving Bregman functions can be used to suggest some iterative schemes for solving the higher order biequilibrium problems. We prove that the convergence of these methods requires either pseudomonotonicity or partially relaxed strongly monotonicity. These are weaker conditions than monotonicity. As special cases, we obtain iterative schemes for solving biequilibrium problems and related optimization problems. The comparison of these methods with other methods is a subject of future

In section 2, we introduce and study some new concepts regarding higher order strongly biconvex functions and their properties are discussed. Main properties of the higher order strongly biconvex are discussed in Section 3. Various parallelograms have been derived in Section 4 as applications of higher order strongly *k*-biconvex affine functions, which can be used to characterize several inner product and Banach spaces. In Section 5, we show that the optimality conditions of the higher order strongly generalize biconvex functions can be characterized by biequilibrium problems. Using the auxiliary technique coupled with Bregman functions, we have suggested and investigated several iterative methods for solving the higher order strongly biequilibrium problems. Convergence of the proposed methods is considered using the pseudomonotone operators, which is weaker condition than monotonicity. Several special cases are discussed as applications of our results. Comparison with other methods and implementation of the suggested methods need further efforts. Several open problems are discussed for future research.

2. Preliminaries

Let
${K}_{k}$ be a nonempty closed set in a normed space *H*. We denote by
$\langle \cdot \mathrm{,}\cdot \rangle $ and
$\Vert \text{\hspace{0.05em}}\cdot \text{\hspace{0.05em}}\Vert $ the inner product and norm, respectively.

Definition 2.1.* *The set*
${K}_{k\beta}$ *is said to be k-biconvex set with respect to arbitrary function* k *and the bifunction*
$\beta \left(\mathrm{.}-\mathrm{.}\right)$ *,* *if

$u+k\left(\lambda \right)\beta \left(v-u\right)\in {K}_{k\beta}\mathrm{,}\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{,}\text{\hspace{1em}}\lambda \in \left[\mathrm{0,1}\right]\mathrm{.}$ (2.1)

Clearly, for $k\left(\lambda \right)=\lambda $, the set ${K}_{k\beta}$ is a biconvex set ${K}_{\beta}$, which was introduced and studied by Noor and Noor [20].

Definition 2.2. [20] [21] The set* K *is said to be biconvex set with respect to arbitrary the bifunction*
$\beta \left(\mathrm{.}-\mathrm{.}\right)$ *, if

$u+\lambda \beta \left(v-u\right)\in {K}_{\beta},\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{\beta}\mathrm{,}\text{\hspace{1em}}\lambda \in \left[\mathrm{0,1}\right]\mathrm{.}$ (2.2)

If $\beta \left(v\mathrm{,}u\right)=v-u$, then (2.1) reduces to:

Definition 2.3. [23] [24] * *The set*
${K}_{k}$ *is said to be k-convex set with respect to arbitrary function* k*, if

$u+k\left(\lambda \right)\left(v-u\right)\in {K}_{k},\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k}\mathrm{,}\text{\hspace{1em}}\lambda \in \left[\mathrm{0,1}\right]\mathrm{.}$ (2.3)

which was introduced by Hazy [13]. Also see [25] for relevant literature.

From now onwards, the set
${K}_{k\beta}$ is a *k*-biconvex set, unless otherwise specified.

We now introduce the concept of higher order generalized biconvex function with respect to an arbitrary function *k* and bifunction
$\beta \left(\mathrm{.}-\mathrm{.}\right)$.

Definition 2.4.* *The function *f* on*
${K}_{k\eta}$ *is called higher order generalized biconvex function, if there exists an arbitrary function *k*, bifunction*
$\beta \left(\mathrm{.}-\mathrm{.}\right)$ * and a constant*
$\nu $ *,* *such that

$\begin{array}{l}f\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\le \left(1-k\left(\lambda \right)\right)f\left(u\right)+k\left(\lambda \right)f\left(v\right)\\ \text{\hspace{0.05em}}-\nu \left\{{\left(k\left(\lambda \right)\right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{,}\\ \forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{,\text{\hspace{0.17em}}}\lambda \in \left[\mathrm{0,1}\right]\mathrm{.}\end{array}$ (2.4)

If $k\left(\lambda \right)=\lambda $, then definition 2.4 reduces to:

Definition 2.5.* *The function* f *on*
${K}_{k\beta}$ *is called higher order generalized biconvex function, if there exists an arbitrary function* k*, bifunction*
$\beta \left(\mathrm{.}-\mathrm{.}\right)$ *and* *a constant
$\nu $, such that

$\begin{array}{l}f\left(u+\lambda \beta \left(v-u\right)\right)\le \left(1-\lambda \right)f\left(u\right)+\lambda f\left(v\right)\\ \text{\hspace{0.05em}}-\nu \left\{{\lambda}^{p}\left(1-\lambda \right)+\lambda {\left(1-\lambda \right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{,}\\ \forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{,}\lambda \in \left[\mathrm{0,1}\right]\mathrm{.}\end{array}$ (2.5)

The higher order biconvex functions were introduced and studied by Noor *et al*. [20] [21]. Obviously every higher order
$k\left(\lambda \right)$ biconvex function with
$k\left(\lambda \right)=\lambda $ is a higher order biconvex function, but the converse may is not true.

If
$\beta \left(v\mathrm{,}u\right)=v-u$, then higher order *k*-biconvex functions collapse to:

Definition 2.6. The function* f *on*
${K}_{k}$ *is called higher order* **k-*convex function, if there exists an arbitrary function* **k*,* *and a constant*
$\nu $ *such that

$\begin{array}{l}f\left(u+k\left(\lambda \right)\left(v-u\right)\right)\le \left(1-k\left(\lambda \right)\right)f\left(u\right)+k\left(t\right)f\left(v\right)\\ \text{\hspace{0.05em}}-\nu \left\{{\left(k\left(\lambda \right)\right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert v-u\Vert}^{p}\mathrm{,}\\ \forall u\mathrm{,}v\in {K}_{k}\mathrm{,}\lambda \in \left[\mathrm{0,1}\right]\mathrm{.}\end{array}$ (2.6)

If $k\left(\lambda \right)=\lambda ,\beta \left(v,u\right)=v-u$, then definition 2.6 reduces to:

Definition 2.7.* *The function* f *on* K *is called higher order strongly convex function, if there exists a constant*
$\nu $ *,* *such that

$\begin{array}{l}f\left(u+\lambda \left(v-u\right)\right)\le \left(1-\lambda \right)f\left(u\right)+\lambda f\left(v\right)\\ \text{\hspace{0.05em}}-\nu \left\{{\lambda}^{p}\left(1-\lambda \right)+\lambda {\left(1-\lambda \right)}^{p}\right\}{\Vert v-u\Vert}^{p},\\ \forall u,v\in K,\lambda \in \left[0,1\right],p>1.\end{array}$ (2.7)

Higher order convex functions and their variant forms were introduced and studied by Mohsin *et al**.* [28] and Noor *et al*. [15] [22] recently. It is worth mentioning that these concepts were also represent significant refinement and improvement of the concepts of higher order functions considered in Alabdali *et al*. [29], Olbrys [30], Hazy [24], Lin *e**t al*. [6] and Mako *et al*. [31]. It has been shown [15] [28] that parallelogram laws for Banach spaces can be obtained as novel applications of the higher order convex affine functions, which have application in various fields of mathematical and engineering sciences including prediction theory.

For
$\lambda =1$, the *k*-convex function reduces to:

Condition *A* The function *f* is said to satisfy

$f\left(u+k\left(1\right)\beta \left(v-u\right)\right)\le f\left(v\right)\mathrm{,}\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{,}$ (2.8)

which is called the condition *A*.

Definition 2.8.* *The function *f* on
${K}_{k\beta}$ is said to be higher order strongly generalized quasi *k*-biconvex function, if there exist a function* k* and a constant
$\nu $ such that

$\begin{array}{l}f\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\le \mathrm{max}\left\{f\left(u\right),f\left(v\right)\right\}\\ \text{\hspace{0.05em}}-\nu \left\{{\left(k\left(\lambda \right)\right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p},\\ \forall u,v\in {K}_{k\beta},\lambda \in \left[0,1\right],p>1.\end{array}$ (2.9)

Definition 2.9. The function* f *on*
${K}_{k\beta}$ *is said to be higher order strongly generalized logarithmic biconvex function, if there exist a function* k *and the bifunction*
$\eta \left(\mathrm{.}-\mathrm{.}\right)$ *, such that

$\begin{array}{l}f\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\le {\left(f\left(u\right)\right)}^{1-k\left(\lambda \right)}{\left(f\left(v\right)\right)}^{k\left(\lambda \right)}\\ \text{\hspace{0.05em}}-\nu \left\{{\left(k\left(\lambda \right)\right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p},\\ \forall u,v\in {K}_{k\beta},\lambda \in \left[0,1\right],p>1.\end{array}$ (2.10)

*where
$f(\cdot )>0$. *

We note that, if $\mu =0$, then definition 2.9 reduces to:

Definition 2.10. The function* f *on*
${K}_{k\beta}$ *is said to be logarithmic biconvex function, if there exist a function* k*, such that

$f\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\le {\left(f\left(u\right)\right)}^{1-k\left(\lambda \right)}{\left(f\left(v\right)\right)}^{k\left(\lambda \right)}\mathrm{,}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{,\text{\hspace{0.17em}}}\lambda \in \left[\mathrm{0,1}\right]\mathrm{,}$ (2.11)

or equivalently, it can be written as:

Definition 2.11 The function* f *on*
${K}_{k\beta}$ *is said to be higher order strongly logarithmic generalized biconvex function, if there exist a function *k*,* *such that

$\begin{array}{l}\mathrm{log}f\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\\ \le \left(1-k\left(\lambda \right)\right)\mathrm{log}f\left(u\right)+k\left(\lambda \right)\mathrm{log}f\left(v\right)\mathrm{,}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{,\text{\hspace{0.17em}}}\lambda \in \left[\mathrm{0,1}\right]\mathrm{,}\end{array}$ (2.12)

Using this idea, we can define the higher order logarithmic generalized biconvex functions as:

Definition 2.12. The function* f *on*
${K}_{k\beta}$ *is said to be higher order strongly logarithmic generalized biconvex function, if there exist a function* k*, such that

$\mathrm{log}f\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\le \left(1-k\left(\lambda \right)\right)\mathrm{log}f\left(u\right)+k\left(\lambda \right)\mathrm{log}f\left(v\right)\mathrm{,}$

$\text{\hspace{0.05em}}-\nu \left\{{\left(k\left(\lambda \right)\right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}$ (2.13)

$\forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{,\text{\hspace{0.17em}}}\lambda \in \left[\mathrm{0,1}\right]\mathrm{,}$ (2.14)

which appears to be a new one.

From the above definitions, we have:

$\begin{array}{l}f\left(u+k\left(\lambda \right)\beta \left(v,u\right)\right)\\ \le {\left(f\left(u\right)\right)}^{1-k\left(\lambda \right)}{\left(f\left(v\right)\right)}^{k\left(\lambda \right)}-\nu \left\{{\left(k\left(\lambda \right)\right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\\ \le \left(1-k\left(\lambda \right)\right)f\left(u\right)+k\left(\lambda \right)f\left(v\right)\\ \text{\hspace{0.05em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\nu \left\{{\left(k\left(\lambda \right)\right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\\ \le \mathrm{max}\left\{f\left(u\right),f\left(v\right)\right\}-\nu \left\{{\left(k\left(\lambda \right)\right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p},\\ \text{\hspace{1em}}\forall u,v\in {K}_{k\beta},\text{\hspace{0.17em}}\lambda \in \left[0,1\right],\end{array}$ (2.15)

This shows that higher order strongly logarithmic generalized biconvex function $\Rightarrow $ higher order strongly generalized biconvex functions and higher order strongly generalized biconvex functions $\Rightarrow $ higher order strongly generalized quasi biconvex functions, but the converse is not true.

We remark that, if *F* is both higher order strongly biconvex function and higher order strongly biconcave function, then we can define the following new concept.

Definition 2.13. The function *f *on*
${K}_{k\beta}$ *is called higher order strongly biconvex affine function, if there* *exist an arbitrary function* k* and a constant*
$\nu $ *such that

$\begin{array}{l}f\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)=\left(1-k\left(\lambda \right)\right)f\left(u\right)+k\left(\lambda \right)f\left(v\right)\\ \text{\hspace{0.05em}}-\nu \left\{{\left(k\left(\lambda \right)\right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p},\\ \forall u,v\in {K}_{k\beta},\mathrm{\text{\hspace{0.17em}}}\lambda \in \left[0,1\right],p>1.\end{array}$ (2.16)

For proper and suitable choice of the functions $k\left(\lambda \right)$ and the bifunction $\beta \left(\mathrm{.}-\mathrm{.}\right)$, one can a wide classes of higher order generalized biconvex functions and their variant forms, see [32] [33] [34] [35]. This shows that the concepts introduced in this paper are quite flexible and general ones.

We also need the following assumption regarding the bifunction $\eta \left(\mathrm{.}-\mathrm{.}\right)$ and the function $k\left(\lambda \right)$.

Condition M. Let $\beta \left(\mathrm{.}-\mathrm{.}\right)\mathrm{:}{K}_{k\beta}\times {K}_{k\beta}\to H$ satisfy assumptions:

1) $\beta \left(k\left(\lambda \right)\beta \left(v-u\right)\right)=-k\left(\lambda \right)\beta \left(v-u\right)$ (2.17)

2) $\beta \left(v-u-k\left(\lambda \right)\beta \left(v-u\right)\right)=\left(1-k\left(\lambda \right)\right)\beta \left(v-u\right)\mathrm{,}\forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{,}\lambda \in \left[\mathrm{0,1}\right]\mathrm{.}$ (2.18)

Remark 2.1. Let

$\beta \left(v-u\right)=\beta \left(v-z\right)+\beta \left(z-u\right)\mathrm{,}\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{.}$

Clearly*
$\beta \left(u-v\right)=0$ *,if and only if,*
$u=v,\forall u,v\in {K}_{k\beta}$ *,and

$\beta \left(v-u\right)+\beta \left(u-v\right)=\mathrm{0,}\forall u\mathrm{,}v\in {k}_{k\beta}\mathrm{.}$

that is, the bifunction*
$\beta \left(\mathrm{,}-\mathrm{,}\right)$ *is skew symmetric*. *

We now introduce the concept of *k*-directional derivative, which is mainly due to Noor [26].

Definition 2.14.* *We define the k-directional derivative of f at a point*
$u\text{\hspace{0.17em}}\in \text{\hspace{0.17em}}{K}_{k\beta}$ *in the direction*
$v\in {K}_{k\beta}$ *by

$\begin{array}{c}Df\left(u\mathrm{,}\beta \left(v-u\right)\right)={{f}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)\\ =\underset{\lambda \to {0}^{+}}{\mathrm{lim}}\left\{\frac{f\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)-f\left(u\right)}{k\left(\lambda \right)}\right\}\mathrm{.}\end{array}$ (2.19)

Note that for
$k\left(\lambda \right)=\lambda $ and
$\beta \left(v-u\right)=v$, the *k*-directional derivative of *f* at
$u\in K$ in the direction
$v\in K$ coincides with the usual directional derivative of *f* at *u* in a direction *v* given by

$Df\left(u,v\right):={f}^{\prime}\left(u,v\right)=\underset{\lambda \to {0}^{+}}{\mathrm{lim}}\frac{f\left(u+tv\right)-f\left(u\right)}{\lambda}.$

It is well known that the function $v\to {{f}^{\prime}}_{k}\left(u\mathrm{,}v-u\right)$ is subadditive, positively homogeneous.

Definition 2.15.* *A differentiable function *f* on the biconvex set*
${K}_{k\beta}$ *is said to be higher order strongly pseudo *k*-biconvex function, iff, if there exists a constant
$\mu >0$ such that

${{f}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)+\mu {\Vert \beta \left(v-u\right)\Vert}^{p}\ge 0\Rightarrow f\left(v\right)-f\left(u\right)\ge \mathrm{0,}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{,}p\ge 1.$

Definition 2.16.* *A differentiable function *f* on
${K}_{k\beta}$ is said to be higher order strongly quasi-*k*-biconvex function, iff, if there exists a constant
$\mu >0$ such that

$f\left(v\right)\le f(\; u\; )$

$\Rightarrow $

${{f}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)+\mu {\Vert \beta \left(v-u\right)\Vert}^{p}\le \mathrm{0,}\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{,}p\ge 1.$

Definition 2.17.* *The function *f* on the set
${K}_{k\beta}$ is said to be pseudo-
$k\beta $ - biconvex, iff,

${{f}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)\ge 0\Rightarrow f\left(v\right)\ge f\left(u\right)\mathrm{,}\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{.}$

Definition 2.18.* *The differentiable function *f* on the
${K}_{k\beta}$ is said to be quasi *k*-biconvex function, iff,

$f\left(v\right)\le f\left(u\right)\Rightarrow {{f}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)\le \mathrm{0,}\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{.}$

I All these new concepts may play important and fundamental part in the developments of mathematical programming and optimization theory.

3. Properties of Biconvex Functions

In this section, we consider some basic properties of higher order strongly generalized biconvex functions on the biconvex set ${K}_{k\beta}$.

Theorem 3.1. Let *F* be a differentiable function on the biconvex set
${K}_{k\beta}$ in *H* and let the condition *M* hold. Then the function *F* is higher order generalized biconvex function, if and only if,

$F\left(v\right)-F\left(u\right)\ge {{F}^{\prime}}_{k\beta}\left(u,\beta \left(v-u\right)\right)+\mu {\Vert \beta \left(v-u\right)\Vert}^{p},\text{\hspace{0.17em}}\forall u,v\in {K}_{k\beta},t\in \left[0,1\right],p>1.$ (3.1)

Proof. Let *f* be a higher order generalized biconvex function on biconvex set
${K}_{k\beta}$. Then,

$\begin{array}{l}F\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\le \left(1-k\left(\lambda \right)\right)F\left(u\right)+\lambda F\left(v\right)\\ \text{\hspace{0.05em}}-\mu \left\{k{\left(\lambda \right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{,}\\ \forall u\mathrm{,}v\in {K}_{k\beta},\lambda \in \left[\mathrm{0,1}\right]\mathrm{,}p\ge \mathrm{1,}\end{array}$

which can be written as:

$\begin{array}{l}F\left(v\right)-F\left(u\right)\ge \left\{\frac{F\left(u+\lambda \beta \left(v-u\right)\right)-F\left(u\right)}{k\left(\lambda \right)}\right\}\\ +\mu \left\{k{\left(\lambda \right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{.}\end{array}$

Taking the limit in the above inequality as $\lambda \to 0$, we have

$F\left(v\right)-F\left(u\right)\ge {{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)+\mu {\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{,}$

which is the required (20).

Conversely, let *F* satisfy (20). Then,
$\forall u\mathrm{,}v\in {K}_{k\beta},\lambda \in \left[\mathrm{0,1}\right]$,
${v}_{\lambda}=u+k\left(\lambda \right)\beta \left(v-u\right)\in {K}_{k\beta}$ and using the condition *M*, we have

$\begin{array}{l}F\left(v\right)-F\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\\ \ge \langle {F}^{\prime}\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\mathrm{,}\eta \left(v\mathrm{,}u+k\left(\lambda \right)\beta \left(v-u\right)\right)\rangle \\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}+\mu {\Vert \eta \left(v-u-k\left(\lambda \right)\beta \left(v-u\right)\right)\Vert}^{p}\\ =\left(1-k\left(\lambda \right)\right)\langle {{F}^{\prime}}_{k\beta}\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\mathrm{,}\beta \left(v-u\right)\rangle \\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}+\mu {\left(1-k\left(\lambda \right)\right)}^{p}{\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{.}\end{array}$ (3.2)

In a similar way, we have

$\begin{array}{l}F\left(u\right)-F\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\\ \ge \langle {{F}^{\prime}}_{k\beta}\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\mathrm{,}\beta \left(u-u-k\left(\lambda \right)\beta \left(v-u\right)\right)\rangle \\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\mu {\Vert \eta \left(u-u-k\left(\lambda \right)\beta \left(v-u\right)\right)\Vert}^{p}\\ =-k\left(\lambda \right)\langle {{F}^{\prime}}_{k\beta}\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\mathrm{,}\beta \left(v-u\right)\rangle +\mu k{\left(\lambda \right)}^{p}{\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{.}\end{array}$ (3.3)

Multiplying (3.2) by $k\left(\lambda \right)$ and (3.3) by $\left(1-k\left(\lambda \right)\right)$ and adding the resultant, we have

$\begin{array}{l}F\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\le \left(1-k\left(\lambda \right)\right)F\left(u\right)+k\left(\lambda \right)F\left(v\right)\\ \text{\hspace{0.05em}}-\nu \left\{{\left(k\left(\lambda \right)\right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{,}\end{array}$

showing that *F* is a higher order strongly biconvex function.

Theorem 3.2.* *Let *F* be differentiable higher order generalized biconvex function on the biconvex set
${K}_{k\beta}$. If *F* is a higher order strongly biconvex function, then,

$\begin{array}{l}{{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)+{{F}^{\prime}}_{k\beta}\left(v\mathrm{,}\beta \left(u-v\right)\right)\\ \le -\mu \left\{{\Vert \beta \left(v-u\right)\Vert}^{p}+{\Vert \beta \left(u-v\right)\Vert}^{p}\right\}\mathrm{,}\forall u\mathrm{,}v\in {K}_{k\beta}.\end{array}$ (3.4)

Proof. Let *F* be a higher order strongly biconvex function on the biconvex set
${K}_{k\beta}$. Then,

$F\left(v\right)-F\left(u\right)\ge {{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)+\mu {\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{,}\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k\beta}.$ (3.5)

Changing the role of *u* and *v* in (24), we have

$F\left(u\right)-F\left(v\right)\ge {{F}^{\prime}}_{k\beta}\left(v\mathrm{,}\beta \left(u-v\right)\right)+\mu {\Vert \beta \left(u-v\right)\Vert}^{p}\mathrm{,}\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k\beta}.$ (3.6)

Adding (24) and (25), we have

$\begin{array}{l}{{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)+{{F}^{\prime}}_{k\beta}\left(v\mathrm{,}\beta \left(u-v\right)\right)\\ \le -\mu \left\{{\Vert \beta \left(v-u\right)\Vert}^{p}+{\Vert \beta \left(u-v\right)\Vert}^{p}\right\}\mathrm{,}\forall u\mathrm{,}v\in {K}_{k\beta}.\end{array}$ (3.7)

which shows that ${{F}^{\prime}}_{k\beta}\left(\mathrm{.}\right)$ is a higher order strongly monotone operator.

Theorem 3.3.* *If the differential
${{F}^{\prime}}_{k\beta}\left(\mathrm{.}\right)$ is a higher order strongly monotone, then,

$F\left(v\right)-F\left(u\right)\ge {{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)+\mu {\Vert \beta \left(v-u\right)\Vert}^{p}{\displaystyle {\int}_{0}^{1}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}k{\left(\lambda \right)}^{p-1}\text{d}\lambda \mathrm{.}$

Proof. Let ${{F}^{\prime}}_{k\beta}\left(\mathrm{.}\right)$ be a higher order strongly $k\beta $ -monotone. From (26), we have

${{F}^{\prime}}_{k\beta}\left(v\mathrm{,}\beta \left(u-v\right)\right)\ge {{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)-\mu \left\{{\Vert \beta \left(v-u\right)\Vert}^{p}+{\Vert \beta \left(u-v\right)\Vert}^{p}\right\}\mathrm{.}$ (3.8)

Since *K* is an biconvex set,
$\forall u\mathrm{,}v\in {K}_{k\beta}$,
$\lambda \in \left[\mathrm{0,1}\right]$,
${v}_{\lambda}\text{\hspace{0.17em}}=u+k\left(\lambda \right)\beta \left(v-u\right)\in {K}_{k\beta}$. Taking
$v={v}_{\lambda}$ in (27) and using Condition *M*, we have

$\begin{array}{l}{{F}^{\prime}}_{k\beta}\left({v}_{\lambda}\mathrm{,}\beta \left(u\mathrm{,}u+k\left(\lambda \right)\beta \left(v-u\right)\right)\right)\\ \le {{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(u+k\left(\lambda \right)\beta \left(v-u\right)\mathrm{,}u\right)\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}-\mu \left\{{\Vert \beta \left(u+k\left(\lambda \right)\beta \left(v-u\right)\mathrm{,}u\right)\Vert}^{p}+{\Vert \beta \left(u\mathrm{,}u+k\left(\lambda \right)\beta \left(v-u\right)\right)\Vert}^{p}\right\}\\ \mathrm{=}-k\left(\lambda \right){{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)-2k{\left(\lambda \right)}^{p}\mu {\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{,}\end{array}$

which implies that

${{F}^{\prime}}_{k\beta}\left({v}_{\lambda}\mathrm{,}\beta \left(v-u\right)\right)\ge {{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)+2\mu k{\left(\lambda \right)}^{p-1}{\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{.}$ (3.9)

Let $\xi \left(\lambda \right)=F\mathrm{(}u+k\left(\lambda \right)\beta \left(v-u\right)$. Then, from (28), we have

$\begin{array}{c}{\xi}^{\prime}\left(\lambda \right)={{F}^{\prime}}_{k}\left(u+k\left(\lambda \right)\beta \left(v-u\right)\mathrm{,}\beta \left(v-u\right)\right)\\ \ge {{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)+2\mu k{\left(\lambda \right)}^{p-1}{\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{.}\end{array}$ (3.10)

Integrating (3.10) between 0 and 1, we have

$\xi \left(1\right)-\xi \left(0\right)\ge {{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)+\mu {\Vert \beta \left(v-u\right)\Vert}^{p}{\displaystyle {\int}_{0}^{1}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}k{\left(\lambda \right)}^{p-1}\text{d}\lambda \mathrm{.}$

that is

$F\left(u+k\left(1\right)\beta \left(v-u\right)\right)-F\left(u\right)\ge {{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)+\mu {\Vert \beta \left(v-u\right)\Vert}^{p}{\displaystyle {\int}_{0}^{1}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}k{\left(\lambda \right)}^{p-1}\text{d}\lambda \mathrm{.}$

By using Condition A, we have

$F\left(v\right)-F\left(u\right)\ge {{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)+\mu {\Vert \beta \left(v-u\right)\Vert}^{p}{\displaystyle {\int}_{0}^{1}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}k{\left(\lambda \right)}^{p-1}\text{d}\lambda \mathrm{.}$

the required result.

We now give a necessary condition for higher order strongly pseudo-biconvex function.

Theorem 3.4. Let
${{F}^{\prime}}_{k\beta}\left(\mathrm{.}\right)$ be a higher order strongly relaxed *k*-pseudomonotone operator and Conditions *A* and *M* hold. Then

$F\left(v\right)-F\left(u\right)\ge \alpha {\Vert \eta \left(v\mathrm{,}u\right)\Vert}^{p}{\displaystyle {\int}_{0}^{1}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}k\left(\lambda \right)\text{d}\lambda \mathrm{,}$ (3.11)

*Proof*. Let
${F}^{\prime}$ be higher order relaxed *k*-pseudomonotone. Then,
$\forall u\mathrm{,}v\in {K}_{k\eta}$,

${{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)\ge \mathrm{0,}$

implies that

$-{{F}^{\prime}}_{k\beta}\left(v\mathrm{,}\beta \left(u-v\right)\right)\ge \alpha {\Vert \beta \left(u-v\right)\Vert}^{p}\mathrm{.}$ (3.12)

Since ${K}_{k\beta}$ is an biconvex set, $\forall u\mathrm{,}v\in {K}_{k\beta}$, $\lambda \in \left[\mathrm{0,1}\right]$, ${v}_{\lambda}=u+k\left(\lambda \right)\beta \left(v-u\right)\in {K}_{k\beta}$.

Taking
$v={v}_{\lambda}$ in (3.12) and using condition Condition *M*, we have

$-{{F}^{\prime}}_{k\beta}\left(u+\lambda \beta \left(v-u\right)\mathrm{,}\beta \left(u-v\right)\right)\ge k\left(\lambda \right)\alpha {\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{.}$ (3.13)

Let

$\xi \left(\lambda \right)=F\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\mathrm{,}\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k\beta},\lambda \in \left[\mathrm{0,1}\right].$

Then, using (3.13), we have

${\xi}^{\prime}\left(\lambda \right)={{F}^{\prime}}_{k\beta}\left(u+k\left(\lambda \right)\beta \left(v-u\right)\mathrm{,}\beta \left(u-v\right)\right)\ge k\left(\lambda \right)\alpha {\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{.}$

Integrating the above relation between 0 to 1, we have

$\xi \left(1\right)-\xi \left(0\right)\ge \alpha {\Vert \beta \left(v-u\right)\Vert}^{p}{\displaystyle {\int}_{0}^{1}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}k\left(\lambda \right)\text{d}\lambda \mathrm{,}$

that is,

$F\left(u+k\left(1\right)\beta \left(v-u\right)\right)-F\left(u\right)\ge \alpha {\Vert \beta \left(v-u\right)\Vert}^{p}{\displaystyle {\int}_{0}^{1}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}k\left(\lambda \right)\text{d}\lambda \mathrm{,}$

which implies, using Condition A,

$F\left(v\right)-F\left(u\right)\ge \alpha {\Vert \beta \left(v-u\right)\Vert}^{p}{\displaystyle {\int}_{0}^{1}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}k\left(\lambda \right)\text{d}\lambda \mathrm{,}$

which is the required (3.11).

Theorem 3.5. Let the differential F’(*u*) of a differentiable biconvex function
$F\left(u\right)$ be Lipschitz continuous on the biconvex set
${K}_{k\beta}$ with a constant
$\beta \text{\hspace{0.17em}}>\text{\hspace{0.17em}}0$. Then,

$F\left(u+\beta \left(v-u\right)\right)-F\left(u\right)\le {{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\eta \left(v\mathrm{,}u\right)\right)+\beta {\Vert \beta \left(v-u\right)\Vert}^{p}{\displaystyle {\int}_{0}^{1}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}k\left(\lambda \right)\text{d}\lambda \mathrm{,}\text{\hspace{1em}}u\mathrm{,}v\in {K}_{k\beta}.$

Proof. Its proof follow from Noor and Noor [22].

Definition 3.1. The function *F* is said to be sharply higher order strongly pseudo biconvex, if there exists a constant
$\mu >0$ such that

${{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)\ge 0$

$\Rightarrow $

$\begin{array}{l}F\left(v\right)\ge F\left(v+k\left(\lambda \right)\beta \left(v-u\right)\right)\text{\hspace{0.05em}}+\mu \left\{k{\left(\lambda \right)}^{p}\left(1-k\left(\lambda \right)\right)+\lambda {\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{,}\\ \forall u\mathrm{,}v\in {K}_{k\beta},\lambda \in \left[\mathrm{0,1}\right]\mathrm{.}\end{array}$

Theorem 3.6. Let *F* be a sharply higher order strongly pseudo biconvex function on
${K}_{k\beta}$ with a constant
$\mu >0$. Then

$-{{F}^{\prime}}_{k\beta}\left(v\mathrm{,}\beta \left(v-u\right)\right)\ge \mu {\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{,}\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k\beta}.$

*Proof*. Let *F* be a sharply higher strongly pseudo generalized biconvex function on
${K}_{\eta}$. Then

$\begin{array}{l}F\left(v\right)\ge F\left(v+k\left(\lambda \right)\beta \left(v-u\right)\right)\\ \text{\hspace{0.05em}}+\mu \left\{k{\left(\lambda \right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{,}\\ \forall u\mathrm{,}v\in {K}_{k\beta},\lambda \in \left[\mathrm{0,1}\right]\mathrm{,}\end{array}$

from which we have

$\begin{array}{l}\frac{F\left(v+k\left(\lambda \right)\beta \left(v-u\right)\right)-F\left(v\right)}{k\left(\lambda \right)}\\ \text{\hspace{0.05em}}+\mu \left\{k{\left(\lambda \right)}^{p-1}\left(1-k\left(\lambda \right)\right)+{\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\le 0.\end{array}$

Taking limit in the above-mentioned inequality, as $\lambda \to 0$, we have

$-{{F}^{\prime}}_{k\beta}\left(v\mathrm{,}\beta \left(v-u\right)\right)\ge \mu {\Vert \eta \left(v\mathrm{,}u\right)\Vert}^{p}\mathrm{,}\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{,}$

the required result.

Definition 3.2.* *A function *F* is said to be a higher order strongly pseudo biconvex function with respect to strictly positive bifunction
$W\left(\mathrm{.,.}\right)$, if

$F\left(v\right)<F(\; u\; )$

$\Rightarrow $

$F\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)<F\left(u\right)+k\left(\lambda \right)\left(k\left(\lambda \right)-1\right)W\left(v,u\right),\forall u,v\in {K}_{k\beta},\lambda \in \left[0,1\right].$

Theorem 3.7.* *If the function *F* is higher order strongly biconvex function such that
$F\left(v\right)<F\left(u\right)$, then the function *F* is higher order strongly pseudo biconvex.

Proof. Since
$F\left(v\right)<F\left(u\right)$ and *F* is higher order strongly biconvex function, then
$\forall u\mathrm{,}v\in {K}_{k\eta}\mathrm{,}\lambda \in \left[\mathrm{0,1}\right]$, we have

$\begin{array}{l}F\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\\ \le F\left(u\right)+k\left(\lambda \right)\left(F\left(v\right)-F\left(u\right)\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\mu \left\{k{\left(\lambda \right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\\ <F\left(u\right)+k\left(\lambda \right)\left(1-k\left(\lambda \right)\right)\left(F\left(v\right)-F\left(u\right)\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\mu \left\{k{\left(\lambda \right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\end{array}$

$\begin{array}{l}=F\left(u\right)+k\left(\lambda \right)\left(k\left(\lambda \right)-1\right)\left(F\left(u\right)-F\left(v\right)\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\mu \left\{k{\left(\lambda \right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\\ <F\left(u\right)+k\left(\lambda \right)\left(k\left(\lambda \right)-1\right)W\left(u,v\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\mu \left\{k{\left(\lambda \right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p},\forall u,v\in {K}_{k\beta}.\end{array}$

where
$W\left(u,v\right)=F\left(u\right)-F\left(v\right)>0$. This shows that the function *F* is higher order strongly pseudo biconvex.

4. Applications

In this section, we derive new parallelogram laws for uniformly Banach spaces as a novel application of higher order strongly biconvex functions.

From definition 2.13, we have

$\begin{array}{l}F\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)=\left(1-k\left(\lambda \right)\right)F\left(u\right)+k\left(\lambda \right)F\left(v\right)\\ \text{\hspace{0.05em}}-\mu \left\{k{\left(\lambda \right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{,}\\ \forall u\mathrm{,}v\in {K}_{k\beta},\lambda \in \left[\mathrm{0,1}\right]\mathrm{,}p\ge 1.\end{array}$ (4.1)

Taking $\lambda =\frac{1}{2}$ in (4.1), we have

$\begin{array}{l}F\left(u+k\left(\frac{1}{2}\right)\beta \left(v-u\right)\right)\\ \text{\hspace{0.05em}}+\mu \left\{{\left(k\left(\frac{1}{2}\right)\right)}^{p}\left(1-k\left(\frac{1}{2}\right)\right)+k\left(\frac{1}{2}\right){\left(1-k\left(\frac{1}{2}\right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\\ =k\left(\frac{1}{2}\right)F\left(u\right)+k\left(\frac{1}{2}\right)F\left(v\right)\mathrm{,}\forall u\mathrm{,}v\in {K}_{\eta}\mathrm{.}\end{array}$ (4.2)

which is known as the generalized parallelogram-like laws for the Banach spaces involving higher order generalized biconvex functions.

We now discuss some special cases of the generalized parallelogram-like laws.

1) If $k\left(\frac{1}{2}\right)=\frac{1}{2}$, $F\left(u\right)={\Vert u\Vert}^{p},p>1$, then, Equation (4.2) becomes

${\Vert u+\frac{1}{2}\beta \left(v-u\right)\Vert}^{p}+\mu \frac{1}{{2}^{p}}{\Vert \beta \left(v-u\right)\Vert}^{p}=\frac{1}{2}{\Vert u\Vert}^{p}+\frac{1}{2}{\Vert v\Vert}^{p},\forall u,v\in {K}_{k\beta}.$ (4.3)

2) If $\beta \left(v,u\right)=v-u$, $k\left(\frac{1}{2}\right)=\frac{1}{2}$, then (4.3) reduces to the parallelogram-like law as:

${\Vert v+u\Vert}^{p}+\mu {\Vert v-u\Vert}^{p}={2}^{p-1}\left\{{\Vert u\Vert}^{p}+{\Vert v\Vert}^{p}\right\},$ (4.4)

which is known as the parallelogram-like law for the uniform Banach spaces involving the biconvex functions. Xu [36] obtained these characterizations of *p*- uniform convexity and *q*-uniform smoothness of a Banach space via the functionals
${\Vert \text{\hspace{0.05em}}\mathrm{.}\text{\hspace{0.05em}}\Vert}^{p}$ and
${\Vert \text{\hspace{0.05em}}\mathrm{.}\text{\hspace{0.05em}}\Vert}^{q}$, respectively. Bynum [37] and Chen *et al*. [38] [39] have studied the properties and applications of the parallelogram laws for the Banach spaces. For the applications of the parallelogram laws in Banach spaces in prediction theory and applied sciences, see [36] [37] [38] [39] and the references therein.

3) If $p=2$, $K\left(\frac{1}{2}\right)=\frac{1}{2}$, then parallelogram law (4.2) reduces to:

${\Vert 2u+\beta \left(v-u\right)\Vert}^{2}+\mu {\Vert \beta \left(v-u\right)\Vert}^{2}=2\left\{{\Vert u\Vert}^{2}+{\Vert v\Vert}^{2}\right\}\mathrm{.}$ (4.5)

This is a new parallelogram law, which characterizes the inner product spaces involving bifunction $\beta \left(\mathrm{.}-\mathrm{.}\right)$.

4) If $\beta \left(v,u\right)=v-u$, $\mu =1$, then parallelogram law (37) becomes

${\Vert v+u\Vert}^{2}+{\Vert v-u\Vert}^{2}=2\left\{{\Vert u\Vert}^{2}+{\Vert v\Vert}^{2}\right\}\mathrm{,}$ (4.6)

which is a well known parallelogram and characterizes the inner product spaces.

Remark 4.1. In this section, we have derived some new parallelogram laws as important applications of the higher order strongly biconvex functions. For suitable and appropriate choice of the function* k*,* *bifunction*
$\beta \left(\mathrm{.}-\mathrm{.}\right)$ *and parameter* p*,* *one can obtain a wide class of parallelogram laws, which can be used to characterize various classes of inner products. Applications of the new parallelogram laws in different areas are an open problem for future research.

5. Biequilibrium Problems

In this section, we introduce the biequilibrium problems. We also use the auxiliary principle techniques coupled with Bregman functions to consider the convergence criteria of the proposed iterative methods. This is the main motivation of this section.

For the readers convenience, we recall some basic properties of the Bregman [27] distance functions on the convex set *K*. For strongly convex function *F*, we define the Bregman distance function as:

$B\left(v\mathrm{,}u\right)=F\left(v\right)-F\left(u\right)-\langle {F}^{\prime}\left(u\right)\mathrm{,}v-u\rangle \ge \alpha {\Vert v-u\Vert}^{2}\mathrm{,}\forall u\mathrm{,}v\in K\mathrm{.}$ (5.1)

It is important to emphasize that various types of function *F* gives different Bregman distance. We give the following important examples of some practical important types of function *F* and their corresponding Bregman distance, see [2] [3].

Examples

1) If
$f\left(v\right)={\Vert v\Vert}^{2}$, then
$B\left(v\mathrm{,}u\right)=\Vert v-u\Vert $, which is the squared Euclidean distance (*SE*).

2) If $f\left(v\right)={\displaystyle {\sum}_{i=1}^{n}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{a}_{i}\mathrm{log}\left({v}_{i}\right)$, which is known as Shannon entropy, then its corresponding Bregman distance is given as:

$B\left(v,u\right)={\displaystyle \underset{i=1}{\overset{n}{\sum}}\left({v}_{i}\mathrm{log}\left(\frac{{v}_{i}}{{u}_{i}}\right)+{u}_{i}-{v}_{i}\right)}$.

This distance is called KullbackLeibler distance (*KL*) and as become a very important tool in several areas of applied mathematics such as machine learning.

3) If $f\left(v\right)=-{\displaystyle {\sum}_{i=1}^{n}}\mathrm{log}\left({v}_{i}\right)$, which is called Burg entropy, then its corresponding Bregman distance is given as:

$B\left(v,u\right)={\displaystyle \underset{i=1}{\overset{n}{\sum}}}\left(\mathrm{log}\left(\frac{{v}_{i}}{{u}_{i}}\right)+\frac{{v}_{i}}{{u}_{i}}-1\right).$

This is called ItakuraSaito distance (*IS*), which is very important in information theory.

It is a challenging problem to explore the applications of Bregman distance for other types of nonconvex functions as biconvex, *k*-convex functions and harmonic functions.

We now discuss the optimality for the differentiable higher order generalized biconvex functions.

Theorem 5.1.* *Let *F* be a differentiable higher order strongly biconvex function with modulus
$\mu >0$. If
$u\in {K}_{k\beta}$ is the minimum of the function *F*, then

$F\left(v\right)-F\left(u\right)\ge \mu {\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{,}\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k\beta}.$ (5.2)

*Proof*. Let
$u\in {K}_{k\beta}$ be a minimum of the function *F*. Then,

$F\left(u\right)\le F\left(v\right)\mathrm{,}\forall v\in {K}_{k\beta}.$ (5.3)

Since ${K}_{k\beta}$ is an biconvex set, so, $\forall u\mathrm{,}v\in {K}_{k\beta}$, $\lambda \in \left[\mathrm{0,1}\right]$, ${v}_{\lambda}=u+k\left(\lambda \right)\beta \left(v-u\right)\in {K}_{k\beta}$.

Taking $v={v}_{\lambda}$ in (5.3), we have

$0\le \underset{\lambda \to 0}{\mathrm{lim}}\left\{\frac{F\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)-F\left(u\right)}{k\left(\lambda \right)}\right\}={{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)\mathrm{.}$ (5.4)

Since *F* is differentiable higher order strongly biconvex function, so

$\begin{array}{l}F\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)\le F\left(u\right)+k\left(\lambda \right)\left(F\left(v\right)-F\left(u\right)\right)\\ \text{\hspace{0.05em}}-\mu \left\{k{\left(\lambda \right)}^{p}\left(1-k\left(\lambda \right)\right)+k\left(\lambda \right){\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{,}\forall u\mathrm{,}v\in {K}_{k\beta}.\end{array}$

from which, using (5.4), we have

$\begin{array}{c}F\left(v\right)-F\left(u\right)\ge \underset{\lambda \to 0}{\mathrm{lim}}\left\{\frac{F\left(u+k\left(\lambda \right)\beta \left(v-u\right)\right)-F\left(u\right)}{k\left(\lambda \right)}\right\}\\ \text{\hspace{0.05em}}\text{\hspace{0.17em}}+\mu \left\{k{\left(\lambda \right)}^{p-1}\left(1-k\left(\lambda \right)\right)+{\left(1-k\left(\lambda \right)\right)}^{p}\right\}{\Vert \beta \left(v-u\right)\Vert}^{p}\\ ={{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)+\mu {\Vert \beta \left(v-u\right)\Vert}^{p}\mathrm{,}\end{array}$

the required result (5.2).

We would like to mention that, if $u\in {K}_{k\beta}$ satisfies the inequality,

${{F}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(v-u\right)\right)+\mu {\Vert \beta \left(v-u\right)\Vert}^{p}\ge \mathrm{0,}\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{,}$ (5.5)

then
$u\in {K}_{k\beta}$ is the minimum of the function *F*. The inequality of the type (5.5) is called the higher order strongly biequilibrium problem and appears to new one.

It is worth mentioning that inequalities of the type (5.5) may not arise as the minimization of the higher order strongly biconvex functions. This motivated us to consider a more general higher order strongly biequilibrium problem of which (5.5) is a special case.

To be more precise, for given operator *D*, bifunction
$\beta \left(\mathrm{.}-\mathrm{.}\right)$ and a constant
$\mu >0$, consider the problem of finding
$u\in {K}_{k\beta}$, such that

$D\left(u\mathrm{,}\beta \left(v-u\right)\right)+\mu {\Vert \beta \left(v-u\right)\Vert}^{p}\ge \mathrm{0,}\forall v\in {K}_{k\beta}\mathrm{,}p\ge \mathrm{1,}$ (5.6)

which is called the higher order strongly biequilibrium problem.

1) If $\mu =0$, then the higher order strongly biequilibrium problem (2.1) reduces to finding $u\in {K}_{k\beta}$ such that

$D\left(u\mathrm{,}\beta \left(v-u\right)\right)\ge \mathrm{0,}\text{\hspace{1em}}\forall v\in {K}_{k\beta}\mathrm{,}$

which is called the biequilibrium problem. For the formulation, numerical methods and other aspects of the bivariational inequalities and biequilibrium problem, see [28] [30].

2) If f
$D\left(u\mathrm{,}\beta \left(v-u\right)\right)=\langle Au\mathrm{,}\beta \left(v-u\right)\rangle $, where *A* is a nonlinear operator, then the higher order strongly biequilibrium problem (2.1) is equivalent to finding
$u\in {K}_{k\beta}$ such that

$\langle Au\mathrm{,}\beta \left(v-u\right)\rangle +\mu {\Vert \beta \left(v-u\right)\Vert}^{p}\ge \mathrm{0,}\text{\hspace{1em}}\forall v\in {K}_{k\eta}\mathrm{,}$

which is the higher order strongly bivariational inequality and appears to be a new one.

It is worth mentioning that for suitable and appropriate choice of the operators, $k\beta $ -biconvex sets and spaces, one can obtain a wide class of equilibrium problem, variational-like inequalities and optimization programming. This shows that the higher order strongly biequilibrium problems are quite flexible and unified ones.

Due to the inherent nonlinearity, the projection method and its variant form can not be used to suggest the iterative methods for solving these higher order strongly biequilibrium problems. To overcome these drawback, one may use the auxiliary principle technique of Glowinski *et al*. [11] as developed by Noor [26] [28] [30] Noor *et al*. [14] [15] [16] to suggest and analyze some iterative methods for solving the higher order strongly biequilibrium problems (5.6). This technique does not involve the concept of the projection, which is the main advantage of this technique.

We again use the auxiliary principle technique coupled with Bergman distance functions. These applications are based on the type of convex functions associated with the Bregman distance function. We now suggest and analyze some iterative methods for higher order strongly biequilibrium problems (5.6) using the auxiliary principle technique coupled with Bregman functions as developed by Noor [28] [30].

For a given $u\in {K}_{k\beta}$ satisfying the biequilibrium problem (5.6), we consider the auxiliary problem of finding a $w\in {K}_{k\beta}$ such that

$\begin{array}{l}\rho D\left(w\mathrm{,}\beta \left(v-w\right)\right)+\langle {{E}^{\prime}}_{k\beta}\left(w\right)-{{E}^{\prime}}_{k\beta}\left(u\right)\mathrm{,}\beta \left(v-w\right)\rangle \\ \text{\hspace{0.05em}}+\rho {\Vert \beta \left(v-w\right)\Vert}^{p}\ge \mathrm{0,}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall v\in {K}_{k\beta}\mathrm{,}\end{array}$ (5.7)

where $\rho >0$ is a constant and ${{E}^{\prime}}_{k\beta}\left(u\right)$ is the differential of a strongly biconvex function $E\left(u\right)$ at $u\in K$.

Remark 3.1: The function $B\left(w\mathrm{,}u\right)=E\left(w\right)-E\left(u\right)-{{E}^{\prime}}_{k\beta}\left(u\mathrm{,}\beta \left(w-u\right)\right)$ associated with the biconvex function $E\left(u\right)$ is called the generalized Bregman function. By the strongly biconvexity of the function $E\left(u\right)$, the Bregman function $B\left(\mathrm{.,.}\right)$ is nonnegative and $B\left(w\mathrm{,}u\right)=0$, if and only if $u=w$, $\forall u,w\in {K}_{k\beta}$. For the applications of the Bregman function in solving variational inequalities and complementarity problems, see [28] [30].

We note that, if
$w=u$, then clearly *w* is solution of the higher order strongly biequilibrium problem (5.6). This observation enables us to suggest and analyze the following iterative method for solving (2.1).

Algorithm 3.1. For a given ${u}_{0}\in H$, compute the approximate solution ${u}_{n+1}$ by the iterative scheme:

$\begin{array}{l}\rho D\left({u}_{n+1}\mathrm{,}\beta \left(v-{u}_{n+1}\right)\right)+\langle {{E}^{\prime}}_{k\beta}\left({u}_{n+1}\right)-{{E}^{\prime}}_{k\beta}\left({u}_{n}\mathrm{,}\beta \left(v-{u}_{n+1}\right)\right)\rangle \\ \text{\hspace{0.05em}}+\rho {\Vert \beta \left(v-u\right)\Vert}^{p}\ge \mathrm{0,}\text{\hspace{1em}}\forall v\in {K}_{k\beta}\mathrm{,}\end{array}$ (5.8)

where $\rho >0$ is a constant. Algorithm 3.1 is called the proximal method for solving higher order strongly biequilibrium-like problem (5.6). In passing we remark that the proximal point method was suggested in the context of convex programming problems as a regularization technique.

If $\mu =0$, then Algorithm 3.1 collapses to:

Algorithm 3.2. For a given ${u}_{0}\in H$, compute the approximate solution ${u}_{n+1}$ by the iterative scheme:

$\rho D\left({u}_{n+1}\mathrm{,}\beta \left(v-{u}_{n+1}\right)\right)+\langle {{E}^{\prime}}_{k\beta}\left({u}_{n+1}\right)-{{E}^{\prime}}_{k\beta}\left({u}_{n}\mathrm{,}\beta \left(v-{u}_{n+1}\right)\right)\rangle \ge \mathrm{0,}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\forall v\in {K}_{k\beta}\mathrm{,}$

for solving the bivariational inequality.

For suitable and appropriate choice of the operators and the spaces, one can obtain a number of known and new algorithms for solving variational inequalities and related problems.

Definition 2.2. The bifunction $D\left(\mathrm{.,.}\right)$ is said to be:

1)* *monotone, iff,

$D\left(u\mathrm{,}\beta \left(v-u\right)\right)+D\left(v\mathrm{,}\beta \left(u-v\right)\right)\le \mathrm{0,}\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{.}$

2)* *pseudomonotone* *with respect to
$\mu {\Vert \beta \left(v-u\right)\Vert}^{p}$, iff,

$D\left(u\mathrm{,}\beta \left(v-u\right)\right)+\mu {\Vert \beta \left(v-u\right)\Vert}^{p}\ge 0$

$\Rightarrow $

$-D\left(v\mathrm{,}\beta \left(u-v\right)\right)-\mu {\Vert \beta \left(v-u\right)\Vert}^{p}\ge \mathrm{0,}\text{\hspace{1em}}\forall u\mathrm{,}v\in {K}_{k\beta}\mathrm{.}$

3)* *higher order strongly partially relaxed monotone, if there exists a constant
$\gamma >0$ such that

$D\left(u\mathrm{,}\beta \left(u-v\right)\right)+D\left(v\mathrm{,}\beta \left(z-v\right)\right)\le \gamma {\Vert \beta \left(z-u\right)\Vert}^{p}\mathrm{,}\text{\hspace{1em}}\forall u\mathrm{,}v\mathrm{,}z\in {K}_{k\beta}\mathrm{.}$

Note that for $z=u$ higher order strongly partially relaxed monotonicity reduces to monotonicity. This shows that higher order strongly partially relaxed monotonicity implies monotonicity, but the converse is not true.

Theorem 5.2. Let the bifunction
$D\left(\mathrm{.,.}\right)$ be pseudomonotone with respect to
$\mu {\Vert \beta \left(v-u\right)\Vert}^{p}$. If *E* be differentiable higher order strongly biconvex function with module
$\beta >0$, then the approximate solution
${u}_{n+1}$ obtained from Algorithm 3.1 converges to a solution
$u\in K$ satisfying the biequilibrium problem (5.6).

Proof. Let $u\in {K}_{k\beta}$ be a solution of (5.6). Then,

$D\left(u\mathrm{,}\beta \left(v-u\right)\right)+\mu {\Vert \beta \left(v-u\right)\Vert}^{p}\ge \mathrm{0,}\text{\hspace{1em}}\forall v\in {K}_{k\beta}\mathrm{,}$

implies that

$-D\left(v\mathrm{,}\beta \left(u-v\right)\right)-\mu {\Vert \beta \left(v-u\right)\Vert}^{p}\ge \mathrm{0,}\text{\hspace{1em}}\forall v\in {K}_{k\beta}\mathrm{,}$ (5.9)

since $D\left(\mathrm{.,.}\right)$ is pseudomonotone with respect to $\mu {\Vert \beta \left(v-u\right)\Vert}^{p}$.

Taking $v=u$ in (5.8) and $v={u}_{n+1}$ in (5.9), we have

$\begin{array}{l}\rho D\left({u}_{n+1}\mathrm{,}\beta \left(u-{u}_{n+1}\right)\right)+\langle {{E}^{\prime}}_{k\beta}\left({u}_{n+1}\right)-{{E}^{\prime}}_{k\beta}\left({u}_{n}\mathrm{,}\beta \left(u-{u}_{n+1}\right)\right)\rangle \\ \ge -\rho \mu {\Vert \beta \left(u-{u}_{n+1}\right)\Vert}^{p}\mathrm{.}\end{array}$ (5.10)

And,

$-D\left({u}_{n+1}\mathrm{,}\beta \left(u-{u}_{n+1}\right)\right)-\mu {\Vert \beta \left(u-{u}_{n+1}\right)\Vert}^{p}\ge 0.$ (5.11)

We now consider the Bregman function

$B\left(u\mathrm{,}w\right)=E\left(u\right)-E\left(w\right)-\langle {{E}^{\prime}}_{k\beta}\left(w\mathrm{,}\beta \left(u-w\right)\right)\rangle \ge \mu {\Vert \beta \left(u-w\right)\Vert}^{p}\mathrm{,}$ (5.12)

using higher order strongly biconvexity of *E*.

Now combining (5.10), (5.11) and (5.12), we have

$\begin{array}{l}B\left(u,{u}_{n}\right)-B\left(u,{u}_{n+1}\right)\\ =E\left({u}_{n+1}\right)-E\left({u}_{n}\right)-\langle {{E}^{\prime}}_{k\beta}\left({u}_{n},\beta \left(u-{u}_{n}\right)\right)\rangle +\langle {{E}^{\prime}}_{k\beta}\left({u}_{n+1},\beta \left(u-{u}_{n+1}\right)\right)\rangle \\ =E\left({u}_{n+1}\right)-E\left({u}_{n}\right)-\langle {{E}^{\prime}}_{k\beta}\left({u}_{n}\right)-{{E}^{\prime}}_{k\beta}\left({u}_{n+1},\beta \left(u-{u}_{n+1}\right)\right)\rangle -{{E}^{\prime}}_{k\beta}\left({u}_{n},\beta \left({u}_{n+1}-{u}_{n}\right)\right)\\ \ge \mu {\Vert \beta \left({u}_{n+1}-{u}_{n}\right)\Vert}^{p}+\langle {{E}^{\prime}}_{k\beta}\left({u}_{n+1}\right)-{{E}^{\prime}}_{k\beta}\left({u}_{n},\beta \left(u-{u}_{n+1}\right)\right)\rangle \\ \ge \mu {\Vert \beta \left({u}_{n+1}-{u}_{n}\right)\Vert}^{p}-\rho D\left({u}_{n+1},\beta \left(u-{u}_{n+1}\right)\right)-\rho \mu {\Vert \beta \left(u-{u}_{n+1}\right)\Vert}^{p}\\ \ge \mu {\Vert \beta \left({u}_{n+1}-{u}_{n}\right)\Vert}^{p}.\end{array}$

If ${u}_{n+1}={u}_{n}$, then clearly ${u}_{n}$ is a solution of the problem (5.6). Otherwise, it follows that $B\left(u\mathrm{,}{u}_{n}\right)-B\left(u\mathrm{,}{u}_{n+1}\right)$ is nonnegative and we must have

$\underset{n\to \infty}{\mathrm{lim}}\left(\beta \left({u}_{n+1}-{u}_{n}\right)\right)=0.$

from which, we have

$\underset{n\to \infty}{\mathrm{lim}}\Vert {u}_{n+1}-{u}_{n}\Vert =0.$

It follows that the sequence $\left\{{u}_{n}\right\}$ is bounded. Let $\stackrel{\xaf}{u}$ be a cluster point of the subsequence $\left\{{u}_{{n}_{i}}\right\}$, and let $\left\{{u}_{{n}_{i}}\right\}$ be a subsequence converging toward $\stackrel{\xaf}{u}$. Now using the technique of Zhu and Marcotte [17], it can be shown that the entire sequence $\left\{{u}_{n}\right\}$ converges to the cluster point $\stackrel{\xaf}{u}$ satisfying the biequilibrium problem (5.6).

It is well-known that to implement the proximal point methods, one has to find the approximate solution implicitly, which is itself a difficult problem. To overcome this drawback, we now consider another method for solving the higher order strongly biequilibrium problem (5.6) using the auxiliary principle technique.

For a given $u\in {K}_{k\beta}$, find $w\in {K}_{k\beta}$, such that

$\begin{array}{l}\rho D\left(u\mathrm{,}\beta \left(v-w\right)\right)+\langle {{E}^{\prime}}_{k\beta}\left(w\right)-{{E}^{\prime}}_{k\beta}\left(u\right)\mathrm{,}\beta \left(v-w\right)\rangle \\ \text{\hspace{0.05em}}+\rho \mu {\Vert \beta \left(v-w\right)\Vert}^{p}\ge \mathrm{0,}\text{\hspace{1em}}\forall v\in {K}_{k\beta}\mathrm{,}\end{array}$ (5.13)

where
${{E}^{\prime}}_{k}\left(u\right)$ is the differential of a strongly *k*-biconvex function
$E\left(u\right)$ at
$u\in {K}_{k\beta}$. Note that problems (5.13) and (5.8) are quite different problems. It is clear that for
$w=u$, *w* is a solution of (5.6). This fact allows us to suggest and analyze another iterative method for solving the higher order strongly biequilibrium problem (5.6).

Algorithm 3.3. For a given ${u}_{0}\in H$, compute the approximate solution ${u}_{n+1}$ by the iterative scheme:

$\begin{array}{l}\rho D\left({u}_{n}\mathrm{,}\beta \left(v-{u}_{n+1}\right)\right)+\langle {{E}^{\prime}}_{k\beta}\left({u}_{n+1}\right)-{{E}^{\prime}}_{k\beta}\left({u}_{n}\right)\mathrm{,}\beta \left(v-{u}_{n+1}\right)\rangle \\ \ge -\rho \mu {\Vert \beta \left(v-{u}_{n+1}\right)\Vert}^{p}\mathrm{,}\text{\hspace{1em}}\forall v\in {K}_{k\beta}\mathrm{,}\end{array}$ (5.14)

for solving the higher order strongly biequilibrium problem (5.6).

If $D\left(u\mathrm{,}\beta \left(v\mathrm{,}u\right)\right)=\langle Au\mathrm{,}\beta \left(v-u\right)\rangle $, Algorithm 3.3 collapses to:

Algorithm 3.4. For a given ${u}_{0}\in H$, compute the approximate solution ${u}_{n+1}$ by the iterative schemes:

$\begin{array}{l}\rho \langle A{u}_{n}\mathrm{,}\beta \left(v-{u}_{n+1}\right)\rangle +\langle {{E}^{\prime}}_{k\beta}\left({u}_{n+1}\right)-{{E}^{\prime}}_{k\beta}\left({u}_{n}\right)\mathrm{,}\beta \left(v-{u}_{n+1}\right)\rangle \\ \ge -\rho \mu {\Vert \beta \left(v-{u}_{n+1}\right)\Vert}^{p}\mathrm{,}\text{\hspace{1em}}\forall v\in {K}_{k\beta}\mathrm{,}\end{array}$

for solving the higher order strongly bivariational inequalities and appears to be a new one.

Discussion of Results: We have introduced and investigated higher order strongly biequilibrium problems. The auxiliary principle technique coupled with Bregman distance functions is used to suggest some iterative methods for solving biequilibrium problems, which is the main contributions of this section. We also mention that for suitable and appropriate choice of the operators and the spaces, one can obtain various known and new algorithms for solving higher order strongly equilibrium problem (5.6) and related optimization problems. It is an interesting problem from both analytically and numerically point of view.

6. Conclusion

In this paper, we have introduced and studied some new classes of biconvex functions, which are called higher order strongly biconvex functions. These concepts are more general and unifying than the previous ones. Several new properties of these higher order strongly biconvex functions are discussed and their relations with previously known results are highlighted. Several parallelogram laws for inner product spaces and Banach spaces are obtained as applications of the higher order strongly *k*-convex functions, which is its can viewed as novel applications of higher order strongly biconvex functions. It is shown that the optimality conditions of the differentiable *k*-biconvex functions can be characterized by a class of higher order strongly biequilibrium problems. This result is used to introduce some new classes of higher order strongly biequlibrium problems. The auxiliary principle technique involving the Bregman functions are used to suggest and investigate some iterative methods for solving higher order strongly biequilibrium problem. These techniques include the entropy method and its variant forms. Consequently, we can conclude that the entropy methods have played an important in the analysis of the convergence criteria of the proposed methods implicitly. It is itself an interesting problem to develop some efficient numerical methods for solving higher order strongly biequilibrium problems along with their applications in pure and applied sciences. Despite the current activity, much clearly remains to be done in these fields. It is expected that the ideas and techniques of this paper may be the starting point for future research activities.

Acknowledgements

We wish to express our deepest gratitude to our colleagues, collaborators, and friends, who have direct or indirect contributions in the process of this paper. We are also like to thank the Rector, COMSATS University Islamabad, Pakistan for the research facilities and support in our research endeavors.

References

[1] Stampacchia, G. (1964) Formes bilineaires coercivites sur les ensembles convexes. Comptes rendus de l’Académie des Sciences Paris, 258, 4413-4416.

[2] Jolaoso, L.O., Aphane, M. and Khan, S.H. (2020) Two Bregman Projection Methods for Solving Variational Inequality Problems in Hilbert Spaces with Applications to Signal Processing. Symmetry, 12, 2007.

https://doi.org/10.3390/sym12122007

[3] Sunthrayuth, P. and Cholamjiak, P. (2020) Modified Extragradient Method with Bregman Distance for Variational Inequalities. Applicable Analysis.

[4] Glowinski, R., Lions, J.L. and Tremolieres, R. (1981) Numerical Analysis of Variational Inequalities. North-Holland, Amsterdam.

[5] Kinderlehrer, D. and Stampaccia, G. (1980) Introduction to Variational Inequalties and Applications. Academic Press, New York.

[6] Lin, G.H. and Fukushima, M. (2003) Some Exact Penalty Results for Nonlinear Programs and Mathematical Programs with Equilibrium Constraints. Journal of Optimization Theory and Applications, 118, 67-80.

https://doi.org/10.1023/A:1024787424532

[7] Lions, J.L. and Stampacchia, G. (1967) Variational Inequalities. Communications on Pure and Applied Mathematics, 20, 493-519.

https://doi.org/10.1002/cpa.3160200302

[8] Noor, M.A. (1988) General Variational Inequalities. Applied Mathematics Letters, 1, 119-121.

https://doi.org/10.1016/0893-9659(88)90054-7

[9] Noor, M.A. (2000) New Approximation Schemes for General Variational Inequalities. Journal of Mathematical Analysis and Applications, 251, 217-229.

https://doi.org/10.1006/jmaa.2000.7042

[10] Noor, M.A. (2006) Fundamental of Equilibrium Problems. Mathematical Inequalities & Applications, 9, 520-566.

https://doi.org/10.7153/mia-09-51

[11] Noor, M.A. (2004) Some Developments in General Variational Inequalities. Applied Mathematics and Computation, 152, 199-277.

https://doi.org/10.1016/S0096-3003(03)00558-7

[12] Noor, M.A. (1994) Variational-Like Inequalities. Optimization, 30, 323-330.

https://doi.org/10.1080/02331939408843995

[13] Noor, M.A. (2005) Invex Equilibrium Problems. Journal of Mathematical Analysis and Applications, 302, 463-475.

https://doi.org/10.1016/j.jmaa.2004.08.014

[14] Noor, M.A., Noor, K.I. and Al-Said, E. (2011) Auxiliary Principle Technique for Solving Bifunction Variational Inequalities. Journal of Optimization Theory and Applications, 149, 441-445.

https://doi.org/10.1007/s10957-010-9785-z

[15] Noor, M.A., Noor, K.I. and Rassias, M.Th. (2020) New Trends in General Variational Inequalities. Acta Applicandae Mathematicae, 170, 981-1064.

https://doi.org/10.1007/s10440-020-00366-2

[16] Noor, M.A., Noor, K.I. and Rassias, Th.M. (1993) Some Aspects of Variational Inequalities. Journal of Applied Mathematics and Computing, 47, 485-512.

https://doi.org/10.1016/0377-0427(93)90058-J

[17] Zhu, D.L. and Marcotte, P. (1966) Co-Coercvity and Its Role in the Convergence of Iterative Schemes for Solving Variational Inequalities. SIAM Journal on Optimization, 6, 714-726.

https://doi.org/10.1137/S1052623494250415

[18] Blum, E. and Oettli, W. (1994) From Optimization and Variational Inequalities to Equilibrium Problems. Mathematics Student, 63, 123-145.

[19] Noor, M.A. and Oettli, W. (1994) On General Nonlinear Complementarity Problems and Quasi-Equilibria. Le Matematiche (Catania), 49, 313-331.

[20] Noor, M.A. and Noor, K.I. (2021) Exponentially Biconvex Functions and Bivariational Inequalties. In: Hazarika, B., Acharjee, S. and Srivastava, H.M., Eds., Advances in Mathematical Analysis and Multidisciplinary Applications, CRC Press, Boca Raton, 23 p.

[21] Noor, M.A., Noor, K.I. and Rassias, M.T. (2021) Strongly Biconvex Functions and Bi-Variational Inequalities. In: Pardalos, P.M. and Rassias, Th.M., Eds., Mathematical Analysis, Optimization, Approximation and Applications, World Scientific Publishing Company, Singapore, 21 p.

[22] Noor, M.A. and Noor, K.I. (2021) Properties of Higher Order Preinvex Functions. Numerical Algebra, Control and Optimization, 11, 431-441.

https://doi.org/10.3934/naco.2020035

[23] Micherda, B. and Rajba, T. (2012) On Some Hermite-Hadamard-Fejer Inequalities for (k,h)-Convex Functions. Mathematical Inequalities & Applications, 12, 931-940.

https://doi.org/10.7153/mia-15-79

[24] Hazy, A. (2012) Bernstein-Doetsch Type Results for (k,h)-Convex Functions. Miskolc Mathematical Notes, 13, 325-336.

https://doi.org/10.18514/MMN.2012.538

[25] Crestescu, G., Gaianu, M. and Awan, M.U. (2015) Regularity Properties and Integral Inequalities Related to (k,h1,h2)-Convexity of Functions. Annals of West University of Timisoara—Mathematics and Computer Science, 53, 19-35.

https://doi.org/10.1515/awutm-2015-0002

[26] Noor, M.A. (2006) Some New Classes of Nonconvex Functions. Nonlinear Functional Analysis and Applications, 11, 165-171.

[27] Bregman, L.M. (1967) The Relaxation Method for Finding Common Points of Convex Sets and Its Application to the Solution of Problems in Convex Programming. USSR Computational Mathematics and Mathematical Physics, 7, 200-217.

https://doi.org/10.1016/0041-5553(67)90040-7

[28] Mohsen, B.B., Noor, M.A., Noor, K.I. and Postolache, M. (2019) Strongly Convex Functions of Higher Order Involving Bifunction. Mathematics, 7, 1028.

https://doi.org/10.3390/math7111028

[29] Alabdali, O., Guessab, A. and Schmeisser, G. (2019) Characterizations of Uniform Convexity for Differentiable Functions. Applicable Analysis and Discrete Mathematics, 13, 721-732.

https://doi.org/10.2298/AADM190322029A

[30] Olbrys, A. (2018) A Support Theorem for Generalized Convexity and Its Applications. Journal of Mathematical Analysis and Applications, 458, 1044-1058.

https://doi.org/10.1016/j.jmaa.2017.09.038

[31] Mako, J. and Pales, Z. (2012) On ø-Convexity. Publicationes Mathematicae Debrecen, 80, 107-126.

[32] Cristescu, G. and Lupsa, L. (2002) Non-Connected Convexities and Applications. Kluwer Academic Publisher, Dordrecht.

[33] Niculescu, C.P. and Persson, L.E. (2018) Convex Functions and Their Applications. Springer-Verlag, New York.

[34] Pecaric, J.E., Proschan, F. and Tong, Y.L. (1992) Convex Functions and Statistical Applications. Academic Press, New York.

[35] Tiel, J.V. (1984) Convex Analysis. John Wiley and Sons, New York.

[36] Xu, H.K. (1991) Inequalities in Banach Spaces with Applications. Nonlinear Analysis, Theory, Methods and Applications, 16, 1127-1138.

https://doi.org/10.1016/0362-546X(91)90200-K

[37] Bynum, W.L. (1976) Weak Parallelogram laws for Banach Spaces. Canadian Mathematical Bulletin, 19, 269-275.

https://doi.org/10.4153/CMB-1976-042-4

[38] Cheng, R. and Harris, C.B. (2013) Duality of the Weak Parallelogram Laws on Banach Spaces. Journal of Mathematical Analysis and Applications, 404, 64-70.

https://doi.org/10.1016/j.jmaa.2013.02.064

[39] Cheng, R. and Ross, W.T. (2015) Weak Parallelogram Laws on Banach Spaces and Applications to Prediction. Periodica Mathematica Hungarica, 71, 45-58.

https://doi.org/10.1007/s10998-014-0078-4