login
Number of permutations p() of 1..n+6 with centered difference p(i+1)-p(i-1) < 0 exactly 5 times
1

%I #3 Mar 31 2012 12:35:48

%S 35,1540,45738,877968,13547688,174550992,1997701992,20732684544,

%T 200430861345,1825340594220,15885688179790,132967715234640,

%U 1079711852655436,8537605886212832,66114368122853232,502516350655667712

%N Number of permutations p() of 1..n+6 with centered difference p(i+1)-p(i-1) < 0 exactly 5 times

%C Column 6 of A180887

%H R. H. Hardin, <a href="/A180883/b180883.txt">Table of n, a(n) for n=1..100</a>

%o /* bc , formula r(n,k) from A000892 */

%o define factorial(n) {

%o auto prod;

%o prod=1;

%o while(n>=2)prod*=n--;

%o return prod;

%o }

%o define binomial(n,i) {

%o if(i<0||i>n)return 0;

%o return factorial(n)/(factorial(i)*factorial(n-i));

%o }

%o define r(n,k) {

%o auto j,sum;

%o sum=0;

%o for(j=0; j<=k+1; j++) {

%o sum+=(-1)^j*(k+1-j)^n*binomial(n+1,j);

%o }

%o return sum;

%o }

%o define t(n,k) {

%o auto sum,i;

%o sum=0;

%o for(i=0; i<=(k-1); i++)sum+=r((n+k)/2,i)*r((n+k)-(n+k)/2,(k-1)-i);

%o return sum*binomial((n+k),(n+k)/2);

%o }

%o for(n=1; n<=100; n++) {

%o print n, " ", t(n,6), "\n";

%o }

%o quit

%K nonn

%O 1,1

%A _R. H. Hardin_ Sep 23 2010