[seqfan] Re: More terms for A185074 (editing).

[Heinz, Alois]

I used a recursive Maple procedure with dynamic programming
to compute more terms:

b:= proc(r, i, n) option remember;
       `if`(r=0, 1, `if`(i>n, 0,
       add (b(r-j/i, i+1, n), j=0..min(n, r*i))))
     end:
a:= n-> b(n, 1, n):
seq (a(n), n=1..13);

1, 2, 4, 16, 36, 447, 1274, 9443, 54094, 995169,
3013040, 79403971, 244277081

Alois

Am 06.03.2012 19:50, schrieb Richard Mathar:
> http://list.seqfan.eu/pipermail/seqfan/2012-March/016529.html jl> 
> Joerg Arndt feels that A185074 <https://oeis.org/A185074> (Number of 
> jl> representations of n in the form sum(i=1..n, c(i)/i ), where each 
> of the jl> c(i)'s is in {0,1,...,n}.), submitted recently, needs more 
> terms before