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

A027193 Number of partitions of n into an odd number of parts.

Original entry on oeis.org

0, 1, 1, 2, 2, 4, 5, 8, 10, 16, 20, 29, 37, 52, 66, 90, 113, 151, 190, 248, 310, 400, 497, 632, 782, 985, 1212, 1512, 1851, 2291, 2793, 3431, 4163, 5084, 6142, 7456, 8972, 10836, 12989, 15613, 18646, 22316, 26561, 31659, 37556, 44601, 52743, 62416, 73593, 86809, 102064, 120025, 140736
Offset: 0

Views

Author

Keywords

Comments

Number of partitions of n in which greatest part is odd.
Number of partitions of n+1 into an even number of parts, the least being 1. Example: a(5)=4 because we have [5,1], [3,1,1,1], [2,1,1] and [1,1,1,1,1,1].
Also number of partitions of n+1 such that the largest part is even and occurs only once. Example: a(5)=4 because we have [6], [4,2], [4,1,1] and [2,1,1,1,1]. - Emeric Deutsch, Apr 05 2006
Also the number of partitions of n such that the number of odd parts and the number of even parts have opposite parities. Example: a(8)=10 is a count of these partitions: 8, 611, 521, 431, 422, 41111, 332, 32111, 22211, 2111111. - Clark Kimberling, Feb 01 2014, corrected Jan 06 2021
In Chaves 2011 see page 38 equation (3.20). - Michael Somos, Dec 28 2014
Suppose that c(0) = 1, that c(1), c(2), ... are indeterminates, that d(0) = 1, and that d(n) = -c(n) - c(n-1)*d(1) - ... - c(0)*d(n-1). When d(n) is expanded as a polynomial in c(1), c(2),..,c(n), the terms are of the form H*c(i_1)*c(i_2)*...*c(i_k). Let P(n) = [c(i_1), c(i_2), ..., c(i_k)], a partition of n. Then H is negative if P has an odd number of parts, and H is positive if P has an even number of parts. That is, d(n) has A027193(n) negative coefficients, A027187(n) positive coefficients, and A000041 terms. The maximal coefficient in d(n), in absolute value, is A102462(n). - Clark Kimberling, Dec 15 2016

Examples

			G.f. = x + x^2 + 2*x^3 + 2*x^4 + 4*x^5 + 5*x^6 + 8*x^7 + 10*x^8 + 16*x^9 + 20*x^10 + ...
From _Gus Wiseman_, Feb 11 2021: (Start)
The a(1) = 1 through a(8) = 10 partitions into an odd number of parts are the following. The Heinz numbers of these partitions are given by A026424.
  (1)  (2)  (3)    (4)    (5)      (6)      (7)        (8)
            (111)  (211)  (221)    (222)    (322)      (332)
                          (311)    (321)    (331)      (422)
                          (11111)  (411)    (421)      (431)
                                   (21111)  (511)      (521)
                                            (22111)    (611)
                                            (31111)    (22211)
                                            (1111111)  (32111)
                                                       (41111)
                                                       (2111111)
The a(1) = 1 through a(8) = 10 partitions whose greatest part is odd are the following. The Heinz numbers of these partitions are given by A244991.
  (1)  (11)  (3)    (31)    (5)      (33)      (7)        (53)
             (111)  (1111)  (32)     (51)      (52)       (71)
                            (311)    (321)     (322)      (332)
                            (11111)  (3111)    (331)      (521)
                                     (111111)  (511)      (3221)
                                               (3211)     (3311)
                                               (31111)    (5111)
                                               (1111111)  (32111)
                                                          (311111)
                                                          (11111111)
(End)
		

References

  • N. J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; p. 39, Example 7.

Crossrefs

