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-5 of 5 results.

A329889 a(n) is the unique integer k such that A108951(k) = A260633(n).

Original entry on oeis.org

1, 3, 6, 12, 5, 10, 20, 15, 30, 60, 28, 45, 21, 56, 90, 42, 180, 84, 63, 168, 70, 126, 140, 252, 189, 280, 504, 210, 378, 264, 1008, 420, 315, 220, 840, 630, 1680, 792, 330, 1260, 1584, 945, 1400, 660, 2520, 495, 1890, 882, 1320, 2100, 990, 1764, 2640, 4200, 1980, 3528, 1485, 2200, 8400, 2646, 3960, 6300, 2970, 5292, 7920, 3300, 5940, 2772
Offset: 1

Views

Author

Antti Karttunen, Dec 28 2019

Keywords

Crossrefs

Formula

a(n) = A329900(A260633(n)).

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

Views

Author

Keywords

Comments

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

Examples

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

Crossrefs

Subsequence of A055932, A191743, and A324583.
Cf. A085089, A101296 (left inverses).
Equals range of values taken by A046523.
Cf. A178799 (first differences), A247451 (squarefree kernel), A146288 (number of divisors).
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.

Programs

  • Haskell
    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
    
  • Maple
    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
  • Mathematica
    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 *)
  • 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))) \\ Charles R Greathouse IV, Jun 10 2011
    
  • PARI
    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 */
    
  • PARI
    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
    
  • PARI
    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
    
  • PARI
    \\ 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
    
  • Sage
    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

Formula

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

Extensions

Offset corrected by Matthew Vandermast, Oct 19 2008
Minor correction by Charles R Greathouse IV, Sep 03 2010

A008480 Number of ordered prime factorizations of n.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

a(n) depends only on the prime signature of n (cf. A025487). So a(24) = a(375) since 24 = 2^3 * 3 and 375 = 3 * 5^3 both have prime signature (3,1).
Multinomial coefficients in prime factorization order. - Max Alekseyev, Nov 07 2006
The Dirichlet inverse is given by A080339, negating all but the A080339(1) element in A080339. - R. J. Mathar, Jul 15 2010
Number of (distinct) permutations of the multiset of prime factors. - Joerg Arndt, Feb 17 2015
Number of not divisible chains in the divisor lattice of n. - Peter Luschny, Jun 15 2013

References

  • A. Knopfmacher, J. Knopfmacher, and R. Warlimont, "Ordered factorizations for integers and arithmetical semigroups", in Advances in Number Theory, (Proc. 3rd Conf. Canadian Number Theory Assoc., 1991), Clarendon Press, Oxford, 1993, pp. 151-165.
  • Steven R. Finch, Mathematical Constants, Cambridge, 2003, pp. 292-295.

Crossrefs

Cf. A124010, record values and where they occur: A260987, A260633.
Absolute values of A355939.

