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 657 results. Next

A035310 Let f(n) = number of ways to factor n = A001055(n); a(n) = sum of f(k) over all terms k in A025487 that have n factors.

Original entry on oeis.org

1, 4, 12, 47, 170, 750, 3255, 16010, 81199, 448156, 2579626, 15913058, 102488024, 698976419, 4976098729, 37195337408, 289517846210, 2352125666883, 19841666995265, 173888579505200, 1577888354510786, 14820132616197925, 143746389756336173, 1438846957477988926
Offset: 1

Views

Author

Keywords

Comments

Ways of partitioning an n-multiset with multiplicities some partition of n.
Number of multiset partitions of strongly normal multisets of size n, where a finite multiset is strongly normal if it covers an initial interval of positive integers with weakly decreasing multiplicities. The (weakly) normal version is A255906. - Gus Wiseman, Dec 31 2019

Examples

			a(3) = 12 because there are 3 terms in A025487 with 3 factors, namely 8, 12, 30; and f(8)=3, f(12)=4, f(30)=5 and 3+4+5 = 12.
From _Gus Wiseman_, Dec 31 2019: (Start)
The a(1) = 1 through a(3) = 12 multiset partitions of strongly normal multisets:
  {{1}}  {{1,1}}    {{1,1,1}}
         {{1,2}}    {{1,1,2}}
         {{1},{1}}  {{1,2,3}}
         {{1},{2}}  {{1},{1,1}}
                    {{1},{1,2}}
                    {{1},{2,3}}
                    {{2},{1,1}}
                    {{2},{1,3}}
                    {{3},{1,2}}
                    {{1},{1},{1}}
                    {{1},{1},{2}}
                    {{1},{2},{3}}
(End)
		

Crossrefs

Sequence A035341 counts the ordered cases. Tables A093936 and A095705 distribute the values; e.g. 81199 = 30 + 536 + 3036 + 6181 + 10726 + 11913 + 14548 + 13082 + 21147.
Row sums of A317449.
The uniform case is A317584.
The case with empty intersection is A317755.
The strict case is A317775.
The constant case is A047968.
The set-system case is A318402.
The case of strict parts is A330783.
Multiset partitions of integer partitions are A001970.
Unlabeled multiset partitions are A007716.

