A061394 -id:A061394 - cp's OEIS frontend

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

1, 1, 1, 1, 2, 1, 2, 1, 1, 3, 1, 3, 2, 1, 4, 1, 4, 3, 1, 3, 3, 1, 1, 5, 1, 4, 4, 1, 5, 4, 1, 4, 5, 2, 1, 6, 1, 5, 6, 1, 6, 5, 1, 5, 7, 3, 1, 7, 1, 6, 8, 1, 5, 8, 4, 1, 7, 6, 1, 4, 6, 4, 1, 1, 6, 9, 1, 6, 9, 4, 1, 8, 1, 7, 10, 1, 6, 11, 6, 1, 8, 7, 1, 5, 9, 7, 2, 1, 7, 12, 1, 7, 11, 5, 1, 9, 1, 8, 12, 1, 7, 14

Offset: 1

Matthew Vandermast, Nov 11 2008

The formula used in obtaining the A025487(n)th row (see below) also gives the number of divisors of the k-th power of A025487(n).
Every row that appears in A146289 appears exactly once in the table. Rows appear in order of first appearance in A146289.
T(n,0)=1.

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

Anonymous?, Polynomial calculator
Eric Weisstein's World of Mathematics, Distinct Prime Factors
G. Xiao, WIMS server, Factoris (both expands and factors polynomials)

For the number of distinct prime factors of n, see A001221.
Row sums equal A146288(n). T(n, 1)=A036041(n) for n>1. T(n, A061394(n))=A052306(n).
Row A098719(n) of this table is identical to row n of A007318.
Cf. A146289. Also cf. A146291, A146292.

If A025487(n)'s canonical factorization 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 + ek).

A025487 Least integer of each prime signature A124832; also products of primorial numbers A002110.

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 16, 24, 30, 32, 36, 48, 60, 64, 72, 96, 120, 128, 144, 180, 192, 210, 216, 240, 256, 288, 360, 384, 420, 432, 480, 512, 576, 720, 768, 840, 864, 900, 960, 1024, 1080, 1152, 1260, 1296, 1440, 1536, 1680, 1728, 1800, 1920, 2048, 2160, 2304, 2310

Offset: 1

David W. Wilson

All numbers of the form 2^k1*3^k2*...*p_n^k_n, where k1 >= k2 >= ... >= k_n, sorted.
A111059 is a subsequence. - Reinhard Zumkeller, Jul 05 2010
Choie et al. (2007) call these "Hardy-Ramanujan integers". - Jean-François Alcover, Aug 14 2014
The exponents k1, k2, ... can be read off Abramowitz & Stegun p. 831, column labeled "pi".
For all such sequences b for which it holds that b(n) = b(A046523(n)), the sequence which gives the indices of records in b is a subsequence of this sequence. For example, A002182 which gives the indices of records for A000005, A002110 which gives them for A001221 and A000079 which gives them for A001222. - Antti Karttunen, Jan 18 2019
The prime signature corresponding to a(n) is given in row n of A124832. - M. F. Hasler, Jul 17 2019

			The first few terms are 1, 2, 2^2, 2*3, 2^3, 2^2*3, 2^4, 2^3*3, 2*3*5, ...

Franklin T. Adams-Watters, Table of n, a(n) for n = 1..10001 (first 291 terms from Will Nicholes)
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972.
Kevin Broughan, Equivalents of the Riemann Hypothesis, Vol. 1: Arithmetic Equivalents, Cambridge University Press, 2017. See section 8.2, "Hardy-Ramanujan Numbers".
YoungJu Choie, Nicolas Lichiardopol, Pieter Moree and Patrick Solé, On Robin's criterion for the Riemann hypothesis, Journal de théorie des nombres de Bordeaux, Vol. 19, No. 2 (2007), pp. 357-372. See section 5, p. 367.
Asaf Cohen Antonir and Asaf Shapira, An Elementary Proof of a Theorem of Hardy and Ramanujan (2022). arXiv:2207.09410 [math.NT]
Michael De Vlieger, Relations of A025487 to A002110, A002182, and A002201.
Steven R. Finch, Errata and Addenda to Mathematical Constants, arXiv:2001.00578 [math.HO], 2020, pp. 9-10.
G. H. Hardy and S. Ramanujan, Asymptotic formulae for the distribution of integers of various types, Proc. London Math. Soc, Ser. 2, Vol. 16 (1917), pp. 112-132. Also published in the book Collected Papers of Srinivasa Ramanujan, Chelsea, 1962, pages 245-261.
Jeffery Kline, On the eigenstructure of sparse matrices related to the prime number theorem, Linear Algebra and its Applications (2020) Vol. 584, 409-430.
L. B. Richmond, Asymptotic results for partitions (I) and the distribution of certain integers, Journal of Number Theory, Vol. 8, No. 4 (1976), pp. 372-389. See page 388.

