cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 10 results.

A146291 Triangle T(n,m) read by rows (n >= 1, 0 <= m <= A001222(n)), giving the number of divisors of n with m prime factors (counted with multiplicity).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 2, 2, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 2, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 3, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 2
Offset: 1

Views

Author

Matthew Vandermast, Nov 11 2008

Keywords

Comments

All rows are palindromic. T(n,0) = T(n,A001222(n)) = 1.
Two numbers have identical rows in the table if and only if they have the same prime signature.
If n is a perfect square then Sum_{even m} T(n,m) = 1 + Sum_{odd m} T(n,m), otherwise Sum_{even m} T(n,m) = Sum_{odd m} T(n,m). - Geoffrey Critzer, Feb 08 2015

Examples

			Rows begin:
1;
1, 1;
1, 1;
1, 1, 1;
1, 1;
1, 2, 1;
1, 1;
1, 1, 1, 1;
1, 1, 1;
1, 2, 1;
...
12 has 1 divisor with 0 total prime factors (1), 2 with 1 (2 and 3), 2 with 2 (4 and 6) and 1 with 3 (12), for a total of 6. The 12th row of the table therefore reads (1, 2, 2, 1). These are the positive coefficients of the polynomial 1 + 2k + 2k^2 + (1)k^3 = (1 + k + k^2)(1 + k), derived from the prime factorization of 12 (namely, 2^2*3^1).
		

Crossrefs

Row sums equal A000005(n). T(n,1) = A001221(n) for n>1.
Row n of A007318 is identical to row A002110(n) of this table and also identical to the row for any squarefree number with n prime factors.

Programs

  • Maple
    with(numtheory):
    T:= n-> (p-> seq(coeff(p, x, i), i=0..degree(p)))(
             add(x^bigomega(d), d=divisors(n))):
    seq(T(n), n=1..100);  # Alois P. Heinz, Feb 25 2015
  • Mathematica
    Join[{{1}},
    Table[nn = DivisorSigma[0, n];
      CoefficientList[
       Series[Product[(1 - x^i)/(1 - x), {i,
    FactorInteger[n][[All, 2]] + 1}], {x, 0, nn}], x], {n, 2, 100}]] (* Geoffrey Critzer, Jan 01 2015 *)

Formula

If the canonical factorization of n into prime powers is the product of p^e(p), then T(n,m) is the coefficient of k^m in the polynomial expansion of Product_p (sum_{i=0..e} k^i).

A146289 Triangle T(n,m) read by rows (n >= 1, 0 <= m <= A001221(n)), giving the number of divisors of n with m distinct prime factors.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 3, 1, 2, 1, 2, 1, 1, 1, 1, 3, 2, 1, 1, 1, 2, 1, 1, 2, 1, 1, 4, 1, 1, 1, 3, 2, 1, 1, 1, 3, 2, 1, 2, 1, 1, 2, 1, 1, 1, 1, 4, 3, 1, 2, 1, 2, 1, 1, 3, 1, 3, 2, 1, 1, 1, 3, 3, 1, 1, 1, 1, 5, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 4, 4, 1, 1, 1, 2, 1, 1, 2, 1, 1, 4, 3, 1, 1, 1, 3, 3
Offset: 1

Views

Author

Matthew Vandermast, Nov 11 2008

Keywords

Comments

The formula used in obtaining the n-th row (see below) also gives the number of divisors of the k-th power of n.
Two numbers have identical rows in the table if and only if they have the same prime signature.
T(n,0)=1.

Examples

			Rows begin:
1;
1,1;
1,1;
1,2;
1,1;
1,2,1;
1,1;
1,3;
1,2;
1,2,1;
...
12 has 1 divisor with 0 distinct prime factors (1); 3 with 1 (2, 3 and 4); and 2 with 2 (6 and 12), for a total of 6. The 12th row of the table therefore reads (1, 3, 2). These are the positive coefficients of the polynomial equation 1 + 3k + 2k^2 = (1 + 2k)(1 + k), derived from the prime factorization of 12 (namely, 2^2*3^1).
		

