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 11-20 of 73 results. Next

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

A059897 Symmetric square array read by antidiagonals: A(n,k) is the product of all factors that occur in one, but not both, of the Fermi-Dirac factorizations of n and k.

Original entry on oeis.org

1, 2, 2, 3, 1, 3, 4, 6, 6, 4, 5, 8, 1, 8, 5, 6, 10, 12, 12, 10, 6, 7, 3, 15, 1, 15, 3, 7, 8, 14, 2, 20, 20, 2, 14, 8, 9, 4, 21, 24, 1, 24, 21, 4, 9, 10, 18, 24, 28, 30, 30, 28, 24, 18, 10, 11, 5, 27, 2, 35, 1, 35, 2, 27, 5, 11, 12, 22, 30, 36, 40, 42, 42, 40, 36, 30, 22, 12, 13, 24, 33
Offset: 1

Views

Author

Marc LeBrun, Feb 06 2001

Keywords

Comments

Old name: Square array read by antidiagonals: T(i,j) = product prime(k)^(Ei(k) XOR Ej(k)) where Ei and Ej are the vectors of exponents in the prime factorizations of i and j; XOR is the bitwise operation on binary representation of the exponents.
Analogous to multiplication, with XOR replacing +.
From Peter Munn, Apr 01 2019: (Start)
(1) Defines an abelian group whose underlying set is the positive integers. (2) Every element is self-inverse. (3) For all n and k, A(n,k) is a divisor of n*k. (4) The terms of A050376, sometimes called Fermi-Dirac primes, form a minimal set of generators. In ordered form, it is the lexicographically earliest such set.
The unique factorization of positive integers into products of distinct terms of the group's lexicographically earliest minimal set of generators seems to follow from (1) (2) and (3).
From (1) and (2), every row and every column of the table is a self-inverse permutation of the positive integers. Rows/columns numbered by nonmembers of A050376 are compositions of earlier rows/columns.
It is a subgroup of the equivalent group over the nonzero integers, which has -1 as an additional generator.
As generated by A050376, the subgroup of even length words is A000379. The complementary set of odd length words is A000028.
The subgroup generated by A000040 (the primes) is A005117 (the squarefree numbers).
(End)
Considered as a binary operation, the result is (the squarefree part of the product of its operands) times the square of (the operation's result when applied to the square roots of the square parts of its operands). - Peter Munn, Mar 21 2022

Examples

			A(864,1944) = A(2^5*3^3,2^3*3^5) = 2^(5 XOR 3) * 3^(3 XOR 5) = 2^6 * 3^6 = 46656.
The top left 12 X 12 corner of the array:
   1,  2,  3,  4,  5,  6,  7,  8,   9,  10,  11,  12
   2,  1,  6,  8, 10,  3, 14,  4,  18,   5,  22,  24
   3,  6,  1, 12, 15,  2, 21, 24,  27,  30,  33,   4
   4,  8, 12,  1, 20, 24, 28,  2,  36,  40,  44,   3
   5, 10, 15, 20,  1, 30, 35, 40,  45,   2,  55,  60
   6,  3,  2, 24, 30,  1, 42, 12,  54,  15,  66,   8
   7, 14, 21, 28, 35, 42,  1, 56,  63,  70,  77,  84
   8,  4, 24,  2, 40, 12, 56,  1,  72,  20,  88,   6
   9, 18, 27, 36, 45, 54, 63, 72,   1,  90,  99, 108
  10,  5, 30, 40,  2, 15, 70, 20,  90,   1, 110, 120
  11, 22, 33, 44, 55, 66, 77, 88,  99, 110,   1, 132
  12, 24,  4,  3, 60,  8, 84,  6, 108, 120, 132,   1
From _Peter Munn_, Apr 04 2019: (Start)
The subgroup generated by {6,8,10}, the first three integers > 1 not in A050376, has the following table:
    1     6     8    10    12    15    20   120
    6     1    12    15     8    10   120    20
    8    12     1    20     6   120    10    15
   10    15    20     1   120     6     8    12
   12     8     6   120     1    20    15    10
   15    10   120     6    20     1    12     8
   20   120    10     8    15    12     1     6
  120    20    15    12    10     8     6     1
(End)
		

Crossrefs

Cf. A284567 (A000142 or A003418-analog for this operation).
Rows/columns: A073675 (2), A120229 (3), A120230 (4), A307151 (5), A307150 (6), A307266 (8), A307267 (24).
Particularly significant subgroups or cosets: A000028, A000379, A003159, A005117, A030229, A252895. See also the lists in A329050, A352273.
Sequences that relate this sequence to multiplication: A000188, A007913, A059895.

Programs

  • Mathematica
    a[i_, i_] = 1;
    a[i_, j_] := Module[{f1 = FactorInteger[i], f2 = FactorInteger[j], e1, e2}, e1[] = 0; Scan[(e1[#[[1]]] = #[[2]])&, f1]; e2[] = 0; Scan[(e2[#[[1]]] = #[[2]])&, f2]; Times @@ (#^BitXor[e1[#], e2[#]]& /@ Union[f1[[All, 1]], f2[[All, 1]]])];
    Table[a[i - j + 1, j], {i, 1, 15}, {j, 1, i}] // Flatten (* Jean-François Alcover, Jun 19 2018 *)
  • PARI
    T(n,k) = {if (n==1, return (k)); if (k==1, return (n)); my(fn=factor(n), fk=factor(k)); vp = setunion(fn[,1]~, fk[,1]~); prod(i=1, #vp, vp[i]^(bitxor(valuation(n, vp[i]), valuation(k, vp[i]))));} \\ Michel Marcus, Apr 03 2019
    
  • PARI
    T(i, j) = {if(gcd(i, j) == 1, return(i * j)); if(i == j, return(1)); my(f = vecsort(concat(factor(i)~, factor(j)~)), t = 1, res = 1); while(t + 1 <= #f, if(f[1, t] == f[1, t+1], res *= f[1, t] ^ bitxor(f[2, t] , f[2, t+1]); t+=2; , res*= f[1, t]^f[2, t]; t++; ) ); if(t == #f, res *= f[1, #f] ^ f[2, #f]); res } \\ David A. Corneth, Apr 03 2019
    
  • PARI
    A059897(n,k) = if(n==k, 1, core(n*k) * A059897(core(n,1)[2],core(k,1)[2])^2) \\ Peter Munn, Mar 21 2022
  • Scheme
    (define (A059897 n) (A059897bi (A002260 n) (A004736 n)))
    (define (A059897bi a b) (let loop ((a a) (b b) (m 1)) (cond ((= 1 a) (* m b)) ((= 1 b) (* m a)) ((equal? (A020639 a) (A020639 b)) (loop (A028234 a) (A028234 b) (* m (expt (A020639 a) (A003987bi (A067029 a) (A067029 b)))))) ((< (A020639 a) (A020639 b)) (loop (/ a (A028233 a)) b (* m (A028233 a)))) (else (loop a (/ b (A028233 b)) (* m (A028233 b)))))))
    ;; Antti Karttunen, Apr 11 2017
    

Formula

For all x, y >= 1, A(x,y) * A059895(x,y)^2 = x*y. - Antti Karttunen, Apr 11 2017
From Peter Munn, Apr 01 2019: (Start)
A(n,1) = A(1,n) = n
A(n, A(m,k)) = A(A(n,m), k)
A(n,n) = 1
A(n,k) = A(k,n)
if i_1 <> i_2 then A(A050376(i_1), A050376(i_2)) = A050376(i_1) * A050376(i_2)
if A(n,k_1) = n * k_1 and A(n,k_2) = n * k_2 then A(n, A(k_1,k_2)) = n * A(k_1,k_2)
(End)
T(k, m) = k*m for coprime k and m. - David A. Corneth, Apr 03 2019
if A(n*m,m) = n, A(n*m,k) = A(n,k) * A(m,k) / k. - Peter Munn, Apr 04 2019
A(n,k) = A007913(n*k) * A(A000188(n), A000188(k))^2. - Peter Munn, Mar 21 2022

Extensions

New name from Peter Munn, Mar 21 2022

A030059 Numbers that are the product of an odd number of distinct primes.

Original entry on oeis.org

2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 30, 31, 37, 41, 42, 43, 47, 53, 59, 61, 66, 67, 70, 71, 73, 78, 79, 83, 89, 97, 101, 102, 103, 105, 107, 109, 110, 113, 114, 127, 130, 131, 137, 138, 139, 149, 151, 154, 157, 163, 165, 167, 170, 173, 174, 179, 181, 182, 186, 190, 191, 193
Offset: 1

Views

Author

Keywords

Comments

From Enrique Pérez Herrero, Jul 06 2012: (Start)
This sequence and A030229 partition the squarefree numbers: A005117.
Also solutions to the equation mu(n) = -1.
Sum_{n>=1} 1/a(n)^s = (zeta(s)^2 - zeta(2*s))/(2*zeta(s)*zeta(2*s)). (End) [See A088245 and the Hardy reference. - Wolfdieter Lang, Oct 18 2016]
The lexicographically least sequence of integers > 1 such that for each entry, the number of proper divisors occurring in the sequence is equal to 0 modulo 3. - Masahiko Shin, Feb 12 2018
The asymptotic density of this sequence is 3/Pi^2 (A104141). - Amiram Eldar, May 22 2020
Solutions to the equation Sum_{d|n} mu(d)*sigma(d) = -n, where sigma(n) is the sum of divisors function (A000203). - Robert D. Rosales, May 20 2024

References

  • B. C. Berndt & R. A. Rankin, Ramanujan: Letters and Commentary, see p. 23; AMS Providence RI 1995.
  • G. H. Hardy, Ramanujan, AMS Chelsea Publishing, 2002, pp. 64 - 65, (misprint on p. 65, line starting with Hence: it should be ... -1/Zeta(s) not ... -Zeta(s)).
  • S. Ramanujan, Collected Papers, Ed. G. H. Hardy et al., Cambridge 1927; Chelsea, NY, 1962, p. xxiv, 21.

Crossrefs

Programs

  • Maple
    a := n -> `if`(numtheory[mobius](n)=-1,n,NULL); seq(a(i),i=1..193); # Peter Luschny, May 04 2009
    # alternative
    A030059 := proc(n)
        option remember;
        local a;
        if n = 1 then
            2;
        else
            for a from procname(n-1)+1 do
                if numtheory[mobius](a) = -1 then
                    return a;
                end if;
            end do:
        end if;
    end proc: # R. J. Mathar, Sep 22 2020
  • Mathematica
    Select[Range[300], MoebiusMu[#] == -1 &] (* Enrique Pérez Herrero, Jul 06 2012 *)
  • PARI
    is(n)=my(f=factor(n)[,2]); #f%2 && vecmax(f)==1 \\ Charles R Greathouse IV, Oct 16 2015
    
  • PARI
    is(n)=moebius(n)==-1 \\ Charles R Greathouse IV, Jan 31 2017
    
  • Python
    from math import isqrt, prod
    from sympy import primerange, integer_nthroot, primepi
    def A030059(n):
        def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b+1,isqrt(x//c)+1),a+1)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b+1,integer_nthroot(x//c,m)[0]+1),a+1) for d in g(x,a2,b2,c*b2,m-1)))
        def f(x): return int(n+x-primepi(x)-sum(sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x,0,1,1,i)) for i in range(3,x.bit_length(),2)))
        kmin, kmax = 0,1
        while f(kmax) > kmax:
            kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax # Chai Wah Wu, Aug 29 2024

Formula

omega(a(n)) = A001221(a(n)) gives A005408. {primes A000040} UNION {sphenic numbers A007304} UNION {numbers that are divisible by exactly 5 different primes A051270} UNION {products of 7 distinct primes (squarefree 7-almost primes) A123321} UNION {products of 9 distinct primes; also n has exactly 9 distinct prime factors and n is squarefree A115343} UNION.... - Jonathan Vos Post, Oct 19 2007
a(n) < n*Pi^2/3 infinitely often; a(n) > n*Pi^2/3 infinitely often. - Charles R Greathouse IV, Sep 07 2017

Extensions

More terms from David W. Wilson

A067661 Number of partitions of n into distinct parts such that number of parts is even.

Original entry on oeis.org

1, 0, 0, 1, 1, 2, 2, 3, 3, 4, 5, 6, 7, 9, 11, 13, 16, 19, 23, 27, 32, 38, 45, 52, 61, 71, 83, 96, 111, 128, 148, 170, 195, 224, 256, 292, 334, 380, 432, 491, 556, 630, 713, 805, 908, 1024, 1152, 1295, 1455, 1632, 1829, 2049, 2291, 2560, 2859, 3189, 3554, 3959, 4404
Offset: 0

Views

Author

Naohiro Nomoto, Feb 23 2002

Keywords

Comments

Ramanujan theta functions: phi(q) (A000122), chi(q) (A000700).

Examples

			G.f. = 1 + x^3 + x^4 + 2*x^5 + 2*x^6 + 3*x^7 + 3*x^8 + 4*x^9 + 5*x^10 + ...
From _Gus Wiseman_, Jan 08 2021: (Start)
The a(3) = 1 through a(14) = 11 partitions (A-D = 10..13):
  21   31   32   42   43   53   54   64     65     75     76     86
            41   51   52   62   63   73     74     84     85     95
                      61   71   72   82     83     93     94     A4
                                81   91     92     A2     A3     B3
                                     4321   A1     B1     B2     C2
                                            5321   5421   C1     D1
                                                   6321   5431   5432
                                                          6421   6431
                                                          7321   6521
                                                                 7421
                                                                 8321
(End)
		

References

  • B. C. Berndt, Ramanujan's Notebooks Part III, Springer-Verlag, see p. 18 Entry 9 Corollary (2).

Crossrefs

Dominates A000009.
Numbers with these strict partitions as binary indices are A001969.
The non-strict case is A027187, ranked by A028260.
The Heinz numbers of these partitions are A030229.
The odd version is A067659, ranked by A030059.
The version for rank is A117192, with positive case A101708.
Other cases of even length:
- A024430 counts set partitions of even length.
- A034008 counts compositions of even length.
- A052841 counts ordered set partitions of even length.
- A174725 counts ordered factorizations of even length.
- A332305 counts strict compositions of even length
- A339846 counts factorizations of even length.
A008289 counts strict partitions by sum and length.
A026805 counts partitions whose least part is even.

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, 1):
    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, 1]; Table[a[n], {n, 0, 80}] (* Jean-François Alcover, Jan 16 2015, after Alois P. Heinz *)
    a[ n_] := SeriesCoefficient[ (QPochhammer[ -x, x] + QPochhammer[ x]) / 2, {x, 0, n}]; (* Michael Somos, May 06 2015 *)
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&EvenQ[Length[#]]&]],{n,0,30}] (* Gus Wiseman, Jan 08 2021 *)
  • PARI
    {a(n) = my(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=0, S, (n%2==0) * q^(n*(n+1)/2) / prod(k=1, n, 1-q^k ) );
    Vec(gf)  \\ Joerg Arndt, Apr 01 2014

Formula

G.f.: A(q) = Sum_{n >= 0} a(n) q^n = 1 + q^3 + q^4 + 2 q^5 + 2 q^6 + 3 q^7 + ... = Sum_{n >= 0} q^(n(2n+1))/(q; q){2n} [_Bill Gosper, Jun 25 2005]
Also, let B(q) = Sum_{n >= 0} A067659(n) q^n = q + q^2 + q^3 + q^4 + q^5 + 2 q^6 + ... Then B(q) = Sum_{n >= 0} q^((n+1)(2n+1))/(q; q)_{2n+1}.
Also we have the following identity involving 2 X 2 matrices:
Prod_{k >= 1} [ 1, q^k; q^k, 1 ] = [ A(q), B(q); B(q), A(q) ] [Bill Gosper, Jun 25 2005]
a(n) = (A000009(n)+A010815(n))/2. - Vladeta Jovovic, Feb 24 2002
Expansion of (1 + phi(-x)) / (2*chi(-x)) in powers of x where phi(), chi() are Ramanujan theta functions. - Michael Somos, Feb 14 2006
a(n) + A067659(n) = A000009(n). - R. J. Mathar, Jun 18 2016
a(n) ~ exp(Pi*sqrt(n/3)) / (8*3^(1/4)*n^(3/4)). - Vaclav Kotesovec, May 24 2018
A000009(n) = a(n) + A067659(n). - Gus Wiseman, Jan 09 2021
From Peter Bala, Feb 05 2021: (Start)
G.f.: A(x) = (1/2)*((Product_{n >= 0} 1 + x^n) + (Product_{n >= 0} 1 - x^n)).
Let B(x) denote the g.f. of A067659. Then
A(x)^2 - B(x)^2 = A(x^2) - B(x^2) = Product_{n >= 1} 1 - x^(2*n) = Sum_{n in Z} (-1)^n*x^(n*(3*n+1)).
A(x) + B(x) is the g.f. of A000009.
1/(A(x) - B(x)) is the g.f. of A000041.
(A(x) + B(x))/(A(x) - B(x)) is the g.f. of A015128.
A(x)/(A(x) + B(x)) = Sum_{n >= 0} (-1)^n*x^n^2 = (1 + theta_3(-x))/2.
B(x)/(A(x) - B(x)) is the g.f. of A014968.
A(x)/(A(x^2) - B(x^2)) is the g.f. of A027187.
B(x)/(A(x^2) - B(x^2)) is the g.f. of A027193. (End)

A236913 Number of partitions of 2n of type EE (see Comments).

Original entry on oeis.org

1, 1, 3, 6, 12, 22, 40, 69, 118, 195, 317, 505, 793, 1224, 1867, 2811, 4186, 6168, 9005, 13026, 18692, 26613, 37619, 52815, 73680, 102162, 140853, 193144, 263490, 357699, 483338, 650196, 870953, 1161916, 1544048, 2044188, 2696627, 3545015, 4644850, 6066425
Offset: 0

Views

Author

Clark Kimberling, Feb 01 2014

Keywords

Comments

The partitions of n are partitioned into four types:
EO, even # of odd parts and odd # of even parts, A236559;
OE, odd # of odd parts and even # of even parts, A160786;
EE, even # of odd parts and even # of even parts, A236913;
OO, odd # of odd parts and odd # of even parts, A236914.
A236559 and A160786 are the bisections of A027193;
A236913 and A236914 are the bisections of A027187.

Examples

			The partitions of 4 of type EE are [3,1], [2,2], [1,1,1,1], so that a(2) = 3.
type/k . 1 .. 2 .. 3 .. 4 .. 5 .. 6 .. 7 .. 8 ... 9 ... 10 .. 11
EO ..... 0 .. 1 .. 0 .. 2 .. 0 .. 5 .. 0 .. 10 .. 0 ... 20 .. 0
OE ..... 1 .. 0 .. 2 .. 0 .. 4 .. 0 .. 8 .. 0 ... 16 .. 0 ... 29
EE ..... 0 .. 1 .. 0 .. 3 .. 0 .. 6 .. 0 .. 12 .. 0 ... 22 .. 0
OO ..... 0 .. 0 .. 1 .. 0 .. 3 .. 0 .. 7 .. 0 ... 14 .. 0 ... 27
From _Gus Wiseman_, Feb 09 2021: (Start)
This sequence counts even-length partitions of even numbers, which have Heinz numbers given by A340784. For example, the a(0) = 1 through a(4) = 12 partitions are:
  ()  (11)  (22)    (33)      (44)
            (31)    (42)      (53)
            (1111)  (51)      (62)
                    (2211)    (71)
                    (3111)    (2222)
                    (111111)  (3221)
                              (3311)
                              (4211)
                              (5111)
                              (221111)
                              (311111)
                              (11111111)
(End)
		

Crossrefs

Note: A-numbers of ranking sequences are in parentheses below.
The ordered version is A000302.
The case of odd-length partitions of odd numbers is A160786 (A340931).
The Heinz numbers of these partitions are (A340784).
A027187 counts partitions of even length/maximum (A028260/A244990).
A034008 counts compositions of even length.
A035363 counts partitions into even parts (A066207).
A047993 counts balanced partitions (A106529).
A058695 counts partitions of odd numbers (A300063).
A058696 counts partitions of even numbers (A300061).
A067661 counts strict partitions of even length (A030229).
A072233 counts partitions by sum and length.
A339846 counts factorizations of even length.
A340601 counts partitions of even rank (A340602).
A340785 counts factorizations into even factors.
A340786 counts even-length factorizations into even factors.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, [1, 0$3],
          `if`(i<1, [0$4], b(n, i-1)+`if`(i>n, [0$4], (p->
          `if`(irem(i, 2)=0, [p[3], p[4], p[1], p[2]],
              [p[2], p[1], p[4], p[3]]))(b(n-i, i)))))
        end:
    a:= n-> b(2*n$2)[1]:
    seq(a(n), n=0..40);  # Alois P. Heinz, Feb 16 2014
  • Mathematica
    z = 25; m1 = Map[Length[Select[Map[{Count[#, True], Count[#, False]} &,  OddQ[IntegerPartitions[2 #]]], EvenQ[(*Odd*)First[#]] && OddQ[(*Even*)Last[#]] &]] &, Range[z]]; m2 = Map[Length[Select[Map[{Count[#, True], Count[#, False]} &,      OddQ[IntegerPartitions[2 # - 1]]], OddQ[(*Odd*)First[#]] && EvenQ[(*Even*)Last[#]] &]] &, Range[z]]; m3 = Map[Length[Select[Map[{Count[#, True], Count[#, False]} &,
    OddQ[IntegerPartitions[2 #]]], EvenQ[(*Odd*)First[#]] && EvenQ[(*Even*)Last[#]] &]] &, Range[z]] ; m4 = Map[Length[Select[Map[{Count[#, True], Count[#, False]} &,
    OddQ[IntegerPartitions[2 # - 1]]], OddQ[(*Odd*)First[#]] && OddQ[(*Even*)Last[#]] &]] &, Range[z]];
    m1 (* A236559, type EO*)
    m2 (* A160786, type OE*)
    m3 (* A236913, type EE*)
    m4 (* A236914, type OO*)
    (* Peter J. C. Moses, Feb 03 2014 *)
    b[n_, i_] := b[n, i] = If[n == 0, {1, 0, 0, 0}, If[i < 1, {0, 0, 0, 0}, b[n, i - 1] + If[i > n, {0, 0, 0, 0}, Function[p, If[Mod[i, 2] == 0, p[[{3, 4, 1, 2}]], p[[{2, 1, 4, 3}]]]][b[n - i, i]]]]]; a[n_] := b[2*n, 2*n][[1]]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Oct 27 2015, after Alois P. Heinz *)
    Table[Length[Select[IntegerPartitions[2n],EvenQ[Length[#]]&]],{n,0,15}] (* Gus Wiseman, Feb 09 2021 *)

Extensions

More terms from Alois P. Heinz, Feb 16 2014

A338899 Concatenated sequence of prime indices of squarefree semiprimes (A006881).

Original entry on oeis.org

1, 2, 1, 3, 1, 4, 2, 3, 2, 4, 1, 5, 1, 6, 2, 5, 1, 7, 3, 4, 1, 8, 2, 6, 1, 9, 2, 7, 3, 5, 2, 8, 1, 10, 1, 11, 3, 6, 2, 9, 1, 12, 4, 5, 1, 13, 3, 7, 1, 14, 2, 10, 4, 6, 2, 11, 1, 15, 3, 8, 1, 16, 2, 12, 3, 9, 1, 17, 4, 7, 1, 18, 2, 13, 2, 14, 4, 8, 1, 19, 2, 15
Offset: 1

Views

Author

Gus Wiseman, Nov 16 2020

Keywords

Comments

This is a triangle with two columns and strictly increasing rows, namely {A270650(n), A270652(n)}.
A squarefree semiprime is a product of any two distinct prime numbers. A prime index of n is a number m such that the m-th prime number divides n. The multiset of prime indices of n is row n of A112798.

Examples

			The sequence of terms together with their prime indices begins:
      6: {1,2}     57: {2,8}     106: {1,16}    155: {3,11}
     10: {1,3}     58: {1,10}    111: {2,12}    158: {1,22}
     14: {1,4}     62: {1,11}    115: {3,9}     159: {2,16}
     15: {2,3}     65: {3,6}     118: {1,17}    161: {4,9}
     21: {2,4}     69: {2,9}     119: {4,7}     166: {1,23}
     22: {1,5}     74: {1,12}    122: {1,18}    177: {2,17}
     26: {1,6}     77: {4,5}     123: {2,13}    178: {1,24}
     33: {2,5}     82: {1,13}    129: {2,14}    183: {2,18}
     34: {1,7}     85: {3,7}     133: {4,8}     185: {3,12}
     35: {3,4}     86: {1,14}    134: {1,19}    187: {5,7}
     38: {1,8}     87: {2,10}    141: {2,15}    194: {1,25}
     39: {2,6}     91: {4,6}     142: {1,20}    201: {2,19}
     46: {1,9}     93: {2,11}    143: {5,6}     202: {1,26}
     51: {2,7}     94: {1,15}    145: {3,10}    203: {4,10}
     55: {3,5}     95: {3,8}     146: {1,21}    205: {3,13}
		

Crossrefs

A270650 is the first column.
A270652 is the second column.
A320656 counts multiset partitions using these rows, or factorizations into squarefree semiprimes.
A338898 is the version including squares, with columns A338912 and A338913.
A338900 gives row differences.
A338901 gives the row numbers for first appearances.
A001221 and A001222 count distinct/all prime indices.
A001358 lists semiprimes.
A004526 counts 2-part partitions, with strict case shifted right once.
A005117 lists squarefree numbers.
A006881 lists squarefree semiprimes.
A046315 and A100484 list odd and even semiprimes.
A046388 lists odd squarefree semiprimes.
A166237 gives first differences of squarefree semiprimes.

Programs

  • Mathematica
    Join@@Cases[Select[Range[100],SquareFreeQ[#]&&PrimeOmega[#]==2&],k_:>PrimePi/@First/@FactorInteger[k]]

A000379 Numbers where total number of 1-bits in the exponents of their prime factorization is even; a 2-way classification of integers: complement of A000028.

Original entry on oeis.org

1, 6, 8, 10, 12, 14, 15, 18, 20, 21, 22, 26, 27, 28, 32, 33, 34, 35, 36, 38, 39, 44, 45, 46, 48, 50, 51, 52, 55, 57, 58, 62, 63, 64, 65, 68, 69, 74, 75, 76, 77, 80, 82, 85, 86, 87, 91, 92, 93, 94, 95, 98, 99, 100, 106, 111, 112, 115, 116, 117, 118, 119, 120, 122, 123, 124, 125, 129
Offset: 1

Views

Author

Keywords

Comments

This sequence and A000028 (its complement) give the unique solution to the problem of splitting the positive integers into two classes in such a way that products of pairs of distinct elements from either class occur with the same multiplicities [Lambek and Moser]. Cf. A000069, A001969.
See A000028 for precise definition, Maple program, etc.
The sequence contains products of even number of distinct terms of A050376. - Vladimir Shevelev, May 04 2010
From Vladimir Shevelev, Oct 28 2013: (Start)
Numbers m such that the infinitary Möbius function (A064179) of m equals 1. (This follows from the definition of A064179.)
A number m is in the sequence iff the number k = k(m) of terms of A050376 that divide m with odd maximal exponent is even (see example).
(End)
Numbers k for which A064547(k) [or equally, A268386(k)] is even. Numbers k for which A010060(A268387(k)) = 0. - Antti Karttunen, Feb 09 2016
The sequence is closed under the commutative binary operation A059897(.,.). As integers are self-inverse under A059897(.,.), it therefore forms a subgroup of the positive integers considered as a group under A059897(.,.). Specifically (expanding on the comment above dated May 04 2010) it is the subgroup of even length words in A050376, which is the group's lexicographically earliest ordered minimal set of generators. A000028, the set of odd length words in A050376, is its complementary coset. - Peter Munn, Nov 01 2019
From Amiram Eldar, Oct 02 2024: (Start)
Numbers whose number of infinitary divisors (A037445) is a square.
Numbers whose exponentially odious part (A367514) has an even number of distinct prime factors, i.e., numbers k such that A092248(A367514(k)) = 0. (End)

Examples

			If m = 120, then the maximal exponent of 2 that divides 120 is 3, for 3 it is 1, for 4 it is 1, for 5 it is 1. Thus k(120) = 4 and 120 is a term. - _Vladimir Shevelev_, Oct 28 2013
		

References

  • Joe Roberts, Lure of the Integers, Math. Assoc. America, 1992, p. 22.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Subsequences: A030229, A238748, A262675, A268390.
Subsequence of A268388 (apart from the initial 1).
Complement: A000028.
Sequences used in definitions of this sequence: A133008, A050376, A059897, A064179, A064547, A124010 (prime exponents), A268386, A268387, A010060.
Other 2-way classifications: A000069/A001969 (to which A000120 and A010060 are relevant), A000201/A001950.
This is different from A123240 (e.g., does not contain 180). The first difference occurs already at n=31, where A123240(31) = 60, a value which does not occur here, as a(31+1) = 62. The same is true with respect to A131181, as A131181(31) = 60.

Programs

  • Haskell
    a000379 n = a000379_list !! (n-1)
    a000379_list = filter (even . sum . map a000120 . a124010_row) [1..]
    -- Reinhard Zumkeller, Oct 05 2011
    
  • Mathematica
    Select[ Range[130], EvenQ[ Count[ Flatten[ IntegerDigits[#, 2]& /@ Transpose[ FactorInteger[#]][[2]]], 1]]&] // Prepend[#, 1]& (* Jean-François Alcover, Apr 11 2013, after Harvey P. Dale *)
  • PARI
    is(n)=my(f=factor(n)[,2]); sum(i=1,#f,hammingweight(f[i]))%2==0 \\ Charles R Greathouse IV, Aug 31 2013
    (Scheme, two variants)
    (define A000379 (MATCHING-POS 1 1 (COMPOSE even? A064547)))
    (define A000379 (MATCHING-POS 1 1 (lambda (n) (even? (A000120 (A268387 n))))))
    ;; Both require also my IntSeq-library. - Antti Karttunen, Feb 09 2016

Extensions

Edited by N. J. A. Sloane, Dec 20 2007, to restore the original definition.

A104141 Decimal expansion of 3/Pi^2.

Original entry on oeis.org

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

Views

Author

Lekraj Beedassy, Mar 07 2005

Keywords

Comments

3/Pi^2 is the limit of (Sum_{k=1..n} phi(k))/n^2, where phi(k) is Euler's totient A000010(k), i.e., of A002088(n)/A000290(n) as n tends to infinity.
The previous comment in the context of Farey series means that the length of the n-th Farey series can be approximated by multiplying this constant by n^2, "and that the approximation gets proportionally better as n gets larger", according to Conway and Guy. - Alonso del Arte, May 28 2011
The asymptotic density of the sequences of squarefree numbers with even number of prime factors (A030229), odd number of prime factors (A030059), and coprime to 6 (A276378). - Amiram Eldar, May 22 2020

Examples

			3/Pi^2 = 0.303963550927013314331638389629...
		

References

  • J. H. Conway and R. K. Guy, The Book of Numbers, New York: Springer-Verlag, 1995, p. 156.
  • L. E. Dickson, History of the Theory of Numbers, Vol. I pp. 126 Chelsea NY 1966.
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 184.

Crossrefs

Programs

Formula

Equals Sum_{n>=1} 1/A039956(n)^2. - Amiram Eldar, May 22 2020
From Terry D. Grant, Oct 31 2020: (Start)
Equals (-1)*zeta(0)/zeta(2).
Equals 1/(zeta(2)/2).
Equals 1/A195055.
Equals (1/2)*Sum_{k>=1} mu(k)/k^2. (End)
From Hugo Pfoertner, Apr 23 2024: (Start)
Equals A059956/2.
Equals A082020/5. (End)

Extensions

More terms from Ryan Propper, Aug 04 2005

A340601 Number of integer partitions of n of even rank.

Original entry on oeis.org

1, 1, 0, 3, 1, 5, 3, 11, 8, 18, 16, 34, 33, 57, 59, 98, 105, 159, 179, 262, 297, 414, 478, 653, 761, 1008, 1184, 1544, 1818, 2327, 2750, 3480, 4113, 5137, 6078, 7527, 8899, 10917, 12897, 15715, 18538, 22431, 26430, 31805, 37403, 44766, 52556, 62620, 73379
Offset: 0

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. For this sequence, the rank of an empty partition is 0.

Examples

			The a(1) = 1 through a(9) = 18 partitions (empty column indicated by dot):
  (1)  .  (3)    (22)  (5)      (42)    (7)        (44)      (9)
          (21)         (41)     (321)   (43)       (62)      (63)
          (111)        (311)    (2211)  (61)       (332)     (81)
                       (2111)           (322)      (521)     (333)
                       (11111)          (331)      (2222)    (522)
                                        (511)      (4211)    (531)
                                        (2221)     (32111)   (711)
                                        (4111)     (221111)  (4221)
                                        (31111)              (4311)
                                        (211111)             (6111)
                                        (1111111)            (32211)
                                                             (33111)
                                                             (51111)
                                                             (222111)
                                                             (411111)
                                                             (3111111)
                                                             (21111111)
                                                             (111111111)
		

Crossrefs

Note: Heinz numbers are given in parentheses below.
The positive case is A101708 (A340605).
The Heinz numbers of these partitions are A340602.
The odd version is A340692 (A340603).
- Rank -
A047993 counts partitions of rank 0 (A106529).
A072233 counts partitions by sum and length.
A101198 counts partitions of rank 1 (A325233).
A101707 counts partitions of odd positive rank (A340604).
A101708 counts partitions of even positive rank (A340605).
A257541 gives the rank of the partition with Heinz number n.
A340653 counts factorizations of rank 0.
- Even -
A024430 counts set partitions of even length.
A027187 counts partitions of even length (A028260).
A027187 (also) counts partitions of even maximum (A244990).
A034008 counts compositions of even length.
A035363 counts partitions into even parts (A066207).
A052841 counts ordered set partitions of even length.
A058696 counts partitions of even numbers (A300061).
A067661 counts strict partitions of even length (A030229).
A236913 counts even-length partitions of even numbers (A340784).
A339846 counts factorizations of even length.

Programs

  • Maple
    b:= proc(n, i, r) option remember; `if`(n=0, 1-max(0, r),
          `if`(i<1, 0, b(n, i-1, r) +b(n-i, min(n-i, i), 1-
          `if`(r<0, irem(i, 2), r))))
        end:
    a:= n-> b(n$2, -1):
    seq(a(n), n=0..55);  # Alois P. Heinz, Jan 22 2021
  • Mathematica
    Table[If[n==0,1,Length[Select[IntegerPartitions[n],EvenQ[Max[#]-Length[#]]&]]],{n,0,30}]
    (* Second program: *)
    b[n_, i_, r_] := b[n, i, r] = If[n == 0, 1 - Max[0, r], If[i < 1, 0, b[n, i - 1, r] + b[n - i, Min[n - i, i], 1 - If[r < 0, Mod[i, 2], r]]]];
    a[n_] := b[n, n, -1];
    a /@ Range[0, 55] (* Jean-François Alcover, May 10 2021, after Alois P. Heinz *)
  • PARI
    p_q(k) = {prod(j=1, k, 1-q^j); }
    GB_q(N, M)= {if(N>=0 && M>=0,  p_q(N+M)/(p_q(M)*p_q(N)), 0 ); }
    A_q(N) = {my(q='q+O('q^N), g=1+sum(i=1,N, sum(j=1,N/i, q^(i*j) * ( ((1/2)*(1+(-1)^(i+j))) + sum(k=1,N-(i*j), ((q^k)*GB_q(k,i-2)) * ((1/2)*(1+(-1)^(i+j+k)))))))); Vec(g)}
    A_q(50) \\ John Tyler Rascoe, Apr 15 2024

Formula

G.f.: 1 + Sum_{i, j>0} q^(i*j) * ( (1+(-1)^(i+j))/2 + Sum_{k>0} q^k * q_binomial(k,i-2) * (1+(-1)^(i+j+k))/2 ). - John Tyler Rascoe, Apr 15 2024
a(n) ~ exp(Pi*sqrt(2*n/3)) / (8*n*sqrt(3)). - Vaclav Kotesovec, Apr 17 2024

A339560 Number of integer partitions of n that can be partitioned into distinct pairs of distinct parts, i.e., into a set of edges.

Original entry on oeis.org

1, 0, 0, 1, 1, 2, 2, 4, 5, 8, 8, 13, 17, 22, 28, 39, 48, 62, 81, 101, 127, 167, 202, 253, 318, 395, 486, 608, 736, 906, 1113, 1353, 1637, 2011, 2409, 2922, 3510, 4227, 5060, 6089, 7242, 8661, 10306, 12251, 14503, 17236, 20345, 24045, 28334, 33374, 39223, 46076
Offset: 0

Views

Author

Gus Wiseman, Dec 10 2020

Keywords

Comments

Naturally, such a partition must have an even number of parts. Its multiplicities form a graphical partition (A000569, A320922), and vice versa.

Examples

			The a(3) = 1 through a(11) = 13 partitions (A = 10):
  (21)  (31)  (32)  (42)  (43)    (53)    (54)    (64)    (65)
              (41)  (51)  (52)    (62)    (63)    (73)    (74)
                          (61)    (71)    (72)    (82)    (83)
                          (3211)  (3221)  (81)    (91)    (92)
                                  (4211)  (3321)  (4321)  (A1)
                                          (4221)  (5221)  (4322)
                                          (4311)  (5311)  (4331)
                                          (5211)  (6211)  (4421)
                                                          (5321)
                                                          (5411)
                                                          (6221)
                                                          (6311)
                                                          (7211)
For example, the partition y = (4,3,3,2,1,1) can be partitioned into a set of edges in two ways:
  {{1,2},{1,3},{3,4}}
  {{1,3},{1,4},{2,3}},
so y is counted under a(14).
		

Crossrefs

A338916 allows equal pairs (x,x).
A339559 counts the complement in even-length partitions.
A339561 gives the Heinz numbers of these partitions.
A339619 counts factorizations of the same type.
A000070 counts non-multigraphical partitions of 2n, ranked by A339620.
A000569 counts graphical partitions, ranked by A320922.
A001358 lists semiprimes, with squarefree case A006881.
A002100 counts partitions into squarefree semiprimes.
A058696 counts partitions of even numbers, ranked by A300061.
A209816 counts multigraphical partitions, ranked by A320924.
A320655 counts factorizations into semiprimes.
A320656 counts factorizations into squarefree semiprimes.
A339617 counts non-graphical partitions of 2n, ranked by A339618.
A339655 counts non-loop-graphical partitions of 2n, ranked by A339657.
A339656 counts loop-graphical partitions, ranked by A339658.
A339659 counts graphical partitions of 2n into k parts.
The following count partitions of even length and give their Heinz numbers:
- A027187 has no additional conditions (A028260).
- A096373 cannot be partitioned into strict pairs (A320891).
- A338914 can be partitioned into strict pairs (A320911).
- A338915 cannot be partitioned into distinct pairs (A320892).
- A338916 can be partitioned into distinct pairs (A320912).
- A339559 cannot be partitioned into distinct strict pairs (A320894).

Programs

  • Mathematica
    strs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],And[SquareFreeQ[#],PrimeOmega[#]==2]&]}]];
    Table[Length[Select[IntegerPartitions[n],strs[Times@@Prime/@#]!={}&]],{n,0,15}]

Formula

A027187(n) = a(n) + A339559(n).

Extensions

More terms from Jinyuan Wang, Feb 14 2025
Previous Showing 11-20 of 73 results. Next