cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A351759 Number of words either empty or beginning with the first letter of the n-ary alphabet, where each letter of the alphabet occurs n times and letters of neighboring word positions are equal or neighbors in the alphabet.

Original entry on oeis.org

1, 1, 3, 37, 3593, 2336376, 12085125703, 492300320300251, 163856834420168837331, 450436951076447377275760495, 10360995769563446558192576106611177, 2010994252878923176109086647579398800496256, 3317129856969862808949985510559979253492213624483355
Offset: 0

Views

Author

Alois P. Heinz, Feb 18 2022

Keywords

Comments

Also the number of (n^2-1)-step walks on n-dimensional cubic lattice from (1,0,...,0) to (n,n,...,n) with positive unit steps in all dimensions such that the absolute difference of the dimension indices used in consecutive steps is <= 1.

Examples

			a(2) = 3:
  +----+   +----+   +----+
  |aabb|   |abab|   |abba|
  +----+   +----+   +----+
  |1222|   |1122|   |1112|
  |0012|   |0112|   |0122|
  +----+   +----+   +----+
  |xx  |   |x x |   |x  x|
  |  xx|   | x x|   | xx |
  +----+   +----+   +----+
a(3) = 37: aaabbbccc, aaabbcbcc, aaabbccbc, aaabbcccb, aaabcbbcc, aaabcbcbc, aaabcbccb, aaabccbbc, aaabccbcb, aaabcccbb, aababbccc, aababcbcc, aababccbc, aababcccb, aabbabccc, aabbcccba, aabcbabcc, aabcbccba, aabccbabc, aabccbcba, aabcccbab, aabcccbba, abaabbccc, abaabcbcc, abaabccbc, abaabcccb, abababccc, ababcccba, abbaabccc, abbcccbaa, abcbaabcc, abcbccbaa, abccbaabc, abccbcbaa, abcccbaab, abcccbaba, abcccbbaa.
		

Crossrefs

Main diagonal of A208673.

Programs

  • Maple
    b:= proc(l, t) option remember; (n-> `if`(l=[0$n], 1,
          add(`if`(l[i]=0, 0, b(subsop(i=l[i]-1, l), i)),
               i=max(1, t-1)..min(n, t+1))))(nops(l))
        end:
    a:= n-> b([n$n], 0):
    seq(a(n), n=0..6);
  • Mathematica
    b[l_, t_] := b[l, t] = Function [n, If[l == Array[0&, n], 1,
         Sum[If[l[[i]] == 0, 0, b[ReplacePart[l, i -> l[[i]] - 1], i]],
         {i, Max[1, t - 1], Min[n, t + 1]}]]][Length[l]];
    a[n_] := b[Array[n&, n], 0];
    Table[Print[n, " ", a[n]]; a[n], {n, 0, 7}] (* Jean-François Alcover, May 16 2022, after Alois P. Heinz *)

Formula

a(n) = A208673(n,n).