login
Row sums of A197653.
8

%I #87 Jul 14 2024 10:10:25

%S 1,5,46,485,5626,69062,882540,11614437,156343330,2142556130,

%T 29791689148,419260001030,5960334608788,85469709312860,

%U 1234797737654296,17955907741675749,262607675818816050,3860239468267647914,57002176852356800700,845159480056345448610

%N Row sums of A197653.

%C Number of meanders of length (n+1)*4 which are composed by arcs of equal length and a central angle of 90 degrees.

%C Definition of a meander:

%C A binary curve C is a triple (m, S, dir) such that

%C (a) S is a list with values in {L,R} which starts with an L,

%C (b) dir is a list of m different values, each value of S being allocated a value of dir,

%C (c) consecutive Ls increment the index of dir,

%C (d) consecutive Rs decrement the index of dir,

%C (e) the integer m > 0 divides the length of S and

%C (f) C is a meander if each value of dir occurs length(S)/m times.

%C For this sequence, m = 4.

%C The terms are proved by brute force for 0 <= n <= 8, but not yet in general. - _Susanne Wienand_, Oct 29 2011

%H Peter Luschny, <a href="http://oeis.org/wiki/User:Peter_Luschny/Meander">Meanders and walks on the circle</a>.

%F a(n) = Sum_{k=0..n} Sum_{j=0..3} Sum_{i=0..3} (-1)^(j+i)*C(i,j)*C(n,k)^4*(n+1)^j*(k+1)^(3-j)/(k+1)^3. - _Peter Luschny_, Nov 02 2011

%F a(n) = Sum_{k=0..n} h(n,k)*binomial(n,k)^4, where h(n,k) = (1+k)*(1-((n-k)/(1+k))^4)/(1+2*k-n) if 1+2*k-n <> 0 else h(n,k) = 4. - _Peter Luschny_, Nov 24 2011

%F From _Peter Bala_, Mar 21 2023: (Start)

%F Conjecture 1: a(n) = Sum_{k = 0..n} binomial(n+1,k)^2*binomial(n,k)^2.

%F If true, then we have the third-order recurrence equation

%F n^2*(n + 1)^3*P(n-1)*a(n) = 2*n^2*(400*n^8 - 1260*n^7 + 20*n^6 + 3020*n^5 - 1646*n^4 - 1951*n^3 + 1142*n^2 + 465*n - 290)*a(n-1) + 4*(n - 1)*(2800*n^9 - 15420*n^8 + 30620*n^7 - 23710*n^6 + 808*n^5 + 6863*n^4 - 1309*n^3 - 1218*n^2 + 496*n - 60)*a(n-2) + 8*(n - 2)^2*(2*n - 3)*(4*n - 5)*(4*n - 7)*P(n)*a(n-3) with a(0) = 1, a(1) = 5 and a(2) = 46 and where P(n) = 100*n^5 - 65*n^4 - 35*n^3 + 25*n^2 + 6*n - 5.

%F Conjecture 2: working with offset 1, that is, a(1) = 1, a(2) = 5, ..., then the supercongruence a(n*p^r) == a(n*p^(r-1)) (mod p^(3*r)) holds for positive integers n and r and all primes p >= 5. (End)

%F a(n) ~ 2^(4*n + 5/2) / (Pi*n)^(3/2). - _Vaclav Kotesovec_, Apr 17 2023

%F Peter Bala's conjecture 1 can equivalently written: a(n) = hypergeom([-n - 1, -n - 1, -n, -n], [1, 1, 1], 1). - _Detlef Meya_, May 28 2024

%F a(n) = Sum_{k=0..n+1} (k/(n+1))^2 * binomial(n+1,k)^4. - _Seiichi Manyama_, Jul 14 2024

%e Some examples of list S and allocated values of dir if n = 4:

%e Length(S) = (4+1)*4 = 20.

%e S: L,L,L,L,L,L,L,L,L,L,L,L,L,L,L,L,L,L,L,L

%e dir: 1,2,3,0,1,2,3,0,1,2,3,0,1,2,3,0,1,2,3,0

%e S: L,L,L,L,R,L,R,R,L,R,R,R,L,R,L,L,R,L,L,L

%e dir: 1,2,3,0,0,0,0,3,3,3,2,1,1,1,1,2,2,2,3,0

%e S: L,L,R,L,L,L,L,R,R,L,R,R,L,R,L,L,L,L,R,R

%e dir: 1,2,2,2,3,0,1,1,0,0,0,3,3,3,3,0,1,2,2,1

%e Each value of dir occurs 20/4 = 5 times.

%t A198256[n_] := Sum[Sum[Sum[(-1)^(j + i)* Binomial[i, j]*Binomial[n, k]^4*(n + 1)^j*(k + 1)^(3 - j)/(k + 1)^3, {i, 0, 3}], {j, 0, 3}], {k, 0, n}]; Table[A198256[n], {n, 0, 16}] (* _Peter Luschny_, Nov 02 2011 *)

%o (Sage)

%o from mpmath import mp, hyper

%o def A198256(n) : return hyper([1-n, 1-n, 1-n, 1-n], [3, 3, 3], 1)*(n^4-n^6)/4 + hyper([-n, -n, -n, -n], [2, 2, 2], 1)*(1+n+n^2+n^3) + hyper([2, 1-n, 1-n, 1-n, 1-n], [1, 3, 3, 3], 1)*(n^4+n^5)/4

%o mp.dps = 32

%o for n in (0..19) : print(int(A198256(n))) # _Peter Luschny_, Oct 24 2011

%o (PARI)

%o A198256(n) = {sum(k=0,n,if(n == 1+2*k,4,(1+k)*(1-((n-k)/(1+k))^4)/(1+2*k-n))*binomial(n,k)^4)} \\ _Peter Luschny_, Nov 24 2011

%Y Cf. A197653, A198060, A198257, A198258.

%Y Cf. A005260, A181069.

%Y Cf. A181067, A374614.

%K nonn

%O 0,2

%A _Susanne Wienand_, Oct 22 2011