4306 days ago by dmanna

#sequence A(l,m) l,m = var('l,m') def A(l,m): sum=0 for k in range(l,m+1): sum = sum + factorial(l)*factorial(m)/2^(m-l)*2^k*binomial(2*(m-k),m-k)*binomial(m+k,m)*binomial(k,l) return sum 
# p-adic valuation of an integer n def nu(p,n): out=oo if n!=0: out=0 while n%p==0: out=out+1 n=n/p return out 
#creating a list-type sequence of the first N p-adic valuations of {A_{l,m} m>=l} for given l def Seq(p,l,N): list=[] for j in range(N): list.append(nu(p,A(l,l+j))) return list 
#creates a plot of the previous sequence. I have the options set to connect the dots and make the plot black. def PlotSeq(p,l,N): return list_plot(Seq(p,l,N),plotjoined=True,rgbcolor=(0,0,0)) 
#On the next three lines, some examples: PlotSeq(2,18,100)