Crossrefs

Row sums equal A000005(n).
T(n, 1) = A001222(n) for n>1. T(n, A001221(n)) = A005361(n).
Row n of A007318 is identical to row A002110(n) of this table and also identical to the row for any squarefree number with n prime factors.
Cf. A146290. Also cf. A146291, A146292.

Programs

  • Maple
    f:= proc(n)
       local F,G,f,t,k;
       F:= ifactors(n)[2];
       G:= mul(1+f[2]*t, f= F);
       seq(coeff(G,t,k),k=0..nops(F));
    end proc:
    seq(f(n),n=1..100); # Robert Israel, Feb 10 2015
  • Mathematica
    Join[{{1}}, Table[nn = DivisorSigma[0, n];CoefficientList[Series[Product[1 + i x, {i, FactorInteger[n][[All, 2]]}], {x, 0,nn}], x], {n, 2, 100}]] // Grid (* Geoffrey Critzer, Feb 09 2015 *)
  • PARI
    tabf(nn) = {for (n=1, nn, vd = divisors(n); vo = vector(#vd, k, omega(vd[k])); for (k=0, vecmax(vo), print1(#select(x->x==k, vo), ", ");); print(););} \\ Michel Marcus, Apr 22 2017

Formula

If the canonical factorization of n into prime powers is Product p^e(p), then T(n, m) is the coefficient of k^m in the polynomial expansion of Product_p (1 + e(p) k).

A181555 a(n) = A002110(n)^n.

Original entry on oeis.org

1, 2, 36, 27000, 1944810000, 65774855015100000, 733384949590939374729000000, 9037114296609938214167920266348510000000, 78354300210436852307898467208663359164858967744100000000
Offset: 0

Views

Author

Matthew Vandermast, Oct 31 2010

Keywords

Comments

For n>0, a(n)= first counting number whose prime signature consists of n repeated n times (cf. A002024). Subsequence of A025487.

Examples

			a(4) = 1944810000 = 210^4 = 2^4 * 3^4 * 5^4 * 7^4.
		

Crossrefs

A061742(n) = A002110(n)^2. See also A006939, A066120, A166475, A167448.
A000005(a(n)) = A000169(n). The divisors of a(n) appear as the first A000169(n) terms of A178479, with A178479(A000169(n)) = a(n).
A071207(n, k) gives the number of divisors of n with (n-k) distinct prime factors, A181567(n, k) gives the number of divisors of n with k prime factors counted with multiplicity.

Programs

  • Mathematica
    a[0] = 1; a[n_] := Product[Prime[i], {i, 1, n}]^n; Array[a, 9, 0] (* Amiram Eldar, Aug 08 2019 *)

Formula

a(n) = A079474(2n,n). - Alois P. Heinz, Aug 22 2019

A168264 For all sufficiently high values of k, d(n^k) > d(m^k) for all m < n. (Let k, m, and n represent positive integers only.)

Original entry on oeis.org

1, 2, 4, 6, 12, 24, 30, 60, 120, 180, 210, 420, 840, 1260, 1680, 2310, 4620, 9240, 13860, 18480, 27720, 30030, 60060, 120120, 180180, 240240, 360360, 510510, 1021020, 2042040, 3063060, 4084080, 6126120, 9699690, 19399380, 38798760, 58198140
Offset: 1

Views

Author

Matthew Vandermast, Nov 23 2009

Keywords

Comments

d(n) is the number of divisors of n (A000005(n)).

Examples

			Since the exponents in 1680's prime factorization are (4,1,1,1), the k-th power of 1680 has (4k+1)(k+1)^3 = 4k^4 + 13k^3 + 15k^2 + 7k + 1 divisors. Comparison with the analogous formulas for all smaller members of A025487 shows the following:
a) No number smaller than 1680 has a positive coefficient in its "power formula" for any exponent larger than k^4.
b) The only power formula with a k^4 coefficient as high as 4 is that for 1260 (4k^4 + 12k^3 + 13k^2 + 6k + 1).
c) The k^3 coefficient for 1680 is higher than for 1260.
So for all sufficiently high values of k, d(1680^k) > d(m^k) for all m < 1680.
		

