|
%I #13 Apr 09 2014 11:59:59
%S 1,3,7,22,43,170,321,1127,2816,9221,16403,90121,157741,475003,1377765,
%T 5194616,8925845,41509031,70873504,314321071,846062365,2349224135,
%U 3973214615,25816158418,53403343055,145093944173,373986049879
%N a(n) = the numerator of the continued fraction [[n/1];[n/2],[n/3],..,[n/n]] = the numerator of [[n/n];[n/(n-1)],[n/(n-2)],..,[n/1]], where [x] is floor(x).
%p A128599 := proc(n) local a,i ; a := 1 ; for i from n-1 to 1 by -1 do a := floor(n/i)+1/a ; od ; RETURN(numer(a)) ; end: for n from 1 to 40 do printf("%d, ",A128599(n)) ; od ; # _R. J. Mathar_, Jun 02 2007
%t Table[Numerator[FromContinuedFraction[Table[Floor[n/i],{i,1,n}]]], {n, 1, 30}] (* _Stefan Steinerberger_, Jun 02 2007 *)
%Y Cf. A128600, A128601.
%K frac,nonn
%O 1,2
%A _Leroy Quet_, Mar 12 2007
%E More terms from _Stefan Steinerberger_ and _R. J. Mathar_, Jun 02 2007
|
