cp's OEIS Frontend

A287231 Number of matchings in the n-triangular graph.

%I A287231 #30 Feb 16 2025 08:33:46
%S A287231 1,4,51,2460,513619,509709696,2590569730617,71972142178289680,
%T A287231 11572569464349559854105,11332749125368045400133079296,
%U A287231 70775590368575601248957366910425851,2939823814188321813975498471683171002746816,844162736935477006294039214093750952242356035727995,1736712038520659436678773853448507425382701807453031820800000
%N A287231 Number of matchings in the n-triangular graph.
%H A287231 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/IndependentEdgeSet.html">Independent Edge Set</a>
%H A287231 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/JohnsonGraph.html">Johnson Graph</a>
%H A287231 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/Matching.html">Matching</a>
%H A287231 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/TriangularGraph.html">Triangular Graph</a>
%o A287231 (PARI)
%o A287231 \\ groups all labeled oriented graphs on n vertices by out degree configuration.
%o A287231 OrientedByOutDegrees(n)={ \\ high memory usage and slow for n > 10
%o A287231 local(M=Map());
%o A287231 my(acc(p,v)=my(z); mapput(M, p, if(mapisdefined(M, p, &z), z+v, v)));
%o A287231 my(recurse(p,i,q,v,e)=if(i<0, for(k=0,e,acc(x^k+q, binomial(e,k)*v)), my(t=polcoeff(p,i));for(k=0,t,self()(p,i-1,(t-k+x*k)*x^i+q,binomial(t,k)*v,e+t-k))));
%o A287231 my(iterate(v,k,f)=for(i=1,k,v=f(v));v);
%o A287231 iterate(Mat([1,1]), n-1, src->M=Map();for(i=1, matsize(src)[1], my(p=src[i,1]); recurse(p,poldegree(p),0,src[i,2],0)); Mat(M))
%o A287231 }
%o A287231 a(n)={
%o A287231 my(v=vector(n\2,n,(2*n)!/(2^n*n!)));
%o A287231 my(c(p)=my(h=(poldegree(p)+1)\2); my(r=n-1-sum(i=1,h,polcoeff(p,2*i-1))); (1+sum(i=1,r\2,binomial(r,2*i)*v[i]))*prod(i=1,h,v[i]^(polcoeff(p,2*i)+polcoeff(p,2*i-1))));
%o A287231 my(M=OrientedByOutDegrees(n-1));
%o A287231 sum(i=1,matsize(M)[1],M[i,2]*c(M[i,1]))
%o A287231 } \\ _Andrew Howroyd_, Aug 25 2017
%K A287231 nonn
%O A287231 2,2
%A A287231 _Eric W. Weisstein_, May 22 2017
%E A287231 a(9)-a(12) from _Andrew Howroyd_, Aug 25 2017
%E A287231 a(13)-a(14) from _Eric W. Weisstein_, Oct 01 2017
%E A287231 a(15) from _Eric W. Weisstein_, Oct 15 2017