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 21-30 of 225 results. Next

A339890 Number of odd-length factorizations of n into factors > 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Dec 28 2020

Keywords

Examples

			The a(n) factorizations for n = 24, 48, 60, 72, 96, 120:
  24      48          60       72          96          120
  2*2*6   2*3*8       2*5*6    2*4*9       2*6*8       3*5*8
  2*3*4   2*4*6       3*4*5    2*6*6       3*4*8       4*5*6
          3*4*4       2*2*15   3*3*8       4*4*6       2*2*30
          2*2*12      2*3*10   3*4*6       2*2*24      2*3*20
          2*2*2*2*3            2*2*18      2*3*16      2*4*15
                               2*3*12      2*4*12      2*5*12
                               2*2*2*3*3   2*2*2*2*6   2*6*10
                                           2*2*2*3*4   3*4*10
                                                       2*2*2*3*5
		

Crossrefs

The case of set partitions (or n squarefree) is A024429.
The case of partitions (or prime powers) is A027193.
The ordered version is A174726 (even: A174725).
The remaining (even-length) factorizations are counted by A339846.
A000009 counts partitions into odd parts, ranked by A066208.
A001055 counts factorizations, with strict case A045778.
A027193 counts partitions of odd length, ranked by A026424.
A058695 counts partitions of odd numbers, ranked by A300063.
A160786 counts odd-length partitions of odd numbers, ranked by A300272.
A316439 counts factorizations by product and length.
A340101 counts factorizations into odd factors.
A340102 counts odd-length factorizations into odd factors.

