Re^2: Error-finding program
Richard Mathar
mathar at strw.leidenuniv.nl
Thu Oct 26 12:52:55 CEST 2006
> From seqfan-owner at ext.jussieu.fr Wed Oct 25 21:07:37 2006
> Return-Path: <seqfan-owner at ext.jussieu.fr>
> To: seqfan at ext.jussieu.fr
> Subject: Re: Error-finding program
> ..
> A026302 Should certainly be monotonic; I don't know what the correct
> value of a(18) is.
The fourth digit of a(18) was missing:
0 1
1 2
2 9
3 44
4 230
5 1242
6 6853
7 38376
8 217242
9 1239980
10 7123765
11 41141916
12 238637282
13 1389206210
14 8112107475
15 47495492400
16 278722764954
17 1638970147188
18 9654874654438
19 56965811111240
20 336590781348276
21 1991357644501170
in PARI:
A026300(n,k)={
if(n<0 || k < 0,
return(0) ;
) ;
if(n<=1,
1,
if(k==0,
1,
sum(i=0,k/2,
binomial(n,2*i+n-k)*(binomial(2*i+n-k,i)-binomial(2*i+n-k,i-1))
) ;
) ;
) ;
}
A026302(n)={
A026300(2*n,n) ;
}
{
for(n=0,21,
print(n," ",A026302(n))
) ;
}
More information about the SeqFan
mailing list