We use the notation (2n−1)!!=1⋅3⋯(2n−1) and (2n)!!=2⋅4⋯(2n)=2nn! for any positive integer n. Observe that (2n)!=(2n)!!(2n−1)!!=2nn!(2n−1)!!.
For any positive integer n we have
(2n4n)=(2n)!2(4n)!=(2n)!222n(2n)!(4n−1)!!=(2n)!22n(4n−1)!!(n2n)=(2n)!1(n!(2n)!)2=(2n)!1(2n(2n−1)!!)2=(2n)!22n(2n−1)!!2
Then expression (1) can be rewritten as follows:
(2k2k+1)−(2k−12k)=(2k)!22k(2k+1−1)!!−(2k)!22k(2k−1)!!2=(2k)!22k(2k−1)!!⋅((2k+1)(2k+3)…(2k+2k−1)−(2k−1)(2k−3)…(2k−2k+1))
We compute the exponent of 2 in the prime decomposition of each factor (the first one is a rational number but not necessarily an integer; it is not important).
First, we show by induction on n that the exponent of 2 in (2n)! is 2n−1. The base case n=1 is trivial. Suppose that (2n)!=22n−1(2d+1) for some integer d. Then we have
(2n+1)!=22n(2n)!(2n+1−1)!!=22n22n−1⋅(2d+1)(2n+1−1)!!=22n+1−1⋅(2q+1)
for some integer q. This finishes the induction step.
Hence, the exponent of 2 in the first factor in (2) is 2k−(2k−1)=1.
The second factor in (2) can be considered as the value of the polynomial
P(x)=(x+1)(x+3)…(x+2k−1)−(x−1)(x−3)…(x−2k+1)
at x=2k. Now we collect some information about P(x).
Observe that P(−x)=−P(x), since k≥2. So P(x) is an odd function, and it has nonzero coefficients only at odd powers of x. Hence P(x)=x3Q(x)+cx, where Q(x) is a polynomial with integer coefficients.
Compute the exponent of 2 in c. We have
c=2(2k−1)!!i=1∑2k−12i−11=(2k−1)!!i=1∑2k−1(2i−11+2k−2i+11)=(2k−1)!!i=1∑2k−1(2i−1)(2k−2i+1)2k=2ki=1∑2k−1(2i−1)(2k−2i+1)(2k−1)!!=2kS
For any integer i=1,…,2k−1, denote by a2i−1 the residue inverse to 2i−1 modulo 2k. Clearly, when 2i−1 runs through all odd residues, so does a2i−1, hence
S=i=1∑2k−1(2i−1)(2k−2i+1)(2k−1)!!≡−i=1∑2k−1(2i−1)2(2k−1)!!≡−i=1∑2k−1(2k−1)!!a2i−12=−(2k−1)!!i=1∑2k−1(2i−1)2=−(2k−1)!!32k−1(22k−1)(mod2k)
Therefore, the exponent of 2 in S is k−1, so c=2kS=22k−1(2t+1) for some integer t.
Finally we obtain that
P(2k)=23kQ(2k)+2kc=23kQ(2k)+23k−1(2t+1)
which is divisible exactly by 23k−1. Thus, the exponent of 2 in (2) is 1+(3k−1)=3k.