Programs

  • Maple
    g:= proc(n, k, t) option remember; `if`(n>k, 0, t)+
          `if`(isprime(n), 0, add(`if`(d>k, 0, g(n/d, d, 1-t)),
              d=numtheory[divisors](n) minus {1, n}))
        end:
    a:= n-> `if`(n<2, 0, g(n$2, 1)):
    seq(a(n), n=1..100);  # Alois P. Heinz, Dec 30 2020
  • 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],OddQ@Length[#]&]],{n,100}]

Formula

a(n) + A339846(n) = A001055(n).

A339846 Number of even-length factorizations of n into factors > 1.

Original entry on oeis.org

1, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 2, 0, 1, 1, 3, 0, 2, 0, 2, 1, 1, 0, 4, 1, 1, 1, 2, 0, 3, 0, 3, 1, 1, 1, 5, 0, 1, 1, 4, 0, 3, 0, 2, 2, 1, 0, 6, 1, 2, 1, 2, 0, 4, 1, 4, 1, 1, 0, 6, 0, 1, 2, 6, 1, 3, 0, 2, 1, 3, 0, 8, 0, 1, 2, 2, 1, 3, 0, 6, 3, 1, 0, 6, 1, 1, 1, 4, 0, 6, 1, 2, 1, 1, 1, 10, 0, 2, 2, 5, 0, 3, 0, 4, 3
Offset: 1

Views

Author

Gus Wiseman, Dec 28 2020

Keywords

Examples

			The a(n) factorizations for n = 12, 16, 24, 36, 48, 72, 96, 120:
  2*6  2*8      3*8      4*9      6*8      8*9      2*48         2*60
  3*4  4*4      4*6      6*6      2*24     2*36     3*32         3*40
       2*2*2*2  2*12     2*18     3*16     3*24     4*24         4*30
                2*2*2*3  3*12     4*12     4*18     6*16         5*24
                         2*2*3*3  2*2*2*6  6*12     8*12         6*20
                                  2*2*3*4  2*2*2*9  2*2*3*8      8*15
                                           2*2*3*6  2*2*4*6      10*12
                                           2*3*3*4  2*3*4*4      2*2*5*6
                                                    2*2*2*12     2*3*4*5
                                                    2*2*2*2*2*3  2*2*2*15
                                                                 2*2*3*10
		

Crossrefs

The case of set partitions (or n squarefree) is A024430.
The case of partitions (or prime powers) is A027187.
The ordered version is A174725, odd: A174726.
The odd-length factorizations are counted by A339890.
A001055 counts factorizations, with strict case A045778.
A001358 lists semiprimes, with squarefree case A006881.
A027187 counts partitions of even length, ranked by A028260.
A058696 counts partitions of even numbers, ranked by A300061.
A316439 counts factorizations by product and length.
A340102 counts odd-length factorizations into odd factors.

Programs

  • Maple
    g:= proc(n, k, t) option remember; `if`(n>k, 0, t)+
          `if`(isprime(n), 0, add(`if`(d>k, 0, g(n/d, d, 1-t)),
              d=numtheory[divisors](n) minus {1, n}))
        end:
    a:= n-> `if`(n=1, 1, g(n$2, 0)):
    seq(a(n), n=1..100);  # Alois P. Heinz, Dec 30 2020
  • 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],EvenQ@Length[#]&]],{n,100}]
  • PARI
    A339846(n, m=n, e=1) = if(1==n, e, sumdiv(n, d, if((d>1)&&(d<=m), A339846(n/d, d, 1-e)))); \\ Antti Karttunen, Oct 22 2023

Formula

a(n) + A339890(n) = A001055(n).

Extensions

Data section extended up to a(105) by Antti Karttunen, Oct 22 2023

A253249 Number of nonempty chains in the divides relation on the divisors of n.

Original entry on oeis.org

1, 3, 3, 7, 3, 11, 3, 15, 7, 11, 3, 31, 3, 11, 11, 31, 3, 31, 3, 31, 11, 11, 3, 79, 7, 11, 15, 31, 3, 51, 3, 63, 11, 11, 11, 103, 3, 11, 11, 79, 3, 51, 3, 31, 31, 11, 3, 191, 7, 31, 11, 31, 3, 79, 11, 79, 11, 11, 3, 175, 3, 11, 31, 127, 11, 51, 3, 31, 11, 51
Offset: 1

Views

Author

Geoffrey Critzer, Jun 04 2015

Keywords

Comments

For prime p, a(p)=3.
a(2^k) = 2^(k+1)-1.
For integers of the form n = p_1*p_2*...*p_k we have a(n) = A007047(k).
The value of a(n) depends only on the exponents in the prime factorization of n.

Examples

			a(10) = 11 because we have: {1}, {2}, {5}, {10}, {1|2}, {1|5}, {1|10}, {2|10}, {5|10}, {1|2|10}, {1|5|10}.
		

Crossrefs

Cf. A002033, A007047, A074206, A107311, A378219 (Dirichlet inverse).

Programs

  • Maple
    with(numtheory):
    b:= proc(n) option remember: 1+ `if`(n=1, 0,
           add(b(d), d=divisors(n) minus {n}))
        end:
    a:= n-> add(b(d), d=divisors(n)):
    seq(a(n), n=1..100);  # Alois P. Heinz, Jun 04 2015
  • Mathematica
    Table[Total[Table[Length[Select[Subsets[Divisors[n], {k}],Apply[And, Map[Apply[Divisible, #] &,Partition[Reverse[#], 2, 1]]] &]], {k, 1,PrimeOmega[n] + 1}]], {n, 1, 100}]

Formula

Dirichlet g.f.: zeta(s)^2*A(s) where A(s) is the Dirichlet g.f. for A074206. - Geoffrey Critzer, May 23 2018
Sum_{k=1..n} a(k) ~ -4*n^r / (r*Zeta'(r)), where r = A107311 = 1.728647238998183618135103... is the root of the equation zeta(r) = 2. - Vaclav Kotesovec, Jan 31 2019
a(n) = 4*A002033(n-1) - 1 for n > 1. - Geoffrey Critzer, Aug 19 2020

A070824 Number of divisors of n which are > 1 and < n (nontrivial divisors).

Original entry on oeis.org

0, 0, 0, 1, 0, 2, 0, 2, 1, 2, 0, 4, 0, 2, 2, 3, 0, 4, 0, 4, 2, 2, 0, 6, 1, 2, 2, 4, 0, 6, 0, 4, 2, 2, 2, 7, 0, 2, 2, 6, 0, 6, 0, 4, 4, 2, 0, 8, 1, 4, 2, 4, 0, 6, 2, 6, 2, 2, 0, 10, 0, 2, 4, 5, 2, 6, 0, 4, 2, 6, 0, 10, 0, 2, 4, 4, 2, 6, 0, 8, 3, 2, 0, 10, 2, 2
Offset: 1

Views

Author

Wolfdieter Lang, May 08 2002

Keywords

Comments

These are sometimes called the proper divisors, but see A032741 for the usual meaning of that term.
a(n) = number of ordered factorizations of n into two factors, n = 2, 3, ... If n has the prime factorization n=Product p^e(j), j=1..r, the number of compositions of the vector (e(1), ..., e(r)) equals the number of ordered factorizations of n. Andrews (1998, page 59) gives a formula for the number of m-compositions of (e(1), ..., e(r)) which equals the number f(n,m) of ordered m-factorizations of n. But with m=2 the formula reduces to f(n,2) = d(n)-2 = a(n). - Augustine O. Munagi, Mar 31 2005
a(n) = 0 if and only if n is 1 or prime. - Jon Perry, Nov 08 2008
For n > 2: number of zeros in n-th row of triangle A051778. - Reinhard Zumkeller, Dec 03 2014
a(n) = number of partitions of n in which largest and least parts occur exactly once and their difference is 2. Example: a(12) = 4 because we have [7,5], [5,4,3], [4,3,3,2], and [3,2,2,2,2,1]. In general, if d is a nontrivial divisor of n, then [d+1,{d}^(n/d-2),d-1] is a partition of n of the prescribed type. - Emeric Deutsch, Nov 03 2015
Absolute values of the inverse Möbius transform of (-1)^prime(n), n >= 2. - Wesley Ivan Hurt, Jun 22 2024

Examples

			a(12) = 4 with the nontrivial divisors 2,3,4,6.
a(24) = 6 = card({{2,12},{3,8},{4,6},{6,4},{8,3},{12,2}}). - _Peter Luschny_, Nov 14 2011
		

References

  • George E. Andrews, The Theory of Partitions, Addison-Wesley, Reading 1976; reprinted, Cambridge University Press, Cambridge, 1984, 1998.

Crossrefs

First column in the matrix power A175992^2.
Row sums of A175992 starting from the second column.
Column k=2 of A251683.

Programs

  • Haskell
    a070824 n = if n == 1 then 0 else length $ tail $ a027751_row n -- Reinhard Zumkeller, Dec 03 2014
    
  • Maple
    0, seq(numtheory[tau](n)-2,n=2..100); # Augustine O. Munagi, Mar 31 2005
  • Mathematica
    Join[{0},Rest[DivisorSigma[0,Range[90]]-2]] (* Harvey P. Dale, Jun 23 2012 *)
    a[ n_] := SeriesCoefficient[ Sum[x^(2 k)/(1 - x^k), {k, 2, n/2}], {x, 0, n}]; (* Michael Somos, Jun 24 2019 *)
  • PARI
    {a(n) = if( n<1, 0, my(v = vector(n, i, i>1)); dirmul(v, v)[n])}; /* Michael Somos, Jun 24 2019 */
    
  • PARI
    apply( A070824(n)=numdiv(n+(n<2))-2, [1..90]) \\ M. F. Hasler, Oct 11 2019
    
  • Python
    from sympy import divisor_count
    def A070824(n): return 0 if n == 1 else divisor_count(n)-2 # Chai Wah Wu, Jun 03 2022

Formula

a(n) = A000005(n)-2, n>=2 (with the number-of-divisors function d(n) = A000005(n)).
a(n) = d(n)-2, for n>=2, where d(n) is the number-of-divisors function. E.g., a(12) = 4 because 12 has 4 ordered factorizations into two factors: 2*6, 6*2, 3*4, 4*3. - Augustine O. Munagi, Mar 31 2005
G.f.: Sum_{k>=2} x^(2k)/(1-x^k). - Jon Perry, Nov 08 2008
Dirichlet generating function: (zeta(s)-1)^2. - Mats Granvik May 25 2013
Sum_{k=1..n} a(k) ~ n*log(n) + (2*gamma - 3)*n, where gamma is Euler's constant (A001620). - Amiram Eldar, Nov 27 2022
a(n) = abs( Sum_{d|n} (-1)^prime(d) ), n >= 2 with a(1) = 0. - Wesley Ivan Hurt, Jun 22 2024
a(n) = Sum_{k=2..n-1} floor(n/k) - floor((n-1)/k), see Chhimpa and Yadav. - Stefano Spezia, Oct 13 2024

Extensions

a(1)=0 added by Peter Luschny, Nov 14 2011
Several minor edits by M. F. Hasler, Oct 14 2019

A153881 1 followed by -1, -1, -1, ... .

Original entry on oeis.org

1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1
Offset: 1

Views

Author

Mats Granvik, Jan 03 2009

Keywords

Comments

Dirichlet inverse of A074206.

Crossrefs

If prefixed by initial 0, we get A134824.
Cf. A074206 (Dirichlet inverse).

Programs

Formula

G.f: x*(1-2*x)/(1-x). - Mats Granvik, Mar 09 2009, rewritten R. J. Mathar, Mar 31 2010
a(n) = (-1)^A000040(n). - Juri-Stepan Gerasimov, Sep 10 2009
G.f.: x / (1 + x / (1 - 2*x)). - Michael Somos, Apr 02 2012
From Wesley Ivan Hurt, Jun 20 2014: (Start)
a(1) = 1; a(n) = -1, n > 1.
a(n) = 1 - 2*sign(n-1) = 1 - 2*A057427(n-1).
a(n) = (-1)^sign(1-n) = (-1)^A057427(1-n).
a(n) = 2*floor(1/n)-1 = 2*A063524(n)-1. (End)
Dirichlet g.f.: 2 - zeta(s). - Álvar Ibeas, Dec 30 2018
a(n) = Sum_{d|n} A033879(d)*A055615(n/d) = Sum_{d|n} A344587(d)*A346234(n/d). - Antti Karttunen, Nov 22 2024

Extensions

Edited by Charles R Greathouse IV, Mar 18 2010
More terms from Antti Karttunen, Nov 22 2024

A069916 Number of log-concave compositions (ordered partitions) of n.

Original entry on oeis.org

1, 1, 2, 4, 6, 9, 14, 20, 26, 36, 47, 60, 80, 102, 127, 159, 194, 236, 291, 355, 425, 514, 611, 718, 856, 1009, 1182, 1381, 1605, 1861, 2156, 2496, 2873, 3299, 3778, 4301, 4902, 5574, 6325, 7176, 8116, 9152, 10317, 11610, 13028, 14611, 16354, 18259, 20365
Offset: 0

Views

Author

Pontus von Brömssen, Apr 24 2002

Keywords

Comments

These are compositions with weakly decreasing first quotients, where the first quotients of a sequence are defined as if the sequence were an increasing divisor chain, so for example the first quotients of (6,3,1) are (1/2,1/3). - Gus Wiseman, Mar 16 2021

Examples

			Out of the 8 compositions of 4, only 2+1+1 and 1+1+2 are not log-concave, so a(4)=6.
From _Gus Wiseman_, Mar 15 2021: (Start)
The a(1) = 1 through a(6) = 14 compositions:
  (1)  (2)   (3)    (4)     (5)      (6)
       (11)  (12)   (13)    (14)     (15)
             (21)   (22)    (23)     (24)
             (111)  (31)    (32)     (33)
                    (121)   (41)     (42)
                    (1111)  (122)    (51)
                            (131)    (123)
                            (221)    (132)
                            (11111)  (141)
                                     (222)
                                     (231)
                                     (321)
                                     (1221)
                                     (111111)
(End)
		

Crossrefs

The version for differences instead of quotients is A070211.
A000005 counts constant compositions.
A000009 counts strictly increasing (or strictly decreasing) compositions.
A000041 counts weakly increasing (or weakly decreasing) compositions.
A001055 counts factorizations.
A002843 counts compositions with adjacent parts x <= 2y.
A003238 counts chains of divisors summing to n-1, with strict case A122651.
A003242 counts anti-run compositions.
A074206 counts ordered factorizations.
A167865 counts strict chains of divisors summing to n.

Programs

  • Mathematica
    (* This program is not suitable for computing a large number of terms *)
    compos[n_] := Permutations /@ IntegerPartitions[n] // Flatten[#, 1]&;
    logConcaveQ[p_] := And @@ Table[p[[i]]^2 >= p[[i-1]]*p[[i+1]], {i, 2, Length[p]-1}]; a[n_] := Count[compos[n], p_?logConcaveQ]; Table[an = a[n]; Print["a(", n, ") = ", an]; a[n], {n, 0, 25}] (* Jean-François Alcover, Feb 29 2016 *)
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],GreaterEqual@@Divide@@@Reverse/@Partition[#,2,1]&]],{n,0,15}] (* Gus Wiseman, Mar 15 2021 *)
  • Sage
    def A069916(n) : return sum(all(p[i]^2 >= p[i-1] * p[i+1] for i in range(1, len(p)-1)) for p in Compositions(n)) # Eric M. Schmidt, Sep 29 2013

A251683 Irregular triangular array: T(n,k) is the number of ordered factorizations of n with exactly k factors, n >= 1, 1 <= k <= A086436(n).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 4, 3, 1, 1, 2, 1, 2, 1, 3, 3, 1, 1, 1, 4, 3, 1, 1, 4, 3, 1, 2, 1, 2, 1, 1, 6, 9, 4, 1, 1, 1, 2, 1, 2, 1, 1, 4, 3, 1, 1, 6, 6, 1, 1, 4, 6, 4, 1, 1, 2, 1, 2, 1, 2, 1, 7, 12, 6, 1, 1, 2, 1, 2, 1, 6, 9, 4
Offset: 1

Views

Author

Geoffrey Critzer, Dec 06 2014

Keywords

Comments

Row sums = A074206.
Row lengths give A086436.
T(n,2) = A070824(n).
T(n,3) = A200221(n).
Sum_{k>=1} k*T(n,k) = A254577.
For all n > 1, Sum_{k=1..A086436(n)} (-1)^k*T(n,k) = A008683(n). - Geoffrey Critzer, May 25 2018
From Gus Wiseman, Aug 21 2020: (Start)
Also the number of strict length k + 1 chains of divisors from n to 1. For example, row n = 24 counts the following 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
(End)

Examples

			Triangle T(n,k) begins:
  1;
  1;
  1;
  1, 1;
  1;
  1, 2;
  1;
  1, 2, 1;
  1, 1;
  1, 2;
  1;
  1, 4, 3;
  1;
  1, 2;
  1, 2;
  ...
There are 8 ordered factorizations of the integer 12: 12, 6*2, 4*3, 3*4, 2*6, 3*2*2, 2*3*2, 2*2*3.  So T(12,1)=1, T(12,2)=4, and T(12,3)=3.
		

Crossrefs

A008480 gives rows ends.
A086436 gives row lengths.
A124433 is the same except for signs and zeros.
A334996 is the same except for zeros.
A337107 is the restriction to factorial numbers (but with zeros).
A000005 counts divisors.
A001055 counts factorizations.
A001222 counts prime factors with multiplicity.
A074206 counts strict chains of divisors from n to 1.
A067824 counts strict chains of divisors starting with n.
A122651 counts strict chains of divisors summing to n.
A167865 counts strict chains of divisors > 1 summing to n.
A253249 counts strict nonempty chains of divisors of n.
A337071 counts strict chains of divisors starting with n!.
A337256 counts strict chains of divisors of n.

Programs

  • Maple
    with(numtheory):
    b:= proc(n) option remember; expand(x*(1+
          add(b(n/d), d=divisors(n) minus {1, n})))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=1..degree(p)))(b(n)):
    seq(T(n), n=1..100);  # Alois P. Heinz, Dec 07 2014
  • Mathematica
    f[1] = {{}};
    f[n_] := f[n] =
      Level[Table[
        Map[Prepend[#, d] &, f[n/d]], {d, Rest[Divisors[n]]}], {2}];
    Prepend[Map[Select[#, # > 0 &] &,
      Drop[Transpose[
        Table[Map[Count[#, k] &,
          Map[Length, Table[f[n], {n, 1, 40}], {2}]], {k, 1, 10}]],
       1]],{1}] // Grid
    (* Second program: *)
    b[n_] := b[n] = x(1+Sum[b[n/d], {d, Divisors[n]~Complement~{1, n}}]);
    T[n_] := CoefficientList[b[n]/x, x];
    Array[T, 100] // Flatten (* Jean-François Alcover, Nov 17 2020, after Alois P. Heinz *)

Formula

Dirichlet g.f.: 1/(1 - y*(zeta(x)-1)).

A340101 Number of factorizations of 2n + 1 into odd factors > 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Dec 28 2020

Keywords

Examples

			The factorizations for 2n + 1 = 27, 45, 135, 225, 315, 405, 1155:
  27      45      135       225       315       405         1155
  3*9     5*9     3*45      3*75      5*63      5*81        15*77
  3*3*3   3*15    5*27      5*45      7*45      9*45        21*55
          3*3*5   9*15      9*25      9*35      15*27       33*35
                  3*5*9     15*15     15*21     3*135       3*385
                  3*3*15    5*5*9     3*105     5*9*9       5*231
                  3*3*3*5   3*3*25    5*7*9     3*3*45      7*165
                            3*5*15    3*3*35    3*5*27      11*105
                            3*3*5*5   3*5*21    3*9*15      3*5*77
                                      3*7*15    3*3*5*9     3*7*55
                                      3*3*5*7   3*3*3*15    5*7*33
                                                3*3*3*3*5   3*11*35
                                                            5*11*21
                                                            7*11*15
                                                            3*5*7*11
		

Crossrefs

The version for partitions is A160786, ranked by A300272.
The even version is A340785.
The odd-length case is A340102.
A000009 counts partitions into odd parts, ranked by A066208.
A001055 counts factorizations, with strict case A045778.
A027193 counts partitions of odd length, ranked by A026424.
A058695 counts partitions of odd numbers, ranked by A300063.
A316439 counts factorizations by product and length.
Odd bisection of A001055, and also of A349907.

Programs

  • Maple
    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=numtheory[divisors](n) minus {1, n}))
        end:
    a:= n-> g(2*n+1$2):
    seq(a(n), n=0..100);  # Alois P. Heinz, Dec 30 2020
  • 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],OddQ[Times@@#]&]],{n,1,100,2}]
  • PARI
    A001055(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A001055(n/d, d))); (s)); \\ After code in A001055
    A340101(n) = A001055(n+n+1); \\ Antti Karttunen, Dec 13 2021

Formula

a(n) = A001055(2n+1).
a(n) = A349907(2n+1). - Antti Karttunen, Dec 13 2021

Extensions

Data section extended up to 105 terms by Antti Karttunen, Dec 13 2021

A347456 Number of factorizations of n with alternating product >= 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Oct 09 2021

Keywords

Comments

We define the alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^(i-1)).
A factorization of n is a weakly increasing sequence of positive integers > 1 with product n.
Also the number of factorizations of n with alternating sum >= 0.

Examples

			The a(n) factorizations for n = 4, 16, 24, 36, 60, 64, 96:
  4     16        24      36        60       64            96
  2*2   4*4       2*2*6   6*6       2*5*6    8*8           2*6*8
        2*2*4     2*3*4   2*2*9     3*4*5    2*4*8         3*4*8
        2*2*2*2           2*3*6     2*2*15   4*4*4         4*4*6
                          3*3*4     2*3*10   2*2*16        2*2*24
                          2*2*3*3            2*2*4*4       2*3*16
                                             2*2*2*2*4     2*4*12
                                             2*2*2*2*2*2   2*2*2*2*6
                                                           2*2*2*3*4
		

Crossrefs

The case of partitions is A000041, reverse A344607.
The reverse version is A001055, strict A347705.
Positions of 3's appear to be A065036.
Positions of 1's are 1 and A167171.
The opposite version (<= instead of >=) is A339846.
The strict version (> instead of >=) is A339890, also the odd-length case.
Allowing any integer alternating product gives A347437.
The case of alternating product 1 is A347438, also the even-length case.
Allowing any integer reciprocal alternating product gives A347439.
The complement (< instead of >=) is A347440.
Allowing any integer reverse-alternating product gives A347442.
A038548 counts factorizations with a wiggly permutation.
A045778 counts strict factorizations.
A074206 counts ordered factorizations.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A119620 counts partitions with alternating product 1.
A347447 counts strict factorizations with alternating product > 1.

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[facs[n],altprod[#]>=1&]],{n,100}]

Formula

a(n) = A347438(n) + A347440(n).

A348379 Number of factorizations of n with an alternating permutation.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 2, 3, 1, 4, 1, 4, 2, 2, 1, 6, 1, 2, 2, 4, 1, 5, 1, 5, 2, 2, 2, 8, 1, 2, 2, 6, 1, 5, 1, 4, 4, 2, 1, 10, 1, 4, 2, 4, 1, 6, 2, 6, 2, 2, 1, 11, 1, 2, 4, 6, 2, 5, 1, 4, 2, 5, 1, 15, 1, 2, 4, 4, 2, 5, 1, 10, 3, 2, 1, 11, 2
Offset: 1

Views

Author

Gus Wiseman, Oct 28 2021

Keywords

Comments

First differs from A335434 at a(216) = 27, A335434(216) = 28. Also differs from A335434 at a(270) = 19, A335434(270) = 20.
A factorization of n is a weakly increasing sequence of positive integers > 1 with product n.
All of the counted factorizations are separable (A335434).
A sequence is alternating if it is alternately strictly increasing and strictly decreasing, starting with either. For example, the partition (3,2,2,2,1) has no alternating permutations, even though it does have the anti-run permutations (2,3,2,1,2) and (2,1,2,3,2). Alternating permutations of multisets are a generalization of alternating or up-down permutations of {1..n}.

Examples

			The a(270) = 19 factorizations:
  (2*3*3*15)  (2*3*45)  (2*135)  (270)
  (2*3*5*9)   (2*5*27)  (3*90)
  (3*3*5*6)   (2*9*15)  (5*54)
              (3*3*30)  (6*45)
              (3*5*18)  (9*30)
              (3*6*15)  (10*27)
              (3*9*10)  (15*18)
              (5*6*9)
		

Crossrefs

Partitions not of this type are counted by A345165, ranked by A345171.
Partitions of this type are counted by A345170, ranked by A345172.
Twins and partitions of this type are counted by A344740, ranked by A344742.
The case with twins is A347050.
The complement is counted by A348380, without twins A347706.
The ordered version is A348610.
A001055 counts factorizations, strict A045778, ordered A074206.
A001250 counts alternating permutations.
A025047 counts alternating or wiggly compositions, ranked by A345167.
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]]}]];
    wigQ[y_]:=Or[Length[y]==0,Length[Split[y]]==Length[y]&&Length[Split[Sign[Differences[y]]]]==Length[y]-1];
    Table[Length[Select[facs[n],Select[Permutations[#],wigQ]!={}&]],{n,100}]

Formula

a(2^n) = A345170(n).
Previous Showing 21-30 of 225 results. Next