Subsequence of A055932, A191743, and A324583.
Cf. A025488, A051282, A036041, A051466, A061394, A124832, A161360, A166469, A181815, A181817, A283980, A306802, A322584, A322585 (characteristic function), A329897, A329898, A329899, A329900, A329904, A330683.
Cf. A085089, A101296 (left inverses).
Equals range of values taken by A046523.
Cf. A178799 (first differences), A247451 (squarefree kernel), A146288 (number of divisors).
Subsequences of this sequence include: A000079, A000142, A000400, A001013, A001813, A002110, A002182, A005179, A006939, A025527, A056836, A061742, A064350, A066120, A087980, A097212, A097213, A111059, A119840, A119845, A126098, A129912, A140999, A166338, A166470, A166472, A166473, A166475, A167448, A168262, A168263, A168264, A179215, A181555, A181804, A181806, A181809, A181818, A181822, A181823, A181824, A181825, A181826, A181827, A182763, A182862, A182863, A212170, A220264, A220423, A250269, A250270, A260633, A266047, A284456, A300357, A304938, A329894, A330687; also A037019 and A330681 (when sorted), possibly also A289132.
Rearrangements of this sequence include A036035, A059901, A063008, A077569, A085988, A086141, A087443, A108951, A181821, A181822, A322827, A329886, A329887.
Cf. also array A124832 (row n = prime signature of a(n)) and A304886, A307056.

