Divisibility of matrices associated with multiplicative functions

Let S={x₁,…,x_n} be a set of n distinct positive integers. For x,y∈S and y<x, we say the y is a greatest-type divisor of x in S if y∣x and it can be deduced that z=y from y∣z,z∣x,z<x and z∈S. For x∈S, let G_S(x) denote the set of all greatest-type divisors of x in S. For any arithmetic function f, let (f(x_i,x_j)) denote the n×n matrix having f evaluated at the greatest common divisor (x_i,x_j) of x_i and x_j as its i,j-entry and let (f[x_i,x_j]) denote the n×n matrix having f evaluated at the least common multiple [x_i,x_j] of x_i and x_j as its i,j-entry. In this paper, we assume that S is a gcd-closed set and . We show that if f is a multiplicative function such that (f∗μ)(d)∈Z whenever and f(a)|f(b) whenever a|b and a,b∈S and (f(x_i,x_j)) is nonsingular, then the matrix (f(x_i,x_j)) divides the matrix (f[x_i,x_j]) in the ring M_n(Z) of n×n matrices over the integers. As a consequence, we show that (f(x_i,x_j)) divides (f[x_i,x_j]) in the ring M_n(Z) if (f∗μ)(d)∈Z whenever and f is a completely multiplicative function such that (f(x_i,x_j)) is nonsingular. This confirms a conjecture of Hong raised in 2004.     

A non-separable solution of the diffusion equation based on the Galerkin’s method using cubic splines

The two dimensional diffusion equation of the form is considered in this paper. We try a bi-cubic spline function of the form as its solution. The initial coefficients C_i,j(0) are computed simply by applying a collocation method; C_i,j = f(x_i, y_j) where f(x, y) = u(x, y, 0) is the given initial condition. Then the coefficients C_i,j(t) are computed by X(t) = e^tQX(0) where X(t) = (C_0,1, C_0,1, C_0,2, … , C_0,N, C_1,0, … , C_N,N) is a one dimensional array and the square matrix Q is derived from applying the Galerkin’s method to the diffusion equation. Note that this expression provides a solution that is not necessarily separable in space coordinates x, y. The results of sample calculations for a few example problems along with the calculation results of approximation errors for a problem with known analytical solution are included.     

Symmetry and specializability in the continued fraction expansions of some infinite products

Let f(x)∈Z[x]. Set f₀(x)=x and, for n?1, define fn(x)=f(fn−1(x)). We describe several infinite families of polynomials for which the infinite product

     

Twisted Hecke L-values and period polynomials

Let f₁,…,fd be an orthogonal basis for the space of cusp forms of even weight 2k on Γ₀(N). Let L(fi,s) and L(fi,χ,s) denote the L-function of fi and its twist by a Dirichlet character χ, respectively. In this note, we obtain a “trace formula” for the values at integers m and n with 0<m,n<2k and proper parity. In the case N=1 or N=2, the formula gives us a convenient way to evaluate precisely the value of the ratio L(f,χ,m)/L(f,n) for a Hecke eigenform f.     

Stochastic matrices and a property of the infinite sequences of linear functionals

Our starting point is the proof of the following property of a particular class of matrices. Let T={T_i,j} be a n×m non-negative matrix such that ∑_jT_i,j=1 for each i. Suppose that for every pair of indices (i,j), there exists an index l such that T_i,l≠T_j,l. Then, there exists a real vector k=(k₁,k₂,…,k_m)^T,k_i≠k_j,i≠j;0<k_i?1, such that, if i≠j.Then, we apply that property of matrices to probability theory. Let us consider an infinite sequence of linear functionals , corresponding to an infinite sequence of probability measures {μ_(·)(i)}_i∈N, on the Borel σ-algebra such that, . The property of matrices described above allows us to construct a real bounded one-to-one piecewise continuous and continuous from the left function f such that

     

Isospectral flows that preserve matrix structure

The matrix A = (a_ij) ∈ S_n is said to lie on a strict undirected graph G if a_ij = 0 (i ≠ j) whenever (i, j) is not in E(G). If S is skew-symmetric, the isospectral flow maintains the spectrum of A. We consider isospectral flows that maintain a matrix A(t) on a given graph G. We review known results for a graph G that is a (generalised) path, and construct isospectral flows for a (generalised) ring, and a star, and show how a flow may be constructed for a general graph. The analysis may be applied to the isospectral problem for a lumped-mass finite element model of an undamped vibrating system. In that context, it is important that the flow maintain other properties such as irreducibility or positivity, and we discuss whether they are maintained.     

An effective iterative method for computing real and complex roots of systems of nonlinear equations

A simple and flexible iterative method is proposed to determine the real or complex roots of any system of nonlinear equations F(x)=0. The idea is based on passing defined functions G_j(x_j),j=1,…,n tangent to F_i(x_j),i,j=1,…,n at an arbitrary starting point. Choosing G_j(x_j) in the form of or or any other reversible function compatible to F_i(x_j), where k is obtained for the best correlation with the function F_i(x_j), gives an added freedom, which in contrast with all existing methods, accelerates the convergence.The method that was first proposed for computing the roots of any single function is now adopted for a system of nonlinear equations. This method is compared to some classical and famous methods such as Newton’s method and Newton-Simpson’s method. The results show the effectiveness and robustness of this new method.     

On almost prime k-tuples

Let k?2 and ai,bi(1?i?k) be integers such that ai>0 and _∏1?i<j?k(aibj−ajbi)≠0. Let Ω(m) denote the total number of prime factors of m. Suppose has no fixed prime divisors. Results of the form where rk is asymptotic to klogk have been obtained by using sieve methods, in particular weighted sieves. In this paper, we use another kind of weighted sieve due to Selberg to obtain improved admissible values for rk.     

