%I #8 May 25 2016 22:01:11
%S 1,10,300,3003,20100,112101,2001000,33020001,200010000,3200120001,
%T 20000100000,320001200001,2000001000000,32000012000001,
%U 200000010000000,3200000120000001,20000000100000000,320000001200000001,2000000001000000000,32000000012000000001,200000000010000000000,3200000000120000000001,20000000000100000000000,320000000001200000000001,2000000000001000000000000,32000000000012000000000001,200000000000010000000000000,3200000000000120000000000001,20000000000000100000000000000,320000000000001200000000000001
%N Largest n-digit triangular number with minimum digit sum.
%C Does the obvious pattern continue? - _Vladeta Jovovic_, Apr 07 2002
%p F:= proc(d)
%p local s, P, nP, S, x, bestx;
%p bestx:= 0;
%p for s in [1,3,6,9] do
%p for P in map(op @combinat:-permute, combinat:-partition(s)) do
%p nP:= nops(P);
%p for S in map(t -> [d-1, op(t)], combinat:-choose([$0..d-2],nP-1)) do
%p x:= add(P[i]*10^S[i],i=1..nP);
%p if x > bestx and issqr(1+8*x) then bestx:= x fi;
%p od;
%p od;
%p if bestx > 0 then return bestx fi;
%p od;
%p end proc:
%p seq(F(d),d=1..30); # _Robert Israel_, May 25 2016
%Y Cf. A069661, A069662, A069663, A069664, A069665, A069666, A069667, A069668, A069669, A069670, A069671.
%K nonn,base
%O 1,2
%A _Amarnath Murthy_, Apr 06 2002
%E More terms from _Vladeta Jovovic_, Apr 07 2002
%E a(15)..a(30) from _Robert Israel_, May 25 2016
|