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.

Previous Showing 51-60 of 225 results. Next

A273013 Number of different arrangements of nonnegative integers on a pair of n-sided dice such that the dice can add to any integer from 0 to n^2-1.

Original entry on oeis.org

1, 1, 1, 3, 1, 7, 1, 10, 3, 7, 1, 42, 1, 7, 7, 35, 1, 42, 1, 42, 7, 7, 1, 230, 3, 7, 10, 42, 1, 115, 1, 126, 7, 7, 7, 393, 1, 7, 7, 230, 1, 115, 1, 42, 42, 7, 1, 1190, 3, 42, 7, 42, 1, 230, 7, 230, 7, 7, 1, 1158, 1, 7, 42, 462, 7, 115, 1, 42, 7, 115, 1, 3030
Offset: 1

Views

Author

Elliott Line, May 13 2016

Keywords

Comments

The set of b values (see formula), and therefore also a(n), depends only on the prime signature of n. So, for example, a(24) will be identical to a(n) of any other n which is also of the form p_1^3*p_2, (e.g., 40, 54, 56).
The value of b_1 will always be 1. When n is prime, the only nonzero b will be b_1, so therefore a(n) will be 1.
In any arrangement, both dice will have a 0, and one will have a 1 (here called the lead die). To determine any one of the actual arrangements to numbers on the dice, select one of the permutations of divisors (for the lead die), then select another permutation that is either the same length as that of the lead die, or one less. For example, if n = 12, we might select 2*3*2 for the lead die, and 3*4 for the second die. These numbers effectively tell you when to "switch track" when numbering the dice, and will uniquely result in the numbering: (0,1,6,7,12,13,72,73,78,79,84,85; 0,2,4,18,20,22,36,38,40,54,56,58).
a(n) is the number of (unordered) pairs of polynomials c(x) = x^c_1 + x^c_2 + ... + x^c_n, d(x) = x^d_1 + x^d_2 + ... + x^d_n with nonnegative integer exponents such that c(x)*d(x) = (x^(n^2)-1)/(x-1). - Alois P. Heinz, May 13 2016
a(n) is also the number of principal reversible squares of order n. - S. Harry White, May 19 2018
From Gus Wiseman, Oct 29 2021: (Start)
Also the number of ordered factorizations of n^2 with alternating product 1. This follows from the author's formula. Taking n instead of n^2 gives a(sqrt(n)) if n is a perfect square, otherwise 0. Here, an ordered factorization of n is a sequence of positive integers > 1 with product n, and the alternating product of a sequence (y_1,...,y_k) is Product_i y_i^((-1)^(i-1)). For example, the a(1) = 1 through a(9) = 3 factorizations are:
() (22) (33) (44) (55) (66) (77) (88) (99)
(242) (263) (284) (393)
(2222) (362) (482) (3333)
(2233) (2244)
(2332) (2442)
(3223) (4224)
(3322) (4422)
(22242)
(24222)
(222222)
The even-length case is A347464.
(End)

Examples

			When n = 4, a(n) = 3; the three arrangements are (0,1,2,3; 0,4,8,12), (0,1,4,5; 0,2,8,10), (0,1,8,9; 0,2,4,6).
When n = 5, a(n) = 1; the sole arrangement is (0,1,2,3,4; 0,5,10,15,20).
		

Crossrefs

Positions of 1's are 1 and A000040.
A000290 lists squares, complement A000037.
A001055 counts factorizations, ordered A074206.
A119620 counts partitions with alternating product 1, ranked by A028982.
A339846 counts even-length factorizations, ordered A174725.
A339890 counts odd-length factorizations, ordered A174726.
A347438 counts factorizations with alternating product 1.
A347460 counts possible alternating products of factorizations.
A347463 counts ordered factorizations with integer alternating product.
A347466 counts factorizations of n^2.