Linearized stability analysis of discrete Volterra equations

Various different types of stability are defined, in a unified framework, for discrete Volterra equations of the type x(n)=f(n)+∑ⁿ_j=0K(n,j,x(n)) (n?0). Under appropriate assumptions, stability results are obtainable from those valid in the linear case (K(n,j,x(n))=B(n,j)x(j)), and a linearized stability theory is studied here by using the fundamental and resolvent matrices. Several necessary and sufficient conditions for stability are obtained for solutions of the linear equation by considering the equations in various choices of Banach space , the elements of which are sequences of vectors (, , n,j?0, etc.). We show that the theory, including a number of new results as well as results already known, can be presented in a systematic framework, in which results parallel corresponding results for classical Volterra integral equations of the second kind.     

Multiple partitions, lattice paths and a Burge-Bressoud-type correspondence

A bijection is presented between (1): partitions with conditions f_j+f_j+1≤k−1 and f₁≤i−1, where f_j is the frequency of the part j in the partition, and (2): sets of k−1 ordered partitions (n⁽¹⁾,n⁽²⁾,…,n^(k−1)) such that and , where m_j is the number of parts in n^(j). This bijection entails an elementary and constructive proof of the Andrews multiple-sum enumerating partitions with frequency conditions. A very natural relation between the k−1 ordered partitions and restricted paths is also presented, which reveals our bijection to be a modification of Bressoud’s version of the Burge correspondence.     

On the parity of the partition function

Let N be the set of all positive integers and D a subset of N. Let p(D,n) be the number of partitions of n with parts in D and let |D(x)| denote the number of elements of D not exceeding x. It is proved that if D is an infinite subset of N such that p(D,n) is even for all n?n₀, then |D(x)|?logx/log2−logn₀/log2. Moreover, if D is an infinite subset of N such that p(D,n) is odd for all n?n₀ and , then |D(x)|?logx/log2−logn₀/log2. These lower bounds are essentially the best possible.     

Congruences involving Bernoulli polynomials

Let {B_n(x)} be the Bernoulli polynomials. In the paper we establish some congruences for , where p is an odd prime and x is a rational p-integer. Such congruences are concerned with the properties of p-regular functions, the congruences for and the sum , where h(d) is the class number of the quadratic field of discriminant d and p-regular functions are those functions f such that are rational p-integers and for n=1,2,3,… . We also establish many congruences for Euler numbers.     

Bernoulli matrix and its algebraic properties

In this paper, we define the generalized Bernoulli polynomial matrix B^(α)(x) and the Bernoulli matrix B. Using some properties of Bernoulli polynomials and numbers, a product formula of B^(α)(x) and the inverse of B were given. It is shown that not only B(x)=P[x]B, where P[x] is the generalized Pascal matrix, but also B(x)=FM(x)=N(x)F, where F is the Fibonacci matrix, M(x) and N(x) are the (n+1)×(n+1) lower triangular matrices whose (i,j)-entries are and , respectively. From these formulas, several interesting identities involving the Fibonacci numbers and the Bernoulli polynomials and numbers are obtained. The relationships are established about Bernoulli, Fibonacci and Vandermonde matrices.     

Transition matrices for well-conditioned Markov chains

Let T∈R^n×n be an irreducible stochastic matrix with stationary distribution vector π. Set A = I − T, and define the quantity , where A_j, j = 1, … , n, are the (n − 1) × (n − 1) principal submatrices of A obtained by deleting the jth row and column of A. Results of Cho and Meyer, and of Kirkland show that κ₃ provides a sensitive measure of the conditioning of π under perturbation of T. Moreover, it is known that .In this paper, we investigate the class of irreducible stochastic matrices T of order n such that , for such matrices correspond to Markov chains with desirable conditioning properties. We identify some restrictions on the zero-nonzero patterns of such matrices, and construct several infinite classes of matrices for which κ₃ is as small as possible.     

Nonsingularity of least common multiple matrices on gcd-closed sets

Let n be a positive integer. Let S={x₁,…,x_n} be a set of n distinct positive integers. The least common multiple (LCM) matrix on S, denoted by [S], is defined to be the n×n matrix whose (i,j)-entry is the least common multiple [x_i,x_j] of x_i and x_j. The set S is said to be gcd-closed if for any x_i,x_j∈S,(x_i,x_j)∈S. For an integer m>1, let ω(m) denote the number of distinct prime factors of m. Define ω(1)=0. In 1997, Qi Sun conjectured that if S is a gcd-closed set satisfying max_x∈S{ω(x)}?2, then the LCM matrix [S] is nonsingular. In this paper, we settle completely Sun's conjecture. We show the following result: (i). If S is a gcd-closed set satisfying max_x∈S{ω(x)}?2, then the LCM matrix [S] is nonsingular. Namely, Sun's conjecture is true; (ii). For each integer r?3, there exists a gcd-closed set S satisfying max_x∈S{ω(x)}=r, such that the LCM matrix [S] is singular.     

A complete determination of Rabinowitsch polynomials

Let m be a positive integer and fm(x) be a polynomial of the form fm(x)=x²+x−m. We call a polynomial fm(x) a Rabinowitsch polynomial if for and consecutive integers x=x₀,x₀+1,…,x₀+s−1, |fm(x)| is either 1 or prime. In this paper, we show that there are exactly 14 Rabinowitsch polynomials fm(x).