Derangements and the p-adic incomplete gamma function ANDREW O’DESKY AND DAVID HARRY RICHMAN Abstract. We introduce a new p-adic analogue of the incomplete gamma function. We also introduce a closely related family of combinatorial sequences counting derangements and arrangements in the wreath product of a cyclic group with a symmetric group. 1. Introduction In 1975, Morita discovered a p-adic analogue of the gamma function. The primary object of this paper is to introduce a p-adic analogue of the incomplete gamma function. The incomplete gamma function Γ(s, z) is the meromorphic function of two variables obtained by integrating the differential zs−1e−z dz over a contour disjoint from the negative real axis which begins at a nonzero complex number z and ends at +∞. It can be shown that for any positive integers n and r r, Γ(n, r) is a rational multiple of 1/e (cf. (15)). Let τp : Q(e) → Qp denote 2 3 the unique field homomorphism satisfying τp(1/e)=1+ p + p /2!+ p /3!+ ··· . Then τpΓ(n, r) is a p-adic number and we may ask for the existence of a p-adic continuous function which interpolates these values. Theorem 1.1. There exists a continuous function Γp : Zp × (1 + pZp) → Zp satisfying Γ (n, r)= τ Γ(n, r) arXiv:2012.04615v3 [math.NT] 1 Apr 2021 p p for any positive integers n and r satisfying r ≡ 1 (mod p). Explicitly, for any s ∈ Zp and r ∈ 1+ pZp we have the formula ∞ s − 1 Γ (s,r) = exp (pr) rs−k−1k! p p k Xk=0 where expp is the p-adic exponential function. Date: April 2, 2021. 2020 Mathematics Subject Classification. 33B20, 11S80, 11B75, 05A05. Key words and phrases. derangements, incomplete gamma function, p-adic analysis. 1 2 ANDREWO’DESKYANDDAVIDHARRYRICHMAN We also prove the following curious formula for Γp involving the floor function. For any n ∈ N and r ∈ Z satisfying r ≡ 1 (mod p), 1/r n 1 −n ⌊e r n!+ 2 ⌋ if r < 0, Γp(n +1, 1/r)= r expp(p/r) · 1/r n (⌊e r n!⌋ if r > 0. This formula is rather surprising as it expresses a p-adically continuous function in terms of a function which crucially uses the linear ordering on the real numbers. 1.1. Combinatorial interpretation using wreath products. The second ob- ject of this paper is to give a combinatorial interpretation for the integral values of the incomplete gamma function. Just as the values of the gamma function at positive integers count permutations, we show that the values of the incomplete gamma function at integral pairs (n, r) ∈ N>0 × Z count more sophisticated com- binatorial gadgets involving the wreath products of symmetric groups with cyclic groups, an interesting and well-studied family of groups. Let Cr denote the cyclic group of order r and let Sn denote the symmetric group of degree n. The permutation action of Sn on [n] := {1, 2,...,n} extends naturally to a permutation action of Cr ≀ Sn on the set Cr × [n] (for details see §3.1). An r-cyclic derangement of degree n is an element of Cr ≀ Sn whose action on Cr × [n] has no fixed point. An r-cyclic arrangement of degree n is a choice of subset A ⊂ [n] and an element of Cr ≀S|A| (r-cyclic derangements were introduced by Assaf [4] while r-cyclic arrangements appear to be a new notion). We assemble these into a single count by defining, for any nonzero integer r, the quantity (−1)n(# of |r|-cyclic derangements of degree n) if r< 0, a(n, r) := (1) (# of r-cyclic arrangements of degree n if r> 0. By use of an inclusion-exclusion argument one can prove the algebraic identity n n a(n, r)= k! rk . k Xk=0 This identity holds for positive and negative r, which serves as motivation for combining these two counts into one. This identity also allows for rapid calcula- tion of a(n, r). For instance, a(1, 2)=3 and a(2, −2) = 5. The values of the incomplete gamma function at integral values are closely related to a(n, r). In analogy to the classical formula Γ(n+1) = n!, one can show that for any nonzero integer r and nonnegative integer n, Γ(n +1, 1/r)= a(n, r)r−ne−1/r. Thus the incomplete gamma function essentially interpolates the counts of cyclic arrangements and cyclic derangements (up to the factor r−ne−1/r). The presence of this factor means that the two problems of p-adic interpolation for the functions DERANGEMENTS AND THE p-ADIC INCOMPLETE GAMMA FUNCTION 3 Γ(n +1, 1/r) and a(n, r) are not the same. Nonetheless we will show that the function a(n, r) also interpolates p-adically in both variables. Theorem 1.2. For any (n, r) ∈ N × Z let a(n, r) ∈ Z be defined by (1). (1) For any fixed nonzero integer r, the function n 7→ a(n, r) extends to a locally analytic metric map on Zp. Furthermore, n 7→ a(n, r) extends to an analytic function on Zp if and only if p divides r. (2) As a function of two variables, (n, r) 7→ a(n, r) extends to a continuous map ap : Zp × Zp → Zp satisfying the inequality |ap(x, r) − ap(y,s)|p ≤ max{|x − y|p, |r − s|p} for any (x, r), (y,s) ∈ Zp × Zp. Connection with pseudo-polynomials. The existence of a p-adic metric map inter- polating a(n, r) means that the values of a(n, r) satisfy rather special congruences. In two prominent cases these congruences have been noticed before. One such case concerns derangements of degree n, which are permutations on n letters without a fixed point. Let dn = a(n, −1) denote the number of derangements of degree n. This defines a sequence whose first few terms are 1, 0, 1, 2, 9, 44, 265, 1854, 14833, 133496, .... The arithmetic properties of this sequence have been the subject of some recent interest, e.g. [30], [29], [18], [28]. It was observed in [14, §2.4] that for any positive integer m, a b (−1) da ≡ (−1) db (mod m) if a ≡ b (mod m). (2) The other prominent case, which is in some sense dual to derangements, concerns arrangements of degree n, which are pairs (A, σ) consisting of a subset A ⊂ [n] and a permutation σ of A. Indeed it was observed in [11] that for any positive integer m, the number an = a(n, 1) of arrangements of degree n satisfies the congruence ac ≡ ad (mod m) if c ≡ d (mod m). (3) In general, a function f : N → Z satisfying the congruence fa ≡ fb (mod m) whenever a, b, m ∈ N satisfy a ≡ b (mod m) is called a pseudo-polynomial. The first half of Theorem 1.2 says that n 7→ a(n, r) is a pseudo-polynomial for any fixed nonzero integer r. In other words, for any fixed nonzero integer r, a(c,r) ≡ a(d,r) if c ≡ d (mod m). The congruences (2) and (3) follow from the first half of Theorem 1.2 by taking r = −1 and r = 1, respectively. Other special cases of Theorem 1.2. The assertion in the first half of Theorem 1.2 that n 7→ a(n, r) is p-adic locally analytic also appears in the literature for the special cases r ∈ {±1}. The fact that n 7→ a(n, −1) (counting derangements) extends to a p-adic locally analytic function was first proven by Barsky [5] using the formal Laplace transform. The fact that n 7→ a(n, 1) (counting arrangements) 4 ANDREWO’DESKYANDDAVIDHARRYRICHMAN extends to a p-adic locally analytic function was proven more recently by Pa¸sol– Zaharescu [23]. 1.2. Permutations with restricted cycle lengths. The third object of this paper is to demonstrate p-adic interpolation properties for permutations with restricted cycle lengths. This gives a complementary generalization of the con- gruence (2) for derangements which is distinct from its generalization to wreath products in Theorem 1.2. A derangement may be thought of as a permutation with no cycle of length one. We generalize this restriction by letting L ⊂ N>0 be any subset, which we consider as the set of ‘permitted’ cycle lengths. Define the count of ‘permissible’ permutations, L dn := #{σ ∈ Sn : every cycle length of σ is in L}. L We prove the following p-adic interpolation result for n 7→ dn . L Theorem 1.3. Let L be a set of positive integers, and let dn denote the number of permutations in Sn whose cycle lengths are in L. L (a) Suppose 1 ∈ L. Then the function n 7→ dn extends to a p-adic continuous function on Zp if and only if p 6∈ L. n L (b) Suppose 1 6∈ L. Then the function n 7→ (−1) dn extends to a p-adic contin- uous function on Zp if and only if p ∈ L. This theorem encompasses a large class of sequences. To illustrate, we record a few examples here; the reader may find more examples in §5. For brevity, we write PL for the set of primes for which the displayed sequence (depending on L) extends to a p-adic continuous function on Zp. Example 1 (Squares). If L1 := {1, 4, 9, 16,...} consists of the perfect squares, then L1 (dn )=(1, 1, 1, 1, 7, 31, 91, 211, 1681, 52417,...) and PL1 = {all primes}.