Sparse geometric graphs with small dilation

Given a set S of n points in , and an integer k such that 0k<n, we show that a geometric graph with vertex set S, at most n−1+k edges, maximum degree five, and dilation O(n/(k+1)) can be computed in time O(nlogn). For any k, we also construct planar n-point sets for which any geometric graph with n−1+k edges has dilation Ω(n/(k+1)); a slightly weaker statement holds if the points of S are required to be in convex position.     

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_π.     

Positive solutions of nonlinear fractional boundary value problems using Adomian decomposition method

We use Adomian decomposition method for solving the fractional nonlinear two-point boundary value problem

where D is Caputo fractional derivative, c is a constant, μ > 0, and F:[0,1]×[0,∞)→[0,∞) a continuous function. The fractional Bratu problem is solved as an illustrative example.     

Fractional Brownian motion with complex variance via random walk in the complex plane and applications

A model of complex-valued fractional Brownian motion has been built up recently as the limit of a random walk in the complex plane, but this model involves radial steps only. It is shown that, by using non-radial steps, this model can be easily extended to define a fractional Brownian motion with complex-valued variance. The relations between complex-valued Brownian motion and the heat equation of order n is clarified and mainly one obtains the general expression of the probability density functions for these processes. One shows that the maximum entropy principle (MPE) provides the probability density of the complex-valued fractional Brownian motion, exactly like for the standard Brownian motion. And lastly, one shows that the heat equation of order 2n (which is the Fokker–Planck equation (FPE) of the complex-valued Brownian motion) has a solution which is similar to that of the FPE of fractional order introduced before by the author, therefore, to some extent, an identification between the complex-valued model via random walk in the complex plane and the model involving a derivative of fractional order.     

Weak solutions for stochastic differential equations with additive fractional noise

In this paper, we show the existence of a weak solution for a stochastic differential equation driven by an additive fractional Brownian motion with Hurst parameter , and a discontinuous drift. The proof of this result is based on the Girsanov theorem for the fractional Brownian motion.     

Covering point sets with two disjoint disks or squares

We study the following problem: Given a set of red points and a set of blue points on the plane, find two unit disks C_R and C_B with disjoint interiors such that the number of red points covered by C_R plus the number of blue points covered by C_B is maximized. We give an algorithm to solve this problem in O(n^8/3log²n) time, where n denotes the total number of points. We also show that the analogous problem of finding two axis-aligned unit squares S_R and S_B instead of unit disks can be solved in O(nlogn) time, which is optimal. If we do not restrict ourselves to axis-aligned squares, but require that both squares have a common orientation, we give a solution using O(n³logn) time.     

Computing homotopic shortest paths efficiently

We give deterministic and randomized algorithms to find shortest paths homotopic to a given collection Π of disjoint paths that wind amongst n point obstacles in the plane. Our deterministic algorithm runs in time , and the randomized algorithm runs in expected time O(k_out+k_inlogn+n(logn)^1+ε). Here k_in is the number of edges in all the paths of Π, and k_out is the number of edges in the output paths.     

The second kind Chebyshev–Newton–Cotes quadrature rule (open type) and its numerical improvement

The weighted Newton–Cotes quadrature rules of open type are denoted by

where w(x) is a positive function and is the step size. Various cases can be selected for the weight function of the above formula. In this paper, we consider as the main weight function and study the general formula:

The precision degree of the above formula is n + 1 for even n’s and is n for odd n’s but if one considers its upper and lower bounds as two additional variables, a nonlinear system will be derived whose solution improves the precision degree of above formula up to degree n + 2 numerically. In this way, some examples are given to show the numerical superiority of our idea.     

Minimal polynomials of the modified de Bruijn sequences

It is well known that a de Bruijn sequence over has the minimal polynomial (x+1)^d, where 2^n-1+nd2ⁿ-1. We study the minimal polynomials of the modified de Bruijn sequences.     

The oscillation of higher-ordernonlinear difference equations

The authors discuss the relation of the oscillation of the following two difference equations,

where m ≥ 2, τ : N → N, N isthe set of integers, |n − τ(n)| ≤ Mfor n N₀, M is a positive integer, is nondecreasing in x, xf(n, x)> 0, as x ≠ 0. Wewill show some relations of the oscillation of the above two equations. Especially, for m to be even, we establish the equivalenceof the oscillation of the above two difference equations.     

Uniform asymptotic expansions of the Pollaczek polynomials

Two uniform asymptotic expansions are obtained for the Pollaczek polynomials P_n(cosθ;a,b). One is for , , in terms of elementary functions and in descending powers of . The other is for , in terms of a special function closely related to the modified parabolic cylinder functions, in descending powers of n. This interval contains a turning point and all possible zeros of P_n(cosθ) in θ(0,π/2].     

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.     

Eigenvalues of association schemes of quadratic forms

The duality and primitivity of the association scheme Qua(n,q) of quadratic forms in n variables and the association scheme Sym(n,q) of symmetric bilinear forms in n variables over the finite field are discussed by Wang et al. [Association schemes of quadratic forms and symmetric bilinear forms, J. Algebraic Combin. 17 (2003) 149–161]. In this paper, eigenvalues of Qua(n,q) are computed, where q is a power of 2. As an application, two fusion schemes of Qua(n,q) are discussed and their dual schemes are constructed.     

