Multiple Periodic Solutions and Fractal Attractors of Differential Equations with n -Valued Impulses

: Ordinary differential equations with n -valued impulses are examined via the associated Poincaré translation operators from three perspectives: (i) the lower estimate of the number of periodic solutions on the compact subsets of Euclidean spaces and, in particular, on tori; (ii) weakly locally stable (i.e., non-ejective in the sense of Browder) invariant sets; (iii) fractal attractors determined implicitly by the generating vector ﬁelds, jointly with Devaney’s chaos on these attractors of the related shift dynamical systems. For (i), the multiplicity criteria can be effectively expressed in terms of the Nielsen numbers of the impulsive maps. For (ii) and (iii), the invariant sets and attractors can be obtained as the ﬁxed points of topologically conjugated operators to induced impulsive maps in the hyperspaces of the compact subsets of the original basic spaces, endowed with the Hausdorff metric. Five illustrative examples of the main theorems are supplied about multiple periodic solutions (Examples 1–3) and fractal attractors (Examples 4 and 5).

The impulsive maps can be deterministic or stochastic (random), crisp or fuzzy, state dependent or independent, time dependent or independent, single-valued or multivalued.Here, we will be exclusively interested in multivalued, deterministic non-fuzzy state and time independent impulses in subsets of Euclidean spaces.A particular attention will be paid to a subclass of n-valued maps (see Definition 1 below), whence the title of our article.For other sorts of multivalued impulses, see e.g., (Chapter 11 in Reference [12]), References [13][14][15][16][17][18].
As far as we know, the differential equations with n-valued impulses have been tendentiously considered only in Reference [19] for multiple periodic solutions.On the other hand, the recent research of n-valued maps is very active in the topological (i.e., mainly Nielsen) fixed point theory (see e.g., References [20,21] and the earlier survey article of Brown in the handbook [22]).This research is far from being trivial, because there is for instance nothing known about the lower estimates for the number of periodic points of such maps, especially because the number of points of their iterates can be quite arbitrary in general.Moreover, it seems to be difficult to find conditions under which the iterates have an exact given number of points.
Relaxing the strict requirement of exactly n-values in Definition 1 below, we can consider the union operators of n single-valued maps, called the Hutchinson-Barnsley operators.These operators play a crucial role in constructing the fractals as attractors of iterated function systems (see References [23,24]).This relaxation makes the study in a certain sense more liberal, but the n-valued Hutchinson-Barnsley operators then become nothing else but split n-valued maps, which might not be so interesting.Nevertheless, the application of the deep results for the iterated function systems, including the chaotic dynamics on the Hutchinson-Barnsley (fractal) attractors, to impulsive differential equations via the Poincaré translation operators along the trajectories is quite original.
Besides these two novelty applications, our research in this field can be justified by a simple argument that the n-valued impulses extend with no doubts the variability in practical applications.For instance, the repeated vaccination need not be always the same, but they can differ each from other just by finitely many possibilities (e.g., when the doctors have at the same time to their disposal vaccines made by n different producers).
Our paper is organized as follows.After the useful definitions in Preliminaries, we will recall the basic properties and results about n-valued maps (in Section 3) and the Hutchinson-Barnsley operators (in Section 4).These results are neither new, but (in case of n-valued maps) nor so well known.New and original are the applications of these results to impulsive differential equations in R n (in Section 5) and R n /Z n (in Section 6), jointly with the obtained theorems about topological fractals and deterministic chaos in the sense of Devaney (in Section 7).Several illustrative examples are supplied in Sections 5 and 6, jointly with concluding remarks in Section 8.

Preliminaries
In the entire text all topological spaces will be metric.A space X is an absolute neighbourhood retract (written X ∈ ANR) if, for every space Y and every closed subset A ⊂ Y, each continuous map f : A → X is extendable over some open neighbourhood U of A in Y.A space X is an absolute retract (written X ∈ AR) if each f : A → X is extendable over Y. Evidently, if X ∈ AR, then X ∈ ANR.
By a polyhedron, we understand as usually a triangulable space.It is well known that every polyhedron is an ANR-space.An important example of a compact polyhedron will be for us a torus.By the n-torus T n , n ≥ 1, we will mean here either the factor space R n /Z n = (R/Z) n or the Cartesian product S 1 × . . .× S 1

