|
%I #24 Mar 13 2014 10:11:19
%S 1,1,2,1,2,6,1,2,4,12,1,2,4,12,60,1,2,4,6,12,60,1,2,4,6,12,60,420,1,2,
%T 4,6,12,24,120,840,1,2,4,6,12,24,72,360,2520,1,2,4,6,12,24,60,120,360,
%U 2520,1,2,4,6,12,24,60,120,360,2520,27720,1,2,4,6,12,12,24,60,120,360
%N Triangle read by rows: T(n,k) is the smallest positive integer having at least k of the first n positive integers as divisors.
%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Table_of_divisors">Table of divisors</a>.
%F T(n,k) = min { lcm(x_1,...,x_k) ; 0 < x_1 < ... < x_k <= n }
%e Triangle begins:
%e 1
%e 1 2
%e 1 2 6
%e 1 2 4 12
%e 1 2 4 12 60
%e 1 2 4 6 12 60
%p with(combstruct):
%p a096179_row := proc(n) local k,L,l,R,LCM,comb;
%p R := NULL; LCM := ilcm(seq(i,i=[$1..n]));
%p for k from 1 to n-1 do
%p L := LCM;
%p comb := iterstructs(Combination(n),size=k):
%p while not finished(comb) do
%p l := nextstruct(comb);
%p L := min(L,ilcm(op(l)));
%p od;
%p R := R,L;
%p od;
%p R,LCM end; # _Peter Luschny_, Dec 06 2010
%t (* Triangular *)
%t A096179[n_,k_]:=Min[LCM@@@Subsets[Range[n],{k}]];
%t A002024[n_]:=Floor[1/2+Sqrt[2*n]];
%t A002260[n_]:=n-Binomial[Floor[1/2+Sqrt[2*n]],2];
%t (* Linear *)
%t A096179[n_]:=A096179[n]=A096179[A002024[n],A002260[n]];
%t (* Enrique Pérez Herrero_, Dec 08 2010 *)
%o (PARI) A096179(n,k)={ my(m=lcm(vector(k,i,i))); forvec(v=vector(k-1,i,[2,n]), m>lcm(v) & m=lcm(v), 2); m } \\ _M. F. Hasler_, Nov 30 2010
%Y Main diagonal is A003418. Minimum in column k is A061799(k). See also A094348, A096180.
%K nonn,tabl
%O 1,3
%A _Matthew Vandermast_, Jun 19 2004
|
