<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>On Global Partial Calmness for Bilevel Programming Problems with Linear Lower Level Problem</article-title>
      </title-group>
      <contrib-group>
        <aff id="aff0">
          <label>0</label>
          <institution>Leonid I. Minchenko Belarusian State University of Informatics and Radioelectronics 6 Brovki str</institution>
          ,
          <addr-line>220013 Minsk</addr-line>
          ,
          <country country="BY">Belarus</country>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2000</year>
      </pub-date>
      <fpage>412</fpage>
      <lpage>419</lpage>
      <abstract>
        <p>Bilevel programming problems occupy an important place in optimization theory and its various applications. We consider bilevel programs such that their lower level problem is linear with respect to the lower level variable. It is known that such programs are partially calm and, therefore, can be reduced to a one level optimization problem by means of local penalization. One of the goals of our paper is to show that these programs can be reduced to a one level optimization problem via a global exact penalization. In the paper we also prove sufficient conditions of partial calmness for bilevel programs and study the pseudoLipschitzian continuity (Aubin property) of the solution mapping of the lower level problem.</p>
      </abstract>
      <kwd-group>
        <kwd>bilevel programming</kwd>
        <kwd>partial calmness</kwd>
        <kwd>optimal value function</kwd>
        <kwd>penalization</kwd>
        <kwd>solution mapping</kwd>
        <kwd>pseudo-Lipschitzian continuity</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>Introduction</title>
      <p>A bilevel programming problem is a hierarchical optimization problem with two players, where the upper level
player (leader) pursues the goal to minimize his objective function while taking into account the reaction of the
lower level player (follower).</p>
      <p>Let x ∈ Rn, y ∈ Rm. Consider the following bilevel programming problem (BLPP):</p>
      <p>G(x, y) → min,
x ∈ X = {x ∈ Rn |gj (x) ≤ 0 j ∈ J },
y ∈ S(x) =∆ Arg min{f (x, y) |y ∈ F (x) },
(1.1)
(1.2)
(1.3)
where F (x) = {y ∈ Rm |hi(x, y) ≤ 0 i ∈ I }, functions G(x, y), f (x, y), gj (x) and hi(x, y) are continuously
differentiable, I = {1, ..., p}, J = {1, ..., s}.</p>
      <p>Inequalities gj (x) ≤ 0 j ∈ J and hi(x, y) ≤ 0 i ∈ I are usually referred to as the upper level and the lower