Crossrefs

Subsequence of A025487, A060735, A116998. Includes A002110, A168262, A168263.
See also A168265, A168266, A168267.

Formula

If the canonical factorization of n into prime powers is Product p^e(p), then the formula for the number of divisors of the k-th power of n is Product_p (ek + 1). (See also A146289, A146290.)
For two positive integers m and n with different prime signatures, let j be the largest exponent of k for which m and n have different coefficients, after the above formula for each integer is expanded as a polynomial. Let m_j and n_j denote the corresponding coefficients. d(n^k) > d(m^k) for all sufficiently high values of k if and only if n_j > m_j.

A071207 Triangular array T(n,k) read by rows, giving number of rooted trees on the vertex set {1..n+1} where k children of the root have a label smaller than the label of the root.

Original entry on oeis.org

1, 1, 1, 4, 4, 1, 27, 27, 9, 1, 256, 256, 96, 16, 1, 3125, 3125, 1250, 250, 25, 1, 46656, 46656, 19440, 4320, 540, 36, 1, 823543, 823543, 352947, 84035, 12005, 1029, 49, 1, 16777216, 16777216, 7340032, 1835008, 286720, 28672, 1792, 64, 1, 387420489
Offset: 0

Views

Author

Cedric Chauve (chauve(AT)lacim.uqam.ca), May 16 2002

Keywords

Comments

The n-th term of the n-th binomial transform of a sequence {b} is given by {d} where d(n) = sum(k=0,n,T(n,k)*b(k)) and T(n,k)=binomial(n,k)*n^(n-k); such diagonals are related to the hyperbinomial transform (A088956). - Paul D. Hanna, Nov 04 2003
T(n,k) gives the number of divisors of A181555(n) with (n-k) distinct prime factors. See also A001221, A146289, A146290, A181567. - Matthew Vandermast, Oct 31 2010
T(n,k) is the number of partial functions on {1,2,...,n} leaving exactly k elements undefined. Row sums = A000169. - Geoffrey Critzer, Jan 08 2012
As a triangular matrix, transforms rows into diagonals in the table of coefficients of successive iterations of x/(1-x). - Paul D. Hanna, Jan 19 2014
Also the number of rooted trees on n+1 labeled vertices in which some specified vertex (say, vertex 1) has k children. - Alan Sokal, Jul 22 2022

Examples

			1
1     1
4     4     1
27    27    9     1
256   256   96    16    1
3125  3125  1250  250   25    1
46656 46656 19440 4320  540   36    1
		

Crossrefs

Programs

  • Maple
    T:= (n, k)-> binomial(n, k)*n^(n-k): seq(seq(T(n, k), k=0..n), n=0..10);
  • Mathematica
    Prepend[Flatten[ Table[Table[Binomial[n, k] n^(n - k), {k, 0, n}], {n, 1, 8}]], 1]  (* Geoffrey Critzer, Jan 08 2012 *)
  • PARI
    T(n,k)=if(k<0 || k>n,0,binomial(n,k)*n^(n-k))
    
  • PARI
    /* Transforms rows into diagonals in the iterations of x/(1-x): */
    {T(n, k)=local(F=x, M, N, P, m=n); M=matrix(m+2, m+2, r, c, F=x; for(i=1, r+c-2, F=subst(F, x, x/(1-x+x*O(x^(m+2))))); polcoeff(F, c)); N=matrix(m+1, m+1, r, c, F=x; for(i=1, r, F=subst(F, x, x/(1-x+x*O(x^(m+2))))); polcoeff(F, c)); P=matrix(m+1, m+1, r, c, M[r+1, c]); (P~*N~^-1)[n+1, k+1]}
    for(n=0, 10, for(k=0, n, print1(T(n, k), ", ")); print("")) \\ Paul D. Hanna, Jan 19 2014

