cp's OEIS Frontend

A218827 Number of indecomposable (by shuffling) alternating n-anagrams.

%I A218827 #12 Aug 03 2019 14:20:03
%S A218827 1,1,3,16,129,1438,20955,384226,8623101,231978454,7359117591,
%T A218827 271673905642,11543742745689,559348370431630,30659822500574739,
%U A218827 1887796293833267746,129757032076160998677,9900820197631733600518,834421415151529202479023,77318409826165250051727514
%N A218827 Number of indecomposable (by shuffling) alternating n-anagrams.
%C A218827 Alternating anagrams enumeration is related to A000366 by a(n) = A000366(n+1).
%C A218827 For n>2, a(n) are alternatively congruent to 3 and -2 modulo 18.
%H A218827 G. Kreweras and J. Barraud, <a href="http://dx.doi.org/10.1006/eujc.1997.0161">Anagrammes alternés</a>, European Journal of Combinatorics,Volume 18, Issue 8, November 1997, Pages 887-891.
%H A218827 G. Kreweras and D. Dumont, <a href="http://dx.doi.org/10.1016/S0012-365X(99)00238-1">Sur les anagrammes alternés</a>, Discrete Mathematics, Volume 211, Issues 1-3, 28 January 2000, Pages 103-110.
%t A218827 m = 20(*terms*); matc = Array[0&, {m, m}];
%t A218827 a366[n_] := (-2^(-1))^(n - 2)*Sum[Binomial[n, k]*(1 - 2^(n + k + 1))* BernoulliB[n + k + 1], {k, 0, n}];
%t A218827 ci[n_, k_] := ci[n, k] = Module[{v}, If[matc[[n, k]] == 0, If[n == k, v = 1, If[k == 1, v = c[n], v = Sum[Binomial[n - 1, i - 1]*c[i]*ci[n - i, k - 1], {i, 1, n - k + 1}]]]; matc[[n, k]] = v]; Return[matc[[n, k]]]];
%t A218827 a[n_] := a366[n + 1] - If[n == 1, 0, Sum[ci[n, i], {i, 2, n}]];
%t A218827 Array[a, m] (* _Jean-François Alcover_, Aug 03 2019, from PARI *)
%o A218827 (PARI) a000366(n)= {return((-1/2)^(n-2)*sum(k=0, n, binomial(n, k)*(1-2^(n+k+1))*bernfrac(n+k+1)));}
%o A218827 ci(n,k) = {if (matc[n,k] == 0, if (n==k, v = 1, if (k==1, v = c(n), v = sum(i=1, n-k+1, binomial(n-1,i-1)*c(i)*ci(n-i,k-1)););); matc[n,k] = v;); return(matc[n,k]);}
%o A218827 c(n) = {if (n==1, return(a000366(n+1)), return(a000366(n+1) - sum(i=2, n, ci(n,i))));}
%o A218827 allc(m) = {matc = matrix(m,m); for (i=1, m, print1(c(i), ", "););}
%Y A218827 Cf. A000366, A218826. First column of A239895.
%K A218827 nonn
%O A218827 1,3
%A A218827 _Michel Marcus_, Nov 07 2012