import Data.Set (singleton, fromList, deleteFindMin, union)
a025487 n = a025487_list !! (n-1)
a025487_list = 1 : h [b] (singleton b) bs where
   (_ : b : bs) = a002110_list
   h cs s xs'@(x:xs)
     | m <= x    = m : h (m:cs) (s' `union` fromList (map (* m) cs)) xs'
     | otherwise = x : h (x:cs) (s  `union` fromList (map (* x) (x:cs))) xs
     where (m, s') = deleteFindMin s
-- Reinhard Zumkeller, Apr 06 2013

isA025487 := proc(n)
    local pset,omega ;
    pset := sort(convert(numtheory[factorset](n),list)) ;
    omega := nops(pset) ;
    if op(-1,pset) <> ithprime(omega) then
        return false;
    end if;
    for i from 1 to omega-1 do
        if padic[ordp](n,ithprime(i)) < padic[ordp](n,ithprime(i+1)) then
            return false;
        end if;
    end do:
    true ;
end proc:
A025487 := proc(n)
    option remember ;
    local a;
    if n = 1 then
        1 ;
    else
        for a from procname(n-1)+1 do
            if isA025487(a) then
                return a;
            end if;
        end do:
    end if;
end proc:
seq(A025487(n),n=1..100) ; # R. J. Mathar, May 25 2017

PrimeExponents[n_] := Last /@ FactorInteger[n]; lpe = {}; ln = {1}; Do[pe = Sort@PrimeExponents@n; If[ FreeQ[lpe, pe], AppendTo[lpe, pe]; AppendTo[ln, n]], {n, 2, 2350}]; ln (* Robert G. Wilson v, Aug 14 2004 *)
(* Second program: generate all terms m <= A002110(n): *)
f[n_] := {{1}}~Join~
  Block[{lim = Product[Prime@ i, {i, n}],
   ww = NestList[Append[#, 1] &, {1}, n - 1], dec},
   dec[x_] := Apply[Times, MapIndexed[Prime[First@ #2]^#1 &, x]];
   Map[Block[{w = #, k = 1},
      Sort@ Prepend[If[Length@ # == 0, #, #[[1]]],
        Product[Prime@ i, {i, Length@ w}] ] &@ Reap[
         Do[
          If[# < lim,
             Sow[#]; k = 1,
             If[k >= Length@ w, Break[], k++]] &@ dec@ Set[w,
             If[k == 1,
               MapAt[# + 1 &, w, k],
               PadLeft[#, Length@ w, First@ #] &@
                 Drop[MapAt[# + Boole[i > 1] &, w, k], k - 1] ]],
           {i, Infinity}] ][[-1]]
] &, ww]]; Sort[Join @@ f@ 13] (* Michael De Vlieger, May 19 2018 *)

isA025487(n)=my(k=valuation(n,2),t);n>>=k;forprime(p=3,default(primelimit),t=valuation(n,p);if(t>k,return(0),k=t);if(k,n/=p^k,return(n==1))) \\ Charles R Greathouse IV, Jun 10 2011

factfollow(n)={local(fm, np, n2);
  fm=factor(n); np=matsize(fm)[1];
  if(np==0,return([2]));
  n2=n*nextprime(fm[np,1]+1);
  if(np==1||fm[np,2]Franklin T. Adams-Watters, Dec 01 2011 */

is(n) = {if(n==1, return(1)); my(f = factor(n));  f[#f~, 1] == prime(#f~) && vecsort(f[, 2],,4) == f[, 2]} \\ David A. Corneth, Feb 14 2019

upto(Nmax)=vecsort(concat(vector(logint(Nmax,2),n,select(t->t<=Nmax,if(n>1,[factorback(primes(#p),Vecrev(p)) || p<-partitions(n)],[1,2]))))) \\ M. F. Hasler, Jul 17 2019

\\ For fast generation of large number of terms, use this program:
A283980(n) = {my(f=factor(n)); prod(i=1, #f~, my(p=f[i, 1], e=f[i, 2]); if(p==2, 6, nextprime(p+1))^e)}; \\ From A283980
A025487list(e) = { my(lista = List([1, 2]), i=2, u = 2^e, t); while(lista[i] != u, if(2*lista[i] <= u, listput(lista,2*lista[i]); t = A283980(lista[i]); if(t <= u, listput(lista,t))); i++); vecsort(Vec(lista)); }; \\ Returns a list of terms up to the term 2^e.
v025487 = A025487list(101);
A025487(n) = v025487[n];
for(n=1,#v025487,print1(A025487(n), ", ")); \\ Antti Karttunen, Dec 24 2019

def sharp_primorial(n): return sloane.A002110(prime_pi(n))
N = 2310
nmax = 2^floor(log(N,2))
sorted([j for j in (prod(sharp_primorial(t[0])^t[1] for k, t in enumerate(factor(n))) for n in (1..nmax)) if j <= N])
# Giuseppe Coppoletta, Jan 26 2015

What can be said about the asymptotic behavior of this sequence? - Franklin T. Adams-Watters, Jan 06 2010
Hardy & Ramanujan prove that there are exp((2 Pi + o(1))/sqrt(3) * sqrt(log x/log log x)) members of this sequence up to x. - Charles R Greathouse IV, Dec 05 2012
From Antti Karttunen, Jan 18 & Dec 24 2019: (Start)
A085089(a(n)) = n.
A101296(a(n)) = n [which is the first occurrence of n in A101296, and thus also a record.]
A001221(a(n)) = A061395(a(n)) = A061394(n).
A007814(a(n)) = A051903(a(n)) = A051282(n).
a(A101296(n)) = A046523(n).
a(A306802(n)) = A002182(n).
a(n) = A108951(A181815(n)) = A329900(A181817(n)).
If A181815(n) is odd, a(n) = A283980(a(A329904(n))), otherwise a(n) = 2*a(A329904(n)).
(End)
Sum_{n>=1} 1/a(n) = Product_{n>=1} 1/(1 - 1/A002110(n)) = A161360. - Amiram Eldar, Oct 20 2020

Offset corrected by Matthew Vandermast, Oct 19 2008

Minor correction by Charles R Greathouse IV, Sep 03 2010

A061395 Let p be the largest prime factor of n; if p is the k-th prime then set a(n) = k; a(1) = 0 by convention.

Original entry on oeis.org

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

Offset: 1

Henry Bottomley, Apr 30 2001

Records occur at the primes. - Robert G. Wilson v, Dec 30 2007
For n > 1: length of n-th row in A067255. - Reinhard Zumkeller, Jun 11 2013
a(n) = the largest part of the partition having Heinz number n. We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product(p_j-th prime, j=1...r) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] we get 2*2*3*7*29 = 2436. Example: a(20) = 3; indeed, the partition having Heinz number 20 = 2*2*5 is [1,1,3]. - Emeric Deutsch, Jun 04 2015

			a(20) = 3 since the largest prime factor of 20 is 5, which is the 3rd prime.

Álvar Ibeas, Table of n, a(n) for n = 1..100000 (first 1000 terms from Harry J. Smith)
Index entries for sequences computed from indices in prime factorization

Cf. A000720, A006530, A055396, A061394, A133674, A243055.

Cf. A029837, A087207.

a061395 = a049084 . a006530  -- Reinhard Zumkeller, Jun 11 2013

with(numtheory):
a:= n-> pi(max(1, factorset(n)[])):
seq(a(n), n=1..100);  # Alois P. Heinz, Aug 03 2013

Insert[Table[PrimePi[FactorInteger[n][[ -1]][[1]]], {n, 2, 120}], 0, 1] (* Stefan Steinerberger, Apr 11 2006 *)
f[n_] := PrimePi[ FactorInteger@n][[ -1, 1]]; Array[f, 94] (* Robert G. Wilson v, Dec 30 2007 *)

a(n) = if (n==1, 0, primepi(vecmax(factor(n)[,1]))); \\ Michel Marcus, Nov 14 2022

from sympy import primepi, primefactors
def a(n): return 0 if n==1 else primepi(primefactors(n)[-1])
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, May 14 2017

A000040(a(n)) = A006530(n); a(n) = A049084(A006530(n)). - Reinhard Zumkeller, May 22 2003
A243055(n) = a(n) - A055396(n). - Antti Karttunen, Mar 07 2017
a(n) = A000720(A006530(n)). - Alois P. Heinz, Mar 05 2020
a(n) = A029837(A087207(n)+1). - Flávio V. Fernandes, Apr 24 2025

Definition reworded by N. J. A. Sloane, Jul 01 2008

A124832 Table of exponents of prime factorizations in A025487.

Original entry on oeis.org

1, 2, 1, 1, 3, 2, 1, 4, 3, 1, 1, 1, 1, 5, 2, 2, 4, 1, 2, 1, 1, 6, 3, 2, 5, 1, 3, 1, 1, 7, 4, 2, 2, 2, 1, 6, 1, 1, 1, 1, 1, 3, 3, 4, 1, 1, 8, 5, 2, 3, 2, 1, 7, 1, 2, 1, 1, 1, 4, 3, 5, 1, 1, 9, 6, 2, 4, 2, 1, 8, 1, 3, 1, 1, 1, 5, 3, 2, 2, 2, 6, 1, 1, 10, 3, 3, 1, 7, 2, 2, 2, 1, 1, 4, 4, 5, 2, 1, 9, 1, 4, 1, 1, 1, 6

Offset: 2

Franklin T. Adams-Watters, Nov 09 2006

This is an enumeration of all partitions.

			From _M. F. Hasler_, Oct 12 2018: (Start)
The table starts as follows:
  n : signature   (A025487(n) = factorization)
  1 : []          (1 = empty product)
  2 : [1]         (2 = 2^1)
  3 : [2]         (4 = 2^2)
  4 : [1, 1]      (6 = 2^1 * 3^1)
  5 : [3]         (8 = 2^3)
  6 : [2, 1]      (12 = 2^2 * 3^1)
  7 : [4]         (16 = 2^4)
  8 : [3,1]       (24 = 2^3 * 3^1)
  9 : [1, 1, 1]   (30 = 2^1 * 3^1 * 5^1)
  etc. (End)

Michael De Vlieger, Table of n, a(n) for n = 2..10820 (rows 2 <= n <= 2500, first 106 terms from Ray Chandler).

Cf. A025487, A036041 (row sums), A061394 (row lengths), A124829, A036036, A080577.

Map[FactorInteger[#][[All, -1]] &, Import["https://oeis.org/A025487/b025487.txt", "Data"][[2 ;; 48, -1]] ] // Flatten (* Michael De Vlieger, Feb 06 2020 *)

A124832_row(n)=factor(A025487(n))[,2] \\ M. F. Hasler, Oct 12 2018

A025487(n) = Product_{k=1..A061394(n)} prime(k)^T(n,k). [Edited by M. F. Hasler, Oct 12 2018]

Erroneous explanations in cross-references corrected by M. F. Hasler, Oct 12 2018

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

Matthew Vandermast, Nov 11 2008

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.

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

Anonymous?, Polynomial calculator
Eric Weisstein's World of Mathematics, Roundness
G. Xiao, WIMS server, Factoris (both expands and factors polynomials)

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.

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

A304886 Irregular triangle where row n contains indices k where the product of A002110(k) = A025487(n).

Original entry on oeis.org

0, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 3, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 3, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 2, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 3, 1, 1, 1, 1, 1, 2, 4, 2, 2, 2, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2

Offset: 1

Michael De Vlieger, May 21 2018

Row n consists of terms k such that A025487(n) = the product of primorials p_k#, the k in row n written least to greatest k.
For m = A025487(n) in A000079 (i.e., m is an integer power of 2), row n contains A000079(m) 1s.
For m = A025487(n) in A002110 (i.e., m is a primorial) row n contains a single term k that is the index of m in A002110.

			Triangle begins as in rightmost column, which lists the terms that occur on row n. Maximum value of each row is given by A061394(n).
   n  A025487(n)   Row n
--------------------------------
   1        1      0
   2        2      1
   3        4      1,1
   4        6      2
   5        8      1,1,1
   6       12      1,2
   7       16      1,1,1,1
   8       24      1,1,2
   9       30      3
  10       32      1,1,1,1,1
  11       36      2,2
  12       48      1,1,1,2
  13       60      1,3
  14       64      1,1,1,1,1,1
  15       72      1,2,2
  16       96      1,1,1,1,2
  17      120      1,1,3
  18      128      1,1,1,1,1,1,1
  19      144      1,1,2,2
  20      180      2,3
  ...

Michael De Vlieger, Table of n, a(n) for n = 1..8600
Michael De Vlieger, Concordance of A025487, A051282, A061394, and A304886
Michael De Vlieger, Indices of primorials whose product is highly composite
Michael De Vlieger, Indices of primorials whose product is superabundant

Cf. A025487, A051282 (row lengths), A061394 (row maximum), A124832, A181815.

Cf. also A307056.

(* Simple (A025487(n) < 10^5): *)
{{0}}~Join~Map[With[{w = #}, Reverse@ Array[Function[k, Count[w, _?(# >= k &)] ], Max@ w]] &, Select[Array[{#, FactorInteger[#][[All, -1]]} &, 400], Times @@ Boole@ {#1 == Times @@ MapIndexed[Prime[First@ #2]^#1 &, #3], #2 == #3} == 1 & @@ {#1, #2, Sort[#2, Greater]} & @@ # &][[All, -1]] ]
(* Efficient (A025487(n) < 10^23): *)
f[n_] := Block[{ww, g, h},
  g[x_] := Apply[Times,
    MapIndexed[Prime[First@ #2]^#1 &, x]];
  h[x_] := Reverse@
    Array[Function[k, Count[x, _?(# >= k &)] ], Max@ x];
  ww = NestList[Append[#, 1] &, {1}, # - 1] &[-2 +
     Length@ NestWhileList[NextPrime@ # &, 1,
     Times @@ {##} <= n &, All] ];
  Map[h, SortBy[Flatten[#, 1], g]] &@
   Map[Block[{w = #, k = 1},
      Apply[
         Join, {{ConstantArray[1, Length@ w]},
           If[Length@ # == 0, #, #[[1]]] }] &@ Reap[
         Do[
          If[# < n,
            Sow[w]; k = 1,
             If[k >= Length@ w, Break[], k++]] &@
               g@ Set[w,
               If[k == 1,
                 MapAt[# + 1 &, w, k],
                 PadLeft[#, Length@ w, First@ #] &@
                   Drop[MapAt[# + Boole[i > 1] &, w, k],
                    k - 1] ]], {i, Infinity}] ][[-1]] ] &, ww]]; {{0}}~Join~f@ 400

For row n > 1, Product_{k=1..A051282(n)} A000040(T(n,k)) = A181815(n). [Product of primes indexed by nonzero terms of row n is equal to A181815(n)] - Antti Karttunen, Dec 28 2019

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

A025487 Least integer of each prime signature A124832; also products of primorial numbers A002110.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

PARI

PARI

Sage

Formula

Extensions

A061395 Let p be the largest prime factor of n; if p is the k-th prime then set a(n) = k; a(1) = 0 by convention.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

Extensions

A124832 Table of exponents of prime factorizations in A025487.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

A304886 Irregular triangle where row n contains indices k where the product of A002110(k) = A025487(n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A124830 Number of distinct prime factors of A055932(n).