A181555 -id:A181555 - cp's OEIS frontend

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

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

A181796 a(n) = number of divisors of n whose canonical prime factorizations contain no repeated positive exponents (cf. A130091).

Original entry on oeis.org

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

Offset: 1

Matthew Vandermast, Nov 22 2010

nonn

The canonical factorization of n into prime powers can be written as Product p(i)^e(i), for example. A host of equivalent notations can also be used (for another example, see Weisstein link). a(n) depends only on prime signature of n (cf. A025487).
a(n) >= A085082(n). (A085082(n) equals the number of members of A025487 that divide A046523(n), and each member of A025487 is divisible by at least one member of A130091 that divides no smaller member of A025487.) a(n) > A085082(n) iff n has in its canonical prime factorization at least two exponents greater than 1.
a(n) = number of such divisors of n that in their prime factorization all exponents are unique. - Antti Karttunen, May 27 2017
First differs from A335549 at a(90) = 7, A335549(90) = 8. First differs from A335516 at a(180) = 9, A335516(180) = 10. - Gus Wiseman, Jun 28 2020

			12 has a total of six divisors (1, 2, 3, 4, 6 and 12). Of those divisors, the number 1 has no prime factors, hence, no positive exponents at all (and no repeated positive exponents) in its canonical prime factorization. The lists of positive exponents for 2, 3, 4, 6 and 12 are (1), (1), (2), (1,1) and (2,1) respectively (cf. A124010). Of all six divisors, only the number 6 (2^1*3^1) has at least one positive exponent repeated (namely, 1). The other five do not; hence, a(12) = 5.
For n = 90 = 2 * 3^2 * 5, the divisors that satisfy the condition are: 1, 2, 3, 3^2, 5, 2 * 3^2, 3^2 * 5, altogether 7, (but for example 90 itself is not included), thus a(90) = 7.

Antti Karttunen, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Prime Factorization.
Gus Wiseman, Sequences counting and encoding certain classes of multisets

Diverges from A088873 at n=24 and from A085082 at n=36. a(36) = 7, while A085082(36) = 6.
Partitions with distinct multiplicities are A098859.
Sorted prime signature is A118914.
Unsorted prime signature is A124010.
a(n) is the number of divisors of n in A130091.
Factorizations with distinct multiplicities are A255231.
The largest of the counted divisors is A327498.
Factorizations using the counted divisors are A327523.
Cf. A000005, A007916, A045778, A056239, A212168, A327498, A335519, A335549.

