OFFSET
0,3
COMMENTS
Conjecture 1: the supercongruence a(p) == a(1) (mod p^5) holds for all primes p >= 7 (checked up to p = 199).
Conjecture 2: for r >= 2, the supercongruence a(p^r) == a(p^(r-1)) (mod p^(3*r+3) holds for all primes p >= 5.
LINKS
Paolo Xausa, Table of n, a(n) for n = 0..650
Peter Bala, Recurrence equation for A361712
FORMULA
a(n) ~ 2^(1/4)*(1 + sqrt(2))^(4*n+1)/(4*Pi^(3/2)*n^(3/2)).
a(n) = hypergeom([-n, -n, n, n + 1], [1, 1, 1], 1) - binomial(2*n, n)*binomial(2*n - 1, n) = A361878(n) - A361877(n). - Peter Luschny, Mar 27 2023
EXAMPLE
a(7) - a(1) = (2^2)*(7^5)*5009 == 0 (mod 7^5)
a(11) - a(1) = (2^5)*(11^5)*45864163 == 0 (mod 11^5)
a(7^2) - a(7) = (2*3)*(7^9)*377052719*240136524699189343838527* 17965610580703155723668147409587 == 0 (mod 7^9)
MAPLE
seq(add(binomial(n, k)^2*binomial(n+k, k)*binomial(n+k-1, k), k = 0..n-1), n = 0..25);
# Alternative:
A361712 := n -> hypergeom([-n, -n, n, n + 1], [1, 1, 1], 1) - binomial(2*n, n)*binomial(2*n-1, n): seq(simplify(A361712(n)), n = 0..18); # Peter Luschny, Mar 27 2023
MATHEMATICA
A361712[n_] := HypergeometricPFQ[{-n, -n, n, n+1}, {1, 1, 1}, 1] - Binomial[2*n, n]*Binomial[2*n-1, n]; Array[A361712, 20, 0] (* Paolo Xausa, Jul 10 2024 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Peter Bala, Mar 21 2023
STATUS
approved