Programs

  • Maple
    with(numtheory):
    g:= proc(n, k) option remember;
          `if`(n>k, 0, 1) +`if`(isprime(n), 0,
          add(`if`(d>k, 0, g(n/d, d)), d=divisors(n) minus {1, n}))
        end:
    b:= proc(n, i, l)
          `if`(n=0, g(mul(ithprime(t)^l[t], t=1..nops(l))$2),
          `if`(i<1, 0, add(b(n-i*j, i-1, [l[], i$j]), j=0..n/i)))
        end:
    a:= n-> b(n$2, []):
    seq(a(n), n=1..10);  # Alois P. Heinz, May 26 2013
  • Mathematica
    g[n_, k_] := g[n, k] = If[n > k, 0, 1] + If[PrimeQ[n], 0, Sum[If[d > k, 0, g[n/d, d]], {d, Divisors[n] ~Complement~ {1, n}}]]; b[n_, i_, l_] := If[n == 0, g[p = Product[Prime[t]^l[[t]], {t, 1, Length[l]}], p], If[i < 1, 0, Sum[b[n - i*j, i-1, Join[l, Array[i&, j]]], {j, 0, n/i}]]]; a[n_] := b[n, n, {}]; Table[Print[an = a[n]]; an, {n, 1, 13}] (* Jean-François Alcover, Dec 12 2013, after Alois P. Heinz *)
  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
    D(p, n)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); my(u=EulerT(v)); Vec(1/prod(k=1, n, 1 - u[k]*x^k + O(x*x^n))-1, -n)/prod(i=1, #v, i^v[i]*v[i]!)}
    seq(n)={my(s=0); forpart(p=n, s+=D(p,n)); s} \\ Andrew Howroyd, Dec 30 2020
  • Python
    from sympy.core.cache import cacheit
    from sympy import divisors, isprime, prime
    from operator import mul
    @cacheit
    def g(n, k):
        return (0 if n > k else 1) + (0 if isprime(n) else sum(g(n//d, d) for d in divisors(n)[1:-1] if d <= k))
    @cacheit
    def b(n, i, l):
        if n==0:
            p = reduce(mul, (prime(t + 1)**l[t] for t in range(len(l))))
            return g(p, p)
        else:
            return 0 if i<1 else sum([b(n - i*j, i - 1, l + [i]*j) for j in range(n//i + 1)])
    def a(n):
        return b(n, n, [])
    for n in range(1, 11): print(a(n)) # Indranil Ghosh, Aug 19 2017, after Maple code
    

Extensions

More terms from Erich Friedman.
81199 from Alford Arnold, Mar 04 2008
a(10) from Alford Arnold, Mar 31 2008
a(10) corrected by Alford Arnold, Aug 07 2008
a(11)-a(13) from Alois P. Heinz, May 26 2013
a(14) from Alois P. Heinz, Sep 27 2014
a(15) from Alois P. Heinz, Jan 10 2015
Terms a(16) and beyond from Andrew Howroyd, Dec 30 2020

A181815 a(n) = largest integer such that, when any of its divisors divides A025487(n), the quotient is a member of A025487.

Original entry on oeis.org

1, 2, 4, 3, 8, 6, 16, 12, 5, 32, 9, 24, 10, 64, 18, 48, 20, 128, 36, 15, 96, 7, 27, 40, 256, 72, 30, 192, 14, 54, 80, 512, 144, 60, 384, 28, 108, 25, 160, 1024, 45, 288, 21, 81, 120, 768, 56, 216, 50, 320, 2048, 90, 576, 11, 42, 162, 240, 1536, 112, 432, 100, 640, 4096, 180, 1152
Offset: 1

Views

Author

Matthew Vandermast, Nov 30 2010

Keywords

Comments

A permutation of the natural numbers.
The number of divisors of a(n) equals the number of ordered factorizations of A025487(n) as A025487(j)*A025487(k). Cf. A182762.
From Antti Karttunen, Dec 28 2019: (Start)
Rearranges terms of A108951 into ascending order, as A108951(a(n)) = A025487(n).
The scatter plot looks like a curtain of fractal spray, which is typical for many of the so-called "entanglement permutations". Indeed, according to the terminology I use in my 2016-2017 paper, this sequence is obtained by entangling the complementary pair (A329898, A330683) with the complementary pair (A005843, A003961), which gives the following implicit recurrence: a(A329898(n)) = 2*a(n) and a(A330683(n)) = A003961(a(n)). An explicit form is given in the formula section.
(End)

Examples

			For any divisor d of 9 (d = 1, 3, 9), 36/d (36, 12, 4) is a member of A025487. 9 is the largest number with this relationship to 36; therefore, since 36 = A025487(11), a(11) = 9.
		

Crossrefs

If this sequence is considered the "primorial deflation" of A025487(n) (see first formula), the primorial inflation of n is A108951(n), and the primorial inflation of A025487(n) is A181817(n).
A181820(n) is another mapping from the members of A025487 to the positive integers.

Programs

  • Mathematica
    (* First, load the program at A025487, then: *)
    Map[If[OddQ@ #, 1, Times @@ Prime@ # &@ Rest@ NestWhile[Append[#1, {#3, Drop[#, -LengthWhile[Reverse@ #, # == 0 &]] &[#2 - PadRight[ConstantArray[1, #3], Length@ #2]]}] & @@ {#1, #2, LengthWhile[#2, # > 0 &]} & @@ {#, #[[-1, -1]]} &, {{0, TakeWhile[If[# == 1, {0}, Function[g, ReplacePart[Table[0, {PrimePi[g[[-1, 1]]]}], #] &@ Map[PrimePi@ First@ # -> Last@ # &, g]]@ FactorInteger@ #], # > 0 &]}}, And[FreeQ[#[[-1, -1]], 0], Length[#[[-1, -1]] ] != 0] &][[All, 1]] ] &, Union@ Flatten@ f@ 6] (* Michael De Vlieger, Dec 28 2019 *)
  • PARI
    A181815(n) = A329900(A025487(n)); \\ Antti Karttunen, Dec 24 2019

Formula

If A025487(n) is considered in its form as Product A002110(i)^e(i), then a(n) = Product p(i)^e(i). If A025487(n) is instead considered as Product p(i)^e(i), then a(n) = Product (p(i)/A008578(i))^e(i).
a(n) = A122111(A181820(n)). - Matthew Vandermast, May 21 2012
From Antti Karttunen, Dec 24-29 2019: (Start)
a(n) = Product_{i=1..A051282(n)} A000040(A304886(i)).
a(n) = A329900(A025487(n)) = A319626(A025487(n)).
a(n) = A163511(A329905(n)).
For n > 1, if A330682(n) = 1, then a(n) = A003961(a(A329904(n))), otherwise a(n) = 2*a(A329904(n)).
A252464(a(n)) = A329907(n).
A330690(a(n)) = A050378(n).
a(A306802(n)) = A329902(n).
(End)

A181822 a(n) = member of A025487 whose prime signature is conjugate to the prime signature of A025487(n).

Original entry on oeis.org

1, 2, 6, 4, 30, 12, 210, 60, 8, 2310, 36, 420, 24, 30030, 180, 4620, 120, 510510, 1260, 72, 60060, 16, 900, 840, 9699690, 13860, 360, 1021020, 48, 6300, 9240, 223092870, 180180, 2520, 19399380, 240, 69300, 216, 120120, 6469693230, 1800, 3063060, 144, 44100, 27720, 446185740, 1680, 900900, 1080, 2042040, 200560490130, 12600, 58198140, 32, 720
Offset: 1

Views

Author

Matthew Vandermast, Dec 07 2010

Keywords

Comments

A permutation of the members of A025487.
If integers m and n have conjugate prime signatures, then A001222(m) = A001222(n), A071625(m) = A071625(n), A085082(m) = A085082(n), and A181796(m) = A181796(n).

Examples

			A025487(5) = 8 = 2^3 has a prime signature of (3). The partition that is conjugate to (3) is (1,1,1), and the member of A025487 with that prime signature is 30 = 2*3*5 (or 2^1*3^1*5^1).  Therefore, a(5) = 30.
		

Crossrefs

Other rearrangements of A025487 include A036035, A059901, A063008, A077569, A085988, A086141, A087443, A108951, A181821.
A181825 lists members of A025487 with self-conjugate prime signatures. See also A181823-A181824, A181826-A181827.

Programs

  • Mathematica
    f[n_] := Block[{ww, dec}, dec[x_] := Apply[Times, MapIndexed[Prime[First@ #2]^#1 &, x]]; ww = NestList[Append[#, 1] &, {1}, # - 1] &[-2 + Length@ NestWhileList[NextPrime@ # &, 1, Times @@ {##} <= n &, All] ]; {{{0}}}~Join~Map[Block[{w = #, k = 1}, Sort@ Apply[Join, {{ConstantArray[1, Length@ w]}, If[Length@ # == 0, #, #[[1]]] }] &@ Reap[Do[If[# <= n, Sow[w]; 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[Map[{Times @@ MapIndexed[Prime[First@ #2]^#1 &, #], Times @@ MapIndexed[Prime[First@ #2]^#1 &, Table[LengthWhile[#1, # >= j &], {j, #2}]] & @@ {#, Max[#]}} &, Join @@ f[2310]]][[All, -1]] (* Michael De Vlieger, Oct 16 2018 *)
  • PARI
    partitionConj(v)=vector(v[1],i,sum(j=1,#v,v[j]>=i))
    primeSignature(n)=vecsort(factor(n)[,2]~,,4)
    f(n)=if(n==1, return(1)); my(e=partitionConj(primeSignature(n))~); factorback(concat(Mat(primes(#e)~),e))
    A025487=[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];
    concat(1, apply(f, A025487)) \\ Charles R Greathouse IV, Jun 02 2016

Formula

If A025487(n) = Product p(i)^e(i), then a(n) = Product A002110(e(i)). I.e., a(n) = A108951(A181819(A025487(n))). a(n) also equals A108951(A181820(n)).

A098719 Position of prime(n)# in A025487.

Original entry on oeis.org

1, 2, 4, 9, 22, 54, 114, 246, 488, 948, 1809, 3327, 6020, 10624, 18246, 30726, 51148, 84074, 135598, 216398, 340886, 529051, 814237, 1240172, 1874464, 2817289, 4195918, 6186286, 9049492, 13121704, 18895821, 27199504, 38892092, 55318849, 78130780, 110028527
Offset: 1

Views

Author

Jeff Burch, Sep 29 2004

Keywords

Crossrefs

Programs

  • Mathematica
    Block[{nn = 24, f, s}, f[n_] := {{1}}~Join~Block[{lim = Product[Prime@ i, {i, n}], ww = NestList[Append[#, 1] &, {1}, n - 1], g}, g[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++]] &@ 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]]; s = Sort[Join @@ f@ nn]; {1}~Join~Array[Position[s, Product[Prime@ i, {i, #}]] &, nn][[All, 1, 1]]] (* Michael De Vlieger, Jul 23 2018 *)

Extensions

Extended by T. D. Noe, Nov 12 2010
More terms from Michael De Vlieger, Jul 23 2018
Name corrected by Amiram Eldar, Jun 05 2022

A064839 List the natural numbers starting a new row only with each new least prime signature (A025487). a(n) is the column position associated with n.

Original entry on oeis.org

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

Views

Author

Alford Arnold, Oct 24 2001

Keywords

Comments

Row 2 records the primes (A000040). Rows 3 and 4 record the semiprimes (A001358). Rows 5, 6 and 9 record the 3-almost primes (A014612) etc. A058933 is a similar sequence based on k-almost primes.
The graph of this sequence is interesting for large n because it shows multiple curves, one for each prime signature. For example, the six highest curves on the graph of a(n) for n up to 10^4 are for the (1,1), (1,1,1), (1), (2,1,1), (2,1), and (1,1,1,1) prime signatures. The (1) curve dominates until n=58; the (1,1) curve dominates until n=1279786, when the (1,1,1) curve intersects the (1,1) curve. Each (1,1,...,1) curve dominates for a finite number of n.
Ordinal transform of A101296. - Antti Karttunen, May 15 2017
a(n) is the number of positive integers up to n with the same prime signature as n. For example, the a(20) = 3 numbers are {12, 18, 20}. - Gus Wiseman, Jul 08 2019
Ordinal transform of A046523. - Alois P. Heinz, May 31 2020

Examples

			The list begins as follows:
1
2 3 5 7 11 13 17 19 23 29 31 37 41 43 47 53 ...
4 9 25 49 ...
6 10 14 15 21 22 26 33 34 35 38 39 46 51 ...
8 27 ...
12 18 20 28 44 45 50 52 ...
16 ...
Note: the above array, without the initial 1, is given by A095904 (and its transpose A179216). - _Antti Karttunen_, May 15 2017
		

Crossrefs

Programs

  • Maple
    p:= proc() 0 end:
    a:= proc(n) option remember; local t; a(n-1);
          t:= (l-> mul(ithprime(i)^l[i], i=1..nops(l)))(
               sort(map(i-> i[2], ifactors(n)[2]), `>`));
          p(t):= p(t)+1
        end: a(0):=0:
    seq(a(n), n=1..100);  # Alois P. Heinz, May 31 2020
  • Mathematica
    prisig[n_]:=If[n==1,{},Sort[Last/@FactorInteger[n]]];
    Table[Count[Array[prisig,n],prisig[n]],{n,30}] (* Gus Wiseman, Jul 08 2019 *)

Extensions

More terms from Naohiro Nomoto, Oct 31 2001

A056808 Members of A055932 which are not least prime signatures (cf. A025487).

Original entry on oeis.org

18, 54, 90, 108, 150, 162, 270, 300, 324, 450, 486, 540, 600, 630, 648, 750, 810, 972, 1050, 1200, 1350, 1458, 1470, 1500, 1620, 1890, 1944, 2100, 2250, 2400, 2430, 2700, 2916, 2940, 3000, 3150, 3240, 3750, 3780, 3888, 4050, 4200, 4374, 4410, 4500, 4800
Offset: 1

Views

Author

Alford Arnold, Aug 22 2000

Keywords

Examples

			18 = 2*3*3 and all prime divisors are consecutive primes but the least prime signature is 12 = 2*2*3; so a(1) = 18.
		

Crossrefs

Programs

  • Mathematica
    With[{nn = 4800}, Select[Range[2, nn], And[#1 != Times @@ MapIndexed[Prime[First@ #2]^#1 &, Sort[#3, Greater]], Last[#2] == Prime@ Length[#2]] & @@ Apply[Join, {{#1}, Transpose@ #2}] & @@ {#, FactorInteger[#]} &] ] (* Michael De Vlieger, Feb 06 2020 *)

Formula

{a(n) : n >= 1} = {A057335(A335485(k)) : k >= 1}. - Peter Munn, Feb 02 2024
Sum_{n>=1} 1/a(n) = A345974 - A161360 = 0.15229524564163275059... . - Amiram Eldar, Jun 26 2025

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Nov 28 2000

A061394 Number of distinct prime factors of n-th least prime signature (A025487); also a(n)-th prime is largest prime factor of n-th least prime signature; also a(n)-th primorial number is largest primorial factor of n-th least product of primorial numbers.

Original entry on oeis.org

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

Views

Author

Henry Bottomley, Apr 30 2001

Keywords

Comments

A002110(a(n)) = A247451(n). - Reinhard Zumkeller, Sep 17 2014
Number of parts of the associated prime signature. - Álvar Ibeas, Nov 01 2014

Crossrefs

Cf. A002110, A247451, A006530, A061395, A025487, A000040, A051903. A001221 by prime signature.

Programs

  • Haskell
    import Data.List (elemIndex); import Data.Maybe (fromJust)
    a061394 = fromJust . (`elemIndex` a002110_list) . a247451
    -- Reinhard Zumkeller, Sep 17 2014
    
  • PARI
    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)))
    [omega(n) | n <- [1..1000], isA025487(n)]
    \\ Or, for older versions:
    apply(omega, select(isA025487, [1..1000])) \\ Charles R Greathouse IV, Nov 07 2014

Formula

a(n) = A061395(A025487(n)) = A001221(A025487(n)) = A051903(A181822(n)).
A000040(a(n)) = A006530(A025487(n)).

Extensions

Offset updated by Matthew Vandermast, Nov 08 2008

A322827 A permutation of A025487: Sequence of least representatives of distinct prime signatures obtained from the run lengths present in the binary expansion of n.

Original entry on oeis.org

1, 2, 6, 4, 36, 30, 12, 8, 216, 180, 210, 900, 72, 60, 24, 16, 1296, 1080, 1260, 5400, 44100, 2310, 6300, 27000, 432, 360, 420, 1800, 144, 120, 48, 32, 7776, 6480, 7560, 32400, 264600, 13860, 37800, 162000, 9261000, 485100, 30030, 5336100, 1323000, 69300, 189000, 810000, 2592, 2160, 2520, 10800, 88200, 4620, 12600
Offset: 0

Views

Author

Antti Karttunen, Jan 16 2019

Keywords

Comments

A101296(a(n)) gives a permutation of natural numbers.

Examples

			The sequence can be represented as a binary tree:
                                      1
                                      |
                   ...................2...................
                  6                                       4
       36......../ \........30                 12......../ \........8
       / \                 / \                 / \                 / \
      /   \               /   \               /   \               /   \
     /     \             /     \             /     \             /     \
   216      180         210    900         72       60         24       16
etc.
Both children are multiples of their common parent, see A323503, A323504 and A323507.
The value of a(n) is computed from the binary expansion of n as follows: Starting from the least significant end of the binary expansion of n (A007088), we record the successive run lengths, subtract one from all lengths except the first one, and use the reversed partial sums of these adjusted values as the exponents of successive primes.
For 11, which is "1011" in base 2, we have run lengths [2, 1, 1] when scanned from the right, and when one is subtracted from all except the first, we have [2, 0, 0], partial sums of which is [2, 2, 2], which stays same when reversed, thus a(11) = 2^2 * 3^2 * 5^2 = 900.
For 13, which is "1101" in base 2, we have run lengths [1, 1, 2] when scanned from the right, and when one is subtracted from all except the first, we have [1, 0, 1], partial sums of which is [1, 1, 2], reversed [2, 1, 1], thus a(13) = 2^2 * 3^1 * 5^1 = 60.
Sequence A227183 is based on the same algorithm.
		

Crossrefs

Cf. A000079 (right edge), A000400 (left edge, apart from 2), A005811, A046523, A101296, A227183, A322585, A322825, A323503, A323504, A323507.
Other rearrangements of A025487 include A036035, A059901, A063008, A077569, A085988, A086141, A087443, A108951, A181821, A181822.
Cf. A005940, A283477, A323505 for other similar trees.

Programs

  • Mathematica
    {1}~Join~Array[Times @@ MapIndexed[Prime[First@ #2]^#1 &, Reverse@ Accumulate@ MapIndexed[Length[#1] - Boole[First@ #2 > 1] &, Split@ Reverse@ IntegerDigits[#, 2]]] &, 54] (* Michael De Vlieger, Feb 05 2020 *)
  • PARI
    A322827(n) = if(!n,1,my(bits = Vecrev(binary(n)), rl=1, o = List([])); for(i=2,#bits,if(bits[i]==bits[i-1], rl++, listput(o,rl))); listput(o,rl); my(es=Vecrev(Vec(o)), m=1); for(i=1,#es,m *= prime(i)^es[i]); (m));

Formula

a(n) = A046523(a(n)) = A046523(A322825(n)).
A001221(a(n)) = A005811(n).
A001222(a(n)) = A227183(n).
A322585(a(n)) = 1.

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

Views

Author

Keywords

Comments

This is an enumeration of all partitions.

Examples

			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)
		

Crossrefs

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

Programs

Formula

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

Extensions

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

A181820 a(1) = 1; for n > 1, if A025487(n) = Product p(i)^e(i), then a(n) = Product p(e(i)).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 10, 8, 11, 9, 14, 12, 13, 15, 22, 20, 17, 21, 18, 26, 16, 25, 28, 19, 33, 30, 34, 24, 35, 44, 23, 39, 42, 38, 40, 55, 27, 52, 29, 50, 51, 36, 49, 66, 46, 56, 65, 45, 68, 31, 70, 57, 32, 60, 77, 78, 58, 88, 85, 63, 76, 37, 110, 69, 48, 84, 91, 75, 102, 62, 54, 98, 104, 95
Offset: 1

Views

Author

Matthew Vandermast, Dec 07 2010

Keywords

Comments

A permutation of the positive integers.
The partition given by the prime signature of A025487(n) has Heinz number a(n). - Pontus von Brömssen, Mar 25 2023

Examples

			A025487(8) = 24 = 2^3*3 has the exponents (3,1) in its canonical prime factorization. Accordingly, a(8) = prime(3)*prime(1) (i.e., A000040(3)*A000040(1)), which equals 5*2=10.
		

Crossrefs

A181815 is another mapping from the members of A025487 to the positive integers. Also see A181819, A181821.
Cf. A000040, A122111, A361808 (inverse), A361809 (fixed points).

Formula

a(n) = A181819(A025487(n)).
a(n) = A122111(A181815(n)).
Showing 1-10 of 657 results. Next