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.
%I A093965 #21 Sep 08 2022 08:45:13
%S A093965 1,4,21,124,825,6186,51961,484968,4988241,56117710,685883121,
%T A093965 9053657196,128397320233,1947359356866,31457343457065,539268744978256,
%U A093965 9778739908939041,187018400758459158,3762370179964296001,79427814910357360020,1755772750650004800441
%N A093965 Number of functions of [n] to [n] that simultaneously avoid the patterns 112 and 221.
%H A093965 Vincenzo Librandi, <a href="/A093965/b093965.txt">Table of n, a(n) for n = 1..200</a>
%F A093965 a(n) = A093966(n, n).
%F A093965 From _Vaclav Kotesovec_, Nov 20 2012: (Start)
%F A093965 E.g.f.: x*(exp(x) - x)/(1-x)^2.
%F A093965 Recurrence: (n-1)*a(n) = n*(n+1)*a(n-1) - (n-1)*n*a(n-2) for n>2.
%F A093965 a(n) ~ n!*n*(e-1). (End)
%t A093965 Rest[CoefficientList[Series[x(E^x-x)/(1-x)^2, {x, 0, 20}], x]* Range[0, 20]!] (* _Vaclav Kotesovec_, Nov 20 2012 *)
%o A093965 (PARI) my(x='x+O('x^66)); Vec(serlaplace(x*(exp(x)-x)/(1-x)^2)) \\ _Joerg Arndt_, May 11 2013
%o A093965 (Magma) [n le 2 select 4^(n-1) else n*((n+1)*Self(n-1) - (n-1)*Self(n-2))/(n-1): n in [1..30]]; // _G. C. Greubel_, Dec 29 2021
%o A093965 (Sage) [factorial(n)*( x*(exp(x) -x)/(1-x)^2 ).series(x,n+1).list()[n] for n in (1..30)] # _G. C. Greubel_, Dec 29 2021
%Y A093965 Cf. A093966.
%K A093965 nonn
%O A093965 1,2
%A A093965 _Ralf Stephan_, Apr 20 2004
%E A093965 Name changed by _Olivier Gérard_, Aug 06 2016