In this paper, we mainly prove the following conjectures of Z.-W. Sun (J. Number Theory 133 (2013), 2914–2928): let $p>2$ be a prime. If $p=x^2+3y^2$ with $x,y\in \mathbb {Z}$ and $x\equiv 1\ ({\rm {mod}}\ 3)$, then

\[ x\equiv\frac14\sum_{k=0}^{p-1}(3k+4)\frac{f_k}{2^k}\equiv\frac12\sum_{k=0}^{p-1}(3k+2)\frac{f_k}{({-}4)^k}\ ({\rm{mod}}\ p^2), \]

$p\equiv 1\pmod 3$

\[ \sum_{k=0}^{p-1}\frac{f_k}{2^k}\equiv\sum_{k=0}^{p-1}\frac{f_k}{({-}4)^k}\ ({\rm{mod}}\ p^3), \]

$f_n=\sum _{k=0}^n\binom {n}k^3$

$n$

We develop multisummability, in the positive real direction, for generalized power series with natural support, and we prove o-minimality of the expansion of the real field by all multisums of these series. This resulting structure expands both $\mathbb {R}_{\mathcal {G}}$ and the reduct of $\mathbb {R}_{\text {an}^*}$ generated by all convergent generalized power series with natural support; in particular, its expansion by the exponential function defines both the gamma function on $(0,\infty )$ and the zeta function on $(1,\infty )$.

We give bounds on the error in the asymptotic approximation of the log-Gamma function $\ln \unicode[STIX]{x1D6E4}(z)$ for complex $z$ in the right half-plane. These improve on earlier bounds by Behnke and Sommer [Theorie der analytischen Funktionen einer komplexen Veränderlichen, 2nd edn (Springer, Berlin, 1962)], Spira [‘Calculation of the Gamma function by Stirling’s formula’, Math. Comp.25 (1971), 317–322], and Hare [‘Computing the principal branch of log-Gamma’, J. Algorithms25 (1997), 221–236]. We show that $|R_{k+1}(z)/T_{k}(z)|<\sqrt{\unicode[STIX]{x1D70B}k}$ for nonzero $z$ in the right half-plane, where $T_{k}(z)$ is the $k$th term in the asymptotic series, and $R_{k+1}(z)$ is the error incurred in truncating the series after $k$ terms. We deduce similar bounds for asymptotic approximation of the Riemann–Siegel theta function $\unicode[STIX]{x1D717}(t)$. We show that the accuracy of a well-known approximation to $\unicode[STIX]{x1D717}(t)$ can be improved by including an exponentially small term in the approximation. This improves the attainable accuracy for real $t>0$ from $O(\exp (-\unicode[STIX]{x1D70B}t))$ to $O(\exp (-2\unicode[STIX]{x1D70B}t))$. We discuss a similar example due to Olver [‘Error bounds for asymptotic expansions, with an application to cylinder functions of large argument’, in: Asymptotic Solutions of Differential Equations and Their Applications (ed. C. H. Wilcox) (Wiley, New York, 1964), 16–18], and a connection with the Stokes phenomenon.

We use properties of the gamma function to estimate the products $\prod _{k=1}^{n}(4k-3)/4k$ and $\prod _{k=1}^{n}(4k-1)/4k$, motivated by the work of Chen and Qi [‘Completely monotonic function associated with the gamma function and proof of Wallis’ inequality’, Tamkang J. Math.36(4) (2005), 303–307] and Mortici et al. [‘Completely monotonic functions and inequalities associated to some ratio of gamma function’, Appl. Math. Comput.240 (2014), 168–174].

This paper concerns the problem of algebraic differential independence of the gamma function and ${\mathcal{L}}$-functions in the extended Selberg class. We prove that the two kinds of functions cannot satisfy a class of algebraic differential equations with functional coefficients that are linked to the zeros of the ${\mathcal{L}}$-function in a domain $D:=\{z:0<\text{Re}\,z<\unicode[STIX]{x1D70E}_{0}\}$ for a positive constant $\unicode[STIX]{x1D70E}_{0}$.

In this paper, we propose a matrix-variate generalization of the Gauss hypergeometric distribution and study several of its properties. We also derive probability density functions of the product of two independent random matrices when one of them is Gauss hypergeometric. These densities are expressed in terms of Appell’s first hypergeometric function F1 and Humbert’s confluent hypergeometric function Φ1of matrix arguments.

