ON THE p-adic VALUATION OF STIRLING NUMBERS OF THE FIRST KIND PAOLO LEONETTI AND CARLO SANNA arxiv:605.07424v2 [math.nt] 4 Jun 206 Abstract. For all integers n k, define H(n, k) := /(i i k ), where the sum is extended over all positive integers i < < i k n. These quantities are closely related to the Stirling numbers of the first kind by the identity H(n, k) = s(n +, k + )/n!. Motivated by the works of Erdős Niven and Chen Tang, we study the p-adic valuation of H(n, k). In particular, for any prime number p, integer k 2, and x (k )p, we prove that ν p(h(n, k)) < (k )(log p (n/(k )) ) for all positive integers n [(k )p, x] whose base p representations start with the base p representation of k, but at most 3x 0.835 exceptions. We also generalize a result of Lengyel by giving a description of ν 2(H(n, 2)) in terms of an infinite binary sequence.. Introduction It is well known that the n-th harmonic number H n := + 2 + + n is not an integer whenever n 2. Indeed, this result has been generalized in several ways (see, e.g., [2, 7, 3]). In particular, given integers n k, Erdős and Niven [8] proved that H(n, k) := i i k i < <i k n is an integer only for finitely many n and k. Precisely, Chen and Tang [4] showed that H(, ) and H(3, 2) are the only integral values. (See also [] for a generalization to arithmetic progressions.) A crucial step in both the proofs of Erdős Niven and Chen Tang s results consists in showing that, when n and k are in an appropriate range, for some prime number p the p-adic valuation of H(n, k) is negative, so that H(n, k) cannot be an integer. Moreover, a study of the p-adic valuation of the harmonic numbers was initiated by Eswarathasan and Levine [9]. They conjectured that for any prime number p the set J p of all positive integers n such that ν p (H n ) > 0 is finite. Although Boyd [3] gave a probabilistic model predicting that #J p = O(p 2 (log log p) 2+ε ), for any ε > 0, and Sanna [2] proved that J p has asymptotic density zero, the conjecture is still open. Another result of Sanna [2] is that ν p (H n ) = log p n for any n in a subset S p of the positive integers with logarithmic density greater than 0.273. In this paper, we study the p-adic valuation of H(n, k). Let s(n, k) denotes an unsigned Stirling number of the first kind [0, 6.], i.e., s(n, k) is the number of permutations of {,..., n} with exactly k disjoint cycles. Then H(n, k) and s(n, k) are related by the following easy identity. 200 Mathematics Subject Classification. Primary: B73; Secondary: B50, A5. Key words and phrases. Stirling number of the first kind, harmonic numbers, p-adic valuation.
2 Paolo Leonetti and Carlo Sanna Lemma.. For all integers n k, we have H(n, k) = s(n +, k + )/n!. In light of Lemma., and since the p-adic valuation of the factorial is given by the formula [0, p. 57, 4.24] ν p (n!) = n s p(n), p where s p (n) is the sum of digits of the base p representation of n, it follows that ν p (H(n, k)) = ν p (s(n +, k + )) n s p(n), () p hence the study of ν p (H(n, k)) is equivalent to the study of ν p (s(n +, k + )). That explains the title of this paper. In this regard, p-adic valuations of sequences with combinatorial meanings have been studied by several authors (see, e.g., [5, 5, 7, 9, 20, 22]). In particular, the p-adic valuation of Stirling numbers of the second kind have been extensively studied [, 6, 2, 4, 6]. On the other hand, very few seems to be known about the p-adic valuation of Stirling numbers of the first kind. Indeed, up to our knowledge, the only systematic work on this topic is due to Lengyel [8]. Among several results, he showed (see the proof of [8, Theorem.2]) that ν p (H(n, k)) > k log p n + O k (), (2) for all prime numbers p and all integers n k. The main aim of this article is to provide an upper bound for ν p (H(n, k)). In this respect, we believe that inequality (2) is nearly optimal, and our Theorem 2.3 confirms this in the special case when the base p representation of n starts with the base p representation of k. We also formulate the following: Conjecture.. For any prime number p and any integer k, there exists a constant c = c(p, k) > 0 such that ν p (H(n, k)) < c log n for all sufficiently large integers n. 2. Notation and Main results Before state our results, we need to introduce some notation and definition. For any prime number p, we write v a 0,..., a v p := a i p v i, where a 0,..., a v {0,..., p }, a 0 0, (3) i=0 to denote a base p representation. In particular, hereafter, the restrictions of (3) on a 0,..., a v will be implicitly assumed any time we will write something like a 0,..., a v p. We call p-tree of root a 0,..., a v p a set of positive integers T such that: (t) a 0,..., a v p T ; (t2) If b 0,..., b u p T then u v and b i = a i for i = 0,..., v; (t3) If b 0,..., b u p T and u > v then b 0,..., b u p T. Moreover, b 0,..., b u p is a leaf of T if u > v and b 0,..., b u p T but b 0,..., b u p / T. The set of all leaves of T is denoted by T. Finally, the girth of T is the greatest integer g such that for all b 0,..., b u p T we have b 0,..., b u, c p T for at least g values of c {0,..., p }. We are ready to state our results about the p-adic valuation of H(n, k).
On the p-adic valuation of Stirling numbers 3 Theorem 2.. Let p be a prime number and k = e 0,..., e t p + 2. Then there exist a p-tree T p (k) of root k and a nonnegative integer W p (k) such that for any integers n = e 0,..., e t, d t+,..., d s p and r [t +, s] we have: (i) If e 0,..., e t, d t+,..., d r p T p (k) then ν p (H(n, k)) W p (k) + r ks + ; (ii) If e 0,..., e t, d t+,..., d r p T p (k) then ν p (H(n, k)) = W p (k) + r ks. Moreover, the girth of T p (k) is less than p 0.835. equal to. In particular, the girth of T 2 (k) is Note that the case k = has been excluded from the statement. (As mentioned in the introduction, see [3, 9, 2] for results on the p-adic valuation of H(n, ) = H n.) Later, in Section 5, we explain a method to effectively compute the elements of T p (k) for given p and k, and we also illustrate some examples of the results of these computations. Lengyel [8, Theorem 2.5] proved that for each integer m 2 it holds ν 2 (s(2 m, 3)) = 2 m 3m + 3 which, in light of identity (), is in turn equivalent to ν 2 (H(2 m, 2)) = 4 2m. (4) As an application of Theorem 2., we give a corollary that generalizes (4) and provides a quite precise description of ν 2 (H(n, 2)). Corollary 2.2. There exists a sequence f 0, f,... {0, } such that for any integer n = d 0,..., d s 2 2 we have: (i) If d 0 = f 0,..., d s = f s, then ν 2 (H(n, 2)) s; (ii) If d 0 = f 0,..., d r = f r, and d r f r, for some positive integer r s, then ν 2 (H(n, 2)) = r 2s. Precisely, the sequence f 0, f,... can be computed recursively by f 0 = and { if ν 2 (H( f 0,..., f s, 2, 2)) s, f s = 0 otherwise, for any positive integer s. In particular, f 0 =, f =, f 2 = 0. Note that (4) is indeed a consequence of Corollary 2.2. In fact, on the one hand, for m = 2 the identity (4) can be checked quickly. On the other hand, for any integer m 3 we have 2 m = d 0,..., d m 2 with d 0 = = d m =, so that d 0 = f 0, d = f, and d 2 f 2, hence (4) follows from Corollary 2.2(ii), with s = m and r = 2. Finally, we obtain the following upper bound for ν p (H(n, k)). Theorem 2.3. Let p be a prime number, k = e 0,..., e t p + 2, and x (k )p. Then the inequality ν p (H(n, k)) < (k )(log p n log p (k ) ) holds for all n = e 0,..., e t, d t+,..., d s p [(k )p, x], but at most 3x 0.835 exceptions. (5)
4 Paolo Leonetti and Carlo Sanna 3. Preliminaries Let us start by proving the identity claimed in Lemma.. Proof of Lemma.. By [0, Eq. 6.] and s(n +, 0) = 0, we have the polynomial identity n n (X + i) = s(n +, k + )X k, hence + n H(n, k)x k = k= and the claim follows. i= n i= k=0 ( ) X i + = n! n (X + i) = i= n k=0 s(n +, k + ) X k n! From here later, let us fix a prime number p and let k = e 0,..., e t p + 2 and n = d 0,..., d s p be positive integers with s t + and d i = e i for i = 0,..., t. For any a 0,..., a v {0,..., p }, define B p (a 0,..., a v ) := a 0,..., a v p a 0,..., a v p, where by convention a 0,..., a v p = 0 if v = 0, and also B p (a 0,..., a v ) := { c p (i) : i =,..., B p (a 0,..., a v ) } where c p () < c p (2) < denotes the sequence of all positive integers not divisible by p. Lastly, put A p (n, v) := { m {,..., n} : ν p (m) = s v }, for each integer v 0. The next lemma relates A p (n, v) and B p (d 0,..., d v ). Lemma 3.. For each nonnegative integer v s, we have A p (n, v) = { jp s v : j B p (d 0,..., d v ) }. In particular, #A p (n, v) = B p (d 0,..., d v ) and A p (n, v) depends only on p, s, d 0,..., d v. Proof. For m {,..., n}, we have m A p (n, v) if and only if p s v n but p s v+ n. Therefore, s s n n #A p (n, v) = p s v p s v+ = d i p v i d i p v i and as claimed. = i=0 i=0 v v d i p v i d i p v i = d 0,..., d v p d 0,..., d v p i=0 i=0 = B p (d 0,..., d v ), A p (n, v) = { c p (i)p s v : i =,..., #A p (n, v) } = { c p (i)p s v : i =,..., B p (d 0,..., d v ) } = { jp s v : j B p (d 0,..., d v ) },
On the p-adic valuation of Stirling numbers 5 Before stating the next lemma, we need to introduce some additional notation. First, we define t C p (n, k) := A p (n, v) and Π p (k) := free p (j), j C p(n,k) where free p (m) := m/p νp(m) for any positive integer m. Note that, since d i = e i for i = 0,..., t, from Lemma 3. it follows easily that indeed Π p (k) depends only on p and k, and not on n. Then we put t U p (k) := B p (e 0,..., e v )v + t +, while, for a 0,..., a t+v+ {0,..., p }, with v 0 and a i = e i for i = 0,..., t, we set H p(a 0,..., a t+v ) := j j k and 0 v,...,v k t+v j /p v < <j k /p v k v + +v k =U p(k)+v j B p(a 0,...,a v ),...,j k B p(a 0,...,a vk ) H p (a 0,..., a t+v+ ) := H p(a 0,..., a t+v ) + Π p (k) j B p(a 0,...,a t+v+ ) Note that ν p (H p (a 0,..., a t+v+ )) 0, this fact will be fundamental later. The following lemma gives a kind of p-adic expansion for H(n, k). We use O(p v ) to denote a rational number with p-adic valuation greater than or equal to v. Lemma 3.2. We have H(n, k) = s t Proof. Clearly, we can write where and j. H p (d 0,..., d t+v+ ) p v ks+up(k) + O ( p s t ks+up(k)). H(n, k) = V p(n,k) J p (n, k, v) p v Vp(n,k), V p (n, k) := max { ν p (i i k ) : i < < i k n }, J p (n, k, v) := i < <i k n ν p(i i k )=V p(n,k) v free p (i i k ), for each nonnegative integer v V p (n, k). We shall prove that V p (n, k) = ks U p (k). On the one hand, we have t B p (e 0,..., e v ) = t ( ) e0,..., e v p e 0,..., e v p = e 0,..., e t p = k. (6)
6 Paolo Leonetti and Carlo Sanna On the other hand, by (6) and thanks to Lemma 3., we obtain t t #C p (n, k) = #A p (n, v) = B p (e 0,..., e v ) = k. (7) Hence, in order to maximize ν p (i i k ) for positive integers i < < i k n, we have to choose i,..., i k by picking all the k elements of C p (n, k) and exactly one element from A p (n, t + ). Therefore, using again (6) and Lemma 3., we get V p (n, k) = = t #A p (n, v)(s v) + (s t ) (8) t B p (e 0,..., e v )(s v) + (s t ) ( t ) = B p (e 0,..., e v ) + s U p (k) = ks U p (k), as desired. Similarly, if ν p (i i k ) = V p (n, k) v, for some positive integers i < < i k n and some nonnegative integer v s t, then only two cases are possible: ν p (i ),..., ν p (i k ) s t v; or i,..., i k consist of all the k elements of C p (n, k) and one element of A p (n, t + v + ). As a consequence, J p (n, k, v) = i < <i k n ν p(i i k )=V p(n,k) v ν p(i ),...,ν p(i k ) s t v free p (i i k ) + Π p(k) i A p(n,t+v+) free p (i), (9) for all nonnegative integers v s t. By putting v l := s ν p (i l ) and j l := free p (i l ) for l =,..., k, the first sum of (9) can be rewritten as free p (i i k ) = 0 v,...,v k t+v i < <i k (s v )+ +(s v k )=V p(n,k) v i A p(n,v ),...,i k A p(n,v k ) 0 v,...,v k t+v j /p v < <j k /p v k v + +v k =U p(k)+v j B p(d 0,...,d v ),...,j k B p(d 0,...,d vk ) where we have also made use of (8) and Lemma 3., hence j j k = H p(d 0,..., d t+v ), J p (n, k, v) = H p (d 0,..., d t+v+ ), (0) for any nonnegative integer v s t. At this point, being s > t, by (8) it follows that V p (n, k) > s t, hence H(n, k) = s t J p (n, k, v) p v ks+up(k) + O ( p s t ks+up(k)), ()
On the p-adic valuation of Stirling numbers 7 since clearly ν p (J p (n, k, v)) 0 for any nonnegative integer v V p (n, k). In conclusion, the claim follows from (0) and (). Finally, we need two lemmas about the number of solutions of some congruences. For rational numbers a and b, we write a b mod p to mean that ν p (a b) > 0. Lemma 3.3. Let r be a rational number and let x, y be positive integers with y < p. Then the number of integers v [x, x+y] such that H v r mod p is less than 3 2 y2/3 +. Proof. The case r = 0 is proved in [2, Lemma 2.2] and the proof works exactly in the same way even for r 0. Lemma 3.4. Let q be a rational number and let a be a positive integer. number of d {0,..., p } such that Then the a+d i=a q mod p (2) c p (i) is less than p 0.835. Proof. It is easy to see that there exists some h {0,..., p 2} such that { c p (a) + i a for i = a,..., a + h, c p (i) = c p (a) + i a + for i = a + h +,..., a + p. Therefore, by putting x := c p (a), y := h, and r := q + H x in Lemma 3.3, we get that the number of d h satisfying (2) is less than 3 2 h2/3 +. Similarly, by putting x := c p (a) + h + 2, y := p h 2, and a+h r := q + H x c p (i) in Lemma 3.3, we get that the number of d [h +, p ] satisfying (2) is less than 3 2 (p h 2)2/3 +. Thus, letting N be the number of d {0,..., p } that satisfy (2), we have i=a N 3 2 h2/3 + + 3 ( ) p 2 2/3 2 (p h 2)2/3 + 3 + 2. 2 Furthermore, it is clear the d and d + cannot both satisfy (2), hence N p/2. Finally, a little computation shows that the maximum of ( ( ) p 2 2/3 p log p (min )) 3 + 2, 2 2 is obtained for p = 59 and is less than 0.835, hence the claim follows.
8 Paolo Leonetti and Carlo Sanna 4. Proof of Theorem 2. Now we are ready to prove Theorem 2.. For any a 0,..., a t+u+ {0,..., p }, with u 0 and a i = e i for i = 0,..., t, let Σ p (a 0,..., a t+u+ ) := u H p (a 0,..., a t+v+ ) p v. Furthermore, define the sequence of sets T (0) p (k), T () p (k),... as follows: T (0) { e 0,..., e t p }, and for any integer u 0 put a 0,..., a t+u+ p T (u+) p if a 0,..., a t+u p T (u) p p (k) := (k) if and only (k) and ν p (Σ p (a 0,..., a t+u+ )) u +. At this point, setting T p (k) := u=0 T p (u) (k), it is straightforward to see that T p (k) is a p-tree of root e 0,..., e t p. Put W p (k) := U p (k) t. If d 0,..., d r p T p (k) then, by the definition of T p (k), we have ν p (Σ p (d 0,..., d r )) r t. Therefore, by Lemma 3.2 it follows that ν p (H(n, k)) W p (k) + r ks +, and this proves (i). If d 0,..., d r p T p (k) then r > t and d 0,..., d r T p (k) but d 0,..., d r / T p (k), so that ν p (Σ p (d 0,..., d r )) r t. (3) Now we distinguish between two cases. If r = t +, then ν p (Σ p (d 0,..., d t+ )) = 0 and by Lemma 3.2 we obtain ν p (H(n, k)) = W p (k) + r ks. If r > t + then by d 0,..., d r T p (k) we get that ν p (Σ p (d 0,..., d r )) r t, which together with (3) and Σ p (d 0,..., d r ) = Σ p (d 0,..., d r ) + H p (d 0,..., d r ) p r t implies that ν p (Σ p (d 0,..., d r )) = r t, hence by Lemma 3.2 we get ν p (H(n, k)) = W p (k) + r ks, and also (ii) is proved. It remains only to bound the girth of T p (k). Let u be a nonnegative integer and pick a 0,..., a t+u p T p (u) (k). By the definition of T p (u+) (k), we have a 0,..., a t+u+ p T (u+) p (k) if and only if ν p (Σ p (a 0,..., a t+u+ )) u +, which in turn is equivalent to u H p (a 0,..., a t+v+ ) p v u + H p(a 0,..., a t+u ) + Π p (k) u H p (a 0,..., a t+v+ ) p v u 0 mod p. Using the definition of B p (a 0,..., a t+u+ ) and the facts that B p (a 0,..., a t+u+ ) = a t+u+ + (p ) j B p(a 0,...,a t+u+ ) u a v p u v, j (4)
On the p-adic valuation of Stirling numbers 9 and ν p (Π p (k)) = 0, we get that (4) is equivalent to where Π p (k) a+a t+u+ i=a ( u a c p (i) i= c p (i) ) H p (a 0,..., a t+v+ ) p v u + H p(a 0,..., a t+u ) mod p, (5) a := (p ) u a v p u v. Note that both a and the right-hand side of (5) do not depend on a t+u+. As a consequence, by Lemma 3.4 we get that a 0,..., a t+u+ p T p (u+) (k) for less than p 0.835 values of a t+u+ {0,..., p }. Thus the girth of T p (k) is less than p 0.835. Finally, consider the case p = 2. Obviously, /c 2 (i) mod 2 for any positive integer i, while the right-hand side of (5) is equal to 0 or (mod 2). Therefore, there exists one and only one choice of a t+u+ {0, } such that (5) is satisfied. This means that the girth of T 2 (k) is exactly. The proof is complete. 5. The computation of T p (k) Given p and k, it might be interesting to effectively compute the elements of T p (k). Clearly, T p (k) could be infinite and this is indeed the case when p = 2, since by Theorem 2. we know that T 2 (k) has girth hence the computation should proceed by first enumerating all the elements of T p (0) (k), then all the elements of T p () (k), and so on. An obvious way to do this is using the recursive definition of the T p (u) (k) s. However, it is easy to see how this method is quite complicated and impractical. A better idea is noting that, taking r = s in Theorem 2., we have T p (u+) (k) = { a 0,..., a t+u, b p : a 0,..., a t+u p T p (u) (k), (6) ν p (H( a 0,..., a t+u, b p, k)) W p (k) (k )s + }, for all integers u 0. Therefore, starting from T p (0) (k) = { e 0,..., e t p }, formula (6) gives a way to compute recursively all the elements of T p (k). In particular, if T p (k) is finite, then after sufficient computation one will get T p (u) (k) = for some positive integer u, so the method actually proves that T p (k) is finite. The authors implemented this algorithm in SageMath, since it allows computations with arbitrary-precision p-adic numbers. In particular, they found that T 3 (2),..., T 3 (6) are all finite sets, with respectively 8, 24, 6, 7, 23 elements, while the cardinality of T 3 (7) is at least 43. Through these numerical experiments, it seems that, in general, T p (k) does not exhibit any trivial structure (see Figures, 2, 3), hence the question of the finiteness of T p (k) is probably a difficult one.
0 Paolo Leonetti and Carlo Sanna 6. Proof of Corollary 2.2 Only for this section, let us focus on the case p = 2 and k = 2, so that t = 0, e 0 =, and W 2 (2) = 0. Thanks to Theorem 2. we know that T 2 (2) has girth. Hence, it follows easily that there exists a sequence f 0, f,... {0, } such that T (u) 2 (2) = { f 0,..., f u 2 } for all integers u 0. In particular, f 0 = e 0 =. At this point, (i) and (ii) are direct consequences of Theorem 2., while the recursive formula (5) is just a special case of (6). 7. Proof of Theorem 2.3 It is easy to see that #T (u) p (k) < p 0.835u, for any positive integer u. On the one hand, the number of n = e 0,..., e t, d t+,..., d s p [(k )p, x] such that d 0,..., d r p / T p (k) for any integer r [t +, s] is less than log p x t u= log p x #T p (u) (k) < u= p 0.835u < 3x 0.835. On the other hand, if n = e 0,..., e t, d t+,..., d s p [(k )p, x] and d 0,..., d r p T p (k) for some integer r [t +, s], then by Theorem 2.(ii) we get that ν p (H(n, k)) = W p (k) + r ks W p (k) (k )s = t B p (e 0,..., e v )v (k )s t B p (e 0,..., e v )t (k )s < (k )(log p n log p (k ) ), where we have made use of (6) and the inequalities s > log p n and t log p (k ). 8. Figures 3, 3, 0 3, 2 3,, 0 3,, 2 3, 2, 0 3, 2, 2 3,, 2, 0 3,, 2, 2 3, 2, 0, 0 3, 2, 0, 2 3,, 2, 2, 3, 2, 0, 0, 3,, 2, 2,, 3 Figure. The 8 elements of T 3 (2) (left tree), and the 7 elements of T 3 (5) (right tree).
On the p-adic valuation of Stirling numbers 2 3 2, 0 3 2, 2 3 2, 0, 3 2, 2, 0 3 2, 2, 2 3 2, 0,, 3 2, 2, 0, 3 2, 2, 2, 0 3 2, 2, 2, 2 3 2, 2, 0,, 0 3 2, 2, 0,, 2 3 2, 2, 2, 0, 0 3 2, 2, 2, 0, 2 3 2, 2, 0,, 0, 3 2, 2, 0,, 2, 3 2, 2, 2, 0, 0, 3 2, 2, 0,, 0,, 0 3 2, 2, 0,, 0,, 2 3 2, 2, 0,, 2,, 3 2, 2, 0,, 0,, 0, 3 2, 2, 0,, 0,, 2, 3 2, 2, 0,, 2,,, 3 2, 2, 0,, 0,, 0,, 3 Figure 2. The 24 elements of T 3 (3)., 0 3, 0, 0 3, 0, 2 3, 0, 0, 0 3, 0, 0, 2 3, 0, 2, 3, 0, 0, 0, 0 3, 0, 0, 0, 2 3, 0, 0, 2, 3, 0, 0, 0, 0, 0 3, 0, 0, 0, 0, 2 3, 0, 0, 0, 2, 3, 0, 0, 2,, 3, 0, 0, 0, 0, 0, 3, 0, 0, 2,,, 3, 0, 0, 0, 0, 0,, 3 Figure 3. The 6 elements of T 3 (4). References [] T. Amdeberhan, D. Manna, and V. H. Moll, The 2-adic valuation of Stirling numbers, Experiment. Math. 7 (2008), no., 69 82. [2] H. Belbachir and A. Khelladi, On a sum involving powers of reciprocals of an arithmetical progression, Ann. Math. Inform. 34 (2007), 29 3. [3] D. W. Boyd, A p-adic study of the partial sums of the harmonic series, Experiment. Math. 3 (994), no. 4, 287 302. [4] Y.-G. Chen and M. Tang, On the elementary symmetric functions of, /2,..., /n, Amer. Math. Monthly 9 (202), no. 0, 862 867.
2 Paolo Leonetti and Carlo Sanna [5] H. Cohn, 2-adic behavior of numbers of domino tilings, Electron. J. Combin. 6 (999), Research Paper 4, 7 pp. (electronic). [6] S. De Wannemacker, On 2-adic orders of Stirling numbers of the second kind, Integers 5 (2005), no., A2, 7. [7] P. Erdős, Verallgemeinerung eines elementar-zahlentheoretischen Satzes von Kürschák., Mat. Fiz. Lapok 39 (932), 7 24 (Hungarian). [8] P. Erdős and I. Niven, Some properties of partial sums of the harmonic series, Bull. Amer. Math. Soc. 52 (946), 248 25. [9] A. Eswarathasan and E. Levine, p-integral harmonic sums, Discrete Math. 9 (99), no. 3, 249 257. [0] R. L. Graham, D. E. Knuth, and O. Patashnik, Concrete mathematics, second ed., Addison-Wesley Publishing Company, Reading, MA, 994, A foundation for computer science. [] S. Hong and C. Wang, The elementary symmetric functions of reciprocals of elements of arithmetic progressions, Acta Math. Hungar. 44 (204), no., 96 2. [2] S. Hong, J. Zhao, and W. Zhao, The 2-adic valuations of Stirling numbers of the second kind, Int. J. Number Theory 8 (202), no. 4, 057 066. [3] J. Kürschák, Über die harmonische Reihe., Mat. Fiz. Lapok 27 (98), 299 300 (Hungarian). [4] T. Lengyel, On the divisibility by 2 of the Stirling numbers of the second kind, Fibonacci Quart. 32 (994), no. 3, 94 20. [5], The order of the Fibonacci and Lucas numbers, Fibonacci Quart. 33 (995), no. 3, 234 239. [6], On the 2-adic order of Stirling numbers of the second kind and their differences, 2st International Conference on Formal Power Series and Algebraic Combinatorics (FPSAC 2009), Discrete Math. Theor. Comput. Sci. Proc., AK, Assoc. Discrete Math. Theor. Comput. Sci., Nancy, 2009, pp. 56 572. [7], Exact p-adic orders for differences of Motzkin numbers, Int. J. Number Theory 0 (204), no. 3, 653 667. [8], On p-adic properties of the Stirling numbers of the first kind, J. Number Theory 48 (205), 73 94. [9] D. Marques and T. Lengyel, The 2-adic order of the Tribonacci numbers and the equation T n = m!, J. Integer Seq. 7 (204), no. 0, Article 4.0., 8. [20] A. Postnikov and B. E. Sagan, What power of two divides a weighted Catalan number?, J. Combin. Theory Ser. A 4 (2007), no. 5, 970 977. [2] C. Sanna, On the p-adic valuation of harmonic numbers, J. Number Theory 66 (206), 4 46. [22], The p-adic valuation of Lucas sequences, Fibonacci Quart. 54 (206), 8 224. Università Luigi Bocconi, Department of Statistics, Milan, Italy E-mail address: leonetti.paolo@gmail.com URL: http://orcid.org/0000-000-789-530 Università degli Studi di Torino, Department of Mathematics, Turin, Italy E-mail address: carlo.sanna.dev@gmail.com URL: http://orcid.org/0000-0002-2-7596