A090623 -id:A090623 - cp's OEIS frontend

A115627 Irregular triangle read by rows: T(n,k) = multiplicity of prime(k) as a divisor of n!.

1, 1, 1, 3, 1, 3, 1, 1, 4, 2, 1, 4, 2, 1, 1, 7, 2, 1, 1, 7, 4, 1, 1, 8, 4, 2, 1, 8, 4, 2, 1, 1, 10, 5, 2, 1, 1, 10, 5, 2, 1, 1, 1, 11, 5, 2, 2, 1, 1, 11, 6, 3, 2, 1, 1, 15, 6, 3, 2, 1, 1, 15, 6, 3, 2, 1, 1, 1, 16, 8, 3, 2, 1, 1, 1, 16, 8, 3, 2, 1, 1, 1, 1

Offset: 2

Franklin T. Adams-Watters, Jan 26 2006

The factorization of n! is n! = 2^T(n,1)*3^T(n,2)*...*p_(pi(n))^T(n,pi(n)) where p_k = k-th prime, pi(n) = A000720(n).
Nonzero terms of A085604; T(n,k) = A085604(n,k), k = 1..A000720(n). - Reinhard Zumkeller, Nov 01 2013
For n=2, 3, 4 and 5, all terms of the n-th row are odd. Are there other such rows? - Michel Marcus, Nov 11 2018
From Gus Wiseman, May 15 2019: (Start)
Differences between successive rows are A067255, so row n is the sum of the first n row-vectors of A067255 (padded with zeros on the right so that all n row-vectors have length A000720(n)). For example, the first 10 rows of A067255 are
  {}
  1
  0 1
  2 0
  0 0 1
  1 1 0
  0 0 0 1
  3 0 0 0
  0 2 0 0
  1 0 1 0
with column sums (8,4,2,1), which is row 10.
(End)
For all prime p > 7, 3*p > 2*nextprime(p), so for any n > 21 there will always be a prime p dividing n! with exponent 2 and there are no further rows with all entries odd. - Charlie Neder, Jun 03 2019

			From _Gus Wiseman_, May 09 2019: (Start)
Triangle begins:
   1
   1  1
   3  1
   3  1  1
   4  2  1
   4  2  1  1
   7  2  1  1
   7  4  1  1
   8  4  2  1
   8  4  2  1  1
  10  5  2  1  1
  10  5  2  1  1  1
  11  5  2  2  1  1
  11  6  3  2  1  1
  15  6  3  2  1  1
  15  6  3  2  1  1  1
  16  8  3  2  1  1  1
  16  8  3  2  1  1  1  1
  18  8  4  2  1  1  1  1
(End)
m such that 5^m||101!: floor(log(101)/log(5)) = 2 terms. floor(101/5) = 20. floor(20/5) = 4. So m = u_1 + u_2 = 20 + 4 = 24. - _David A. Corneth_, Jun 22 2014

T. D. Noe, Rows n = 2..300, flattened
H. T. Davis, Tables of the Mathematical Functions, Vols. 1 and 2, 2nd ed., 1963, Vol. 3 (with V. J. Fisher), 1962; Principia Press of Trinity Univ., San Antonio, TX. [Annotated scan of pages 204-208 of Volume 2.] See Table 2 on page 206.
Wenguang Zhai, On the prime power factorization of n!, Journal of Number Theory, Volume 129, Issue 8, August 2009, pages 1820-1836.
Index entries for sequences related to factorial numbers

Row lengths are A000720.
Row-sums are A022559.
Row-products are A135291.
Row maxima are A011371.
Columns include A011371, A054861, A027868, A054896, A090617, A064458, A090620.
Cf. A090622, A090623, A000142, A115628.
Cf. A085604, A141809.
Cf. A034876, A067255, A071626, A076934, A322583, A325272, A325273, A325276, A325508, A325509.

a115627 n k = a115627_tabf !! (n-2) !! (k-1)
a115627_row = map a100995 . a141809_row . a000142
a115627_tabf = map a115627_row [2..]
-- Reinhard Zumkeller, Nov 01 2013

A115627 := proc(n,k) local d,p; p := ithprime(k) ; n-add(d,d=convert(n,base,p)) ; %/(p-1) ; end proc: # R. J. Mathar, Oct 29 2010

Flatten[Table[Transpose[FactorInteger[n!]][[2]], {n, 2, 20}]] (* T. D. Noe, Apr 10 2012 *)
T[n_, k_] := Module[{p, jm}, p = Prime[k]; jm = Floor[Log[p, n]]; Sum[Floor[n/p^j], {j, 1, jm}]]; Table[Table[T[n, k], {k, 1, PrimePi[n]}], {n, 2, 20}] // Flatten (* Jean-François Alcover, Feb 23 2015 *)

a(n)=my(i=2);while(n-primepi(i)>1,n-=primepi(i);i++);p=prime(n-1);sum(j=1,log(i)\log(p),i\=p) \\ David A. Corneth, Jun 21 2014

T(n,k) = Sum_{i=1..inf} floor(n/(p_k)^i). (Although stated as an infinite sum, only finitely many terms are nonzero.)

