[seqfan] Re: funny recursion producing a known table

Paul D Hanna pauldhanna at juno.com
Sat Apr 27 02:15:21 CEST 2013


Wouter (and SeqFans), 
   For the Mma-less (such as myself), your triangle may be expressed also by: 
 
(PARI) {T(n,k)=if(n>=k,if(n==k,1,if((n-k+1)%2==0,0,if(k==0,sum(m=0, n,T(n\2,m)),T(n-1,k-1)+T(n-2*k,k)))))}  
for(n=0,20,for(k=0,n,print1(T(n,k),", "));print("")) 
 
Another formula to add to your triangle: 
 
A066897(n) = Sum_{k=0..n} k*T(n,k), where A066897(n) = total number of odd parts in all partitions of n. . 
 
I'm sure there are other such latent formulas. 
  Paul 
  
---------- Original Message ----------
From: "Wouter Meeussen" <wouter.meeussen at telenet.be>
To: "Sequence Fanatics Discussion list" <seqfan at list.seqfan.eu>
Subject: [seqfan] funny recursion producing a known table
Date: Fri, 26 Apr 2013 23:52:24 +0200

hi All,

try this one out:

Clear[T];
T[0, 0]:= 1;
T[n_, k_]:= 0 /; ( n<= 0|| OddQ[n +k] ); (* aerated triangular table *)
T[n_, k_]:=T[n, k]=  If[ k==0,   Sum[ T[ n/2, m],{m, 0, n}],  T[n -1, k -1]+ 
T[ n- 2*k ,k] ];

(* and then show it using *)  Table[T[n, k], {n, 0, 12}, {k, 0, n}]

The fun is in the surprising row sums and the familiar columns down to the 
diagonal :
why does such a silly  "T[n -1, k -1]+ T[ n- 2*k ,k] " produce this cute 
result (from Wolfdieter Lang)?

And what on earth does it have to do with its GF? (from Emeric Deutsch)
GF :  1/ Product[ (1-t*x^(2*j-1)) * (1-x^(2*j) )

hope the Mathematica code is readable to all,
enjoy,

Wouter. 


_______________________________________________

Seqfan Mailing list - http://list.seqfan.eu/



More information about the SeqFan mailing list