Formula

T(n,k) = binomial(n, k)*n^(n-k).
E.g.f.: (-LambertW(-y)/y)^x/(1+LambertW(-y)). - Vladeta Jovovic

Extensions

Name edited by Alan Sokal, Jul 22 2022

A146288 Number of divisors of the n-th prime signature number (A025487(n)).

Original entry on oeis.org

1, 2, 3, 4, 4, 6, 5, 8, 8, 6, 9, 10, 12, 7, 12, 12, 16, 8, 15, 18, 14, 16, 16, 20, 9, 18, 24, 16, 24, 20, 24, 10, 21, 30, 18, 32, 24, 27, 28, 11, 32, 24, 36, 25, 36, 20, 40, 28, 36, 32, 12, 40, 27, 32, 48, 30, 42, 22, 48, 32, 45, 36, 13, 48, 30, 48, 60, 35, 48, 48, 24, 54, 50, 56
Offset: 1

Views

Author

Matthew Vandermast, Nov 11 2008

Keywords

Examples

			a(4) = 4 because 4 positive integers divide evenly into A025487(4) = 6: 1, 2, 3 and 6.
		

Crossrefs

a(n) = sum of the n-th row of A146290, A146292.
A rearrangement of A080444.
Cf. A000005 (number of divisors), A025487.

Programs

  • Haskell
    a146288 = a000005 . a025487  -- Reinhard Zumkeller, Sep 17 2014
  • Mathematica
    s = {1}; Do[If[GreaterEqual @@ (f = FactorInteger[n])[[;; , 2]] && PrimePi[f[[-1, 1]]] == Length[f], AppendTo[s, DivisorSigma[0, n]]], {n, 2, 10000}]; s (* Amiram Eldar, Aug 05 2024 *)

Formula

a(n) = A000005(A025487(n)).

A168263 For any m < n, and for all values of k, d(n^k) > d(m^k). (Let k, m, and n represent positive integers only.)

Original entry on oeis.org

1, 2, 4, 6, 12, 24, 60, 120, 180, 840, 1260, 1680, 27720
Offset: 1

Views

Author

Matthew Vandermast, Nov 23 2009

Keywords

Comments

d(n) is the number of divisors of n (A000005(n)).
All members must be highly composite numbers (A002182) with at least as many distinct prime factors as any smaller positive integer (A116998). (See Formula and Example sections.) It turns out that these two conditions are jointly sufficient.
Ramanujan proved that a) for any prime p, there exist a finite number of highly composite numbers with p as its largest prime factor; and b) in the canonical prime factorization of a highly composite number with largest prime factor p, the exponents for all primes > p are never smaller than they are in the factorization of A003418(p). (See formula 54 of the Ramanujan paper.)
It follows that, if the intersection of A003418 and A116998 is finite, so is the intersection of A002182 and A116998. For proof that the former intersection is finite, see A168262.
By using the given formula for the number of divisors, it is possible to define a canonical polynomial p_n(k) for every natural number n. For example, because 60 = (2^2)(3^1)(5^1), we define p_60(k) = (1+2k)(1+k)(1+k). The present sequence is defined only by examining whether p_n(k) achieves a record for natural numbers k, but the question could also be asked whether p_n(k) achieves a record for all k > 0. This stricter requirement does not hold for a(7)-a(13) at various positive values of k < 1, but it does hold for a(1)-a(6). The present sequence is "full", so a(1)-a(6) are the only numbers to satisfy the stronger property. - Hal M. Switkay, Aug 17 2025

Examples

			1) 1680 has more divisors than any smaller positive integer; thus for all m < n, d(1680^1) > d(m^1).
2) Since the exponents in 1680's prime factorization are (4,1,1,1), the k-th power of 1680 has (4k+1)(k+1)^3 = 4k^4 + 13k^3 + 15k^2 + 7k + 1 divisors. Comparison with the analogous formulas for all smaller members of A025487 shows the following:
a) No number smaller than 1680 has a positive coefficient in its "power formula" for any exponent larger than k^4.
b) The only power formula with a k^4 coefficient as high as 4 is that for 1260 (4k^4 + 12k^3 + 13k^2 + 6k + 1).
c) The k^3 coefficient for 1680 is higher than for 1260.
So for all sufficiently high values of k, d(1680^k) > d(m^k) for all m < 1680.
3) Careful comparison of 1680's "power formula" with the analogous formulas for smaller members of A025487 shows that no intermediate value of k can exist for which d(m^k) >= d(1680^k) if m < 1680.
		

References

  • S. Ramanujan, Highly composite numbers, Proc. Lond. Math. Soc. 14 (1915), 347-409; reprinted in Collected Papers, Ed. G. H. Hardy et al., Cambridge 1927; Chelsea, NY, 1962.

Crossrefs

Intersection of A002182 and A116998. Also, intersection of A002182 and A060735, and of A002182 and A168264. (A168264 is a subsequence of A060735, which is a subsequence of A116998.)

Formula

If the canonical factorization of n into prime powers is Product p^e(p), then the formula for the number of divisors of the k-th power of n is Product_p (ek + 1). (See also A146289, A146290.)
For two positive integers m and n with different prime signatures, let j be the largest exponent of k for which m and n have different coefficients, after the above formula for each integer is expanded as a polynomial. Let m_j and n_j denote the corresponding coefficients. d(n^k) > d(m^k) for all sufficiently high values of k if and only if n_j > m_j.

A182860 Number of distinct prime signatures represented among the unitary divisors of n.

Original entry on oeis.org

1, 2, 2, 2, 2, 3, 2, 2, 2, 3, 2, 4, 2, 3, 3, 2, 2, 4, 2, 4, 3, 3, 2, 4, 2, 3, 2, 4, 2, 4, 2, 2, 3, 3, 3, 3, 2, 3, 3, 4, 2, 4, 2, 4, 4, 3, 2, 4, 2, 4, 3, 4, 2, 4, 3, 4, 3, 3, 2, 6, 2, 3, 4, 2, 3, 4, 2, 4, 3, 4, 2, 4, 2, 3, 4, 4, 3, 4, 2, 4, 2, 3, 2, 6, 3, 3, 3, 4, 2, 6, 3, 4, 3, 3, 3, 4, 2, 4, 4, 3, 2, 4, 2, 4, 4
Offset: 1

Views

Author

Matthew Vandermast, Jan 14 2011

Keywords

Comments

a(n) = number of members m of A025487 such that d(m^k) divides d(n^k) for all values of k. (Here d(n) represents the number of divisors of n, or A000005(n).)
a(n) depends only on prime signature of n (cf. A025487).

Examples

			60 has 8 unitary divisors (1, 3, 4, 5, 12, 15, 20 and 60). Primes 3 and 5 have the same prime signature, as do 12 (2^2*3) and 20 (2^2*5); each of the other four numbers listed is the only unitary divisor of 60 with its particular prime signature.  Since a total of 6 distinct prime signatures appear among the unitary divisors of 60, a(60) = 6.
		

Crossrefs