T(n,k) = Sum_{i=1..floor(log(n)/log(p_k))} floor(u_i) where u_0 = n and u_(i+1) = floor((u_i)/p_k). - David A. Corneth, Jun 22 2014

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

Original entry on oeis.org

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).

A240236 Triangle read by rows: sum of digits of n in base k, for 2<=k<=n.

Original entry on oeis.org

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

Offset: 2

Franklin T. Adams-Watters, Apr 02 2014

			Triangle starts:
  1
  2 1
  1 2 1
  2 3 2 1
  2 2 3 2 1
  3 3 4 3 2 1

Reinhard Zumkeller, Rows n = 2..100 of triangle, flattened

Row sums give A043306.
Cf. A000120, A053735, A053737, A007953.
See A138530 for another version.
Cf. A002260, A090623.

a240236 n k = a240236_tabl !! (n-1) !! (k-1)
a240236_row n = a240236_tabl !! (n-1)
a240236_tabl = zipWith (map . flip q)
                       [2..] (map tail $ tail a002260_tabl) where
   q b n = if n < b then n else q b n' + d where (n', d) = divMod n b
-- Reinhard Zumkeller, Apr 29 2015

Table[Total[Flatten[IntegerDigits[n,k]]],{n,20},{k,2,n}]//Flatten (* Harvey P. Dale, Jan 13 2025 *)

T(n,k) = local(r=0);if(k<2,-1,while(n>0,r+=n%k;n\=k);r)

T(n, k) = sumdigits(n, k) \\ Zhuorui He, Aug 25 2025

T(n,k) = n - (k - 1) * Sum_{i=1..floor(log_k(n))} floor(n/k^i). - Ridouane Oudra, Sep 27 2024

T(n,k) = n - (k - 1) * A090623(n,k). - Zhuorui He, Aug 25 2025

A090620 Highest power of 13 dividing n!.

Original entry on oeis.org

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

Offset: 0

Henry Bottomley, Dec 06 2003

Alois P. Heinz, Table of n, a(n) for n = 0..10000 (first 1001 terms from T. D. Noe)

Cf. A011371, A054861.

a:= proc(n) option remember; `if`(n=0, 0,
      a(n-1)+padic[ordp](n, 13))
    end:
seq(a(n), n=0..120);  # Alois P. Heinz, Jun 20 2020

IntegerExponent[Range[0,110]!,13] (* Harvey P. Dale, Aug 22 2011 *)
FoldList[Plus, 0, IntegerExponent[Range[100], 13]] (* T. D. Noe, Apr 10 2012 *)

a(n)=my(t);while(n,t+=n\=13);t \\ Charles R Greathouse IV, Aug 06 2012

a(n) = A090622(n, 13) = A090623(n, 13) = [n/13]+[n/169]+[n/2197]+...

a(n) = n/12 + O(log n). - Charles R Greathouse IV, Aug 06 2012

A078632 Number of geometric subsequences of [1,...,n] with integral successive-term ratio and length > 1.

Original entry on oeis.org

0, 1, 2, 5, 6, 9, 10, 15, 18, 21, 22, 28, 29, 32, 35, 43, 44, 50, 51, 57, 60, 63, 64, 73, 76, 79, 84, 90, 91, 98, 99, 109, 112, 115, 118, 129, 130, 133, 136, 145, 146, 153, 154, 160, 166, 169, 170, 183, 186, 192, 195, 201, 202, 211, 214, 223, 226, 229, 230, 242

Offset: 1

Robert E. Sawyer (rs.1(AT)mindspring.com)

The number of geometric subsequences of [1,...,n] with integral successive-term ratio r and length k is floor(n/r^(k-1))(n > 0, r > 1, k > 0).

			a(2): [1,2]; a(3): [1,2],[1,3]; a(4): [1,2],[1,3],[1,4],[2,4],[1,2,4].

Zhuorui He, Table of n, a(n) for n = 1..10000

Cf. A078651.
Row sums of triangle A090623.
Partial sums of A309891.

g := (n, b) -> local i; add(iquo(n, b^i), i = 1..floor(log(n, b))):
a := n -> local b; add(g(n, b), b = 2..n):
seq(a(n), n = 1..60);  # Peter Luschny, Apr 03 2025

Accumulate[Table[Total[IntegerExponent[n, Rest[Divisors[n]]]], {n, 100}]] (* Paolo Xausa, Aug 27 2025 *)

A078632(n) = {my(s=0, k=2); while(k<=n, s+=(n - sumdigits(n, k))/(k-1); k=k+1); s} \\ Zhuorui He, Aug 26 2025

a(n) = Sum_{r > 1, j > 0} floor(n/r^j).

cp's OEIS Frontend

A115627 Irregular triangle read by rows: T(n,k) = multiplicity of prime(k) as a divisor of n!.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

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

A240236 Triangle read by rows: sum of digits of n in base k, for 2<=k<=n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

PARI

Formula

A090620 Highest power of 13 dividing n!.

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A078632 Number of geometric subsequences of [1,...,n] with integral successive-term ratio and length > 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula