|
%I #26 Oct 30 2022 18:19:59
%S 1,0,0,0,0,0,0,0,1,16,436,6984,114124,1799688,29125117,486980182,
%T 8490078104,154750897552,2951968964768,58917663227568,
%U 1229367602071416,26787823838035750,608794318333169289,14411810690642972432
%N Number of permutations s of {1,2,...,n} such that |s(i)-i|>3 for each i=1,2,...,n.
%C a(n) equals the permanent of the n X n matrix with 0's along the main diagonal, the subdiagonal, the superdiagonal, the sub-subdiagonal, the super-superdiagonal, the sub-sub-subdiagonal, the super-super-superdiagonal, and 1's everywhere else. - _John M. Campbell_, Jul 09 2011
%H Alois P. Heinz, <a href="/A075852/b075852.txt">Table of n, a(n) for n = 0..450</a>
%H George Spahn and Doron Zeilberger, <a href="https://arxiv.org/abs/2108.11285">Automatic Counting of Generalized Latin Rectangles and Trapezoids</a>, arXiv:2108.11285 [math.CO], 2021.
%p b:= proc(s) option remember; (n-> `if`(n=0, 1, add(
%p `if`(abs(n-i)>3, b(s minus {i}), 0), i=s)))(nops(s))
%p end:
%p a:= n-> b({$1..n}):
%p seq(a(n), n=0..15); # _Alois P. Heinz_, Jan 25 2019
%t a[0] = 1; a[n_] := a[n] = If[n < 8, 0, SparseArray[{Band[{1, 1}] -> 0, Band[{2, 1}] -> 0, Band[{3, 1}] -> 0, Band[{4, 1}] -> 0, Band[{1, 2}] -> 0, Band[{1, 3}] -> 0, Band[{1, 4}] -> 0}, {n, n}, 1] // Permanent];
%t Table[Print[n, " ", a[n]]; a[n], {n, 0, 23}] (* _Jean-François Alcover_, May 01 2019 *)
%Y Cf. A001883, A075851.
%K nonn
%O 0,10
%A _Reiner Martin_, Oct 15 2002
%E More terms from _Vladimir Baltic_, _Vladeta Jovovic_, Jan 04 2003
%E a(21) from _Alois P. Heinz_, Jul 04 2015
%E a(22)-a(23) from _Alois P. Heinz_, Jan 22 2019
%E a(0)=1 prepended by _Alois P. Heinz_, Jan 25 2019
|