Programs

  • Mathematica
    Table[Length@ Union@ Map[Sort[FactorInteger[#] /. {p_, e_} /; p > 0 :> If[p == 1, 0, e]] &, Select[Divisors@ n, CoprimeQ[#, n/#] &]], {n, 105}] (* Michael De Vlieger, Jul 19 2017 *)
  • PARI
    A181819(n) = {my(f=factor(n)); prod(k=1, #f~, prime(f[k, 2])); }; \\ From A181819
    A182860(n) = numdiv(A181819(n)); \\ Antti Karttunen, Jul 19 2017

Formula

a(n) = A000005(A181819(n)).
If the canonical factorization of n into prime powers is Product p^e(p), then the formula for d(n^k) is Product_p (ek + 1). (See also A146289, A146290.)
a(n) = A064553(A328830(n)). - Antti Karttunen, Apr 29 2022

A146292 Triangle T(n,m) read by rows (n >= 1, 0 <= m <= A036041(n)), giving the number of divisors of A025487(n) with m prime factors (counted with multiplicity).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 1, 1, 3, 3, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 2, 1, 1, 2, 2, 2, 2, 1, 1, 3, 4, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 3, 2, 1, 1, 2, 2, 2, 2, 2, 1, 1, 3, 4, 4, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 3, 3, 2, 1, 1, 3, 5, 5, 3, 1, 1, 2, 2, 2, 2
Offset: 1

Views

Author

Matthew Vandermast, Nov 11 2008

Keywords

Comments

All rows are palindromic. T(n, 0) = T(n, A036041(n)) = 1.
Every row that appears in A146291 appears exactly once in the table. Rows appear in order of first appearance in A146291.

Examples

			Rows begin:
  1;
  1,1;
  1,1,1;
  1,2,1;
  1,1,1,1;
  1,2,2,1;
  1,1,1,1,1;...
36's 9 divisors include 1 divisor with 0 total prime factors (1);, 2 with 1 (2 and 3); 3 with 2 (4, 6 and 9); 2 with 3 (12 and 18); and 1 with 4 (36). Since 36 = A025487(11), the 11th row of the table therefore reads (1, 2, 3, 2, 1). These are the positive coefficients of the polynomial 1 + 2k + 3k^2 + 2k^3 + (1)k^4 = (1 + k + k^2)(1 + k + k^2), derived from the prime factorization of 36 (namely, 2^2*3^2).
		

Crossrefs

For the number of prime factors of n counted with multiplicity, see A001222.
Row sums equal A146288(n). T(n, 1) = A061394(n) for n>1.
Row A098719(n) of this table is identical to row n of A007318.
Cf. A146291. Also cf. A146289, A146290.

Formula

If A025487(n)'s canonical factorization into prime powers is the product of p^e(p), then T(n, m) is the coefficient of k^m in the polynomial expansion of Product_p (sum_{i=0..e} k^i).

A182861 Number of distinct prime signatures represented among the unitary divisors of A025487(n).

Original entry on oeis.org

1, 2, 2, 3, 2, 4, 2, 4, 4, 2, 3, 4, 6, 2, 4, 4, 6, 2, 4, 6, 4, 5, 3, 6, 2, 4, 8, 4, 8, 4, 6, 2, 4, 8, 4, 8, 4, 4, 6, 2, 6, 4, 9, 3, 8, 4, 8, 4, 6, 6, 2, 8, 4, 6, 12, 4, 8, 4, 8, 4, 6, 6, 2, 8, 4, 10, 12, 4, 6, 8, 4, 8, 6, 8, 4, 6, 9, 6, 3, 2, 8, 4, 10, 12, 4
Offset: 1

Views

Author

Matthew Vandermast, Jan 14 2011

Keywords

Comments

a(n) = number of members m of A025487 such that d(m^k) divides d(A025487(n)^k) for all values of k. (Here d(n) represents the number of divisors of n, or A000005(n).)

Examples

			60 has 8 unitary divisors (1, 3, 4, 5, 12, 15, 20 and 60). Primes 3 and 5 have the same prime signature, as do 12 (2^2*3) and 20 (2^2*5); each of the other four numbers listed is the only unitary divisor of 60 with its particular prime signature.  Since a total of 6 distinct prime signatures appear among the unitary divisors of 60, and since 60 = A025487(13), a(13) = 6.
		

Crossrefs

Formula

a(n) = A000005(A181820(n)).
If the canonical factorization of n into prime powers is Product p^e(p), then the formula for d(n^k) is Product_p (ek + 1). (See also A146289, A146290.)

Extensions

More terms from Amiram Eldar, Jun 20 2019
Showing 1-10 of 10 results.