cp's OEIS Frontend

A234616 Numbers of undirected cycles in the complete tripartite graph K_{n,n,n}.

%I A234616 #37 Feb 16 2025 08:33:21
%S A234616 1,63,6705,1960804,1271288295,1541975757831,3135880743480163,
%T A234616 9904953891455450640,45915662047529291081589,
%U A234616 299038026557168514632822455,2642895689915240835222121682301,30814273315381549790551229559722628
%N A234616 Numbers of undirected cycles in the complete tripartite graph K_{n,n,n}.
%H A234616 Vaclav Kotesovec, <a href="/A234616/b234616.txt">Table of n, a(n) for n = 1..100</a> (terms 1..50 from Andrew Howroyd)
%H A234616 Andrew Howroyd, <a href="/A234616/a234616.txt">Formula for the number of cycles</a>
%H A234616 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/CompleteTripartiteGraph.html">Complete Tripartite Graph</a>
%H A234616 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/GraphCycle.html">Graph Cycle</a>
%F A234616 Row sums of A296546.
%F A234616 a(n) ~ sqrt(3*Pi) * 2^(3*n - 1/2) * n^(3*n - 1/2) / exp(3*n - 3/2). - _Vaclav Kotesovec_, Feb 17 2024
%t A234616 Table[(Sum[Binomial[n, k] Binomial[n, i + p] Binomial[n, j + p] Binomial[k, i] Binomial[k - i, j] (k - 1)! (i + p)! (j + p)! 2^(k - i - j) Binomial[p + i + j - 1, k - 1], {k, n}, {i, 0, k}, {j, 0, k - i}, {p, k - i - j, n}] + Sum[Binomial[n, k]^2 k! (k - 1)!, {k, 2, n}])/2 - n^2, {n, 10}] (* _Eric W. Weisstein_, May 26 2017 *)
%t A234616 Table[(n^2 (HypergeometricPFQ[{1, 1, 1 - n, 1 - n}, {2}, 1] - 3) + Sum[2^(k - i - j) Binomial[k, i] Binomial[k - i, j] Binomial[n, k] Binomial[n, i + p] Binomial[n, j + p] Binomial[i + j + p - 1, k - 1] (k - 1)! (i + p)! (j + p)!, {k, n}, {i, 0, k}, {j, 0, k - i}, {p, k - i - j, n}])/2, {n, 10}] (* _Eric W. Weisstein_, May 25 2023 *)
%o A234616 (PARI)
%o A234616 c(n,k,i,j,p) = {binomial(n,k)*binomial(n,i+p)*binomial(n,j+p)*binomial(k,i)*binomial(k-i,j)*(k-1)!*(i+p)!*(j+p)!*2^(k-i-j)*binomial(p+i+j-1,k-1)}
%o A234616 a(n)={(sum(k=1,n,sum(i=0,k,sum(j=0,k-i,sum(p=k-i-j,n, c(n,k,i,j,p) )))) + sum(k=2,n,binomial(n,k)^2*k!*(k-1)!))/2 - n^2} \\ _Andrew Howroyd_, May 25 2017
%o A234616 (Python)
%o A234616 from sympy import binomial, factorial
%o A234616 def c(n, k, i, j, p): return binomial(n, k)*binomial(n, i + p)*binomial(n, j + p)*binomial(k, i)*binomial(k - i, j)*factorial(k - 1)*factorial(i + p)*factorial(j + p)*2**(k - i - j)*binomial(p + i + j - 1, k - 1)
%o A234616 def a(n): return (sum([sum([sum([sum([c(n, k, i, j, p) for p in range(k - i - j, n + 1)]) for j in range(k - i + 1)]) for i in range(k + 1)]) for k in range(1, n + 1)]) + sum(binomial(n, k)**2*factorial(k)*factorial(k - 1) for k in range(2, n + 1)))/2 - n**2
%o A234616 print([a(k) for k in range(1, 13)]) # _Indranil Ghosh_, Aug 14 2017, after PARI
%Y A234616 Cf. A234365, A234633, A070968.
%Y A234616 Cf. A296546 (cycle polynomial coefficients of K_n,n,n).
%K A234616 nonn
%O A234616 1,2
%A A234616 _Eric W. Weisstein_, Dec 28 2013
%E A234616 a(7)-a(12) from _Andrew Howroyd_, May 25 2017