n-times
, where R denotes the set of reals, Z denotes the set of integers, and In particular, for n = 1, T 1 = S 1 becomes a circle.
If not explicitly specified, we will not distinguish between the additive and multiplicative notations, because the logarithm map e 2πs i → s, s ∈ [0, 1], establishes an isomorphism between these two representations.
Let us also note that the relation between the Euclidean space R n and its factorization R n /Z n can be realized by means of the natural projection, sometimes also called a canonical mapping, By a multivalued map ϕ : X Y, we understand ϕ : X → 2 Y \ {∅}.In the entire text, we will still assume that ϕ has closed values.
A multivalued map ϕ : Obviously, in the single-valued case, if f : X → Y is continuous in a multivalued sense, then it is continuous in the usual (single-valued) sense.Furthermore, every continuous map ϕ : X Y has a closed graph Y is continuous with compact values and A ⊂ X is compact, then ϕ(A) is compact, too.The composition ψ • ϕ : X Z of two continuous maps with compact values, ϕ : X Y and ψ : Y Z, is again continuous with compact values.For more details, see for example, References [25,26].
For the single-valued compact continuous maps f : X → X, where X ∈ ANR, we can define the global topological invariants, namely the Lefschetz number L( f ) ∈ Z and the Nielsen number For the single-valued continuous maps on tori, f : T n → T n , the Anosov-type equality N( f ) = |L( f )| holds.For the single-valued continuous maps on the circle (n = 1), f : S 1 → S 1 , we can also define their degree deg( f Besides their existence property, when the existence of a fixed point The Nielsen number N( f ) of a continuous map f : X → X, where X ∈ ANR, gives still the lower estimate of the number of fixed points, that is, , where the symbol # stands for the cardinality of the fixed point set {x ∈ For the definitions and more details, see e.g., References [27,28].Since in Sections 4 and 7 we will also consider hyperspaces endowed with the Hausdorff metric d H , it will be convenient to recall finally their definitions.Hence, if (X, d) is a metric space endowed with the metric d, then the induced hyperspace (K(X), d H ) is defined as For more details, see for example, References [24,25,29,30].
Let us recall their definition and some basic properties.
Definition 1.An n-valued map ϕ : X Y is a continuous multivalued mapping that associates to each x ∈ X an unordered subset of exactly n points of Y.We say that an n-valued map ϕ is split if there are single-valued continuous maps f 1 , . . ., f n : X → Y such that ϕ(x) = { f 1 (x), . . ., f n (x)}, for all x ∈ X.
One can readily check that unlike to admissible maps, where the sets of values are compact and connected (i.e., continua), those of n-valued maps are disconnected.Moreover, unlike to the above definition of general multivalued maps, for the continuity of n-valued maps is sufficient if, for every closed V ⊂ Y, the set {x ∈ X | ϕ(x) ⊂ V} is closed in X.
Let us recall that X is simply connected if and only if it is path-connected and its fundamental (first homotopy) group is trivial.
If X is a compact ANR-space and f : X → X is a single-valued continuous self-map, then the Nielsen number N( f ) of f is well defined (see e.g., References [27,28]).Since a (compact) ANR-space is (uniformly) locally contractible, and so locally path-connected, if X is still simply connected, then all fixed points of f belong to a single (same) fixed point class, whose index is L( f ), where L( f ) stands for the Lefschetz number of f .For its definition and more details, see for example, Reference [27].Therefore, if X is still simply connected, then Summing up, if X is a simply connected compact ANR-space, then the Nielsen relation of any self-map g Subsequently, the Nielsen numbers N(ψ) and N(ϕ) of homotopic n-valued maps ψ = {g 1 , . . . ,g n } : X X and ϕ = { f 1 , . . . ,f n } : X X can be simply defined and calculated by the formula: where the symbol # denotes the cardinality of a given set.If, in particular, X is a compact AR-space, then formula (1) takes the form because L( f k ) = 1, for all k = 1, . . ., n.

Remark 1.
Let us note that (compact) simply connected ANR-spaces are not necessarily (compact) AR-spaces.For example, every n-dimensional unit sphere S n , where n ≥ 2, is a simply connected ANR-space, but not contractible, and so an AR-space.Moreover, L(id S 2 ) = 2 but, according to (1), N(id S 2 ) = 1.
Remark 2. According to the example due to Jezierski [43], there exists a continuous map ϕ : K 2 K 2 , where K 2 ⊂ C is a two-dimensional closed ball in the complex plane C, whose values consist of 1 or 2 or 3 points, which is fixed point free.It justifies the assumption of exactly n-valued maps in Definition 1.
Of course, the assumption of a simple connectedness of X is not necessary for the splitting of n-valued self maps ϕ : X X.For instance, if X = S 1 = R/Z, then an n-valued map ϕ : S 1 S 1 of degree Deg(ϕ) (for its definition, see for example, Theorem 2.1 in Reference [33]) is split if and only if Deg(ϕ) is a multiple of n (see Corollary 5.1 in Reference [33]).
In the case of splitting, we have to our disposal the following lemma.
Lemma 3 (cf.Theorems 2.2 and 2.3 in Reference [33]).If ϕ : S 1 S 1 is a split n-valued map, then its degree Deg(ϕ) equals n-times the classical degree of the maps in the splitting.Furthermore, if ψ : S 1 → S 1 is homotopic in an n-valued way to ϕ (ψ ∼ ϕ), then where f 1 : R → R is the lift of f 1 : R/Z → R/Z.
In the non-split case, the situation becomes obviously more delicate.Nevertheless, for X = S 1 = R/Z, we have even the Wecken property.

Remark 3.
Observe that the equalities (3) in Lemma 3 can be expressed in terms of the Nielsen numbers as follows: Now, we will briefly sketch the definition and basic properties of the Nielsen number for n-valued maps on compact polyhedra, which is essentially due to Schirmer [46] (see also the recent Reference [20]).The restriction to compact polyhedra is caused by the application of such a fixed point index in Reference [46].For more general indices, see for example, References [40,48,49].On the other hand, it will be quite sufficient for our needs in applications.
Using the fact that n-valued maps are locally (but not globally) equivalent to n single-valued continuous functions, which is important for the choice of a suitable index of isolated fixed points, Schirmer proceeded in Reference [46] analogously to a single-valued case (n = 1).Of course, all the technicalities (especially those related to the fixed point index) must have been appropriately elaborated there.For an alternative approach via the lifting classes, see Reference [20].
Hence, let X be a compact polyhedron and ϕ : X X be an n-valued self-map.To obtain the Nielsen number N(ϕ) of ϕ, the fixed points of ϕ were at first divided into finitely many equivalent classes (see Theorem 5.2 in Reference [46]), called fixed point classes or Nielsen classes.Then a suitable fixed point index was associated with each fixed point class.The (Nielsen) classes with non-zero index are called essential.The Nielsen number N(ϕ) of ϕ is the number of essential (Nielsen) fixed point classes.
The Nielsen number N(ϕ) gives the lower estimate of the number of fixed points of ϕ, that is, that any n-valued self-map ϕ : X X has at least N(ϕ) fixed points (cf.Theorem 5.4 in Reference [46]).Furthermore, ϕ satisfies the homotopy invariance, that is, if X is an n-valued homotopy, then N(ϕ 0 ) = N(ϕ 1 ) (cf.Theorem 6.5 in Reference [46]).
Remark 5. Observe that, unlike to equalities (1) and (2), compact polyhedra in Remark 4 need not be simply connected.Let us note that every compact ANR-space is homotopically equivalent to some polyhedron (see e.g., Reference [26]).On the other hand, since any two continuous maps f , g : X → X, where X is a compact absolute retract, are homotopic and L( f ) = 1 (see again e.g., Reference [26]), we get that N( f ) = N(g) = 1, and subsequently we can put N(ϕ) = ∑ n k=1 N( f k ) = n, as in (2).

Hutchinson-Barnsley's Operators
As far as we know, there are no nontrivial results about the periodic point theory, or more precisely periodic orbit theory, for n-valued maps with n > 1.The reason consists in an almost uncontrollable enormous variability of their iterates, by which the existence (and even worse multiplicity) problems of nontrivial periodic orbits seem to be a difficult task.
On the other hand, the theory of iterated function systems (IFS), originated by Hutchinson [23] and extended and popularized by Barnsley [24], allows us to construct compact invariant subsets A ⊂ X of a complete metric space (X, d) of the Hutchinson-Barnsley operators Since a unique A can be obtained as the limit lim and d H stands for the Hausdorff metric (see Section 2), A is called the (global) attractor of the iterated function system (IFS) {X; f 1 , . . . ,f n }.Moreover, the inequality where L = max k=1,...,n L k , holds for the m-th iterate ϕ m (A 0 ) = A m of ϕ, for all m = 1, 2, . . .The attractor A has usually a fractal structure, whose fractal dimension dim A can be estimated from above (i.e., to get its upper bound dim A ≤ D) by means of the Moran-Hutchinson formula (cf.References [23,24]) provided the sets f k (A), k = 1, . . ., n, are either totally disconnected, that is, or just touching (i.e., neither (7), nor with overlaps).If, in particular, f k are similitudes, that is, Observe that condition (7) is much stronger than the one required in Definition 1 for n-valued maps, because ϕ A : A A is certainly continuous and every image set ϕ A (x) must contain exactly n-points, for each x ∈ A. In other words, ϕ A is in particular a split n-valued map on a compact invariant set A ⊂ X.
Moreover, one can prove that every f k A , k = 1, . . ., n, must be injective, because otherwise if there are two points, say a 1 , a 2 ∈ A, such that f k (a 1 ) = f k (a 2 ) = a ∈ A, then a ∈ A would have two addresses, which is impossible (see e.g., Reference [24]).Therefore, we can assume without any loss of generality that f k A are under (7) invertible, for all k = 1, . . ., n.Furthermore, every inversion Hence, we can define the discrete dynamical system (A, S) by the mapping S : A → A, where which can be uniquely defined by S(a) The system (A, S) is called the shift dynamical system, associated with the iterated function system {A; f 1 , . . . ,f n }.It can be proved that it is chaotic in the sense of Devaney (see e.g., References [24,30]), that is, (i) S is sensitive to initial conditions, (ii) S is transitive (if U, V ⊂ A are open subsets, then there exists an integer n * such that U ∩ S n * (V) = ∅), (iii) the set of periodic points of S is dense in A.
We can extend the definition of S to the hyperspace K(A) := {B ⊂ A | B is a compact subset of A}, endowed with the Hausdorff metric d H . Hence, let us define the hypermap S * : K(A) → K(A), where It is well known (see e.g., References [23,24]) that (K(A), d H ) is also compact and S * is continuous in the metric d H .
Haase Theorem 1 in Reference [52] has proved that the system (K(A), S * ) is also chaotic in the sense of Devaney with respect to the Hausdorff metric d H , provided (5) holds.
Summing up, we can formulate the following two well known propositions.
Proposition 1 (cf.References [23,24]).The iterated function system {X; f 1 , . . . ,f n }, where (X, d) is a complete metric space and f k : X → X are contractions, for all k = 1, . . ., n, admits a unique global attractor A ∈ K(X), which can be obtained as the limit lim m→∞ ϕ m (A 0 ) = A, that is, lim m→∞ d H (A m , A) = 0, where A 0 ⊂ X is an arbitrary compact subset and A = ϕ(A), A m = ϕ m (A 0 ), ϕ := ∪ n k=1 f k : X X.Moreover, the inequality (5) holds for the successive approximations A m , m = 1, 2, . . ., of the attractor A, whose fractal dimension dim A can be estimated from above (i.e., dim A ≤ D) by means of the Moran-Hutchinson formula (6).
Proposition 2 (cf.References [24,52]).Let {X; f 1 , . . . ,f n } be an iterated function system as in Proposition 1. Assume, furthermore, that the sets f k (A) are totally disconnected, for all k = 1, . . ., n (see (7)).Then the shift dynamical system (A, S), where S : A → A is defined in (8), which is associated with the iterated function system {A; f 1 , . . . ,f n }, is chaotic in the sense of Devaney.The same is true for the system (K(A), S * ), where the hypermap S * : K(A) → K(A) is defined in (9), with respect to the Hausdorff metric d H . Remark 6.Since every contraction f k on a complete metric space X has, for every k = 1, . . ., n, exactly one fixed point, the Hutchinson-Barnsley operator ϕ := n k=1 f k : X X must have at most n fixed points.Since the attractor A ⊂ X is compact, and so complete, all the fixed points of the restricted contractions f k A : A → A, k = 1, . . ., n, as well as of ϕ A : A A must belong to A.
The existence of a compact invariant subset of the Hutchinson-Barnsley operator can be also obtained in the frame of topological fixed point theory (unlike to "metric" Proposition 1) as follows.
Let us recall here that a metric space X is locally continuum-connected if, for each neighbourhood U of each point x ∈ X, there is a neighbourhood V ⊂ U of x such that each point of V can be connected with x by a subcontinuum (i.e., a compact, connected subset) of U.
Proposition 3 (cf.References [29,53,54]).Let X be a connected, locally continuum-connected metric space and f k : X → X be compact continuous maps, for all k = 1, 2, . . ., n.Then the Hutchinson-Barnsley operator possesses at least one compact invariant set, say X 0 ∈ K(X), such that X 0 = ϕ(X 0 ).If, in particular, X is a Peano continuum (i.e., compact, connected and locally connected) and f k : X → X are continuous maps, for all k = 1, 2, . . ., n, then ϕ possesses at least one compact invariant set ϕ(X 0 ) = X 0 ⊂ X which is non-ejective in the sense of Browder, that is, Remark 7. The Hutchinson-Barnsley operator in Propositions 1 and 3 need not (but can) be exactly n-valued.Furthermore, the contractions f k in Proposition 1 can be more generally replaced, for the existence of a unique attractor A ∈ K(X), by multivalued contractions with compact values.In Proposition 3, compact continuous maps f k can be also replaced, for the existence of a compact invariant set A ∈ K(X) of ϕ (without uniqueness, but including (10)), by compact continuous multivalued maps with compact values.For more details and further possibilities (see References [29,[53][54][55]).
Remark 8. Let us emphasize that the maps f k as well as the Hutchinson-Barnsley operator ϕ can be, under the assumptions of Proposition 3, fixed point free.For instance, rotations on the circle X = S 1 can be so.On the other hand, the non-ejectivity (10) can be regarded as a weak local stability of an invariant set A ∈ K(X) of ϕ.

Application to Impulsive Differential Equations in R n
In this section, the presented results for n-valued maps will be applied to impulsive differential equations in R n .
Consider the vector differential equation where F : R × R n → R n is the Carathéodory mapping such that F(t, x) ≡ F(t + ω, x), for some given ω > 0, that is, Let, furthermore (11) satisfy a uniqueness condition (e.g., a locally Lipschitz condition) and all solutions of (11) entirely exist on the whole line (−∞, ∞).
We can associate to (11) the Poincaré translation operator T ω : R n → R n along its trajectories as follows: T ω (x 0 ) := {x(ω) : x(•) is a solution of (11) such that x(0 It is well known (see e.g., Chapter 1.1 in Reference [56]) that T ω is a homeomorphism such that T k ω = T kω , (i.e., the semi-group property), for every k ∈ N.
One can easily detect the one-to-one correspondence between the kω-periodic solutions of (11), that is, x(t) ≡ x(t + kω) but x(t) ≡ x(t + jω) for j < k, and k-periodic points of T ω , i.e., x 0 = T k ω (x 0 ) but x 0 = T j ω (x 0 ) for j < k, where x 0 = x(0) and j, k are positive integers.Consider also the vector impulsive differential equation    x = F(t, x), t = t j := jω, for some given ω > 0, where F : R × R n → R n is the Carathéodory mapping such that F(t, x) ≡ F(t + ω, x), Equation ( 11) satisfies a uniqueness condition and a global existence of all its solutions on (−∞, ∞).Let, furthermore, I : R n R n be a continuous impulsive mapping.The solutions of impulsive differential Equation ( 13) will be also understood in the Carathéodory sense, that is, x ∈ AC t j , t j+1 , j ∈ Z.For multivalued impulses I, there need not be any longer one-to-one correspondence between the kω-periodic solutions of ( 13) and k-periodic orbits of the composition I • T ω .Nevertheless, every k-periodic orbit of I • T ω implies the existence of a related kω-periodic solution of ( 13), and vice versa.
The first application deals with compact m-valued impulses I in (13).
Theorem 1. Assume that I : R n R n is a compact m-valued map such that I(K 0 ) = K 0 , where K 0 := I(R n ) is a simply connected ANR-space such that K 0 ⊂ K 1 := T ω (K 0 ).Then I K 0 is an m-valued split map of the form I K 0 = {i 1 , . . . ,i m }, and subsequently the number of ω-periodic solutions of (13) is at least equal to # {k = 1, . . ., m | L(i k ) = 0}, where L is the ordinary Lefschetz number.
Proof.Since I is compact and T ω is a homeomorphism, K 0 := I(R n ) and K 1 := T ω (K 0 ) must be compact sets.Since every ANR-space is locally path-connected, K 0 is a compact simply connected and locally path-connected ANR-space, and I K 0 : K 0 K 0 is (according to Lemma 1) a split m-valued map.
Furthermore, since T µω : R n → R n is also a homeomorphism, for every µ ∈ [0, 1], the composition Letting I K 0 = {i 1 , . . . ,I m } and applying formula (1), we get that N(I It means that the mapping I • T ω K 0 has at least such a number of fixed points, which determine the same number of ω-periodic solutions of (13), as claimed.This completes the proof.

Corollary 1.
Let the assumptions of Theorem 1 be satisfied.Assume additionally that K 0 is a (compact) AR-space.Then Equation (13) admits at least m ω-periodic solutions.
We will supply Corollary 1 by two simple illustrative examples.The first example concerns the scalar case (n = 1).
Example 2. Consider (13), where F is as above, and assume that the inequalities f j (t, . . ., x j , . . . ) > 0 holds for all x j ≥ b j , j = 1, . . ., n, f j (t, . . ., x j , . . . ) < 0 holds for all x j ≤ a j , j = 1, . . ., n, hold uniformly for a.a.t ∈ [0, ω] and all the remaining components of x = (x 1 , . . ., x n ), where F(t, x) = ( f 1 (t, x), . . ., f n (t, x)) T .Let I : R n R n be a compact m-valued map such that Since, in view of (16), the inequalities x j (ω, a j ) ≤ a j and x j (ω, b j ) ≥ b j , j = 1, . . ., n, hold for all the components of the solutions x(•, a) and x(•, b) such that x(0, a) = a and x(0, b) = b, where a = (a 1 , . . ., Therefore, the nonlinear impulsive Equation (13) admits, according to Corollary 1, at least m ω-periodic solutions, provided (16) holds, jointly with I : R n R n being a compact m-valued map such that In the non-splitting case (see Remark 5), we can state the following rather theoretical result.
Theorem 2. Assume that I : R n R n is a compact m-valued map such that I(K 0 ) = K 0 , where K 0 := I(R n ) is a (compact) polyhedron such that K 0 ⊂ K 1 := T ω (K 0 ).Then Equation (13) has at least N(I K 0 ) ω-periodic solutions, where the Nielsen number N(I K 0 ) of the restriction I K 0 : K 0 K 0 was defined in Reference [46] by Schirmer (for the sketch, see Section 3).
Proof.We can proceed quite analogously as in the proof of Theorem 1, when avoiding the parts guaranteeing the split arguments, and use the homotopy invariance to get Remark 9.As already pointed out in Remark 4, to calculate the Nielsen number N(I K 0 ) is not an easy task in general.That is why Theorem 2 is rather theoretical than practical.In the splitting case, its calculations can be much easier (see again Remark 4).On the other hand, although compact polyhedra are only special ANR-spaces, they need not be simply connected, as required in Lemma 1, Theorem 1 and Corollary 1 (cf.Remark 5).

Application to Impulsive Differential Equations on R n /Z n
We will also apply the presented results for m-valued maps to impulsive differential equations on tori.
Hence, consider the Equation ( 11) with the same assumptions as in Section 5.
Assuming still that F(t, . . ., x j , . . . ) ≡ F(t, . . ., x j + 1, . . .), j = 1, . . ., n, where x = (x 1 , . . ., x n ), we can also consider (11) on the torus (factor space) R n /Z n , which can be endowed with the metric , for all a, b ∈ R n .The solutions of (11), considered on R n /Z n , will be also understood in the same Carathéodory sense.The associated Poincaré translation operator T ω : R n /Z n → R n /Z n along the trajectories of ( 11), considered on R n /Z n , takes the form T ω := τ • T ω , where T ω was defined in (12), and It is well known (see e.g., Chapter XVII in Reference [57]) that T ω is also a homeomorphism such that T k ω = T kω (i.e., the semi-group property), for every k ∈ N. In particular, for n = 1, T ω is an orientation-preserving homeomorphism.
The following application is for n > 1, like in Theorem 2, rather theoretical than practical in the non-splitting case.Theorem 3. Assume that I : R n R n is an m-valued (mod 1) map, satisfying (18).Then Equation (13) admits, under (17), at least N( I) ω-periodic (mod 1) solutions, where the Nielsen number N( I) of I = τ • I : R n /Z n R n /Z n was defined in Reference [46] by Schirmer (for the sketch, see Section 3).
Proof.Since the torus R n /Z n is a special compact polyhedron and , of a homeomorphism T µω with I must be also a (compact) m-valued map on a compact polyhedron, for every µ ∈ [0, 1], that is, an m-valued homotopy on R n /Z n .Thus, N( I • T ω ) = N( I) holds for the Nielsen numbers, because of the invariance under homotopy The composition I • T ω has therefore at least N( I) fixed points, and each of them determines the existence of an ω-periodic solution on R n /Z n , that is, an ω-periodic (mod 1) solution of ( 13), as claimed.
Remark 10.In the splitting case, the Nielsen number N( I) in Theorem 3 is much easier for calculation by the formula where I := i 1 , . . ., i m , i k : R n /Z n → R n /Z n are endomorphisms defined by integer matrices A k , k = 1, . . ., m, with eigenvalues λ (k) j , j = 1, . . ., n; k = 1, . . ., m. If, in particular, I is an m-valued constant (mod 1), that is, I is an m-valued constant on R n /Z n , then N( I) = m (cf.Remark 4).
For n = 1, the difficult calculation of the Nielsen number N( I) in Theorem 3 can be simplified (even in the non-splitting case) by means of Lemmas 3 and 4 as follows.
If I is still split, then by means of (3) in Lemma 3 we have that Deg( I) = m [i 1 (1) − i 1 (0)], and so we arrive at Corollary 2 can be illustrated by the following simple example.

Fractals and Chaos Determined by Impulsive Differential Equations
In this section, we would like to apply Propositions 1-3 to fractals and chaos, determined implicitly by impulsive differential Equations (13).
Hence, consider (13) and assume the same as in Section 5. We will suppose that I takes the form of the Hutchinson-Barnsley operator, namely where i k : R n → R n , k = 1, . . ., m, are at least continuous functions.The continuous impulsive mapping I : R n R n need not (but can) be any longer m-valued (cf.Remark 7).We start with the application of Proposition 3. Theorem 4. Let I : R n R n be still a compact map and X ⊂ R n be a connected and locally connected subset containing I(R n ), that is, I(R n ) ⊂ X.Then the composition I • T ω X : X X, where T ω is the Poincaré translation operator along the trajectories of (11), defined in (12), possesses at least one compact invariant set X 0 ⊂ X, that is, I • T ω X (X 0 ) = X 0 .If X is still compact, then I • T ω X possesses at least one compact invariant set X 0 ⊂ X which is non-ejective in the sense of Browder, that is, The same is true for the hypersystem (K(X), (I Proof.Since I is by the hypothesis a compact map, I(R n ) must be compact, and there certainly exists a connected and locally connected subset X ⊂ R n such I(R n ) ⊂ X.Since R n is locally compact, X is also locally continuum-connected, that is, for each neighbourhood U of each point x ∈ X, there is a neighbourhood V ⊂ U of x such that each point of V can be connected with x by a subcontinuum of U.
If X is still compact, then it is a Peano's continuum, that is, compact connected and locally connected.
Then the composition I • T ω X : X X is a compact continuous map, which possesses according to Proposition 3 at least one compact invariant set X 0 ∈ K(X).
If X is still compact, and so a Peano's continuum, then I • T ω X possesses, again according to Proposition 3, at least one compact invariant set X 0 ∈ K(X) which is non-ejective in the sense of Browder, that is, (21), which completes the proof.
The proof for the hypersystem (K(X), (I • T ω X ) * ) follows directly by the arguments in References [29,54].Remark 11.The set X 0 can be called a topological fractal in the sense of References [29,53], because it was obtained by means of the Lefschetz-type fixed point theorem as a fixed point in the hyperspace (K(X), d H ), that is, in the frame of the topological fixed point theory (cf.e.g., Reference [22]).Let us note that our terminology differs from the one in e.g., References [58][59][60] where the notion of a topological fractal is understood in a different way.
Theorem 5.There exists a unique global compact attractor B := T ω (A) ∈ K(R n ) of the system R n ; i 1 , . . ., i m , that is, B = I (B) such that B = lim j→∞ I j (B 0 ) holds for any B 0 ∈ K(R n ).
There also exists a unique global compact attractor Proof.The iterated function system {R n ; i 1 , . . . ,i m } has, according to Proposition 1 a unique global compact attractor A ∈ K(R n ) such that A = I(A) (cf.(20)).It can be obtained as the limit lim j→∞ I j (A 0 ), for A 0 ∈ K(R n ).The inequality (cf.( 5)) ..,m L k , holds for the upper estimate of the Hausdorff distance between A and its successive approximations I j (A 0 ).The fractal (Hausdorff) dimension dim A of A can be estimated from above, that is, dim A < D, by means of a unique solution D of the Moran-Hutchinson equation ∑ m k=1 L D k = 1.In case of similitudes i k , for all k = 1, . . ., m, we have the precise value dim A = D. Now, consider the associated (via the Poincaré operator T ω ) systems R n ; i 1 , . . ., i m and R n ; i 1 , . . ., i m .It is known (see e.g., Lemma 2.8 in Reference [61]) that these associated systems have unique global compact attractor B := T ω (A) and C := T ω −1 (A), respectively, that is, B = I (B) and C = I (C), where B = lim j→∞ I j (B 0 ) holds for any Furthermore, the global attractivity of B, C is preserved from the global attractivity of A, jointly with their fractal (Hausdorff) dimension.
Since A = I(A), the existence of unique compact invariant sets B under I and C under I follows easily from the commutative diagrams: Since being an IFS-attractor is known (see e.g., Reference [59]) to be a bi-Lipschitz invariant and even an equi-Hölder invariant, we can simplify the proof of the attractivity of B and C by showing that the associated Poincaré operator T ω is bi-Lipschitz.Since T ω and T ω −1 are homeomorphisms, it is enough to show that the restrictions T ω K 1 as well as T ω −1 K 2 are Lipschitz on any compact subsets K 1 , K 2 ⊂ R n .It is true, provided the right-hand side F in ( 11) is locally Lipschitz, which is the most common uniqueness condition.The same is then true for T ω (see e.g., Reference [56]).
Since the Hausdorff dimension is also a bi-Lipschitz invariant (see e.g., Corollary 2.4 on p. 32 in Reference [62], Chapter 8.3 in Reference [63]), the proof is complete.
Remark 13.The IFS-attractors B, C in Theorem 5 can be called topological fractals in the sense of References [58][59][60] and, jointly with a metric fractal A in the sense of References [29,53], they can be called more precisely Banach fractals in the sense of References [58,59].Remark 14.Although the fixed points x ∈ A of I, described in Remark 6, correspond to the fixed points T ω ( x) ∈ B of I and T ω −1 ( x) ∈ C of I , they need not be anyhow related to ω-periodic solutions of (13), as in Sections 5 and 6.Nevertheless, if i k in (20) are constants, then they determine ω-periodic solutions of (13) (cf.Remark 4) and their images under T ω and T ω −1 are the fixed points of I and I .
The application of Theorem 5 can be illustrated by the following example.
The Sierpi ński-like attractor A = I(A) of the system R 2 ; i 1 , i 2 , i 3 , the attractor B = T 1 (A) of the system R 2 ; i 1 , i 2 , i 3 , and the attractor C = T −1 1 (A) of the system R 2 ; i 1 , i 2 , i 3 are depicted in Figure 1 Now, consider (13) with the same assumptions as above and suppose additionally that i k : R n → R n in (20) are contractions such that the sets i k (A) are totally disconnected, for all k = 1, . . ., m (see (7)), where A ∈ K(R n ) is a unique global attractor of the iterated function system {R n ; i 1 , . . . ,i m } guaranteed by Proposition 1.Let T ω : R n → R n be the Poincaré translation operator along the trajectories of (11), defined in (12).Then we can define the shift dynamical systems (B, S 1 ) and (C, S 2 ), associated with the systems B; i 1 B , . . ., i m B and C; i 1 C , . . ., i m C , defined in Theorem 5, where B := T ω (A) and C := T ω −1 (A).
Concretely, S 1 : B → B is a single-valued continuous mapping such that Since Devaney's chaos is well known to be a topological invariant (i.e., an invariant under a topological conjugacy) on a compact metric space (see e.g., Reference [64]), the both systems (B, S 1 ) and (C, S 2 ) must be chaotic in the sense of Devaney, because the same was recalled to be true for the system (A, i −1 1 , . . ., i −1 m ) in Section 4. Since the induced maps S

Concluding Remarks
As already pointed out in the Introduction, no Nielsen periodic point theory so far exists for n-valued maps.In case of any progress for at least n = 2, we could think about multiple subharmonic (i.e., kω-periodic with k > 1) solutions of impulsive differential equations.
If a uniqueness condition is omitted at given (impulsive) differential equations or, more generally, at (impulsive) differential inclusions, then the associated Poincaré translation operators become multivalued and, in particular, admissible in the sense of Górniewicz (see Reference [26]).Then, however, the needed information about their compositions with n-valued impulsive maps would be lost or, even worse, useless.On the other hand, in the splitting case, their compositions with single-valued selections might be useful for at least partial answers.As a starting point, one can therefore study the same problems as above for differential inclusions with single-valued impulses, as indicated e.g., in Reference [19].
Another problem is to randomize at least some obtained results or to extend some random results like those in References [15,65] along the lines discussed above.Unfortunately, the randomization technique which we have to our disposal (see e.g., Reference [15], and the references therein) is based on the existence of measurable selections whose multiplicity can be lost.In other words, no Nielsen theory so far exists for random fixed points.
A more promissible situation concerning a possible extension and application of the obtained results seems to be for the Hutchinson-Barnsley operators (see Remark 7).The appropriate application can be possible, for instance, for a fractal image compression of pictures.The results might be also randomized.
X} , and the Hausdorff metric d H (•, •) := K(X) × K(X) → [0, ∞) is induced by d as follows: d H (A, B) := max sup a∈A d(a, B), sup b∈B d(A, b) = max sup a∈A ( inf b∈B d(a, b)), sup b∈B ( inf a∈A d(a, b)) , where d(a, B) and d(A, b) stand for the distances between points a, b and sets A, B, respectively.
The fractal (Hausdorff) dimensions dim B and dim C of both B and C can be estimated in the same way from above, that is, dim B ≤ D and dim C ≤ D, by means of a unique solution D of the equation ∑ m k=1 L D k = 1.If i k are similitudes, for all k = 1, . . ., m, then dim B = dim C = D.