An efficient algorithm for the parametric resource allocation problem

The parametric resource allocation problem asks to minimize the sum of separable single-variable convex functions containing a parameter λ, Σ_{i = 1}ⁿ (ƒ_i(x_i + λg_i(x_i)), under simple constraints Σ_{i = 1}ⁿ x_i = M, l_i≤x_i≤u_i and x_i: nonnegative integers for i = 1, 2, …, n, where M is a given positive integer, and l_i and u_i are given lower and upper bounds on x_i. This paper presents an efficient algorithm for computing the sequence of all optimal solutions when λ is continuously changed from 0 to ∞. The required time is O(G√Mlog² n + n log n + n log(M/n)), where G = Σ_{i = 1}ⁿ u_i − Σ_{i = 1n} l_i and an evaluation of ƒ_i(·) or g_i(·) is assumed to be done in constant time.     

Extremal properties of strong quadrature weights and maximal mass results for truncated strong moment problems

A bisequence of complex numbers {μ_n}_−∞^∞ determines a strong moment functional satisfying L[xⁿ] = μ_n. If is positive-definite on a bounded interval (a,b) R{0}, then has an integral representation , n=0, ±1, ±2,…, and quadrature rules {w_ni,x_ni} exist such that μ_k = ∑_i=i^n_ns_ni^kw_ni. This paper is concerned with establishing certain extremal properties of the weights w_ni and using these properties to obtain maximal mass results satisfied by distributions ψ(x) representing when only a finite bisequence of moments {μ_k}_k=−nⁿ⁻¹ is given.     

A parallel computing scheme for minimizing a class of large scale functions

This paper gives a parallel computing scheme for minimizing a twice continuously differentiable function with the form

where x = (x^T₁,…,x^T_m)^T and x_i R^n_i, ∑^m_{i = 1}n_i = n, and n a very big number. It is proved that we may use m parallel processors and an iterative procedure to find a minimizer of ƒ(x). The convergence and convergence rate are given under some conditions. The conditions for finding a global minimizer of ƒ(x by using this scheme are given, too. A similar scheme can also be used parallelly to solve a large scale system of nonlinear equations in the similar way. A more general case is also investigated.     

Coupling and harmonic functions in the case of continuous time Markov processes

Consider two transient Markov processes (X^v_t)_tεR·, (X^μ_t)_tεR· with the same transition semigroup and initial distributions v and μ. The probability spaces supporting the processes each are also assumed to support an exponentially distributed random variable independent of the process.

We show that there exist (randomized) stopping times S for (X^v_t), T for (X^μ_t) with common final distribution, L(X^v_S|S < ∞) = L(X^μ_T|T < ∞), and the property that for t < S, resp. t < T, the processes move in disjoint portions of the state space. For such a coupling (S, T) it is shown

where denotes the bounded harmonic functions of the Markov transition semigroup. Extensions, consequences and applications of this result are discussed.     

On the degree of regularity of some equations

In this paper we investigate the behaviour of the solutions of equations Σ_I=1ⁿ a_ix_i = b, where Σ_i=1ⁿ, a_i = 0 and b ≠ 0, with respect to colorings of the set N of positive integers. It turns out that for any b ≠ 0 there exists an 8-coloring of N, admitting no monochromatic solution of x₃ − x₂ = x₂ − x₁ + b. For this equation, for b odd and 2-colorings, only an odd-even coloring prevents a monochromatic solution. For b even and 2-colorings, always monochromatic solutions can be found, and bounds for the corresponding Rado numbers are given. If one imposes the ordering x₁ < x₂ < x₃, then there exists already a 4-coloring of N, which prevents a monochromatic solution of x₃ − x₂ = x₂ − x₁ + b, where b ε N.     

On the Glivenko–Cantelli theorem for generalized empirical processes based on strong mixing sequences

Given \s{X_i, i 1\s} as non-stationary strong mixing (n.s.s.m.) sequence of random variables (r.v.'s) let, for 1 i n and some γ ε [0, 1],

F₁(x)=γP(X_i<x)+(1-γ)P(X_ix)

and

I_i(x)=γI(X_i<x)+(1-γ)I(X_ix)

. For any real sequence \s{C_i\s} satisfying certain conditions, let

.

In this paper an exponential type of bound for P(D_n ), for any >0, and a rate for the almost sure convergence of D_n are obtained under strong mixing. These results generalize those of Singh (1975) for the independent and non-identically distributed sequence of r.v.'s to the case of strong mixing.     

A sufficient and necessary condition for an OWA bag mapping having the self-identity

A mapping ƒ : _n=1^∞Iⁿ → I is called a bag mapping having the self-identity if for every (x₁,…,x_n) ε _i=1^∞Iⁿ we have (1) ƒ(x₁,…,x_n) = ƒ(x_i1,…,x_in) for any arrangement (i₁,…,i_n) of {1,…,n}; monotonic; (3) ƒ(x₁,…,x_n, ƒ(x₁,…,x_n)) = ƒ(x₁,…,x_n). Let {ω_i,n : I = 1,…,n;n = 1,2,…} be a family of non-negative real numbers satisfying Σ_i=1ⁿω_i,n = 1 for every n. Then one calls the mapping ƒ : _i=1^∞Iⁿ → I defined as follows an OWA bag mapping: for every (x₁,…,x_n) ε _i=1^∞Iⁿ, ƒ(x₁,…,x_n) = Σ_i=1ⁿω_i,ny_i, where y_i is the it largest element in the set {x₁,…,x_n}. In this paper, we give a sufficient and necessary condition for an OWA bag mapping having the self-identity.     

Sometimes slow growing is fast growing

The slow growing hierarchy is commonly defined as follows: G₀(x) = 0, G_x−1(x) := G_x(x) + 1 and G_λ(x) := G_λ[x](x) where λ<₀ is a limit and ·[·]:₀∩ Lim × ω → ₀ is a given assignment of fundamental sequences for the limits below ₀. The first obvious question which is encountered when one looks at this definition is: How does this hierarchy depend on the choice of the underlying system of fundamental sequences? Of course, it is well known and easy to prove that for the standard assignment of fundamental sequence the hierarchy (G_x)_x<0 is slow growing, i.e. each G_x is majorized by a Kalmar elementary recursive function.

It is shown in this paper that the slow growing hierarchy (G_x)_x<0 — when it is defined with respect to the norm-based assignment of fundamental sequences which is defined in the article by Cichon (1992, pp. 173–193) — is actually fast growing, i.e. each PA-provably recursive function is eventually dominated by G_x for some <₀. The exact classification of this hierarchy, i.e. the problem whether it is slow or fast growing, has been unsolved since 1992. The somewhat unexpected result of this paper shows that the slow growing hierarchy is extremely sensitive with respect to the choice of the underlying system of fundamental sequences.

The paper is essentially self-contained. Only little knowledge about ordinals less than ₀ — like the existence of Cantor normal forms, etc. and the beginnings of subrecursive hierarchy theory as presented, for example, in the 1984 textbook of Rose — is assumed.     

Existence results for doubly near resolvable (υ, 3, 2)-BIBDs

Let V be a set of υ elements. A (1, 2; 3, υ, 1)-frame F is a square array of side v which satisfies the following properties. We index the rows and columns of F with the elements of V, V={x₁,x₂,…,x_υ}. (1) Each cell is either empty or contains a 3-subset of V. (2) Cell (x_i, x_i) is empty for i=1, 2,…, υ. (3) Row x_i of F contains each element of V−{x_i} once and column x_i of F contains each element of V−{x_i} once. (4) The collection of blocks obtained from the nonempty cells of F is a (υ, 3, 2)-BIBD. A (1, 2; 3, υ, 1)-frame is a doubly near resolvable (υ, 3, 2)-BIBD. In this paper, we first present a survey of existence results on doubly near resolvable (υ, 3, 2)-BIBDs and (1, 2; 3, υ, 1)-frames. We then use frame constructions to provide a new infinite class of doubly near resolvable (υ, 3, 2)-BIBDs by constructing (1, 2; 3, υ, 1)-frames.     

Profile minimization on triangulated triangles

Given graph G=(V,E) on n vertices, the profile minimization problem is to find a one-to-one function f:V→{1,2,…,n} such that ∑_vV(G){f(v)−min_xN[v] f(x)} is as small as possible, where N[v]={v}{x: x is adjacent to v} is the closed neighborhood of v in G. The trangulated triangle T_l is the graph whose vertices are the triples of non-negative integers summing to l, with an edge connecting two triples if they agree in one coordinate and differ by 1 in the other two coordinates. This paper provides a polynomial time algorithm to solve the profile minimization problem for trangulated triangles T_l with side-length l.     

A chromatic partition polynomial

A polynomial in two variables is defined by C_n(x,t)=Σ_πΠnx(Gπ,x)t^|π|, where Π_n is the lattice of partitions of the set {1, 2, …, n}, G_π is a certain interval graph defined in terms of the partition gp, χ(G_π, x) is the chromatic polynomial of G_π and |π| is the number of blocks in π. It is shown that , where S(n, i) is the Stirling number of the second kind and (x)_i = x(x − 1) ··· (x − i + 1). As a special case, C_n(−1, −t) = A_n(t), where A_n(t) is the nth Eulerian polynomial. Moreover, A_n(t)=Σ_πΠna_πt^|π| where a_π is the number of acyclic orientations of G_π.     

First passage time for some stationary processes

For a 1-dependent stationary sequence {X_n} we first show that if u satisfies p₁=p₁(u)=P(X₁>u)0.025 and n>3 is such that 88np₁³1, then

P{max(X₁,…,X_n)u}=ν·μⁿ+O{p₁³(88n(1+124np₁³)+561)}, n>3,

where

ν=1−p₂+2p₃−3p₄+p₁²+6p₂²−6p₁p₂,μ=(1+p₁−p₂+p₃−p₄+2p₁²+3p₂²−5p₁p₂)⁻¹

with

p_k=p_k(u)=P{min(X₁,…,X_k)>u}, k1

and

|O(x)||x|.

From this result we deduce, for a stationary T-dependent process with a.s. continuous path {Y_s}, a similar, in terms of P{max_0skTY_s<u}, k=1,2 formula for P{max_0stY_su}, t>3T and apply this formula to the process Y_s=W(s+1)−W(s), s0, where {W(s)} is the Wiener process. We then obtain numerical estimations of the above probabilities.     

Krylov subspace methods for eigenvalues with special properties and their analysis for normal matrices

In this paper we propose a general approach by which eigenvalues with a special property of a given matrix A can be obtained. In this approach we first determine a scalar function ψ: C → C whose modulus is maximized by the eigenvalues that have the special property. Next, we compute the generalized power iterations uinj + 1 = ψ(A)u_j, j = 0, 1,…, where u₀ is an arbitrary initial vector. Finally, we apply known Krylov subspace methods, such as the Arnoldi and Lanczos methods, to the vector u_n for some sufficiently large n. We can also apply the simultaneous iteration method to the subspace span{x⁽ⁿ⁾₁,…,x⁽ⁿ⁾_k} with some sufficiently large n, where x^(j+1)_m = ψ(A)x^(j)_m, j = 0, 1,…, m = 1,…, k. In all cases the resulting Ritz pairs are approximations to the eigenpairs of A with the special property. We provide a rather thorough convergence analysis of the approach involving all three methods as n → ∞ for the case in which A is a normal matrix. We also discuss the connections and similarities of our approach with the existing methods and approaches in the literature.     

A note on a property of a weak self-similar perfect set

Let S be a compact, weak self-similar perfect set based on a system of weak contractions f_j, j=1,…,m each of which is characterized by a variable contraction coefficient _j(l) as d(f_j(x),f_j(y)) _j(l)d(x,y), d(x,y)<l, l>0. If the relation ∑^m_j=1j(l₀)<1 holds at at least one point l₀, then every nonempty compact metric space is a continuous image of the set S.     

Pancyclic orderings of in-tournaments

An in-tournament is an oriented graph such that the negative neighborhood of every vertex induces a tournament. In this paper, pancyclic orderings of a strong in-tournament D are investigated. This is a labeling, say x₁,x₂,…,x_n, of the vertex set of D such that D[{x₁,x₂,…,x_t}] is Hamiltonian for t=3,4,…,n. Moreover, we consider the related problem on weak pancyclic orderings, where the same holds for t4 and x₁ belongs to an arbitrary oriented 3-cycle. We present sharp lower bounds for the minimum indegree ensuring the existence of a pancyclic or a weak pancyclic ordering in strong in-tournaments.     

Invariant record processes and applications to best choice modelling

Let X₁, X₂,…be identically distributed random variables from an unknown continuous distribution. Further let I_r(1), I_r(2),…be a sequence of indicator functions defined on X₁, X₂,…by I_r(k) = 0 if k < r, I_r(k) = 1 if X_k is a r-record AND = 0 otherwise. Suppose that we observe X₁, X₂,… at times T₁ < T₂ <… where the T_k's are realisations of some regular counting process (N(τ)) defined on the positive half-line. Having observed [0, τ], say, the problem is to predict the future behaviour of the counting processes (R_r(τ, s))_sτ = # r-records in [τ, s]. More specifically the objective of this paper is to show that these processes can be (inhomogeneous) Poisson processes even if (N(τ))_τ0 has dependent increments.

The strong link between optimal selection and optimal stopping of record sequences or record processes, perhaps not fully recognized so far, is pointed out in this paper. It is shown to lead to a unification of the treatment of problems which, at first sight, are rather different. Moreover the stopping of record processes in continuous time can lead to rigorous and elegant solutions in cases where dynamic programming is bound to fail. Several examples will be given to facilitate a comparison with other methods.     

“Drawings since hit tables in lotteries” and a new multivariate geometric distribution

In number lotteries people choose r numbers out of s. Weekly published “drawings since hit tables” indicate how many drawings have taken place since each of the s numbers was last selected as a winning number. Among many lotto players, they enhance the widespread belief that numbers should be “due” if they have not come up for a long time. Under the assumptions of independence of the drawings and equiprobability of all possible combinations, the random s-vectors Y_n, n 1, of entries in a drawings since hit table after n drawings form a Markov chain. The limit distribution of Y_n as n → ∞ is a new multivariate generalization of the geometric distribution. The determination of the distribution of the maximum entry in a drawings since hit table within the first n draws of a lottery seems to be an open problem.     

The connectivity of a graph on uniform points on [0,1]

A random graph G_n(x) is constructed on independent random points U₁,…,U_n distributed uniformly on [0,1]^d, d1, in which two distinct such points are joined by an edge if the l_∞-distance between them is at most some prescribed value 0<x<1. The connectivity distance c_n, the smallest x for which G_n(x) is connected, is shown to satisfy

(1)

For d2, the random graph G_n(x) behaves like a d-dimensional version of the random graphs of Erdös and Rényi, despite the fact that its edges are not independent: c_n/d_n→1, a.s., as n→∞, where d_n is the largest nearest-neighbor link, the smallest x for which G_n(x) has no isolated vertices.     

A Helly theorem for geodesic convexity in strongly dismantlable graphs

A (finite or infinite) graph G is strongly dismantlable if its vertices can be linearly ordered x₀,…, x so that, for each ordinal β < , there exists a strictly increasing finite sequence (i_j)_{0 j n} of ordinals such that i₀ = β, i_n = and x_ij+1 is adjacent with x_ij and with all neighbors of x_ij in the subgraph ofG induced by {x_y: β γ }. We show that the Helly number for the geodesic convexity of such a graph equals its clique number. This generalizes a result of Bandelt and Mulder (1990) for dismantlable graphs. We also get an analogous equality dealing with infinite families of convex sets.     

Full quadrature sums for pth powers of polynomials with Freud weights

In this paper, we provide a solution of the quadrature sum problem of R. Askey for a class of Freud weights. Let r> 0, b (− ∞, 2]. We establish a full quadrature sum estimate

1 p < ∞, for every polynomial P of degree at most n + rn^1/3, where W² is a Freud weight such as exp(−¦x¦), > 1, λ_jn are the Christoffel numbers, x_jn are the zeros of the orthonormal polynomials for the weight W², and C is independent of n and P. We also prove a generalisation, and that such an estimate is not possible for polynomials P of degree M = m(n) if m(n) = n + ξ_nn^1/3, where ξ_n → ∞ as n → ∞. Previous estimates could sum only over those x_jn with ¦x_jn¦ σx_1n, some fixed 0 < σ < 1.