OFFSET
0,1
COMMENTS
LINKS
FORMULA
a(n) = 5*Sum_{k = 0..n} binomial(n,k)^2*binomial(n+k,k)^2 - 2*Sum_{k = 0..n} binomial(n,k)^2*binomial(n+k,k).
a(n*p^r - 1) == a(n*p^(r-1) - 1) ( mod p^(3*r) ) for positive integers n and r and for all primes p >= 5.
a(n) = 5*hypergeom([-n, -n, 1 + n, 1 + n], [1, 1, 1], 1) - 2*hypergeom([1 + n, -n, -n], [1, 1], 1). - Peter Luschny, Nov 01 2022
MAPLE
seq(add(5*binomial(n, k)^2*binomial(n+k, k)^2 - 2*binomial(n, k)^2* binomial(n+k, k), k = 0..n), n = 0..20);
# Alternatively:
a := n -> 5*hypergeom([-n, -n, 1 + n, 1 + n], [1, 1, 1], 1) - 2*hypergeom([1 + n, -n, -n], [1, 1], 1): seq(simplify(a(n)), n = 0..18); # Peter Luschny, Nov 01 2022
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Peter Bala, Oct 24 2022
STATUS
approved