Table[DivisorSum[n, 1 &, Length@ Union@ # == Length@ # &@ FactorInteger[#][[All, -1]] &], {n, 105}] (* Michael De Vlieger, May 28 2017 *)

no_repeated_exponents(n) = { my(es = factor(n)[, 2]); if(length(Set(es)) == length(es),1,0); }
A181796(n) = sumdiv(n,d,no_repeated_exponents(d)); \\ Antti Karttunen, May 27 2017

from sympy import factorint, divisors
def ok(n):
    f=factorint(n)
    ex=[f[i] for i in f]
    for i in ex:
        if ex.count(i)>1: return 0
    return 1
def a(n): return sum([1 for i in divisors(n) if ok(i)]) # Indranil Ghosh, May 27 2017

a(A000079(n)) = a(A002110(n)) = n+1.
a(A006939(n)) = A000110(n+1).
a(A181555(n)) = A002720(n).

A212166 Numbers k such that the maximum exponent in its prime factorization equals the number of positive exponents (A051903(k) = A001221(k)).

Original entry on oeis.org

1, 2, 3, 5, 7, 11, 12, 13, 17, 18, 19, 20, 23, 28, 29, 31, 36, 37, 41, 43, 44, 45, 47, 50, 52, 53, 59, 61, 63, 67, 68, 71, 73, 75, 76, 79, 83, 89, 92, 97, 98, 99, 100, 101, 103, 107, 109, 113, 116, 117, 120, 124, 127, 131, 137, 139, 147, 148, 149, 151, 153

Offset: 1

Matthew Vandermast, May 22 2012

			36 = 2^2*3^2 has 2 positive exponents in its prime factorization. The maximal exponent in its prime factorization is also 2. Therefore, 36 belongs to this sequence.

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 844.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
Primefan, The First 2500 Integers Factored (first of 5 pages).

Includes subsequences A000040, A006939, A138534, A181555, A181825.
See also A212164, A212165, A212167, A212168.
Cf. A001221, A050326, A051903, A188654 (complement), A225230.

import Data.List (elemIndices)
a212166 n = a212166_list !! (n-1)
a212166_list = map (+ 1) $ elemIndices 0 a225230_list
-- Reinhard Zumkeller, May 03 2013

okQ[n_] := Module[{f = Transpose[FactorInteger[n]][[2]]}, Max[f] == Length[f]]; Select[Range[424], okQ] (* T. D. Noe, May 24 2012 *)

is(k) = {my(e = factor(k)[, 2]); !(#e) || vecmax(e) == #e;} \\ Amiram Eldar, Sep 08 2024

A225230(a(n)) = 0; A050326(a(n)) = 1. - Reinhard Zumkeller, May 03 2013

A181825 Members of A025487 whose prime signature is self-conjugate (as a partition).

Original entry on oeis.org

1, 2, 12, 36, 120, 360, 1680, 5040, 5400, 27000, 36960, 75600, 110880, 378000, 960960, 1587600, 1663200, 2882880, 7938000, 8316000, 32672640, 34927200, 43243200, 98017920, 174636000, 216216000, 277830000, 908107200, 1152597600, 1241560320, 1470268800, 1944810000

Offset: 1

Matthew Vandermast, Dec 08 2010

A025487(n) is included iff A025487(n) = A181822(n).

Closed under the binary operations of GCD and LCM, since a self-conjugate partition of Omega(a(n)) (which the prime signature of these numbers is) is the concatenation of self-conjugate hooks of decreasing size while moving downward and to the right in the Ferrers diagram, and the GCD (or LCM) of two terms a(i) and a(j) is obtained by taking the smaller (or larger, respectively) of the corresponding hooks. For example, GCD(a(8),a(11)) = GCD(5040,36960) = 1680 = a(7), and LCM(a(8),a(11)) = 110880 = a(13). The two binary operations make the set {a(n)} into a lattice order. - Richard Peterson, May 29 2020

			A025487(11) = 36 = 2^2*3^2 has a prime signature of (2,2), which is a self-conjugate partition; hence, 36 is included in the sequence.

David A. Corneth, Table of n, a(n) for n = 1..11879 (first 578 terms from Amiram Eldar, terms <= 10^70)
David A. Corneth, PARI-program
Eric Weisstein's World of Mathematics, Self-Conjugate Partition

Cf. A025487, A181822, A303557.
Includes subsequences A006939 and A181555.
See also A181823, A181824, A181826, A181827.

\\ See Corneth link \\ David A. Corneth, Jun 03 2020

a(18)-a(32) from Amiram Eldar, Jan 19 2019

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

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

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

			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

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened
C. Chauve, S. Dulucq and O. Guibert, Enumeration of some labeled trees, Proceedings of FPSAC/SFCA 2000 (Moscow), Springer, pp. 146-157.
Alan D. Sokal, A remark on the enumeration of rooted labeled trees, arXiv:1910.14519 [math.CO], 2019 and Discrete Math. 343, 111865 (2020).

Cf. A000169, A000312, A089466, A088956, A166900.

T:= (n, k)-> binomial(n, k)*n^(n-k): seq(seq(T(n, k), k=0..n), n=0..10);

Prepend[Flatten[ Table[Table[Binomial[n, k] n^(n - k), {k, 0, n}], {n, 1, 8}]], 1]  (* Geoffrey Critzer, Jan 08 2012 *)

T(n,k)=if(k<0 || k>n,0,binomial(n,k)*n^(n-k))

/* 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

T(n,k) = binomial(n, k)*n^(n-k).

E.g.f.: (-LambertW(-y)/y)^x/(1+LambertW(-y)). - Vladeta Jovovic

Name edited by Alan Sokal, Jul 22 2022

A181567 Triangle read by rows: T(n,k) is coefficient of k-th power in expansion of ((x^(n+1)-1)/(x-1))^n.

Original entry on oeis.org

1, 1, 1, 1, 2, 3, 2, 1, 1, 3, 6, 10, 12, 12, 10, 6, 3, 1, 1, 4, 10, 20, 35, 52, 68, 80, 85, 80, 68, 52, 35, 20, 10, 4, 1, 1, 5, 15, 35, 70, 126, 205, 305, 420, 540, 651, 735, 780, 780, 735, 651, 540, 420, 305, 205, 126, 70, 35, 15, 5, 1, 1, 6, 21, 56, 126, 252, 462, 786, 1251

Offset: 0

Matthew Vandermast, Oct 31 2010

In each row n>=0, k takes values from 0 to n^2 inclusive. Row sums equal A000169(n+1). All rows are palindromic. Row n is also row n of the (n+1)-nomial array (e.g., row 1 is also row 1 of A007318).
T(n,k) gives the number of divisors of A181555(n) with k prime factors counted with multiplicity. See also A001222, A071207, A146291, A146292.
T(n,k) is the number of size k submultisets of the so-called regular multiset {1_1,1_2,...,1_(n-1),1_n, ... ,i_1,i_2,...,i_(n-1),i_n, ... ,n_1,n_2,...,n_(n-1),n_n} (which contains n copies of i for 0 < i < n). - Thomas Wieder, Dec 28 2013

			Rows begin:
1;
1,1;
1,2,3,2,1;
1,3,6,10,12,12,10,6,3,1;...
T(n=3,k=4) = 12 because we have 12 submultisets (without regard of the order of elements) of size k=4 for the regular multiset (n=3) {1, 1, 1, 2, 2, 2, 3, 3, 3}: {1, 1, 1, 2}, {1, 1, 1, 3}, {1, 1, 2, 2}, {1, 1, 2, 3}, {1, 1, 3, 3}, {1, 2, 2, 2}, {1, 2, 2, 3}, {1, 2, 3, 3}, {1, 3, 3, 3}, {2, 2, 2, 3}, {2, 2, 3, 3}, {2, 3, 3, 3}.

Alois P. Heinz, Rows n = 0..32, flattened
Anonymous, Polynomial calculator
Thomas Wieder, A181567 as Excel table
G. Xiao, WIMS server, Factoris (both expands and factors polynomials)

A163181 gives row n of n-nomial array. See also A000012, A007318, A027907, A008287, A035343, A063260, A063265, A171890.

b:= proc(n, k, i) option remember; `if`(k=0, 1,
     `if`(i<1, 0, add(b(n, k-j, i-1), j=0..n)))
    end:
T:= (n, k)-> b(n, k, n):
seq(seq(T(n, k), k=0..n^2), n=0..8); # Alois P. Heinz, Jul 04 2016

row[n_] := CoefficientList[((x^(n+1) - 1)/(x-1))^n + O[x]^(n^2+1), x]; Table[row[n], {n, 0, 6}] // Flatten (* Jean-François Alcover, Apr 06 2017 *)

A079474 Triangular array: for s=0 to r-1, a(r,s) = p(s)^(r-s), where p(s) is the s-th primorial number. (p(0)=1, p(1)=2, p(2)=23, p(3)=23*5,...).

Original entry on oeis.org

1, 1, 2, 1, 4, 6, 1, 8, 36, 30, 1, 16, 216, 900, 210, 1, 32, 1296, 27000, 44100, 2310, 1, 64, 7776, 810000, 9261000, 5336100, 30030, 1, 128, 46656, 24300000, 1944810000, 12326391000, 901800900, 510510, 1, 256, 279936, 729000000, 408410100000

Offset: 1

Alford Arnold, Jan 15 2003

In the expansion of [1+x+x^2+...+x^(r-s)]^s, the x^n coefficient states how many factors of a(r,s) have n prime factors.

As a square array A(n,k) n>=0 k>=1 read by descending antidiagonals, A(n,k) when n>=1 is the least common period over the positive integers of the occurrence of the first n prime numbers as the k-th least operand in the respective integers' prime factorizations (written without exponents). - Peter Munn, Jan 25 2017

			Triangle starts
  1;
  1,  2;
  1,  4,    6;
  1,  8,   36,    30;
  1, 16,  216,   900,   210;
  1, 32, 1296, 27000, 44100, 2310;
  ...

Cf. A002110, A007318, A008287, A027907, A036035, A061742, A074139.

a(2n,n) gives A181555.

p:= proc(n) option remember; `if`(n=0, 1, ithprime(n)*p(n-1)) end:
a:= (r, s)-> p(s)^(r-s):
seq(seq(a(r, s), s=0..r-1), r=0..10);  # Alois P. Heinz, Aug 22 2019

p[0] = 1; p[s_] := p[s] = Prime[s] p[s-1];
a[r_, s_] := p[s]^(r-s);
Table[a[r, s], {r, 0, 10}, {s, 0, r-1}] // Flatten (* Jean-François Alcover, Dec 07 2019 *)

Edited by Don Reble, Nov 02 2005

A263925 a(n) = least m > 1 such that m + (prime(n)#)^n is prime.

Original entry on oeis.org

3, 5, 11, 19, 89, 323, 29, 61, 79, 199, 563, 181, 353, 1307, 257, 709, 1237, 1277, 1609, 1237, 4157, 2017, 577, 157, 191, 1063, 239, 823, 1607, 4159, 139, 11527, 2339, 18457, 4079, 463, 1861, 1123, 8699, 16561, 719, 4327, 9311, 1693, 3067, 4243, 22397, 4079, 3989, 24071

Offset: 1

Alexei Kourbatov, Oct 30 2015

nonn

Here prime(n)# denotes the primorial A002110(n), i.e., the product of the first n primes. Terms a(n) are often (but not always) prime; out of the first fifty terms, only one (a(6)=323) is composite.

The definition is similar to Fortunate numbers (A005235); however, in A005235 the primorial is not raised to the n-th power. Unlike this sequence, all known Fortunate numbers are prime.

			(prime(2)#)^2=36. a(2)=5 because 5 is the minimal m>1 such that m+36 is prime.

Cf. A002110, A005235, A181555, A264050.

Table[m = 2; While[! PrimeQ[m + Product[Prime@ i, {i, n}]^n], m++]; m, {n, 30}] (* Michael De Vlieger, Nov 11 2015 *)

a(n)=my(s=prod(i=1,n,prime(i))^n); nextprime(s+2)-s

A342455 The fifth powers of primorials: a(n) = A002110(n)^5.

Original entry on oeis.org

1, 32, 7776, 24300000, 408410100000, 65774855015100000, 24421743243121524300000, 34675383095948798128025100000, 85859681408495723096004822084900000, 552622359415801587878908964592391520700000, 11334919554709059323420895730190266747414284300000, 324509123504618420438174660414872405442002404781629300000

Offset: 0

Antti Karttunen, Mar 12 2021

nonn

The ratio G(n) = sigma(n) / (exp(gamma)*n*log(log(n))), where gamma is the Euler-Mascheroni constant (A001620), as applied to these numbers from a(1)=32 onward, develops as:
0.8893323133
0.7551575418
0.7303870617
0.7347890824
0.7263701246
0.7298051649
0.7304358358
0.7354921494
0.7389343933
0.7391912616
0.7416291350
0.7424159544
...
Notably, after its minimum at term a(5) = 65774855015100000, it starts increasing again, albeit rather slowly. At n=10000 the ratio is 0.8632750..., and at n=40000, it is 0.87545260... Question: Does this trend continue indefinitely? In contrast, for primorials, A002110, the ratio appears to be monotonically decreasing, see comments in A342000.

Young Ju Choie, Nicolas Lichiardopol, Pieter Moree and Patrick Solé, On Robin's criterion for the Riemann hypothesis, Journal de théorie des nombres de Bordeaux, 19 no. 2 (2007), pp. 357-372.
Index entries for sequences related to primorial numbers

Cf. A000584, A001620, A002110, A073004, A181555, A342000.

Diagonal in A079474. After the initial term, also the leftmost branch in that subtree of A329886 whose root is 32.

FoldList[Times, 1, Prime@ Range[11]]^5 (* Michael De Vlieger, Mar 14 2021 *)

PARI
```
A342455(n) = prod(i=1,n,prime(i))^5;
```

from sympy.ntheory.generate import primorial
def A342455(n): return primorial(n)**5 if n >= 1 else 1 # Chai Wah Wu, Mar 13 2021

a(n) = A000584(A002110(n)) = A002110(n)^5.

cp's OEIS Frontend

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

A181796 a(n) = number of divisors of n whose canonical prime factorizations contain no repeated positive exponents (cf. A130091).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A212166 Numbers k such that the maximum exponent in its prime factorization equals the number of positive exponents (A051903(k) = A001221(k)).

Views

Author

Keywords

Examples

References

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Formula

A181825 Members of A025487 whose prime signature is self-conjugate (as a partition).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Extensions

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.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula

Extensions

A181567 Triangle read by rows: T(n,k) is coefficient of k-th power in expansion of ((x^(n+1)-1)/(x-1))^n.

Views

Author

Keywords

Comments

Examples

Links

A079474 Triangular array: for s=0 to r-1, a(r,s) = p(s)^(r-s), where p(s) is the s-th primorial number. (p(0)=1, p(1)=2, p(2)=23, p(3)=23*5,...).