%I #27 Aug 01 2023 12:13:59
%S 1,1,1,1,1,1,1,6,27,83,209,461,923,10284,80991,414961,1671853,5699149,
%T 17116009,278723178,3135810159,22493048843,124606826189,574688719793,
%U 2301250545971,49308397822776,721175428306971,6650954153090521,46893517738791361
%N Number of permutations p of [n] such that p(i) > p(i+1) iff i == 0 (mod 6).
%H Alois P. Heinz, <a href="/A250283/b250283.txt">Table of n, a(n) for n = 0..500</a>
%H Ludwig Seidel, <a href="https://babel.hathitrust.org/cgi/pt?id=hvd.32044092897461&view=1up&seq=176">Über eine einfache Entstehungsweise der Bernoulli'schen Zahlen und einiger verwandten Reihen</a>, Sitzungsberichte der mathematisch-physikalischen Classe der königlich bayerischen Akademie der Wissenschaften zu München, volume 7 (1877), 157-187. [USA access only through the <a href="https://www.hathitrust.org/accessibility">HATHI TRUST Digital Library</a>]
%H Ludwig Seidel, <a href="https://www.zobodat.at/pdf/Sitz-Ber-Akad-Muenchen-math-Kl_1877_0157-0187.pdf">Über eine einfache Entstehungsweise der Bernoulli'schen Zahlen und einiger verwandten Reihen</a>, Sitzungsberichte der mathematisch-physikalischen Classe der königlich bayerischen Akademie der Wissenschaften zu München, volume 7 (1877), 157-187. [Access through <a href="https://de.wikipedia.org/wiki/ZOBODAT">ZOBODAT</a>]
%e a(6) = 1: 123456.
%e a(7) = 6: 1234576, 1234675, 1235674, 1245673, 1345672, 2345671.
%p b:= proc(u, o, t) option remember; `if`(u+o=0, 1,
%p `if`(t=0, add(b(u-j, o+j-1, irem(t+1, 6)), j=1..u),
%p add(b(u+j-1, o-j, irem(t+1, 6)), j=1..o)))
%p end:
%p a:= n-> b(n, 0$2):
%p seq(a(n), n=0..35);
%t nmax = 30; CoefficientList[Series[1 + Sum[(x^(6 - k) * HypergeometricPFQ[{1}, {7/6 - k/6, 4/3 - k/6, 3/2 - k/6, 5/3 - k/6, 11/6 - k/6, 2 - k/6}, -x^6/46656])/(6 - k)!, {k, 0, 5}] / HypergeometricPFQ[{}, {1/6, 1/3, 1/2, 2/3, 5/6}, -x^6/46656], {x, 0, nmax}], x] * Range[0, nmax]! (* _Vaclav Kotesovec_, Apr 21 2021 *)
%o (Sage)
%o # From _Peter Luschny_, Feb 06 2017 (Start)
%o @cached_function
%o def b(u, o, t):
%o if u ==-o: return 1
%o if t == 0: return sum(b(u-j, o+j-1, (t+1) % 6) for j in (1..u))
%o return sum(b(u+j-1, o-j, (t+1) % 6) for j in (1..o))
%o a = lambda n: b(n, 0, 0)
%o print([a(n) for n in (0..28)]) # after Maple program
%o # Alternatively:
%o @cached_function
%o def A(m, n):
%o if n == 0: return 1
%o s = -1 if m.divides(n) else 1
%o t = [m*k for k in (0..(n-1)//m)]
%o return s*add(binomial(n, k)*A(m, k) for k in t)
%o A250283 = lambda n: (-1)^int(is_odd(n//6))*A(6, n)
%o print([A250283(n) for n in (0..28)])
%o # (End)
%Y Row n=6 of A181937.
%K nonn
%O 0,8
%A _Alois P. Heinz_, Nov 16 2014