Programs

  • Haskell
    a008480 n = foldl div (a000142 $ sum es) (map a000142 es)
                where es = a124010_row n
    -- Reinhard Zumkeller, Nov 18 2015
    
  • Maple
    a:= n-> (l-> add(i, i=l)!/mul(i!, i=l))(map(i-> i[2], ifactors(n)[2])):
    seq(a(n), n=1..100);  # Alois P. Heinz, May 26 2018
  • Mathematica
    Prepend[ Array[ Multinomial @@ Last[ Transpose[ FactorInteger[ # ] ] ]&, 100, 2 ], 1 ]
    (* Second program: *)
    a[n_] := With[{ee = FactorInteger[n][[All, 2]]}, Total[ee]!/Times @@ (ee!)]; Array[a, 101] (* Jean-François Alcover, Sep 15 2019 *)
  • PARI
    a(n)={my(sig=factor(n)[,2]); vecsum(sig)!/vecprod(apply(k->k!, sig))} \\ Andrew Howroyd, Nov 17 2018
    
  • Python
    from math import prod, factorial
    from sympy import factorint
    def A008480(n): return factorial(sum(f:=factorint(n).values()))//prod(map(factorial,f)) # Chai Wah Wu, Aug 05 2023
  • Sage
    def A008480(n):
        S = [s[1] for s in factor(n)]
        return factorial(sum(S)) // prod(factorial(s) for s in S)
    [A008480(n) for n in (1..101)]  # Peter Luschny, Jun 15 2013
    

Formula

If n = Product (p_j^k_j) then a(n) = ( Sum (k_j) )! / Product (k_j !).
Dirichlet g.f.: 1/(1-B(s)) where B(s) is D.g.f. of characteristic function of primes.
a(p^k) = 1 if p is a prime.
a(A002110(n)) = A000142(n) = n!.
a(n) = A050382(A101296(n)). - R. J. Mathar, May 26 2017
a(n) = 1 <=> n in { A000961 }. - Alois P. Heinz, May 26 2018
G.f. A(x) satisfies: A(x) = x + A(x^2) + A(x^3) + A(x^5) + ... + A(x^prime(k)) + ... - Ilya Gutkovskiy, May 10 2019
a(n) = C(k, n) for k = A001222(n) where C(k, n) is defined as the k-fold Dirichlet convolution of A001221(n) with itself, and where C(0, n) is the multiplicative identity with respect to Dirichlet convolution.
The average order of a(n) is asymptotic (up to an absolute constant) to 2A sqrt(2*Pi) log(n) / sqrt(log(log(n))) for some absolute constant A > 0. - Maxie D. Schmidt, May 28 2021
The sums of a(n) for n <= x and k >= 1 such that A001222(n)=k have asymptotic order of the form x*(log(log(x)))^(k+1/2) / ((2k+1) * (k-1)!). - Maxie D. Schmidt, Feb 12 2021
Other DGFs include: (1+P(s))^(-1) in terms of the prime zeta function for Re(s) > 1 where the + version weights the sequence by A008836(n), see the reference by Fröberg on P(s). - Maxie D. Schmidt, Feb 12 2021
The bivariate DGF (1+zP(s))^(-1) has coefficients a(n) / n^s (-1)^(A001221(n)) z^(A001222(n)) for Re(s) > 1 and 0 < |z| < 2 - Maxie D. Schmidt, Feb 12 2021
The distribution of the distinct values of the sequence for n<=x as x->infinity satisfy a CLT-type Erdős-Kac theorem analog proved by M. D. Schmidt, 2021. - Maxie D. Schmidt, Feb 12 2021
a(n) = abs(A355939(n)). - Antti Karttunen and Vaclav Kotesovec, Jul 22 2022
a(n) = A130675(n)/A112624(n). - Amiram Eldar, Mar 08 2024

Extensions

Edited by N. J. A. Sloane at the suggestion of Andrew S. Plewe, Jun 17 2007

A329900 Primorial deflation of n: starting from x = n, repeatedly divide x by the largest primorial A002110(k) that divides it, until x is an odd number. Then a(n) = Product prime(k_i), for primorial indices k_1 >= k_2 >= ..., encountered in the process.

Original entry on oeis.org

1, 2, 1, 4, 1, 3, 1, 8, 1, 2, 1, 6, 1, 2, 1, 16, 1, 3, 1, 4, 1, 2, 1, 12, 1, 2, 1, 4, 1, 5, 1, 32, 1, 2, 1, 9, 1, 2, 1, 8, 1, 3, 1, 4, 1, 2, 1, 24, 1, 2, 1, 4, 1, 3, 1, 8, 1, 2, 1, 10, 1, 2, 1, 64, 1, 3, 1, 4, 1, 2, 1, 18, 1, 2, 1, 4, 1, 3, 1, 16, 1, 2, 1, 6, 1, 2, 1, 8, 1, 5, 1, 4, 1, 2, 1, 48, 1, 2, 1, 4, 1, 3, 1, 8, 1
Offset: 1

Views

Author

Antti Karttunen, Dec 22 2019

Keywords

Comments

When applied to arbitrary n, the "primorial deflation" (term coined by Matthew Vandermast in A181815) induces the splitting of n to two factors A328478(n)*A328479(n) = n, where we call A328478(n) the non-deflatable component of n (which is essentially discarded), while A328479(n) is the deflatable component. Only if n is in A025487, then the entire n is deflatable, i.e., A328478(n) = 1 and A328479(n) = n.
According to Daniel Suteu, also the ratio (A319626(n) / A319627(n)) can be viewed as a "primorial deflation". That definition coincides with this one when restricted to terms of A025487, as for all k in A025487, A319626(k) = a(k), and A319627(k) = 1. - Antti Karttunen, Dec 29 2019

Crossrefs

Programs

  • Mathematica
    Array[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[f, ReplacePart[Table[0, {PrimePi[f[[-1, 1]]]}], #] &@ Map[PrimePi@ First@ # -> Last@ # &, f]]@ FactorInteger@ #], # > 0 &]}}, And[FreeQ[#[[-1, -1]], 0], Length[#[[-1, -1]] ] != 0] &][[All, 1]] ] &, 105] (* Michael De Vlieger, Dec 28 2019 *)
    Array[Times @@ Prime@(TakeWhile[Reap[FixedPointList[Block[{k = 1}, While[Mod[#, Prime@ k] == 0, k++]; Sow[k - 1]; #/Product[Prime@ i, {i, k - 1}]] &, #]][[-1, 1]], # > 0 &]) &, 105] (* Michael De Vlieger, Jan 11 2020 *)
  • PARI
    A329900(n) = { my(m=1, pp=1); while(1, forprime(p=2, ,if(n%p, if(2==p, return(m), break), n /= p; pp = p)); m *= pp); (m); };
    
  • PARI
    A111701(n) = forprime(p=2, , if(n%p, return(n), n /= p));
    A276084(n) = { for(i=1,oo,if(n%prime(i),return(i-1))); }
    A329900(n) = if(n%2,1,prime(A276084(n))*A329900(A111701(n)));

Formula

For odd n, a(n) = 1, for even n, a(n) = A000040(A276084(n)) * a(A111701(n)).
For even n, a(n) = A000040(A276084(n)) * a(n/A002110(A276084(n))).
A108951(a(n)) = A328479(n), for n >= 1.
a(A108951(n)) = n, for n >= 1.
a(A328479(n)) = a(n), for n >= 1.
a(A328478(n)) = 1, for n >= 1.
a(A002110(n)) = A000040(n), for n >= 1.
a(A000142(n)) = A307035(n), for n >= 0.
a(A283477(n)) = A019565(n), for n >= 0.
a(A329886(n)) = A005940(1+n), for n >= 0.
a(A329887(n)) = A163511(n), for n >= 0.
a(A329602(n)) = A329888(n), for n >= 1.
a(A025487(n)) = A181815(n), for n >= 1.
a(A124859(n)) = A181819(n), for n >= 1.
a(A181817(n)) = A025487(n), for n >= 1.
a(A181821(n)) = A122111(n), for n >= 1.
a(A002182(n)) = A329902(n), for n >= 1.
a(A260633(n)) = A329889(n), for n >= 1.
a(A033833(n)) = A330685(n), for n >= 1.
a(A307866(1+n)) = A330686(n), for n >= 1.
a(A330687(n)) = A330689(n), for n >= 1.

A260987 Record values in A008480.

Original entry on oeis.org

1, 2, 3, 4, 6, 12, 20, 30, 60, 105, 120, 140, 180, 210, 280, 420, 504, 840, 1120, 1512, 1680, 2520, 3780, 5040, 6300, 7560, 9240, 12600, 13860, 15120, 15840, 27720, 34650, 37800, 55440, 83160, 102960, 110880, 138600, 180180, 205920, 216216, 240240, 332640
Offset: 1

Views

Author

Reinhard Zumkeller, Nov 18 2015

Keywords

Examples

			.   n | x = A260633(n)        | y = a(n)             | [x/y] | x mod y
. ----+-----------------------+----------------------+-------+---------
.   1 |     1 | 1             |    1 | 1             |    1  |      0
.   2 |     6 | 2*3           |    2 | 2             |    3  |      0
.   3 |    12 | 2^2*3         |    3 | 3             |    4  |      0
.   4 |    24 | 2^3*3         |    4 | 2^2           |    6  |      0
.   5 |    30 | 2*3*5         |    6 | 2*3           |    5  |      0
.   6 |    60 | 2^2*3*5       |   12 | 2^2*3         |    5  |      0
.   7 |   120 | 2^3*3*5       |   20 | 2^2*5         |    6  |      0
.   8 |   180 | 2^2*3^2*5     |   30 | 2*3*5         |    6  |      0
.   9 |   360 | 2^3*3^2*5     |   60 | 2^2*3*5       |    6  |      0
.  10 |   720 | 2^4*3^2*5     |  105 | 3*5*7         |    6  |     90
.  11 |   840 | 2^3*3*5*7     |  120 | 2^3*3*5       |    7  |      0
.  12 |  1080 | 2^3*3^3*5     |  140 | 2^2*5*7       |    7  |    100
.  13 |  1260 | 2^2*3^2*5*7   |  180 | 2^2*3^2*5     |    7  |      0
.  14 |  1680 | 2^4*3*5*7     |  210 | 2*3*5*7       |    8  |      0
.  15 |  2160 | 2^4*3^3*5     |  280 | 2^3*5*7       |    7  |    200
.  16 |  2520 | 2^3*3^2*5*7   |  420 | 2^2*3*5*7     |    6  |      0
.  17 |  4320 | 2^5*3^3*5     |  504 | 2^3*3^2*7     |    8  |    288
.  18 |  5040 | 2^4*3^2*5*7   |  840 | 2^3*3*5*7     |    6  |      0
.  19 |  7560 | 2^3*3^3*5*7   | 1120 | 2^5*5*7       |    6  |    840
.  20 | 10080 | 2^5*3^2*5*7   | 1512 | 2^3*3^3*7     |    6  |   1008
.  21 | 12600 | 2^3*3^2*5^2*7 | 1680 | 2^4*3*5*7     |    7  |    840
.  22 | 15120 | 2^4*3^3*5*7   | 2520 | 2^3*3^2*5*7   |    6  |      0
.  23 | 25200 | 2^4*3^2*5^2*7 | 3780 | 2^2*3^3*5*7   |    6  |   2520
.  24 | 30240 | 2^5*3^3*5*7   | 5040 | 2^4*3^2*5*7   |    6  |      0
.  25 | 45360 | 2^4*3^4*5*7   | 6300 | 2^2*3^2*5^2*7 |    7  |   1260 .
		

Crossrefs

Programs

  • Haskell
    a260987 n = a260987_list !! (n-1)
    (a260987_list, a260633_list) = unzip $ f 1 0 where
       f x r = if y > r then (y, x) : f (x + 1) y else f (x + 1) r
               where y = a008480 x

Formula

a(n) = A008480(A260633(n)).
Showing 1-5 of 5 results.