The Heinz numbers of these partitions are A026424 or A244991.
The even-length version is A027187.
The case of odd sum as well as length is A160786, ranked by A340931.
The case of odd maximum as well as length is A340385.
Other cases of odd length:
- A024429 counts set partitions of odd length.
- A067659 counts strict partitions of odd length.
- A089677 counts ordered set partitions of odd length.
- A166444 counts compositions of odd length.
- A174726 counts ordered factorizations of odd length.
- A332304 counts strict compositions of odd length.
- A339890 counts factorizations of odd length.
A000009 counts partitions into odd parts, ranked by A066208.
A026804 counts partitions whose least part is odd.
A058695 counts partitions of odd numbers, ranked by A300063.
A072233 counts partitions by sum and length.
A101707 counts partitions of odd positive rank.

Programs

  • Maple
    g:=sum(x^(2*k)/product(1-x^j,j=1..2*k-1),k=1..40): gser:=series(g,x=0,50): seq(coeff(gser,x,n),n=1..45); # Emeric Deutsch, Apr 05 2006
  • Mathematica
    nn=40;CoefficientList[Series[ Sum[x^(2j+1)Product[1/(1- x^i),{i,1,2j+1}],{j,0,nn}],{x,0,nn}],x]  (* Geoffrey Critzer, Dec 01 2012 *)
    a[ n_] := If[ n < 0, 0, Length@Select[ IntegerPartitions[ n], OddQ[ Length@#] &]]; (* Michael Somos, Dec 28 2014 *)
    a[ n_] := If[ n < 1, 0, Length@Select[ IntegerPartitions[ n], OddQ[ First@#] &]]; (* Michael Somos, Dec 28 2014 *)
    a[ n_] := If[ n < 0, 0, Length@Select[ IntegerPartitions[ n + 1], #[[-1]] == 1 && EvenQ[ Length@#] &]]; (* Michael Somos, Dec 28 2014 *)
    a[ n_] := If[ n < 1, 0, Length@Select[ IntegerPartitions[ n + 1], EvenQ[ First@#] && (Length[#] < 2 || #[[1]] != #[[2]]) &]]; (* Michael Somos, Dec 28 2014 *)
  • PARI
    {a(n) = if( n<1, 0, polcoeff( sum( k=1, n, if( k%2, x^k / prod( j=1, k, 1 - x^j, 1 + x * O(x^(n-k)) ))), n))}; /* Michael Somos, Jul 24 2012 */
    
  • PARI
    q='q+O('q^66); concat([0], Vec( (1/eta(q)-eta(q)/eta(q^2))/2 ) ) \\ Joerg Arndt, Mar 23 2014

Formula

a(n) = (A000041(n) - (-1)^n*A000700(n)) / 2.
For g.f. see under A027187.
G.f.: Sum(k>=1, x^(2*k-1)/Product(j=1..2*k-1, 1-x^j ) ). - Emeric Deutsch, Apr 05 2006
G.f.: - Sum(k>=1, (-x)^(k^2)) / Product(k>=1, 1-x^k ). - Joerg Arndt, Feb 02 2014
G.f.: Sum(k>=1, x^(k*(2*k-1)) / Product(j=1..2*k, 1-x^j)). - Michael Somos, Dec 28 2014
a(2*n) = A000701(2*n), a(2*n-1) = A046682(2*n-1); a(n) = A000041(n)-A027187(n). - Reinhard Zumkeller, Apr 22 2006

A067659 Number of partitions of n into distinct parts such that number of parts is odd.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 2, 2, 3, 4, 5, 6, 8, 9, 11, 14, 16, 19, 23, 27, 32, 38, 44, 52, 61, 71, 82, 96, 111, 128, 148, 170, 195, 224, 256, 293, 334, 380, 432, 491, 557, 630, 713, 805, 908, 1024, 1152, 1295, 1455, 1632, 1829, 2048, 2291, 2560, 2859, 3189, 3554, 3958, 4404
Offset: 0

Views

Author

Naohiro Nomoto, Feb 23 2002

Keywords

Comments

Ramanujan theta functions: phi(q) := Sum_{k=-oo..oo} q^(k^2) (A000122), chi(q) := Prod_{k>=0} (1+q^(2k+1)) (A000700).

Examples

			From _Gus Wiseman_, Jan 09 2021: (Start)
The a(5) = 1 through a(15) = 14 partitions (A-F = 10..15):
  5   6     7     8     9     A     B     C     D     E     F
      321   421   431   432   532   542   543   643   653   654
                  521   531   541   632   642   652   743   753
                        621   631   641   651   742   752   762
                              721   731   732   751   761   843
                                    821   741   832   842   852
                                          831   841   851   861
                                          921   931   932   942
                                                A21   941   951
                                                      A31   A32
                                                      B21   A41
                                                            B31
                                                            C21
                                                            54321
(End)
		

Crossrefs

Dominates A000009.
Numbers with these strict partitions as binary indices are A000069.
The non-strict version is A027193.
The Heinz numbers of these partitions are A030059.
The even version is A067661.
The version for rank is A117193, with non-strict version A101707.
The ordered version is A332304, with non-strict version A166444.
Other cases of odd length:
- A024429 counts set partitions of odd length.
- A089677 counts ordered set partitions of odd length.
- A174726 counts ordered factorizations of odd length.
- A339890 counts factorizations of odd length.
A008289 counts strict partitions by sum and length.
A026804 counts partitions whose least part is odd, with strict case A026832.

Programs

  • Maple
    b:= proc(n, i, t) option remember; `if`(n>i*(i+1)/2, 0,
          `if`(n=0, t, add(b(n-i*j, i-1, abs(t-j)), j=0..min(n/i, 1))))
        end:
    a:= n-> b(n$2, 0):
    seq(a(n), n=0..80);  # Alois P. Heinz, Apr 01 2014
  • Mathematica
    b[n_, i_, t_] := b[n, i, t] = If[n > i*(i + 1)/2, 0, If[n == 0, t, Sum[b[n - i*j, i - 1, Abs[t - j]], {j, 0, Min[n/i, 1]}]]]; a[n_] := b[n, n, 0]; Table[a[n], {n, 0, 80}] (* Jean-François Alcover, Jan 16 2015, after Alois P. Heinz *)
    CoefficientList[Normal[Series[(QPochhammer[-x, x]-QPochhammer[x])/2, {x, 0, 100}]], x] (* Andrey Zabolotskiy, Apr 12 2017 *)
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&OddQ[Length[#]]&]],{n,0,30}] (* Gus Wiseman, Jan 09 2021 *)
  • PARI
    {a(n)=local(A); if(n<0, 0, A=x*O(x^n); polcoeff( (eta(x^2+A)/eta(x+A) - eta(x+A))/2, n))} /* Michael Somos, Feb 14 2006 */
    
  • PARI
    N=66;  q='q+O('q^N);  S=1+2*sqrtint(N);
    gf=sum(n=1,S, (n%2!=0) * q^(n*(n+1)/2) / prod(k=1,n, 1-q^k ) );
    concat( [0], Vec(gf) )  /* Joerg Arndt, Oct 20 2012 */
    
  • PARI
    N=66;  q='q+O('q^N);  S=1+sqrtint(N);
    gf=sum(n=1, S, q^(2*n^2-n) / prod(k=1, 2*n-1, 1-q^k ) );
    concat( [0], Vec(gf) )  \\ Joerg Arndt, Apr 01 2014

Formula

For g.f. see under A067661.
a(n) = (A000009(n)-A010815(n))/2. - Vladeta Jovovic, Feb 24 2002
Expansion of (1-phi(-q))/(2*chi(-q)) in powers of q where phi(),chi() are Ramanujan theta functions. - Michael Somos, Feb 14 2006
G.f.: sum(n>=1, q^(2*n^2-n) / prod(k=1..2*n-1, 1-q^k ) ). [Joerg Arndt, Apr 01 2014]
a(n) = A067661(n) - A010815(n). - Andrey Zabolotskiy, Apr 12 2017
A000009(n) = a(n) + A067661(n). - Gus Wiseman, Jan 09 2021

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

A028983 Numbers whose sum of divisors is even.

Original entry on oeis.org

3, 5, 6, 7, 10, 11, 12, 13, 14, 15, 17, 19, 20, 21, 22, 23, 24, 26, 27, 28, 29, 30, 31, 33, 34, 35, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 80, 82
Offset: 1

Views

Author

Keywords

Comments

The even terms of this sequence are the even terms appearing in A178910. [Edited by M. F. Hasler, Oct 02 2014]
A071324(a(n)) is even. - Reinhard Zumkeller, Jul 03 2008
Sigma(a(n)) = A000203(a(n)) = A152678(n). - Jaroslav Krizek, Oct 06 2009
A083207 is a subsequence. - Reinhard Zumkeller, Jul 19 2010
Numbers k such that the number of odd divisors of k (A001227) is even. - Omar E. Pol, Apr 04 2016
Numbers k such that the sum of odd divisors of k (A000593) is even. - Omar E. Pol, Jul 05 2016
Numbers with a squarefree part greater than 2. - Peter Munn, Apr 26 2020
Equivalently, numbers whose odd part is nonsquare. Compare with the numbers whose square part is even (i.e., nonodd): these are the positive multiples of 4, A008586\{0}, and A225546 provides a self-inverse bijection between the two sets. - Peter Munn, Jul 19 2020
Also numbers whose reversed prime indices have alternating product > 1, where we define the alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^(i-1)). Also Heinz numbers of the partitions counted by A347448. - Gus Wiseman, Oct 29 2021
Numbers whose number of middle divisors is not odd (cf. A067742). - Omar E. Pol, Aug 02 2022

Crossrefs

The complement is A028982 = A000290 U A001105.
Subsequences: A083207, A091067, A145204\{0}, A225838, A225858.
Cf. A334748 (a permutation).
Related to A008586 via A225546.
Ranks the partitions counted by A347448, complement A119620.

Programs

  • Mathematica
    Select[Range[82],EvenQ[DivisorSigma[1,#]]&] (* Jayanta Basu, Jun 05 2013 *)
  • PARI
    is(n)=!issquare(n)&&!issquare(n/2) \\ Charles R Greathouse IV, Jan 11 2013
    
  • Python
    from math import isqrt
    def A028983(n):
        def f(x): return n-1+isqrt(x)+isqrt(x>>1)
        kmin, kmax = 1,2
        while f(kmax) >= kmax:
            kmax <<= 1
        while True:
            kmid = kmax+kmin>>1
            if f(kmid) < kmid:
                kmax = kmid
            else:
                kmin = kmid
            if kmax-kmin <= 1:
                break
        return kmax # Chai Wah Wu, Aug 22 2024

Formula

a(n) ~ n. - Charles R Greathouse IV, Jan 11 2013
a(n) = n + (1 + sqrt(2)/2)*sqrt(n) + O(1). - Charles R Greathouse IV, Sep 01 2015
A007913(a(n)) > 2. - Peter Munn, May 05 2020

A340653 Number of balanced factorizations of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jan 15 2021

Keywords

Comments

A factorization into factors > 1 is balanced if it is empty or its length is equal to its maximum Omega (A001222).

Examples

			The balanced factorizations for n = 120, 144, 192, 288, 432, 768:
  3*5*8    2*8*9    3*8*8      4*8*9      6*8*9      8*8*12
  2*2*30   3*6*8    4*6*8      6*6*8      2*8*27     2*2*8*24
  2*3*20   2*4*18   2*8*12     2*8*18     3*8*18     2*3*8*16
  2*5*12   2*6*12   4*4*12     3*8*12     4*4*27     2*4*4*24
           3*4*12   2*2*2*24   4*4*18     4*6*18     2*4*6*16
                    2*2*3*16   4*6*12     4*9*12     3*4*4*16
                               2*12*12    6*6*12     2*2*12*16
                               2*2*2*36   2*12*18    2*2*2*2*48
                               2*2*3*24   3*12*12    2*2*2*3*32
                               2*3*3*16   2*2*2*54
                                          2*2*3*36
                                          2*3*3*24
                                          3*3*3*16
		

Crossrefs

Positions of zeros are A001358.
Positions of nonzero terms are A100959.
The co-balanced version is A340596.
Taking maximum factor instead of maximum Omega gives A340599.
The cross-balanced version is A340654.
The twice-balanced version is A340655.
A001055 counts factorizations.
A045778 counts strict factorizations.
A316439 counts factorizations by product and length.
A320655 counts factorizations into semiprimes.
Other balance-related sequences:
- A010054 counts balanced strict partitions.
- A047993 counts balanced partitions.
- A098124 counts balanced compositions.
- A106529 lists Heinz numbers of balanced partitions.
- A340597 have an alt-balanced factorization.
- A340598 counts balanced set partitions.
- A340600 counts unlabeled balanced multiset partitions.
- A340656 have no twice-balanced factorizations.
- A340657 have a twice-balanced factorization.

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],#=={}||Length[#]==Max[PrimeOmega/@#]&]],{n,100}]
  • PARI
    A340653(n, m=n, mbo=0, e=0) = if(1==n, mbo==e, sumdiv(n, d, if((d>1)&&(d<=m), A340653(n/d, d, max(mbo,bigomega(d)), 1+e)))); \\ Antti Karttunen, Oct 22 2023

Extensions

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

A347437 Number of factorizations of n with integer alternating product.

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, 2, 2, 1, 2, 2, 1, 1, 1, 4, 1, 1, 1, 6, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 5, 2, 2, 1, 2, 1, 3, 1, 2, 1, 1, 1, 2, 1, 1, 2, 8, 1, 1, 1, 2, 1, 1, 1, 6, 1, 1, 2, 2, 1, 1, 1, 5, 4, 1, 1, 2, 1, 1, 1, 2, 1, 3, 1, 2, 1, 1, 1, 6, 1, 2, 2, 6, 1, 1, 1, 2, 1, 1, 1, 7
Offset: 1

Views

Author

Gus Wiseman, Sep 06 2021

Keywords

Comments

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

Examples

			The factorizations for n = 4, 16, 36, 48, 54, 64, 108:
  (4)   (16)      (36)      (48)        (54)    (64)          (108)
  (2*2) (4*4)     (6*6)     (2*4*6)     (2*3*9) (8*8)         (2*6*9)
        (2*2*4)   (2*2*9)   (3*4*4)     (3*3*6) (2*4*8)       (3*6*6)
        (2*2*2*2) (2*3*6)   (2*2*12)            (4*4*4)       (2*2*27)
                  (3*3*4)   (2*2*2*2*3)         (2*2*16)      (2*3*18)
                  (2*2*3*3)                     (2*2*4*4)     (3*3*12)
                                                (2*2*2*2*4)   (2*2*3*3*3)
                                                (2*2*2*2*2*2)
		

Crossrefs

Positions of 1's are A005117, complement A013929.
Allowing any alternating product <= 1 gives A339846.
Allowing any alternating product > 1 gives A339890.
The restriction to powers of 2 is A344607.
The even-length case is A347438, also the case of alternating product 1.
The reciprocal version is A347439.
Allowing any alternating product < 1 gives A347440.
The odd-length case is A347441.
The reverse version is A347442.
The additive version is A347446, ranked by A347457.
Allowing any alternating product >= 1 gives A347456.
The restriction to perfect squares is A347458, reciprocal A347459.
The ordered version is A347463.
A001055 counts factorizations.
A046099 counts factorizations with no alternating permutations.
A071321 gives the alternating sum of prime factors of n (reverse: A071322).
A273013 counts ordered factorizations of n^2 with alternating product 1.
A347460 counts possible alternating products of factorizations.

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],IntegerQ@*altprod]],{n,100}]
  • PARI
    A347437(n, m=n, ap=1, e=0) = if(1==n, if(e%2, 1==denominator(ap), 1==numerator(ap)), sumdiv(n, d, if((d>1)&&(d<=m), A347437(n/d, d, ap * d^((-1)^e), 1-e)))); \\ Antti Karttunen, Oct 22 2023

Formula

a(2^n) = A344607(n).
a(n^2) = A347458(n).

Extensions

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

A347438 Number of unordered factorizations of n with alternating product 1.

Original entry on oeis.org

1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Sep 06 2021

Keywords

Comments

Also the number of unordered factorizations of n with alternating sum 0.
Also the number of unordered factorizations of n with all even multiplicities.
This is the even-length case of A347437, the odd-length case being A347441.
An unordered factorization of n is a weakly increasing sequence of positive integers > 1 with product n.
We define the alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^(i-1)).

Examples

			The a(n) factorizations for n = 16, 64, 144, 256, 576:
  4*4      8*8          12*12        16*16            24*24
  2*2*2*2  2*2*4*4      2*2*6*6      2*2*8*8          3*3*8*8
           2*2*2*2*2*2  3*3*4*4      4*4*4*4          4*4*6*6
                        2*2*2*2*3*3  2*2*2*2*4*4      2*2*12*12
                                     2*2*2*2*2*2*2*2  2*2*2*2*6*6
                                                      2*2*3*3*4*4
                                                      2*2*2*2*2*2*3*3
		

Crossrefs

Positions of zeros are A000037.
Positions of nonzero terms are A000290.
The restriction to perfect squares is A001055 (ordered: A273013).
The restriction to powers of 2 is A035363.
The additive version is A119620, ranked by A028982.
Positions of non-1's are A213367 \ {1}.
Positions of 1's are A280076 = {1} \/ A001248.
Sorted first positions are 1, 2, and all terms of A330972 squared.
Allowing any alternating product <= 1 gives A339846.
Allowing any alternating product > 1 gives A339890.
Allowing any integer alternating product gives A347437.
Allowing any integer reciprocal alternating product gives A347439.
Allowing any alternating product < 1 gives A347440.
Allowing any alternating product >= 1 gives A347456.
A046099 counts factorizations with no alternating permutations.
A071321 gives the alternating sum of prime factors (reverse: A071322).
A316524 gives the alternating sum of prime indices (reverse: A344616).
A344606 counts alternating permutations of prime factors.
A347441 counts odd-length factorizations with integer alternating product.
A347460 counts possible alternating products of factorizations.

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}]
  • PARI
    A347438(n, m=n, k=0, t=1) = if(1==n, (1==t), my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A347438(n/d, d, 1-k, t*(d^((-1)^k))))); (s)); \\ Antti Karttunen, Oct 30 2021

Formula

a(2^n) = A035363(n).
a(n^2) = A001055(n).

Extensions

Name and comments clarified (with unordered) by Jacob Sprittulla, Oct 05 2021

A340596 Number of co-balanced factorizations of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jan 15 2021

Keywords

Comments

We define a factorization of n into factors > 1 to be co-balanced if it has exactly A001221(n) factors.

Examples

			The a(n) co-balanced factorizations for n = 12, 24, 36, 72, 120, 144, 180:
  2*6    3*8     4*9     8*9     3*5*8     2*72     4*5*9
  3*4    4*6     6*6     2*36    4*5*6     3*48     5*6*6
         2*12    2*18    3*24    2*2*30    4*36     2*2*45
                 3*12    4*18    2*3*20    6*24     2*3*30
                         6*12    2*4*15    8*18     2*5*18
                                 2*5*12    9*16     2*6*15
                                 2*6*10    12*12    2*9*10
                                 3*4*10             3*3*20
                                                    3*4*15
                                                    3*5*12
                                                    3*6*10
		

Crossrefs

Positions of terms > 1 are A126706.
Positions of 1's are A303554.
The version for unlabeled multiset partitions is A319616.
The alt-balanced version is A340599.
The balanced version is A340653.
The cross-balanced version is A340654.
The twice-balanced version is A340655.
A001055 counts factorizations.
A045778 counts strict factorizations.
A316439 counts factorizations by product and length.
Other balance-related sequences:
- A010054 counts balanced strict partitions.
- A047993 counts balanced partitions.
- A098124 counts balanced compositions.
- A106529 lists Heinz numbers of balanced partitions.
- A340597 lists numbers with an alt-balanced factorization.
- A340598 counts balanced set partitions.
- A340600 counts unlabeled balanced multiset partitions.

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],Length[#]==PrimeNu[n]&]],{n,100}]
  • PARI
    A340596(n, m=n, om=omega(n)) = if(1==n,(0==om), sumdiv(n, d, if((d>1)&&(d<=m), A340596(n/d, d, om-1)))); \\ Antti Karttunen, Jun 10 2024

Extensions

Data section extended up to a(120) by Antti Karttunen, Jun 10 2024

A340604 Heinz numbers of integer partitions of odd positive rank.

Original entry on oeis.org

3, 7, 10, 13, 15, 19, 22, 25, 28, 29, 33, 34, 37, 42, 43, 46, 51, 52, 53, 55, 61, 62, 63, 69, 70, 71, 76, 77, 78, 79, 82, 85, 88, 89, 93, 94, 98, 101, 105, 107, 113, 114, 115, 116, 117, 118, 119, 121, 123, 130, 131, 132, 134, 136, 139, 141, 146, 147, 148, 151
Offset: 1

Views

Author

Gus Wiseman, Jan 21 2021

Keywords

Comments

The Dyson rank of a nonempty partition is its maximum part minus its number of parts. The rank of an empty partition is 0.
The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.

Examples

			The sequence of partitions with their Heinz numbers begins:
      3: (2)         46: (9,1)       82: (13,1)
      7: (4)         51: (7,2)       85: (7,3)
     10: (3,1)       52: (6,1,1)     88: (5,1,1,1)
     13: (6)         53: (16)        89: (24)
     15: (3,2)       55: (5,3)       93: (11,2)
     19: (8)         61: (18)        94: (15,1)
     22: (5,1)       62: (11,1)      98: (4,4,1)
     25: (3,3)       63: (4,2,2)    101: (26)
     28: (4,1,1)     69: (9,2)      105: (4,3,2)
     29: (10)        70: (4,3,1)    107: (28)
     33: (5,2)       71: (20)       113: (30)
     34: (7,1)       76: (8,1,1)    114: (8,2,1)
     37: (12)        77: (5,4)      115: (9,3)
     42: (4,2,1)     78: (6,2,1)    116: (10,1,1)
     43: (14)        79: (22)       117: (6,2,2)
		

Crossrefs

Note: Heinz numbers are given in parentheses below.
These partitions are counted by A101707.
Allowing negative ranks gives A340692, counted by A340603.
The even version is A340605, counted by A101708.
The not necessarily odd case is A340787, counted by A064173.
A001222 gives number of prime indices.
A061395 gives maximum prime index.
- Rank -
A047993 counts partitions of rank 0 (A106529).
A064173 counts partitions of negative rank (A340788).
A064174 counts partitions of nonnegative rank (A324562).
A064174 (also) counts partitions of nonpositive rank (A324521).
A101198 counts partitions of rank 1 (A325233).
A257541 gives the rank of the partition with Heinz number n.
A340653 counts balanced factorizations.
- Odd -
A000009 counts partitions into odd parts (A066208).
A027193 counts partitions of odd length (A026424).
A027193 (also) counts partitions of odd maximum (A244991).
A058695 counts partitions of odd numbers (A300063).
A067659 counts strict partitions of odd length (A030059).
A160786 counts odd-length partitions of odd numbers (A300272).
A339890 counts factorizations of odd length.
A340101 counts factorizations into odd factors.
A340102 counts odd-length factorizations into odd factors.
A340385 counts partitions of odd length and maximum (A340386).

Programs

  • Mathematica
    rk[n_]:=PrimePi[FactorInteger[n][[-1,1]]]-PrimeOmega[n];
    Select[Range[100],OddQ[rk[#]]&&rk[#]>0&]

Formula

A061395(a(n)) - A001222(a(n)) is odd and positive.

A347439 Number of factorizations of n with integer reciprocal alternating product.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Sep 07 2021

Keywords

Comments

All of these factorizations have an even number of factors, so their reverse-alternating product is also an integer.
A factorization of n is a weakly increasing sequence of positive integers > 1 with product n.
We define the reciprocal alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^i).
The value of a(n) does not depend solely on the prime signature of n. See the example comparing a(144) and a(400). - Antti Karttunen, Jul 28 2024

Examples

			The a(n) factorizations for
n    = 16,       36,       64,           72,       128,          144:
a(n) = 3,        4,        6,            5,        7,            11
--------------------------------------------------------------------------------
       2*8       6*6       8*8           2*36      2*64          2*72
       4*4       2*18      2*32          3*24      4*32          3*48
       2*2*2*2   3*12      4*16          6*12      8*16          4*36
                 2*2*3*3   2*2*2*8       2*2*3*6   2*2*4*8       6*24
                           2*2*4*4       2*3*3*4   2*4*4*4       12*12
                           2*2*2*2*2*2             2*2*2*16      2*2*6*6
                                                   2*2*2*2*2*4   2*3*3*8
                                                                 3*3*4*4
                                                                 2*2*2*18
                                                                 2*2*3*12
                                                                 2*2*2*2*3*3
From _Antti Karttunen_, Jul 28 2024 (Start)
For n=400, there are 12 such factorizations:
  2*200
  4*100
  5*80
  10*40
  20*20
  2*2*2*50
  2*2*5*20
  2*2*10*10
  2*4*5*10
  2*5*5*8
  4*4*5*5
  2*2*2*2*5*5.
Note that 400 = 2^4 * 5^2 has the same prime signature as 144 = 2^4 * 3^2. 400 = 2*4*5*10 is the factorization for which there is no analogous factorization of 144, as 2*3*4*6 doesn't satisfy the condition of having an integer reciprocal alternating product.
(End)
		

Crossrefs

Positions of 0's are A005117 \ {1}.
Positions of non-0's are 1 and A013929.
The restriction to powers of 2 is A027187, reverse A035363.
Positions of 1's are 1 and A082293.
The additive version is A119620, ranked by A347451 and A028982.
Allowing any alternating product <= 1 gives A339846.
Allowing any alternating product > 1 gives A339890.
The non-reciprocal version is A347437.
The reverse version is A347438.
Allowing any alternating product < 1 gives A347440.
The non-reciprocal reverse version is A347442.
Allowing any alternating product >= 1 gives A347456.
The restriction to perfect squares is A347459, non-reciprocal A347458.
A038548 counts possible reverse-alternating products of factorizations.
A046099 counts factorizations with no alternating permutations.
A071321 gives the alternating sum of prime factors (reverse: A071322).
A316524 gives the alternating sum of prime indices (reverse: A344616).
A273013 counts ordered factorizations of n^2 with alternating product 1.
A347441 counts odd-length factorizations with integer alternating product.
A347460 counts possible alternating products of factorizations.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    recaltprod[q_]:=Product[q[[i]]^(-1)^i,{i,Length[q]}];
    Table[Length[Select[facs[n],IntegerQ[recaltprod[#]]&]],{n,100}]
  • PARI
    A347439(n, m=n, ap=1, e=0) = if(1==n, !(e%2) && 1==denominator(ap), sumdiv(n, d, if(d>1 && d<=m, A347439(n/d, d, ap * d^((-1)^e), 1-e)))); \\ Antti Karttunen, Jul 28 2024
    
  • PARI
    A347439(n, m=0, ap=1, e=1) = if(1==n, 1==denominator(ap), sumdiv(n, d, if(d>1 && d>=m, A347439(n/d, d, ap * d^((-1)^e), 1-e)))); \\ Antti Karttunen, Jul 28 2024

Formula

a(2^n) = A027187(n).
a(n^2) = A347459(n).

Extensions

Data section extended up to a(108) by Antti Karttunen, Jul 28 2024
Showing 1-10 of 77 results. Next