level constraints, variables x and y are called the upper level and the lower level variables, respectively.</p>
      <p>Thus, the problem (BLPP) consists in minimizing the upper level (leader's) objective function G(x, y) subject
to the upper level constraints and to solutions y(x) ∈ S(x) of the lower level (follower's ) problem (1.3). The
solution y(x) ∈ S(x) of the lower level problem (1.3) is called the rational reaction of the follower on the leader’s
choice x. The point (x, y) is said to be a feasible point in the problem (BLPP) if x ∈ X, y ∈ S(x). A feasible
point (x0, y0) is called a solution (local solution) of the problem (BLPP) if G(x0, y0) ≤ G(x, y) for all feasible
points (x, y) (for all feasible points from some neighborhood of (x0, y0)).</p>
      <p>Consider the multivalued mappings F : x 7→ F (x) and S : x 7→ S(x). Their graphs and domains are denoted
by grF, domF and grS, domS.</p>
      <p>Note that the problem (BLPP) has been mostly investigated in [Dempe, 2002]-[Henrion, 2011], sometimes it
is referred to as a classical bilevel programming problem or simply a classical bilevel program.</p>
      <p>In [Outrata, 1988] Outrata proposed the following reformulation of the problem (BLPP) into the equivalent
one level problem:</p>
      <p>G(x, y) → mx;iyn, x ∈ X, y ∈ S(x) = {y ∈ F (x) |f (x, y) ≤ φ(x) },
(1.4)
where φ(x) is the optimal value function of the lower level problem (1.3), that is,</p>
      <p>φ(x) = min{f (x, y) |y ∈ F (x) }.</p>
      <p>The problem (BLPP) in the form (1.4) was studied in numerous works [Dempe, 2002], [Ye,
1997][Dempe, 2013]. The main difficulty in solving (1.4) is a nonsmooth constraint involving the value function
φ(x). In [Ye, 1995] Ye and Zhu introduced the now well-known concept of partial calmness and a new method
which allowed to move the nonsmooth constraint from the feasible set to the objective function in (1.4).</p>
      <p>Let (x0, y0) be a feasible point of the problem (BLPP). Recall [Ye, 1995] that the problem (BLPP) in the
form (1.4) is called partial calm at (x0, y0) if there exist a number µ &gt; 0 and a neighborhood V of the point
(x0, y0, 0) in Rn × Rm × R such that G(x, y) − G(x0, y0) + µ |u| ≥ 0 for all (x, y, u) ∈ V such that x ∈ X, y ∈ S(x),
f (x, y) − φ(x) + u = 0.</p>
      <p>In [Ye, 1995] Ye and Zhu proved that the problem (BLPP) in the form (1.4) is partial calm at its local solution
(x0, y0) if and only if there exists a number µ &gt; 0 such that (x0, y0) is a local solution of the partially penalized
problem</p>
      <p>G(x, y) + µ(f (x, y) − φ(x)) → min, x ∈ X, y ∈ F (x).</p>
      <p>
        In their seminal paper [Ye, 1995] where Ye and Zhu introduced the partial calmness, they proved that the
problem (BLPP) with a lower level problem linear in (x, y) is partially calm. Dempe and Zemkoho showed
in [Dempe, 2013] that this proof can be adapted to the case where the lower level problem is linear only in y.
Since the partial exact penalization via partial calmness allows to obtain necessary optimality conditions for
(BLPP) and calculate stationary points for partially calm bilevel programs, partial calmness has drawn a lot
of attention in literature
        <xref ref-type="bibr" rid="ref5 ref6 ref8">(see e.g. [Ye, 2010, Dempe, 2007, Henrion, 2011])</xref>
        . However, an approach involving a
global penalty function can be also interesting.
      </p>
      <p>In Section 2 we consider a bilevel program (BLPP) with the lower level problem linear in y and prove that its
solutions (global) are global solutions of the penalized problem (1.5). Section 3 is devoted to so-called extended
error bound property for multivalued mappings. In Section 4 we prove some sufficient conditions of the
pseudoLipschitz continuity for the solution mapping S and the partial calmness for the problem (BLPP).</p>
      <p>In our paper we consider the bilevel program (BLPP) under the following assumption:</p>
      <p>X ∩ domS = X ∩ domF.</p>
      <p>In other words, (H1) is a quite natural assumption that there exists a rational reaction y(x) ∈ S(x) of the
follower on every leader’s choice x ∈ X ∩ domF .</p>
      <p>Denote by d(v, C) the Euclidean distance between a point v ∈ Rm and a set C ⊂ Rm. Let |v| be the Euclidean
norm of a vector v, B be an open unit ball centered at 0 in Rm or Rn, V (x), V (y) and V"(y) = y + εB, V"(x) =
x+εB (ε &gt; 0) be neighborhoods of x and y. Also denote I(x, y) = {i ∈ I |hi(x, y) = 0 }, h0(x, y) = f (x, y)−φ(x).
(1.5)
(H1)</p>
    </sec>
    <sec id="sec-2">
      <title>Global Penalization in Bilevel Programs</title>
      <p>In this section we consider the problem (BLPP) with an additional assumption
f (x, y) = ⟨a0, y⟩, hi(x, y) = ⟨ai, y⟩ + bi(x),
(H2)
where a0, ai ∈ Rm, bi(x) are scalar functions.</p>
      <p>It is known that under assumption (H2) the tangent cone TF (x)(y) to the set F (x) at a point y ∈ F (x) and a
normal cone NS(x)(y) to the set S(x) at a point y ∈ S(x) can be defined as
and</p>
      <p>TF (x)(y) = clcon(F (x) − y) = {y¯ ∈ Rm |⟨ai, y¯⟩ ≤ 0 i ∈ I(x, y) }</p>
      <p>NS(x)(y) = {
∑</p>
      <p>λiai |λi ≥ 0 i ∈ {0} ∪ I(x, y) },
i2f0g[I(x;y)
respectively.</p>
      <p>Lemma 2.1. Let assumptions (H1) and (H2) hold. Then there exists a number M &gt; 0 such that
d(v, S(x)) ≤ M max{0, ⟨a0, v⟩ − φ(x)}
for all v ∈ F (x) and all x ∈ domS.</p>
      <p>Proof. Set h(x, y) = max{hi(x, y) |i = 0, 1, ..., p } where h0(x, y) = ⟨a0, y⟩ + b0(x), b0(x) = −φ(x). Then
S(x) = {y ∈ Rm |h(x, y) ≤ 0 }. Consider any point x ∈ domS. If v ∈ S(x), the inequality (2.1) holds for all
v ∈ F (x) provided that M &gt; 0.</p>
      <p>In case v ∈ F (x)\S(x) denote as y = y(x, v) the point from S(x) closest to v. Then d(v, S(x)) = |v − y|.
Let l = (v − y) |v − y| 1, v = y(t) = y + tl, where t = |v − y|. In this case l ∈ N (x, y) = {l ∈
Rm l ∈ NS(x)(y) ∩ TF (x)(y), |l| = 1 }.</p>
      <p>Taking into consideration that ⟨ai, v − y⟩ ≤ 0 for i ∈ I(x, y), obtain
h(x, v) = h(x, y(t)) ≥ i2f0mg[aIx(x;y) hi(x, y(t)) = i2f0mg[aIx(x;y){⟨ai, y⟩ + t⟨ai, l⟩ + bi(x)} =</p>
      <p>= t i2f0mg[aIx(x;y)⟨ai, l⟩ = t⟨a0, l⟩ ≥ t l2mN(ixn;y)⟨a0, l⟩ = tδ(x, y)
for all l ∈ N (x, y).</p>
      <p>Suppose that δ(x, y) = ⟨a0, l⟩ ≤ 0. Since hi(x, y) &lt; 0 for i ∈/ I(x, y) and hi(x, y) are continuous in y, there
exists a number ε0 &gt; 0 such that hi(x, y + εl) = ⟨ai, y⟩ + bi(x) + ε⟨ai, l⟩ &lt; 0 for all i ∈/ I(x, y) and all positive
numbers ε ≤ ε0.</p>
      <p>On the other hand, y + εl ∈/ S(x) for all ε0 ≥ ε &gt; 0 and, therefore, the inequality
0 &lt; h(x, y + εl) = i2f0mg[aIx(x;y) hi(x, y + εl) = i2f0mg[aIx(x;y){⟨ai, y⟩ + ε⟨ai, l0⟩ + bi(x)} =</p>
      <p>= ε i2f0mg[aIx(x;y)⟨ai, l⟩ = ε⟨a0, l⟩ = εδ(x, y)
holds for all positive ε ≤ ε0. This contradicts the assumption that δ(x, y) ≤ 0. Therefore, δ(x, y) = ⟨a0(x), (v−
y) |v − y| 1⟩ &gt; 0 for all v ∈ F (x)\S(x).</p>
      <p>The vectors ai don’t depend on x, consequently, the set N (x, y) and the value of δ(x, y) are fully determined
by a choice of subsets I(x, y) in the set I = {1, ..., p}. Since there exists only a finite number of subsets I(x, y)
in the set I, positivity of δ(x, y) implies that δ(x, y) ≥ δ &gt; 0 as y ranges over all boundary points in S(x) and x
ranges over all points in domS.</p>
      <p>Since any point v ∈ F (x)\S(x) can be represented as v = y + tl where l ∈ NS(x)(y) ∩ TF (x)(y), |l| = 1, t &gt; 0,
y is a boundary point in S(x), then from (2.2) follows h(x, v) ≥ δt = δd(v, S(x)) and, therefore,
1
d(v, S(x)) ≤ δ max{0, ⟨ai, v⟩ + bi(x) |i = 0, 1, ..., p } =
1
δ max{0, ⟨a0, v⟩ − φ(x)}.</p>
      <p>Introduce the sets
(2.1)
(2.2)
and their corresponding multivalued mappings</p>
      <p>D(·) : x 7→ D(x) = {y ∈ Rm |(x, y) ∈ D }, C(·) : x 7→ C(x) = {y ∈ Rm |(x, y) ∈ C }.
(2.3)</p>
      <p>Theorem 2.1. Suppose that assumptions (H1) and (H2) hold. Let a feasible point (x0, y0) be a solution
(global) of the problem (BLPP) and the function G be Lipschitz continuous on the set D with Lipschitz constant
l0 &gt; 0. Then there exists a number µ0 &gt; 0 such that for any µ &gt; µ0 the point (x0, y0) is a global solution of the
problem G(x, y) + µ(f (x, y) − φ(x)) → min, (x, y) ∈ D.</p>
      <p>Proof. First of all, note that domD(·) = domC(·) due to (H1). Moreover, from (H1) follows d(y, C(x)) =
d((x, y), {x} × C(x)) ≥ d((x, y), C) for all (x, y) ∈ D.</p>
      <p>In virtue of Proposition 2.4.3 [Clarke, 1983] for all α &gt; l0 the point (x0, y0) is a solution of the problem
G(x, y) + αd((x, y), C) → min, (x, y) ∈ D. Then d((x, y), C) ≤ d(y, C(x)) for all (x, y) ∈ D and, therefore,
G(x, y) + αd(y, C(x)) ≥ G(x0, y0) + αd(y0, C(x0)) = G(x0, y0).</p>
      <p>Then, applying Lemma 2.1 to the set C(x), obtain G(x, y)+αM max{0, h0(x, y)} ≥ G(x0, y0) for all (x, y) ∈ D.
The last inequality is equivalent to the assertion of the theorem with µ0 = l0M .
3</p>
      <p>
        Extended Error Bound Property and Its Application to Bilevel Programming
Let Φ(x) = {y ∈ Rm |qi(x, y) ≤ 0 i ∈ K } where K = {1, ..., r}. In this section we consider a multivalued mapping
Φ : x 7→ Φ(x) assuming that the functions qi : Rn × Rm → R are continuous together with their partial gradients
∇yqi(x, y). Following [Fedorov, 1979, Luderer et al., 2002] introduce the so-called extended error bound property
        <xref ref-type="bibr" rid="ref12 ref13 ref14 ref8">(named R-regularity in [Luderer et al., 2002, Minchenko, 2015, Minchenko, 2011])</xref>
        .
      </p>
      <p>We say the multivalued mapping Φ has the extended error bound property (EEBP) at a point (x0, y0) ∈ grΦ
(relative to X ⊂ Rn) if there exist a number M &gt; 0 and neighborhoods V (x0) and V (y0) such that
d(y, Φ(x)) ≤ M max{0, qi(x, y) i ∈ K}
for all y ∈ V (y0) and all x ∈ V (x0) (x ∈ V (x0) ∩ X).</p>
      <p>We say that the multivalued mapping Φ has the global extended error bound property (GEEBP) if there exists
a number M &gt; 0 such that (3.1) is valid for all y ∈ Rm and all x ∈ domF .</p>
      <p>It is known [Luderer et al., 2002, Borwein, 1986] that in the case of continuously differentiable constraints
qi : Rn ×Rm → R EEBP holds at (x0, y0) ∈ grΦ if y0 satisfies the Mangasarian-Fromovitz constraint qualification
(MFCQ) [Mangasarian, 1967] for the set Φ(x0) and EEBP also holds at (x0, y0) ∈ grΦ relative to domΦ if (x0, y0)
satisfies the constant rank constraint qualification (CRCQ) [Luderer et al., 2002, Rockafellar, 1998].</p>
      <p>
        Recall some definitions
        <xref ref-type="bibr" rid="ref17 ref18">(see, e.g. [Rockafellar, 1998, Klatte, 2015])</xref>
        .
      </p>
      <p>We say that a mapping Φ is pseudo-Lipschitzian or it has the Aubin property relative to X ⊂ Rn at a point
(x0, y0) ∈ grΦ (where x0 ∈ X) if there exist a number lΦ &gt; 0 and neighborhoods V (x0) and V"(y0) such that
Φ(x1) ∩ V"(y0) ⊂ Φ(x2) + lΦ x2 − x1 B
for all x1, x2 ∈ V (x0) (all x1, x2 ∈ V (x0) ∩ X).</p>
      <p>Note that in (3.2) Φ(x) ∩ V"(y0) ̸= ∅ if δ &lt; ε/lΦ.</p>
      <p>A mapping Φ is uniformly bounded at x0 ∈ domΦ if there exist a neighborhood V (x0) and a bounded set Y0
such that Φ(x) ⊂ Y0 for all x ∈ V (x0).</p>
      <p>A mapping Φ is called lower semicontinuous (lsc) at a point (x0, y0) ∈ grΦ (relative to X ⊂ Rn) if for
any neighborhood V (y0) there is a neighborhood V (x0) such that Φ(x) ∩ V (y0) ̸= ∅ for all x ∈ V (x0) (for all
x ∈ V (x0) ∩ X).</p>
      <p>A mapping Φ is called Lipschitz lower semicontinuous at a point (x0, y0) ∈ grΦ if there exist positive numbers
l and δ such that d(y0, Φ(x)) ≤ l x − x0 ∀x ∈ V (x0).</p>
      <p>Let v ∈ Rm, v ∈/ Φ(x). Denote by ΠΦ(x)(v) the set of points from Φ(x) closest to v. Let y ∈ ΠΦ(x)(v). Then
(v − y) |v − y| 1 is a proximal normal [Rockafellar, 1998] to Φ(x) at y.</p>
      <p>Set
(3.1)
(3.2)
y − v + ∑r λi∇yqi(x, y) = 0, λi ≥ 0 and λiqi(x, y) = 0 i ∈ K }.</p>
      <p>|y − v| i=1
The next assertion generalizes results [Minchenko, 2011, Guo, 2013].</p>
      <p>Lemma 3.1. Assume that Φ is lsc at (x0, y0) ∈ grΦ relative to domΦ. Then the following assertions are
equivalent:
(a) Φ has EEBP at (x0, y0) relative to domΦ;
(b) there exists a number M &gt; 0 such that for any sequences
the condition
xk → x0 (xk ∈ domΦ), vk → y0 (vk ∈/ Φ(xk)), yk ∈ ΠΦ(xk)(vk)</p>
      <p>r
ΛvMk (xk, yk) = {λ ∈ Λvk (xk, yk) ∑ |λi| ≤ M } ̸= ∅
i=1
holds starting with some k = k0;
(c) for every pair of sequences xk → x0 (xk ∈ domΦ), vk → y0 (vk ∈/ Φ(xk)) there exists a number M &gt; 0
such that, for all k starting with some k = k0, the condition (3.3) holds for some sequence yk ∈ ΠΦ(xk)(vk).</p>
      <p>Proof. The implication (a) ⇒ (b) follows immediately from the proof of Theorem 2 [Minchenko, 2011]. The
implication (b) ⇒ (c) is obvious. Prove that (c) ⇒ (a). Suppose to the contrary that EEBP does not hold at
(x0, y0) relative to domΦ. This means that there exist sequences xk → x0, xk ∈ domΦ and vk → y0, vk ∈/ Φ(xk)
such that for all k = 1, 2...</p>
      <p>d(vk, Φ(xk)) &gt; k max{0, hi(xk, vk) i ∈ K}.</p>
      <p>Let yk ∈ ΠΦ(xk)(vk). Then there exists a vector λk ∈ Rr such that ∑r λik ≤ M and
i=1
yk − vk
|yk − vk|</p>
      <p>r
+ ∑ λik∇yqi(xk, yk) = 0, λik ≥ 0 i ∈ K(xk, yk), λik = 0 i ∈ K\K(xk, yk)</p>
      <p>i=1
where K(x, y) = {i ∈ K |qi(x, y) = 0 }.</p>
      <p>Since Φ is lsc at (x0, y0), there exists a sequence qk ∈ Φ(xk) such that qk → y0 and vk − yk
This means that yk → y0.</p>
      <p>Then from (3.5) follows
≤ vk − qk .</p>
      <p>r
yk − vk = ⟨∑ λik∇yqi(xk, yk), vk − yk⟩ ≤
i=1
r
= ∑ λikqi(xk, vk) +
i=1</p>
      <p>r
∑ λiko( vk − yk ) ≤
i=1</p>
      <p>r
∑ λik(qi(xk, vk) − qi(xk, yk) + o( vk − yk )) =
i=1</p>
      <p>r
∑ λikqi(xk, vk) +
i=1</p>
      <p>Therefore, d(vk, Φ(xk)) = yk − vk ≤ 2M max{0, qi(xk, vk) i ∈ K}.</p>
      <p>This contradicts (3.4). Thus (c) ⇒ (a).</p>
      <p>Lemma 3.2. Let (x0, y0) ∈ grΦ and |qi(x, y) − qi(x˜, y)| ≤ li |x − x˜| for all x, x˜ ∈ V (x0) and y ∈ V (y0) where
li = const &gt; 0 for all i ∈ K. Assume that the extended error bound property holds at the point (x0, y0) relative
to domΦ. Then Φ is pseudo-Lipschitzian at this point relative to domΦ.</p>
      <p>Proof. If Φ has EEBP at (x0, y0) relative to domΦ, it means that there exist numbers M &gt; 0, δ &gt; 0, ε &gt; 0
such that
d(y, Φ(x)) ≤ M max{0, qi(x, y) i ∈ K}
(3.6)
for all y ∈ V"(y0) and all x ∈ V (x0) ∩ domΦ.</p>
      <p>Denote l = max{li |i = 1, ..., r }. Choose numbers δ and ε such that V (x0) ⊂ V (x0), V"(y0) ⊂ V (y0). Then
d(y0, Φ(x)) ≤ M max{0, qi(x, y0) i ∈ K} ≤
for all x ∈ V (x0) ∩ domΦ.</p>
      <p>This means that Φ is lower Lipschitz continuous at (x0, y0) relative to domΦ and, consequently, Φ(x)∩V"(y0) ̸=
∅ for all x ∈ V 0 (x0) ∩ domΦ and any δ0 &lt; min{δ, ε/l}. Let x, x˜ ∈ V 0 (x0) ∩ domΦ and let y˜ ∈ Φ(x˜) ∩ V"(y0).
Then from (3.6) follows</p>
      <p>d(y˜, Φ(x)) ≤ M max{0, qi(x, y˜) i ∈ K} ≤
≤ M max{0, qi(x, y˜) − qi(x˜, y˜) i ∈ K} ≤ l |x − x˜| .</p>
      <p>This is equivalent to Φ(x˜) ∩ V"(y0) ⊂ Φ(x) + l |x − x˜| B for all x, x˜ ∈ V 0 (x0) ∩ domΦ.
4</p>
      <p>Pseudo-Lipschitz Continuity of Solution Mapping and Partial Calmness in BLPP
In this section we consider the problem BLPP which was stated in Section 1.</p>
      <p>It is well known that the pseudo-Lipschitz continuity of a multivalued mapping at some point of its graph
is closely related to constraint qualifications that hold at that point. However, it is easy to see that if the
Kuhn-Tucker necessary condition holds for a solution of the lower level problem, then such classical constraint
qualifications as the linear independence of the gradients of all active constraints and the Mangasarian-Fromovitz
condition [Mangasarian, 1967] can’t be fulfilled for the mapping S. It means that we need to involve some other
constraint qualifications.</p>
      <p>Following [Janin, 1984] we say that the constant rank constraint quali cation condition (CRCQ) holds for the
mapping S at a point (x0, y0) ∈ C if there exist neighborhoods V (x0) and V (y0) such that rank{∇yhi(x, y) i ∈
K} = const for all x ∈ V (x0), y ∈ V (y0) and any index set K ⊂ {0} ∪ I(x0, y0).</p>
      <p>It is easy to see that C ⊂ grS and CRCQ for the mapping C(·) at a point (x0, y0) ∈ C is equivalent to CRCQ
for the mapping S at this point. Before we prove the next lemma we also note that domC(·) ⊂ domS and the
lower semicontinuity of S relative to domS implies the lower semicontinuity of C(·) relative to domC(·).</p>
      <p>Lemma 4.1. Let (x0, y0) ∈ C. Assume that φ is continuous at x0 relative to domS, S is lsc at (x0, y0)
relative to domS and CRCQ holds for the mapping S at the point (x0, y0). Then S and C(·) have EEBP at
(x0, y0).</p>
      <p>Proof. Suppose that C(·) does not have EEBP at (x0, y0). Then, due to Lemma 3.1 there exist sequences
xk → x0, vk → y0 and {yk} such that C(xk) ̸= ∅, vk ∈/ C(xk), yk ∈ ΠC(xk)(vk), and either Λvk (xk, yk) = ∅ or
d(0, Λvk (xk, yk)) → ∞ where</p>
      <p>p
Λv(x, y) = {λ ∈ Rp+1 y − v + ∑ λi∇yhi(x, y) = 0, λi ≥ 0 and λihi(x, y) = 0 i ∈ {0} ∪ I }.</p>
      <p>|y − v| i=0</p>
      <p>Lower semicontinuity of C(·) implies that yk → y0. Then, without loss of generality, one can assume that
CRCQ holds at all points (xk, yk). This means that Λvk (xk, yk) ̸= ∅ and, therefore, λk → ∅ for each λk ∈
Λvk (xk, yk).</p>
      <p>Since there is only a finite number of possible index sets I(xk, yk), by working with a subsequence if necessary,
we may assume that these index sets are the same for all k = 1, 2, .... In other words, I(xk, yk) = I ⊂ I(x0, y0),
where I doesn’t depend on k = 1, 2, ..., and hi(xk, yk) = 0 i ∈ {0} ∪ I , hi(xk, yk) &lt; 0 i ∈ (I\I ).</p>
      <p>It is known that if Λvk (xk, yk) ̸= ∅, then there exists a maximal linearly independent subfamily
{∇yhi(xk, yk) i ∈ I (xk, yk) ⊂ {0} ∪ I } in the family {∇yhi(xk, yk) i ∈ {0} ∪ I } and a vector λk ∈ Λvk (zk)
such that λik = 0 for all i ∈/ I (xk, yk). By choosing a subsequence if necessary, we may assume that the index
set I (xk, yk) is constant and denote it as I (xk, yk) = I0. Then, for all k = 1, 2, ... {∇yhi(xk, yk) i ∈ I0} is a
maximal linearly independent subfamily in {∇yhi(xk, yk) i ∈ {0} ∪ I } and there exists a vector λk such that
p
|
yk − vk k
yk − vk| + ∑ λi ∇yhi(xk, yk) = 0 λik ≥ 0 i ∈ {0} ∪ I, λik = 0 i ∈/ I0.</p>
      <p>i=0
(4.1)</p>
      <p>Since λk → ∞ in (4.1), without loss of generality one may assume that λk λk 1 → λ¯. Then from (4.2)
follows</p>
      <p>This means that the vectors {∇yhi(x0, y0) i ∈ I0} are linearly dependent and we arrive to the contradiction
with CRCQ for the mapping C(·) at the point (x0, y0).</p>
      <p>Lemma 4.2. Let the solution mapping S be uniformly bounded at x0 and S be lsc relative to domS at some
point (x0, y0) where y0 ∈ S(x0). Then φ is continuous at x0 relative to domS.</p>
      <p>Proof. Consider an arbitrary sequence {xk} ⊂ domS such that xk → x0. Lower semicontinuity of S provides
that there exists yk ∈ S(xk) such that yk → y0. Then
lim sup φ(xk) = lim f (xk, yk) = f (x0, y0) = φ(x0).</p>
      <p>k!1
On the other hand, lim inf φ(xk) ≥ φ(x0) due to Lemma 3.18 [Luderer et al., 2002].</p>
      <p>k!1
Lemma 4.3. Assume that:
1) the solution mapping S is uniformly bounded at x0 and S is lsc at (x0, y0) ∈ grS relative to domS;
2) CRCQ holds for the mapping S at (x0, y0).</p>
      <p>Then the mapping S has EEBP at (x0, y0).</p>
      <p>Proof. According to Lemma 4.2, the function φ is continuous at x0 relative to domS. Then due to Lemma
4.1 EEBP holds at (x0, y0).</p>
      <p>Theorem 4.1. Let (x0, y0) ∈ grS. Assume that:
1) the solution mapping S is uniformly bounded at x0 and S is lsc at (x0, y0) relative to domS;
2) CRCQ holds for the mapping S at (x0, y0).</p>
      <p>Then the mapping S is pseudo-Lipschitzian at (x0, y0) relative to domS.</p>
      <p>Proof. Lemma 4.2 implies that CRCQ holds for the mapping S at (x0, y0). Then due to Lemma 3.2 S is
pseudo-Lipschitzian at (x0, y0).</p>
      <p>Example 4.1. Let y1 − y2 → min, y22 − y1 + x ≤ 0, (y1 − x)2 + y2 ≤ 0, 0 ≤ y3 ≤ 1, 0 ≤ x ≤ 1.</p>
      <p>Let x0 = 0, y0 = (0, 0, 0). CRCQ and other conditions of Theorem 4.1 hold at (x0, y0), therefore, S is
pseudo-Lipschitzian at (x0, y0). On the other hand, it is easy to check that φ(x) = x and S(x) = {y ∈
R3 |y1 = x, y2 = 0, x ≤ y3 ≤ 1 }.</p>
      <p>Theorem 4.2. Let a feasible point (x0, y0) be a local solution of the problem BLPP. Assume that:
1) the function G is Lipschitz continuous on the set D with Lipschitz constant l0 &gt; 0;
2) the solution mapping S is uniformly bounded at x0 and S is lsc at (x0, y0) relative to domS;
3) CRCQ holds for the mapping S at (x0, y0).</p>
      <p>Then there exists a number µ0 &gt; 0 such that for any µ &gt; µ0 the point (x0, y0) is a local solution of the problem
G(x, y) + µ(f (x, y) − φ(x)) → min, (x, y) ∈ D.</p>
      <p>The proof is similar to the proof of Theorem 2.1.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          <source>[Dempe</source>
          , 2002]
          <article-title>Dempe S. Foundations of bilevel programming</article-title>
          .
          <source>Kluwer Acad</source>
          . Publishers, Dordrecht,
          <year>2002</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          <source>[Bard</source>
          , 1998]
          <string-name>
            <surname>Bard J.F.</surname>
          </string-name>
          <article-title>Practical bilevel optimization</article-title>
          .
          <source>Kluwer Acad</source>
          . Publishers, Dordrecht,
          <year>1998</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          <source>[Ye</source>
          , 1995]
          <string-name>
            <given-names>Ye J.J.</given-names>
            ,
            <surname>Zhu</surname>
          </string-name>
          <string-name>
            <surname>D.L.</surname>
          </string-name>
          <article-title>Optimality conditions for bilevel programming problems</article-title>
          . Optimization,
          <volume>33</volume>
          ,
          <year>1995</year>
          .
          <fpage>9</fpage>
          -
          <lpage>27</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          <source>[Ye</source>
          , 1997]
          <string-name>
            <given-names>Ye J.J.</given-names>
            ,
            <surname>Zhu</surname>
          </string-name>
          <string-name>
            <given-names>D.L.</given-names>
            ,
            <surname>Zhu</surname>
          </string-name>
          <string-name>
            <surname>Q.J.</surname>
          </string-name>
          <article-title>Exact penalization and necessary optimality conditions for generalized bilevel programming problems</article-title>
          .
          <source>SIAM J. Optim</source>
          .
          <volume>7</volume>
          (
          <issue>1997</issue>
          ),
          <source>no. 2</source>
          .
          <fpage>481</fpage>
          -
          <lpage>507</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          <source>[Ye</source>
          , 2010]
          <string-name>
            <given-names>Ye J.J.</given-names>
            ,
            <surname>Zhu D</surname>
          </string-name>
          .L.
          <article-title>New necessary optimality conditions for bilevel programs by combining the MPEC and value function approaches</article-title>
          .
          <source>SIAM J. Optim</source>
          .
          <volume>20</volume>
          (
          <year>2010</year>
          ),
          <source>no. 4</source>
          .
          <fpage>1885</fpage>
          <string-name>
            <surname>-</surname>
          </string-name>
          <fpage>1905</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          <source>[Dempe</source>
          , 2007]
          <string-name>
            <surname>Dempe</surname>
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Dutta</surname>
            <given-names>J.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Mordukhovich</surname>
            <given-names>B.S.</given-names>
          </string-name>
          <article-title>New necessary optimality conditions in optimistic bilevel programming</article-title>
          .
          <source>Optimization</source>
          <volume>56</volume>
          (
          <year>2007</year>
          ), no.
          <issue>5-6</issue>
          .
          <fpage>577</fpage>
          -
          <lpage>604</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          <source>[Dempe</source>
          , 2013]
          <string-name>
            <surname>Dempe</surname>
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Zemkoho</surname>
            <given-names>A</given-names>
          </string-name>
          .B.
          <article-title>Bilevel programming: reformulations, constraint qualifications and optimality conditions</article-title>
          .
          <source>Math. Programming 138</source>
          ,
          <year>2013</year>
          .
          <fpage>447</fpage>
          -
          <lpage>473</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          <source>[Henrion</source>
          , 2011]
          <string-name>
            <surname>Henrion</surname>
            <given-names>R.</given-names>
          </string-name>
          , Surowiec T.M.
          <article-title>On calmness conditions in convex bilevel programming</article-title>
          .
          <source>Appl. Anal</source>
          .
          <volume>90</volume>
          (
          <year>2011</year>
          ),
          <source>no. 6</source>
          .
          <fpage>951</fpage>
          -
          <lpage>970</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          <source>[Outrata</source>
          , 1988]
          <string-name>
            <surname>Outrata J.V.</surname>
          </string-name>
          <article-title>A note on the usage of nondifferentiable exact penalties in some special optimization problems</article-title>
          .
          <source>Kybernetika</source>
          <volume>24</volume>
          (
          <year>1988</year>
          ),
          <source>no. 4</source>
          .
          <fpage>251</fpage>
          -
          <lpage>258</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          <source>[Clarke</source>
          , 1983]
          <article-title>Clarke F. Optimization and nonsmooth analysis</article-title>
          . Wiley, New York,
          <year>1983</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          <source>[Fedorov</source>
          , 1979]
          <string-name>
            <surname>Fedorov</surname>
            <given-names>V.V.</given-names>
          </string-name>
          <article-title>Numerical methods of maxmin</article-title>
          . Nauka, Moscow,
          <year>1979</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          [Luderer et al.,
          <year>2002</year>
          ]
          <string-name>
            <surname>Luderer</surname>
            <given-names>B.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Minchenko</surname>
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Satsura</surname>
            <given-names>T</given-names>
          </string-name>
          .
          <article-title>Multivalued analysis and nonlinear programming problems with perturbations</article-title>
          .
          <source>Kluwer Acad</source>
          . Publishers, Dordrecht,
          <year>2002</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          <source>[Minchenko</source>
          , 2015] Minchenko
          <string-name>
            <given-names>L.I.</given-names>
            ,
            <surname>Tarakanov</surname>
          </string-name>
          <string-name>
            <surname>A.N.</surname>
          </string-name>
          <article-title>On second order derivatives of value functions</article-title>
          .
          <source>Optimization</source>
          <volume>64</volume>
          (
          <year>2015</year>
          ).
          <fpage>389</fpage>
          -
          <lpage>407</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          <source>[Minchenko</source>
          , 2011] Minchenko L.,
          <string-name>
            <surname>Stakhovski</surname>
            <given-names>S.</given-names>
          </string-name>
          <article-title>Parametric nonlinear programming problems under the relaxed constant rank condition</article-title>
          .
          <source>SIAM J.Optim</source>
          .
          <volume>21</volume>
          (
          <year>2011</year>
          ),
          <source>no. 1</source>
          .
          <fpage>314</fpage>
          -
          <lpage>332</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          <source>[Borwein</source>
          , 1986]
          <string-name>
            <given-names>Borwein</given-names>
            <surname>J.M.</surname>
          </string-name>
          <article-title>Stability and regular points of inequality systems</article-title>
          .
          <source>J. Optimiz. Theory and Appl</source>
          .
          <volume>48</volume>
          (
          <year>1986</year>
          ).
          <fpage>9</fpage>
          -
          <lpage>52</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref16">
        <mixed-citation>
          <source>[Mangasarian</source>
          , 1967]
          <string-name>
            <given-names>Mangasarian O. L.</given-names>
            ,
            <surname>Fromovitz</surname>
          </string-name>
          <string-name>
            <surname>S.</surname>
          </string-name>
          <article-title>The Fritz John necessary optimality conditions in the presence of equality and inequality constraints</article-title>
          .
          <source>Journal of Mathem. Analysis and Applications</source>
          <volume>17</volume>
          (
          <year>1967</year>
          ).
          <fpage>37</fpage>
          -
          <lpage>47</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref17">
        <mixed-citation>
          <source>[Rockafellar</source>
          , 1998] Rockafellar R.T.,
          <string-name>
            <surname>Wets R.J.-B. Variational</surname>
          </string-name>
          analysis. Springer, Berlin,
          <year>1998</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref18">
        <mixed-citation>
          <source>[Klatte</source>
          , 2015]
          <string-name>
            <given-names>Klatte D.</given-names>
            ,
            <surname>Kummer</surname>
          </string-name>
          <string-name>
            <surname>B</surname>
          </string-name>
          .
          <article-title>On calmness of the argmin mapping in parametric optimization problems</article-title>
          .
          <source>J. Optimiz. Theory and Appl</source>
          .
          <volume>165</volume>
          (
          <year>2015</year>
          ).
          <fpage>708</fpage>
          -
          <lpage>719</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref19">
        <mixed-citation>
          <source>[Guo</source>
          , 2013] Guo
          <string-name>
            <given-names>L.</given-names>
            ,
            <surname>Lin</surname>
          </string-name>
          <string-name>
            <surname>G.H.</surname>
          </string-name>
          <article-title>Notes on some constraint qualifications for mathematical programs with equilibrium</article-title>
          .
          <source>J. Optimiz. Theory and Appl</source>
          .
          <volume>156</volume>
          (
          <year>2013</year>
          ).
          <fpage>600</fpage>
          -
          <lpage>616</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref20">
        <mixed-citation>
          <source>[Janin</source>
          , 1984] Janin R.
          <article-title>Directional derivative of the marginal function in nonlinear programming</article-title>
          .
          <source>Mathematical Programming Studies</source>
          <volume>21</volume>
          (
          <year>1984</year>
          ).
          <fpage>110</fpage>
          -
          <lpage>126</lpage>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>