|
%I #7 Oct 27 2014 19:55:30
%S 1,2,1,6,2,1,6,6,2,1,30,3,6,2,1,5,15,24,6,2,1,35,90,120,24,6,2,1,70,
%T 630,45,120,24,6,2,1,70,630,315,720,120,24,6,2,1,7,210,2520,5040,720,
%U 120,24,6,2,1,77,2100,280,1260,5040,720,120,24,6,2,1,231
%N Rectangular array, by antidiagonals: T(m,n) = greatest (m+1)-th-power-free divisor of n!.
%C Row 1: A055204, greatest squarefree divisor of n!
%C Row 2: A145642, greatest cubefree divisor of n!
%C Row 3: A248766, greatest 4th-power-free divisor of n!
%C Rows 4 to 7: A248769, A248772, A248775, A248778.
%C (The divisors are here called "greatest" rather than "largest" because the name refers to ">", called "greater than".)
%H Clark Kimberling, <a href="/A248779/b248779.txt">Antidiagonals n = 1..60, flattened</a>
%e Northwest corner:
%e 1 2 6 6 30 5 35 70
%e 1 2 6 3 15 90 630 630
%e 1 2 6 24 120 45 315 2520
%e 1 2 6 24 120 720 5040 1260
%t f[n_] := f[n] = FactorInteger[n!]; r[m_, x_] := r[m, x] = m*Floor[x/m];
%t u[n_] := Table[f[n][[i, 1]], {i, 1, Length[f[n]]}];
%t v[n_] := Table[f[n][[i, 2]], {i, 1, Length[f[n]]}];
%t p[m_, n_] := p[m, n] = Product[u[n][[i]]^r[m, v[n]][[i]], {i, 1, Length[f[n]]}]
%t t = Table[n!/p[m, n], {m, 2, 16}, {n, 1, 16}]; TableForm[t] (* A248779 array *)
%t f = Table[t[[n - k + 1, k]], {n, 12}, {k, n, 1, -1}] // Flatten (* A248779 seq. *)
%Y Cf. A055204, A145642, A248766, A248769, A248772, A248775, A248778, A000142.
%K nonn,tabl,easy
%O 1,2
%A _Clark Kimberling_, Oct 14 2014
|