Writing $s\,=\,\sigma \,+\,it$ for a complex variable, it is proved that the modulus of the gamma function, $\left| \Gamma (s) \right|$, is strictly monotone increasing with respect to $\sigma $ whenever $\left| t \right|\,>\,5/4$. It is also shown that this result is false for $\left| t \right|\,\le \,1$.

In this paper we state and prove a new general Hilbert-type inequality in $\mathbb{R}^{n}$ with $k\geq2$ non-conjugate exponents. Using Selberg's integral formula, this result is then applied to obtain explicit upper bounds for the doubly weighted Hardy–Littlewood–Sobolev inequality and some further Hilbert-type inequalities for $k$ non-negative functions and non-conjugate exponents.

The function [Γ(x + 1)]1/x(1 + 1/x)x/x is strictly logarithmically completely monotonic in (0, ∞). The function ψ″ (x + 2) + (1 + x2)/x2(1 + x)2 is strictly completely monotonic in (0, ∞).

In this paper we study different types of planar random motions (performed with constant velocity) with three directions, defined by the vectors dj = (cos(2πj/3), sin(2πj/3)) for j = 0, 1, 2, changing at Poisson-paced times. We examine the cyclic motion (where the change of direction is deterministic), the completely uniform motion (where at each Poisson event each direction can be taken with probability ) and the symmetrically deviating case (where the particle can choose all directions except that taken before the Poisson event). For each of the above random motions we derive the explicit distribution of the position of the particle, by using an approach based on order statistics. We prove that the densities obtained are solutions of the partial differential equations governing the processes. We are also able to give the explicit distributions on the boundary and, for the case of the symmetrically deviating motion, we can write it as the distribution of a telegraph process. For the symmetrically deviating motion we use a generalization of the Bose-Einstein statistics in order to determine the distribution of the triple (N0, N1, N2) (conditional on N(t) = k, with N0 + N1 + N2 = N(t) + 1, where N(t) is the number of Poisson events in [0, t]), where Nj denotes the number of times the direction dj (j = 0, 1, 2) is taken. Possible extensions to four directions or more are briefly considered.

We prove the following.

(1) The inequalities

$$ \biggl(2-\frac{1}{\varGamma(x)}\biggr)^{\!a}+\biggl(2-\frac{1}{\varGamma(1/x)}\biggr)^{\!a}\leq 2\leq\biggl(2-\frac{1}{\varGamma(x)}\biggr)^{\!b}+\biggl(2-\frac{1}{\varGamma(1/x)}\biggr)^{\!b} $$

hold for all $x>0$ if and only if

$$ -1.204\,64\ldots=2+\frac{1}{\gamma}-\frac{1}{6}\biggl(\frac{\pi}{\gamma}\biggr)^{\!2}\leq a\leq0\leq b. $$

(2) For all real numbers $x\in(0,1]$ we have

$$ x^{\alpha}\leq\frac{1}{2}\biggl(\frac{1}{\varGamma(x)}+\frac{1}{\varGamma(1/x)}\biggr)\leq x^{\beta}, $$

with the best possible constants

$$ \alpha=1.321\,76\dots\link{and}\beta=0. $$

These theorems extend and complement a result of Gautschi (from 1974), who proved that for all $x>0$ the harmonic mean of $\varGamma(x)$ and $\varGamma(1/x)$ is greater than or equal to $1$.

AMS 2000 Mathematics subject classification: Primary 33B15; 26D15

Nous prouvons un encadrement optimal pour la quantité $H\left( x,\,s \right)\,=\,\sum{_{n\ge 1}\,\frac{1}{{{\left( x+n \right)}^{s}}}}$ pour $x\,\ge \,0$ et $s\,>\,1$, qui améliore l'encadrement standard par des intégrales. Cet encadrement entraîne des inégalités sur la fonction $\zeta $ de Riemann, et amène à conjecturer la monotonie de la fonction $s\,\mapsto \,{{[(s\,-\,1)\text{ }\!\!\zeta\!\!\text{ (}s\text{)}]}^{\frac{1}{S-1}}}$. On donne des applications à l'étude de la convexité de fonctions liées à la fonction $\Gamma $ d'Euler et à la majoration optimale des fonctions élémentaires intervenant dans les opérateurs de Baskakov. Puis, nous étendons aux fonctions complètement monotones sur ]$0,\,+\infty $[ les résultats établis pour la fonction $x\,\mapsto \,{{x}^{-s}}$, et nous en déduisons des preuves élémentaires du comportement, quand $z$ tend vers 1, des séries génératrices de certaines fonctions arithmétiques. Enfin, nous prouvons qu'une partie du résultat se généralise à une classe de fonctions convexes positives décroissantes.