A090619 -id:A090619 - cp's OEIS frontend

A090622 Square array read by antidiagonals of highest power of k dividing n! (with n,k>1).

1, 0, 1, 0, 1, 3, 0, 0, 1, 3, 0, 0, 1, 1, 4, 0, 1, 0, 1, 2, 4, 0, 0, 1, 1, 2, 2, 7, 0, 0, 0, 1, 1, 2, 2, 7, 0, 0, 1, 0, 2, 1, 3, 4, 8, 0, 0, 0, 1, 0, 2, 1, 3, 4, 8, 0, 0, 0, 0, 1, 1, 2, 1, 4, 4, 10, 0, 0, 0, 1, 1, 1, 1, 4, 2, 4, 5, 10, 0, 0, 1, 0, 1, 1, 2, 1, 4, 2, 5, 5, 11, 0, 0, 0, 1, 0, 1, 1, 2, 1, 4, 2, 5, 5, 11

Offset: 2

Henry Bottomley, Dec 06 2003

			Square array starts:
1, 0, 0, 0, 0, 0, 0, ...
1, 1, 0, 0, 1, 0, 0, ...
3, 1, 1, 0, 1, 0, 1, ...
3, 1, 1, 1, 1, 0, 1, ...
4, 2, 2, 1, 2, 0, 1, ...
4, 2, 2, 1, 2, 1, 1, ...
7, 2, 3, 1, 2, 1, 2, ...

Alois P. Heinz, Antidiagonals n = 2..142, flattened

Columns include A011371, A054861, A090616, A027868, A054861, A054896, A090617, A090618, A027868, A064458, A090619, A090620, A054896, A027868, A090621.
Cf. A090623, A090624, A115627.
Diagonal gives A011776. - Alois P. Heinz, Oct 04 2012

f:= proc(n, p) local c, k; c, k:= 0, p;
       while n>=k do c:= c+iquo(n, k); k:= k*p od; c
    end:
T:= (n, k)-> min(seq(iquo(f(n, i[1]), i[2]), i=ifactors(k)[2])):
seq(seq(T(n, 2+d-n), n=2..d), d=2..20);  # Alois P. Heinz, Oct 04 2012

f[n_, p_] := Module[{c = 0, k = p}, While[n >= k , c = c + Quotient[n, k]; k = k*p ]; c ]; t[n_, k_] := Min[ Table[ Quotient[f[n, i[[1]]], i[[2]]], {i, FactorInteger[k]}]]; Table[ Table[t[n, 2 + d - n], {n, 2, d}], {d, 2, 20}] // Flatten (* Jean-François Alcover, Oct 03 2013, translated from Alois P. Heinz's Maple program *)

For k=p prime: T(n,p) = [n/p] + [n/p^2] + [n/p^3] + .... For k = p^m a prime power: T(n,p^m) = [T(n,p)/m]. For k = b*c with b and c coprime: T(n,a*b) = min(T(n,a), T(n,b)). T(n,k) is close to, but below, n/A090624(k).

A064459 a(n) = Sum_{k>=1} floor(n/12^k).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8

Offset: 0

Robert G. Wilson v, Oct 03 2001

Original incorrect name was: "Highest power of 12 dividing n!": that sequence is A090619. If p is prime, Legendre's formula says the highest power of p dividing n! is Sum_{k>=1} floor(n/p^k), but of course 12 is not prime. - Robert Israel, Mar 23 2018

Harry J. Smith, Table of n, a(n) for n=0..1000
Wikipedia, Legendre's formula.

Cf. A011371, A054861, A090619.

List([0..110],n->Sum([1..n],k-?Int(n/(12^k)))); # Muniru A Asiru, Mar 24 2018

f:= proc(n) add(floor(n/12^k), k=1..floor(log[12](n))) end proc:
f(0):= 0:
map(f, [$0..100]); # Robert Israel, Mar 23 2018

Table[t = 0; p = 12; While[s = Floor[n/p]; t = t + s; s > 0, p *= 12]; t, {n, 0, 100} ]
Join[{0},Accumulate[Table[If[Divisible[n,12],1,0],{n,110}]]] (* Harvey P. Dale, Feb 14 2016 *)

{ for (n=0, 1000, a=0; p=12; while (s = n\p, a+=s; p*=12); write("b064459.txt", n, " ", a) ) } \\ Harry J. Smith, Sep 15 2009

a(n) = floor[n/12] + floor[n/144] + floor[n/1728] + floor[n/20736] + ....

Corrected by Robert Israel, Mar 23 2018

A090621 Exponent of highest power of 16 dividing n!.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 6, 6, 6, 6, 7, 7, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 10, 10, 10, 10, 11, 11, 11, 11, 12, 12, 12, 12, 13, 13, 13, 13, 14, 14, 14, 14, 15, 15, 16, 16, 16, 16, 16, 16, 17, 17, 17, 17, 18, 18, 18, 18, 19, 19, 19, 19

Offset: 0

Henry Bottomley, Dec 06 2003

nonn

			a(10)=2 since 10! = 3628800 = 16^2 * 14175.

Cf. A011371, A054861, A090616, A027868, A054896, A090617, A090618, A090619, A090620.

a(n) = valuation(n!, 16); \\ Michel Marcus, Jul 10 2022

a(n) = A090622(n, 16) = floor(A011371(n)/4) = floor(A090616(n)/2) = floor((floor(n/2) + floor(n/4) + floor(n/8) + floor(n/16) + ...)/4). Almost n/4.

cp's OEIS Frontend

A090622 Square array read by antidiagonals of highest power of k dividing n! (with n,k>1).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A064459 a(n) = Sum_{k>=1} floor(n/12^k).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Maple

Mathematica

PARI

Formula

Extensions

A090621 Exponent of highest power of 16 dividing n!.

Views

Author

Keywords

Examples

Crossrefs

Programs

PARI

Formula