Programs

  • Mathematica
    facs[n_] := If[n <= 1, {{}}, Join@@Table[Map[Prepend[#, d]&, Select[facs[n/d], Min@@# >= d&]], {d, Rest[Divisors[n]]}]];
    altprod[q_] := Product[q[[i]]^(-1)^(i-1), {i, Length[q]}];
    Table[Length[Select[Join@@Permutations/@facs[n^2], altprod[#] == 1&]],{n, 30}]
    (* Gus Wiseman, Oct 29 2021 *)
    (* or *)
    ofc[n_,k_] := If[k > PrimeOmega[n], 0, If[k == 0 && n == 1, 1, Sum[ofc[n/d, k-1],{d, Rest[Divisors[n]]}]]];
    Table[If[n == 1, 1, Sum[ofc[n, x]^2 + ofc[n, x]*ofc[n, x+1], {x, n}]],{n, 30}]
    (* Gus Wiseman, Oct 29 2021, based on author's formula *)
  • PARI
    A273013aux(n, k=0, t=1) = if(1==n, (1==t), my(s=0); fordiv(n, d, if((d>1), s += A273013aux(n/d, 1-k, t*(d^((-1)^k))))); (s));
    A273013(n) = A273013aux(n^2); \\ Antti Karttunen, Oct 30 2021
    
  • SageMath
    @cached_function
    def r(m,n):
        if n==1:
            return(1)
        divList = divisors(m)[:-1]
        return(sum(r(n,d) for d in divList))
    def A273013(n):
        return(r(n,n)) # William P. Orrick, Feb 19 2023

Formula

a(n) = b_1^2 + b_2^2 + b_3^2 + ... + b_1*b_2 + b_2*b_3 + b_3*b_4 + ..., where b_k is the number of different permutations of k divisors of n to achieve a product of n. For example, for n=24, b_3 = 9 (6 permutations of 2*3*4 and 3 permutations of 2*2*6).
a(n) = r(n,n) where r(m,1) = 1 and r(m,n) = Sum_{d|m,dWilliam P. Orrick, Feb 19 2023

A323910 Dirichlet inverse of the deficiency of n, A033879.

Original entry on oeis.org

1, -1, -2, 0, -4, 4, -6, 0, -1, 6, -10, 2, -12, 8, 10, 0, -16, 1, -18, 2, 14, 12, -22, 4, -3, 14, -2, 2, -28, -16, -30, 0, 22, 18, 26, 4, -36, 20, 26, 4, -40, -24, -42, 2, 4, 24, -46, 8, -5, -1, 34, 2, -52, 0, 42, 4, 38, 30, -58, 2, -60, 32, 6, 0, 50, -40, -66, 2, 46, -40, -70, 12, -72, 38, 2, 2, 62, -48, -78, 8, -4, 42, -82, -2, 66, 44, 58, 4, -88, 2, 74, 2
Offset: 1

Views

Author

Antti Karttunen, Feb 12 2019

Keywords

Crossrefs

Cf. A033879, A323911, A323912, A359549 (parity of terms).
Sequences that appear in the convolution formulas: A002033, A008683, A023900, A055615, A046692, A067824, A074206, A174725, A191161, A327960, A328722, A330575, A345182, A349341, A346246, A349387.

Programs

  • Mathematica
    b[n_] := 2 n - DivisorSigma[1, n];
    a[n_] := a[n] = If[n == 1, 1, -Sum[b[n/d] a[d], {d, Most@ Divisors[n]}]];
    Array[a, 100] (* Jean-François Alcover, Feb 17 2020 *)
  • PARI
    up_to = 16384;
    DirInverse(v) = { my(u=vector(#v)); u[1] = (1/v[1]); for(n=2, #v, u[n] = -sumdiv(n, d, if(dA033879(n) = (2*n-sigma(n));
    v323910 = DirInverse(vector(up_to,n,A033879(n)));
    A323910(n) = v323910[n];

Formula

a(1) = 1, and for n > 1, a(n) = -Sum_{d|n, dA033879(n/d) * a(d).
From Antti Karttunen, Nov 14 2024: (Start)
Following convolution formulas have been conjectured for this sequence by Sequence Machine, with each one giving the first 10000 terms correctly:
a(n) = Sum_{d|n} A046692(d)*A067824(n/d).
a(n) = Sum_{d|n} A055615(d)*A074206(n/d).
a(n) = Sum_{d|n} A023900(d)*A174725(n/d).
a(n) = Sum_{d|n} A008683(d)*A323912(n/d).
a(n) = Sum_{d|n} A191161(d)*A327960(n/d).
a(n) = Sum_{d|n} A328722(d)*A330575(n/d).
a(n) = Sum_{d|n} A345182(d)*A349341(n/d).
a(n) = Sum_{d|n} A346246(d)*A349387(n/d).
a(n) = Sum_{d|n} A002033(d-1)*A055615(n/d).
(End)

A336423 Number of strict chains of divisors from n to 1 using terms of A130091 (numbers with distinct prime multiplicities).

Original entry on oeis.org

1, 1, 1, 2, 1, 0, 1, 4, 2, 0, 1, 5, 1, 0, 0, 8, 1, 5, 1, 5, 0, 0, 1, 14, 2, 0, 4, 5, 1, 0, 1, 16, 0, 0, 0, 0, 1, 0, 0, 14, 1, 0, 1, 5, 5, 0, 1, 36, 2, 5, 0, 5, 1, 14, 0, 14, 0, 0, 1, 0, 1, 0, 5, 32, 0, 0, 1, 5, 0, 0, 1, 35, 1, 0, 5, 5, 0, 0, 1, 36, 8, 0, 1, 0
Offset: 1

Views

Author

Gus Wiseman, Jul 27 2020

Keywords

Comments

A number's prime signature (row n of A124010) is the sequence of positive exponents in its prime factorization, so a number has distinct prime multiplicities iff all the exponents in its prime signature are distinct.

Examples

			The a(n) chains for n = 4, 8, 12, 16, 24, 32:
  4/1    8/1      12/1      16/1        24/1         32/1
  4/2/1  8/2/1    12/2/1    16/2/1      24/2/1       32/2/1
         8/4/1    12/3/1    16/4/1      24/3/1       32/4/1
         8/4/2/1  12/4/1    16/8/1      24/4/1       32/8/1
                  12/4/2/1  16/4/2/1    24/8/1       32/16/1
                            16/8/2/1    24/12/1      32/4/2/1
                            16/8/4/1    24/4/2/1     32/8/2/1
                            16/8/4/2/1  24/8/2/1     32/8/4/1
                                        24/8/4/1     32/16/2/1
                                        24/12/2/1    32/16/4/1
                                        24/12/3/1    32/16/8/1
                                        24/12/4/1    32/8/4/2/1
                                        24/8/4/2/1   32/16/4/2/1
                                        24/12/4/2/1  32/16/8/2/1
                                                     32/16/8/4/1
                                                     32/16/8/4/2/1
		

Crossrefs

A336569 is the maximal case.
A336571 does not require n itself to have distinct prime multiplicities.
A000005 counts divisors.
A007425 counts divisors of divisors.
A074206 counts strict chains of divisors from n to 1.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A253249 counts nonempty strict chains of divisors.
A327498 gives the maximum divisor with distinct prime multiplicities.
A336422 counts divisible pairs of divisors, both in A130091.
A336424 counts factorizations using A130091.
A336500 counts divisors of n in A130091 with quotient also in A130091.
A337256 counts strict chains of divisors.

Programs

  • Mathematica
    strchns[n_]:=If[n==1,1,If[!UnsameQ@@Last/@FactorInteger[n],0,Sum[strchns[d],{d,Select[Most[Divisors[n]],UnsameQ@@Last/@FactorInteger[#]&]}]]];
    Table[strchns[n],{n,100}]

A337105 Number of strict chains of divisors from n! to 1.

Original entry on oeis.org

1, 1, 1, 3, 20, 132, 1888, 20128, 584000, 17102016, 553895936, 11616690176, 743337949184, 19467186157568, 999551845713920, 66437400489711616, 10253161206302064640, 388089999627661557760, 53727789519052432998400, 2325767421950553303285760, 365546030278816140131041280
Offset: 0

Views

Author

Gus Wiseman, Aug 17 2020

Keywords

Examples

			The a(4) = 20 chains:
  24/1  24/2/1   24/4/2/1   24/8/4/2/1
        24/3/1   24/6/2/1   24/12/4/2/1
        24/4/1   24/6/3/1   24/12/6/2/1
        24/6/1   24/8/2/1   24/12/6/3/1
        24/8/1   24/8/4/1
        24/12/1  24/12/2/1
                 24/12/3/1
                 24/12/4/1
                 24/12/6/1
		

Crossrefs

A325617 is the maximal case.
A336941 is the version for superprimorials.
A337104 counts the case with distinct prime multiplicities.
A337071 is the case not necessarily ending with 1.
A000005 counts divisors.
A000142 lists factorial numbers.
A001055 counts factorizations.
A027423 counts divisors of factorial numbers.
A067824 counts chains of divisors starting with n.
A074206 counts chains of divisors from n to 1.
A076716 counts factorizations of factorial numbers.
A253249 counts chains of divisors.
A336423 counts chains using A130091, with maximal case A336569.
A336942 counts chains using A130091 from A006939(n) to 1.

Programs

  • Maple
    b:= proc(n) option remember; 1 +
          add(b(d), d=numtheory[divisors](n) minus {n})
        end:
    a:= n-> ceil(b(n!)/2):
    seq(a(n), n=0..14);  # Alois P. Heinz, Aug 23 2020
  • Mathematica
    chnsc[n_]:=Prepend[Join@@Table[Prepend[#,n]&/@chnsc[d],{d,DeleteCases[Divisors[n],1|n]}],{n}];
    Table[Length[chnsc[n!]],{n,0,5}]

Formula

a(n) = A337071(n)/2 for n > 1.
a(n) = A074206(n!).

Extensions

a(19)-a(20) from Alois P. Heinz, Aug 22 2020

A347050 Number of factorizations of n that are a twin (x*x) or have an alternating permutation.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 2, 2, 2, 1, 4, 1, 2, 2, 4, 1, 4, 1, 4, 2, 2, 1, 6, 2, 2, 2, 4, 1, 5, 1, 5, 2, 2, 2, 9, 1, 2, 2, 6, 1, 5, 1, 4, 4, 2, 1, 10, 2, 4, 2, 4, 1, 6, 2, 6, 2, 2, 1, 11, 1, 2, 4, 7, 2, 5, 1, 4, 2, 5, 1, 15, 1, 2, 4, 4, 2, 5, 1, 10, 4, 2, 1, 11, 2
Offset: 1

Views

Author

Gus Wiseman, Oct 15 2021

Keywords

Comments

First differs from A348383 at a(216) = 27, A348383(216) = 28.
A factorization of n is a weakly increasing sequence of positive integers > 1 with product n.
These permutations are ordered factorizations of n with no adjacent triples (..., x, y, z, ...) where x <= y <= z or x >= y >= z.
The version without twins for n > 0 is a(n) + 1 if n is a perfect square; otherwise a(n).

Examples

			The factorizations for n = 4, 12, 24, 30, 36, 48, 60, 64, 72:
  4    12     24     30     36       48       60       64       72
  2*2  2*6    3*8    5*6    4*9      6*8      2*30     8*8      8*9
       3*4    4*6    2*15   6*6      2*24     3*20     2*32     2*36
       2*2*3  2*12   3*10   2*18     3*16     4*15     4*16     3*24
              2*2*6  2*3*5  3*12     4*12     5*12     2*4*8    4*18
              2*3*4         2*2*9    2*3*8    6*10     2*2*16   6*12
                            2*3*6    2*4*6    2*5*6    2*2*4*4  2*4*9
                            3*3*4    3*4*4    3*4*5             2*6*6
                            2*2*3*3  2*2*12   2*2*15            3*3*8
                                     2*2*3*4  2*3*10            3*4*6
                                              2*2*3*5           2*2*18
                                                                2*3*12
                                                                2*2*3*6
                                                                2*3*3*4
                                                                2*2*2*3*3
The a(270) = 19 factorizations:
  (2*3*5*9)   (5*6*9)   (3*90)   (270)
  (3*3*5*6)   (2*3*45)  (5*54)
  (2*3*3*15)  (2*5*27)  (6*45)
              (2*9*15)  (9*30)
              (3*3*30)  (10*27)
              (3*5*18)  (15*18)
              (3*6*15)  (2*135)
              (3*9*10)
Note that (2*3*3*3*5) is separable but has no alternating permutations.
		

Crossrefs

Partitions not of this type are counted by A344654, ranked by A344653.
Partitions of this type are counted by A344740, ranked by A344742.
The complement is counted by A347706, without twins A348380.
The case without twins is A348379.
Dominates A348383, the separable case.
A001055 counts factorizations, strict A045778, ordered A074206.
A001250 counts alternating permutations.
A008480 counts permutations of prime indices, strict A335489.
A025047 counts alternating or wiggly compositions, ranked by A345167.
A056239 adds up prime indices, row sums of A112798.
A325534 counts separable partitions, ranked by A335433.
A325535 counts inseparable partitions, ranked by A335448.
A335452 counts anti-run permutations of prime indices, complement A336107.
A339846 counts even-length factorizations.
A339890 counts odd-length factorizations.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],Function[f,Select[Permutations[f],!MatchQ[#,{_,x_,y_,z_,_}/;x<=y<=z||x>=y>=z]&]!={}]]],{n,100}]

Formula

For n > 1, a(n) = A335434(n) + A010052(n).

A129374 G.f. satisfies: A(x) = 1/(1-x) * A(x^2)*A(x^3)*A(x^4)*...*A(x^n)*...

Original entry on oeis.org

1, 1, 2, 3, 6, 8, 15, 20, 35, 48, 76, 103, 166, 221, 333, 451, 671, 894, 1303, 1730, 2479, 3288, 4615, 6086, 8502, 11142, 15299, 20034, 27285, 35514, 47937, 62168, 83259, 107650, 142929, 184090, 243207, 312041, 409210, 523709, 683261, 871239, 1130703
Offset: 0

Views

Author

Paul D. Hanna, Apr 12 2007

Keywords

Crossrefs

Programs

  • PARI
    {a(n)=local(A=1+x);for(i=2,n,A=1/(1-x)*prod(n=2,i,subst(A,x,x^n+x*O(x^i)))); polcoeff(A,n)}

Formula

G.f.: A(x) = Product_{n>=1} 1/(1 - x^n)^A074206(n) where A074206(n) equals the number of ordered factorizations of n.
a(n) ~ exp((1 + 1/r) * (-Gamma(1+r) * Zeta(1+r) / Zeta'(r))^(1/(1+r)) * n^(r/(1+r))) * (-Gamma(1+r) * Zeta(1+r) / Zeta'(r))^(1/(10*(1+r))) / ((2*Pi)^(29/50) * sqrt(1+r) * n^((6 + 5*r)/(10*(1+r)))), where r = A107311 = 1.7286472389981836181351... is the root of the equation Zeta(r) = 2, Zeta'(r) = -1/A247667. - Vaclav Kotesovec, Nov 04 2018

A212171 Prime signature of n (nonincreasing version): row n of table lists positive exponents in canonical prime factorization of n, in nonincreasing order.

Original entry on oeis.org

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

Views

Author

Matthew Vandermast, Jun 03 2012

Keywords

Comments

Length of row n equals A001221(n).
The multiset of positive exponents in n's prime factorization completely determines a(n) for a host of OEIS sequences, including several "core" sequences. Of those not cross-referenced here or in A212172, many can be found by searching the database for A025487.
(Note: Differing opinions may exist about whether the prime signature of n should be defined as this multiset itself, or as a symbol or collection of symbols that identify or "signify" this multiset. The definition of this sequence is designed to be compatible with either view, as are the original comments. When n >= 2, the customary ways to signify the multiset of exponents in n's prime factorization are to list the constituent exponents in either nonincreasing or nondecreasing order; this table gives the nonincreasing version.)
Table lists exponents in the order in which they appear in the prime factorization of a member of A025487. This ordering is common in database comments (e.g., A008966).
Each possible multiset of an integer's positive prime factorization exponents corresponds to a unique partition that contains the same elements (cf. A000041). This includes the multiset of 1's positive exponents, { } (the empty multiset), which corresponds to the partition of 0.
Differs from A124010 from a(23) on, corresponding to the factorization of 18 = 2^1*3^2 which is here listed as row 18 = [2, 1], but as [1, 2] (in the order of the prime factors) in A124010 and also in A118914 which lists the prime signatures in nondecreasing order (so that row 12 = 2^2*3^1 is also [1, 2]). - M. F. Hasler, Apr 08 2022

Examples

			First rows of table read:
  1;
  1;
  2;
  1;
  1,1;
  1;
  3;
  2;
  1,1;
  1;
  2,1;
  ...
The multiset of positive exponents in the prime factorization of 6 = 2*3 is {1,1} (1s are often left implicit as exponents). The prime signature of 6 is therefore {1,1}.
12 = 2^2*3 has positive exponents 2 and 1 in its prime factorization, as does 18 = 2*3^2. Rows 12 and 18 of the table both read {2,1}.
		

Crossrefs

Cf. A025487, A001221 (row lengths), A001222 (row sums). A118914 gives the nondecreasing version. A124010 lists exponents in n's prime factorization in natural order, with A124010(1) = 0.
A212172 cross-references over 20 sequences that depend solely on n's prime exponents >= 2, including the "core" sequence A000688. Other sequences determined by the exponents in the prime factorization of n include:
Additive: A001221, A001222, A056169.
A highly incomplete selection of sequences, each definable by the set of prime signatures possessed by its members: A000040, A000290, A000578, A000583, A000961, A001248, A001358, A001597, A001694, A002808, A004709, A005117, A006881, A013929, A030059, A030229, A052486.

Programs

  • Magma
    &cat[Reverse(Sort([pe[2]:pe in Factorisation(n)])):n in[1..76]]; // Jason Kimberley, Jun 13 2012
    
  • PARI
    apply( {A212171_row(n)=vecsort(factor(n)[,2]~,,4)}, [1..40])\\ M. F. Hasler, Apr 19 2022

Formula

Row n of A118914, reversed.
Row n of A124010 for n > 1, with exponents sorted in nonincreasing order. Equivalently, row A046523(n) of A124010 for n > 1.

A349059 Number of weakly alternating ordered factorizations of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 3, 1, 4, 2, 3, 1, 8, 1, 3, 3, 8, 1, 8, 1, 8, 3, 3, 1, 18, 2, 3, 4, 8, 1, 11, 1, 16, 3, 3, 3, 22, 1, 3, 3, 18, 1, 11, 1, 8, 8, 3, 1, 38, 2, 8, 3, 8, 1, 18, 3, 18, 3, 3, 1, 32, 1, 3, 8, 28, 3, 11, 1, 8, 3, 11, 1, 56, 1, 3, 8, 8, 3, 11, 1, 38, 8, 3
Offset: 1

Views

Author

Gus Wiseman, Dec 04 2021

Keywords

Comments

An ordered factorization of n is a finite sequence of positive integers > 1 with product n.
We define a sequence to be weakly alternating if it is alternately weakly increasing and weakly decreasing, starting with either.

Examples

			The ordered factorizations for n = 2, 4, 6, 8, 12, 24, 30:
  (2)  (4)    (6)    (8)      (12)     (24)       (30)
       (2*2)  (2*3)  (2*4)    (2*6)    (3*8)      (5*6)
              (3*2)  (4*2)    (3*4)    (4*6)      (6*5)
                     (2*2*2)  (4*3)    (6*4)      (10*3)
                              (6*2)    (8*3)      (15*2)
                              (2*2*3)  (12*2)     (2*15)
                              (2*3*2)  (2*12)     (3*10)
                              (3*2*2)  (2*2*6)    (2*5*3)
                                       (2*4*3)    (3*2*5)
                                       (2*6*2)    (3*5*2)
                                       (3*2*4)    (5*2*3)
                                       (3*4*2)
                                       (4*2*3)
                                       (6*2*2)
                                       (2*2*2*3)
                                       (2*2*3*2)
                                       (2*3*2*2)
                                       (3*2*2*2)
		

Crossrefs

The strong version for compositions is A025047, also A025048, A025049.
The strong case is A348610, complement A348613.
The version for compositions is A349052, complement A349053.
As compositions these are ranked by the complement of A349057.
A001055 counts factorizations, strict A045778, ordered A074206.
A001250 counts alternating permutations, complement A348615.
A335434 counts separable factorizations, complement A333487.
A345164 counts alternating permutations of prime factors, w/ twins A344606.
A345170 counts partitions with an alternating permutation.
A348379 = factorizations w/ alternating permutation, complement A348380.
A348611 counts anti-run ordered factorizations, complement A348616.
A349060 counts weakly alternating partitions, complement A349061.
A349800 = weakly but not strongly alternating compositions, ranked A349799.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&, Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    whkQ[y_]:=And@@Table[If[EvenQ[m],y[[m]]<=y[[m+1]],y[[m]]>=y[[m+1]]], {m,1,Length[y]-1}];
    Table[Length[Select[Join@@Permutations/@facs[n], whkQ[#]||whkQ[-#]&]],{n,100}]

Formula

a(2^n) = A349052(n).

A337135 a(1) = 1; for n > 1, a(n) = Sum_{d|n, d <= sqrt(n)} a(d).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 2, 2, 2, 1, 3, 1, 2, 2, 4, 1, 3, 1, 4, 2, 2, 1, 5, 2, 2, 2, 4, 1, 4, 1, 4, 2, 2, 2, 7, 1, 2, 2, 5, 1, 5, 1, 4, 3, 2, 1, 7, 2, 3, 2, 4, 1, 5, 2, 5, 2, 2, 1, 8, 1, 2, 3, 6, 2, 5, 1, 4, 2, 4, 1, 9, 1, 2, 3, 4, 2, 5, 1, 7, 4, 2, 1, 8, 2, 2, 2, 6, 1, 8, 2, 4, 2, 2, 2
Offset: 1

Views

Author

Ilya Gutkovskiy, Nov 21 2020

Keywords

Comments

From Gus Wiseman, Mar 05 2021: (Start)
This sequence counts all of the following essentially equivalent things:
1. Chains of distinct inferior divisors from n to 1, where a divisor d|n is inferior if d <= n/d. Inferior divisors are counted by A038548 and listed by A161906.
2. Chains of divisors from n to 1 whose first-quotients (in analogy with first-differences) are term-wise greater than or equal to their decapitation (maximum element removed). For example, the divisor chain q = 60/4/2/1 has first-quotients (15,2,2), which are >= (4,2,1), so q is counted under a(60).
3. Chains of divisors from n to 1 such that x >= y^2 for all adjacent x, y.
4. Factorizations of n where each factor is greater than or equal to the product of all previous factors.
(End)

Examples

			From _Gus Wiseman_, Mar 05 2021: (Start)
The a(n) chains for n = 1, 2, 4, 12, 16, 24, 36, 60:
  1  2/1  4/1    12/1    16/1      24/1      36/1      60/1
          4/2/1  12/2/1  16/2/1    24/2/1    36/2/1    60/2/1
                 12/3/1  16/4/1    24/3/1    36/3/1    60/3/1
                         16/4/2/1  24/4/1    36/4/1    60/4/1
                                   24/4/2/1  36/6/1    60/5/1
                                             36/4/2/1  60/6/1
                                             36/6/2/1  60/4/2/1
                                                       60/6/2/1
The a(n) factorizations for n = 2, 4, 12, 16, 24, 36, 60:
    2  4    12   16     24     36     60
       2*2  2*6  2*8    3*8    4*9    2*30
            3*4  4*4    4*6    6*6    3*20
                 2*2*4  2*12   2*18   4*15
                        2*2*6  3*12   5*12
                               2*2*9  6*10
                               2*3*6  2*2*15
                                      2*3*10
(End)
		

Crossrefs

Cf. A002033, A008578 (positions of 1's), A068108.
The restriction to powers of 2 is A018819.
Not requiring inferiority gives A074206 (ordered factorizations).
The strictly inferior version is A342083.
The strictly superior version is A342084.
The weakly superior version is A342085.
The additive version is A000929, or A342098 forbidding equality.
A000005 counts divisors, with sum A000203.
A001055 counts factorizations.
A003238 counts chains of divisors summing to n-1, with strict case A122651.
A038548 counts inferior (or superior) divisors.
A056924 counts strictly inferior (or strictly superior) divisors.
A067824 counts strict chains of divisors starting with n.
A167865 counts strict chains of divisors > 1 summing to n.
A207375 lists central divisors.
A253249 counts strict chains of divisors.
A334996 counts ordered factorizations by product and length.
A334997 counts chains of divisors of n by length.
A342086 counts strict factorizations of divisors.
- Inferior: A033676, A066839, A072499, A161906.
- Superior: A033677, A070038, A161908.
- Strictly Inferior: A060775, A070039, A333806, A341674.
- Strictly Superior: A048098, A064052, A140271, A238535, A341673.

Programs

  • Maple
    a:= proc(n) option remember; `if`(n=1, 1, add(
          `if`(d<=n/d, a(d), 0), d=numtheory[divisors](n)))
        end:
    seq(a(n), n=1..128);  # Alois P. Heinz, Jun 24 2021
  • Mathematica
    a[1] = 1; a[n_] := a[n] = DivisorSum[n, a[#] &, # <= Sqrt[n] &]; Table[a[n], {n, 95}]
    (* second program *)
    asc[n_]:=Prepend[#,n]&/@Prepend[Join@@Table[asc[d],{d,Select[Divisors[n],#Gus Wiseman, Mar 05 2021 *)

Formula

G.f.: Sum_{k>=1} a(k) * x^(k^2) / (1 - x^k).
a(2^n) = A018819(n). - Gus Wiseman, Mar 08 2021

A342087 Number of chains of divisors starting with n and having no adjacent parts x <= y^2.

Original entry on oeis.org

1, 2, 2, 2, 2, 4, 2, 4, 2, 4, 2, 6, 2, 4, 4, 4, 2, 6, 2, 6, 4, 4, 2, 8, 2, 4, 4, 6, 2, 8, 2, 6, 4, 4, 4, 8, 2, 4, 4, 8, 2, 10, 2, 6, 6, 4, 2, 12, 2, 6, 4, 6, 2, 10, 4, 8, 4, 4, 2, 14, 2, 4, 6, 6, 4, 10, 2, 6, 4, 8, 2, 16, 2, 4, 6, 6, 4, 10, 2, 12, 4, 4, 2, 14
Offset: 1

Views

Author

Gus Wiseman, Mar 05 2021

Keywords

Comments

An alternative wording: Number of chains of divisors starting with n and having all adjacent parts x > y^2.

Examples

			The chains for n = 1, 2, 6, 12, 24, 42, 48:
   1    2      6        12        24        42          48
        2/1    6/1      12/1      24/1      42/1        48/1
               6/2      12/2      24/2      42/2        48/2
               6/2/1    12/3      24/3      42/3        48/3
                        12/2/1    24/4      42/6        48/4
                        12/3/1    24/2/1    42/2/1      48/6
                                  24/3/1    42/3/1      48/2/1
                                  24/4/1    42/6/1      48/3/1
                                            42/6/2      48/4/1
                                            42/6/2/1    48/6/1
                                                        48/6/2
                                                        48/6/2/1
		

Crossrefs

The restriction to powers of 2 is A018819.
Not requiring strict inferiority gives A067824.
The weakly inferior version is twice A337135.
The case ending with 1 is counted by A342083.
The strictly superior version is A342084.
The weakly superior version is A342085.
The additive version is A342098, or A000929 allowing equality.
A000005 counts divisors, with sum A000203.
A001055 counts factorizations.
A003238 counts chains of divisors summing to n-1, with strict case A122651.
A038548 counts inferior (or superior) divisors.
A056924 counts strictly inferior (or strictly superior) divisors.
A067824 counts strict chains of divisors starting with n.
A074206 counts ordered factorizations.
A167865 counts strict chains of divisors > 1 summing to n.
A253249 counts strict chains of divisors.
A334997 counts chains of divisors of n by length.

Programs

  • Mathematica
    cem[n_]:=Prepend[Prepend[#,n]&/@Join@@cem/@Most[Divisors[n]],{n}];
    Table[Length[Select[cem[n],And@@Thread[Divide@@@Partition[#,2,1]>Rest[#]]&]],{n,30}]

Formula

For n > 1, a(n) = 2*A342083(n).
Previous Showing 51-60 of 225 results. Next