Area-preserving approximations of polygonal paths

Let P be an x-monotone polygonal path in the plane. For a path Q that approximates P let W_A(Q) be the area above P and below Q, and let W_B(Q) be the area above Q and below P. Given P and an integer k, we show how to compute a path Q with at most k edges that minimizes W_A(Q)+W_B(Q). Given P and a cost C, we show how to find a path Q with the smallest possible number of edges such that W_A(Q)+W_B(Q)C. However, given P, an integer k, and a cost C, it is NP-hard to determine if a path Q with at most k edges exists such that max{W_A(Q),W_B(Q)}C. We describe an approximation algorithm for this setting. Finally, it is also NP-hard to decide whether a path Q exists such that |W_A(Q)−W_B(Q)|=0. Nevertheless, in this error measure we provide an algorithm for computing an optimal approximation up to an additive error.     

Weak convergence to the multiple Stratonovich integral

We have considered the problem of the weak convergence, as tends to zero, of the multiple integral processes

in the space , where fL²([0,T]ⁿ) is a given function, and {η(t)}_>0 is a family of stochastic processes with absolutely continuous paths that converges weakly to the Brownian motion. In view of the known results when n2 and f(t₁,…,t_n)=1_{t1<t2<<tn}, we cannot expect that these multiple integrals converge to the multiple Itô–Wiener integral of f, because the quadratic variations of the η are null. We have obtained the existence of the limit for any {η}, when f is given by a multimeasure, and under some conditions on {η} when f is a continuous function and when f(t₁,…,t_n)=f₁(t₁)f_n(t_n)1_{t1<t2<<tn}, with f_iL²([0,T]) for any i=1,…,n. In all these cases the limit process is the multiple Stratonovich integral of the function f.     

Space-efficient geometric divide-and-conquer algorithms

We develop a number of space-efficient tools including an approach to simulate divide-and-conquer space-efficiently, stably selecting and unselecting a subset from a sorted set, and computing the kth smallest element in one dimension from a multi-dimensional set that is sorted in another dimension. We then apply these tools to solve several geometric problems that have solutions using some form of divide-and-conquer. Specifically, we present a deterministic algorithm running in time using extra memory given inputs of size n for the closest pair problem and a randomized solution running in expected time and using extra space for the bichromatic closest pair problem. For the orthogonal line segment intersection problem, we solve the problem in time using extra space where n is the number of horizontal and vertical line segments and k is the number of intersections.     

Proximity problems on line segments spanned by points

Finding the closest or farthest line segment (line) from a point are fundamental proximity problems. Given a set S of n points in the plane and another point q, we present optimal O(nlogn) time, O(n) space algorithms for finding the closest and farthest line segments (lines) from q among those spanned by the points in S. We further show how to apply our techniques to find the minimum (maximum) area triangle with a vertex at q and the other two vertices in S{q} in optimal O(nlogn) time and O(n) space. Finally, we give an O(nlogn) time, O(n) space algorithm to find the kth closest line from q and show how to find the k closest lines from q in O(nlogn+k) time and O(n+k) space.     

Complete sets of solutions to partition–distribution equations

In the analysis of algebraic equations in discrete statistical systems, it is often necessary to find the complete sets of non-negative integer solutions k₁, k₂, k₃, … to the partition–distribution equations: k₁+k₂+k₃+…=m and k₁+2k₂+3k₃+…=n. In this paper, two efficient methods, indirect and recurrent, are derived to find the desired results for any given positive integers n and m. An example is provided to demonstrate the merits of these two methods.     

An algorithm for calculating the independence and vertex-cover polynomials of a graph

The independence polynomial, ω(G,x)=∑w_kx^k, of a graph, G, has coefficients, w_k, that enumerate the ways of selecting k vertices from G so that no two selected vertices share an edge. The independence number of G is the largest value of k for which w_k≠0. Little is known of less straightforward relationships between graph structure and the properties of ω(G,x), in part because of the difficulty of calculating values of w_k for specific graphs. This study presents a new algorithm for these calculations which is both faster than existing ones and easily adaptable to high-level computer languages.     

Visibility maps of segments and triangles in 3D

Let T be a set of n triangles in three-dimensional space, let s be a line segment, and let t be a triangle, both disjoint from T. We consider the subdivision of T based on (in)visibility from s; this is the visibility map of the segment s with respect to T. The visibility map of the triangle t is defined analogously. We look at two different notions of visibility: strong (complete) visibility, and weak (partial) visibility. The trivial Ω(n²) lower bound for the combinatorial complexity of the strong visibility map of both s and t is almost tight: we prove an O(n²(n)) upper bound for both structures, where (n) is the extremely slowly increasing inverse Ackermann function. Furthermore, we prove that the weak visibility map of s has complexity Θ(n⁵), and the weak visibility map of t has complexity Θ(n⁷). If T is a polyhedral terrain, the complexity of the weak visibility map is Ω(n⁴) and O(n⁵), both for a segment and a triangle. We also present efficient algorithms to compute all discussed structures.