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

A089633 Numbers having no more than one 0 in their binary representation.

Original entry on oeis.org

0, 1, 2, 3, 5, 6, 7, 11, 13, 14, 15, 23, 27, 29, 30, 31, 47, 55, 59, 61, 62, 63, 95, 111, 119, 123, 125, 126, 127, 191, 223, 239, 247, 251, 253, 254, 255, 383, 447, 479, 495, 503, 507, 509, 510, 511, 767, 895, 959, 991, 1007, 1015, 1019, 1021, 1022, 1023
Offset: 0

Views

Author

Reinhard Zumkeller, Jan 01 2004

Keywords

Comments

Complement of A158582. - Reinhard Zumkeller, Apr 16 2009
Also union of A168604 and A030130. - Douglas Latimer, Jul 19 2012
Numbers of the form 2^t - 2^k - 1, 0 <= k < t.
n is in the sequence if and only if 2*n+1 is in the sequence. - Robert Israel, Dec 14 2018
Also the least binary rank of a strict integer partition of n, where the binary rank of a partition y is given by Sum_i 2^(y_i-1). - Gus Wiseman, May 24 2024

Examples

			From _Tilman Piesk_, May 09 2012: (Start)
This may also be viewed as a triangle:             In binary:
                  0                                         0
               1     2                                 01       10
             3    5    6                          011      101      110
           7   11   13   14                  0111     1011     1101     1110
        15   23   27   29   30          01111    10111    11011    11101    11110
      31  47   55   59   61   62
   63   95  111  119  123  125  126
Left three diagonals are A000225,  A055010, A086224. Right diagonal is A000918. Central column is A129868. Numbers in row n (counted from 0) have n binary 1s. (End)
From _Gus Wiseman_, May 24 2024: (Start)
The terms together with their binary expansions and binary indices begin:
   0:      0 ~ {}
   1:      1 ~ {1}
   2:     10 ~ {2}
   3:     11 ~ {1,2}
   5:    101 ~ {1,3}
   6:    110 ~ {2,3}
   7:    111 ~ {1,2,3}
  11:   1011 ~ {1,2,4}
  13:   1101 ~ {1,3,4}
  14:   1110 ~ {2,3,4}
  15:   1111 ~ {1,2,3,4}
  23:  10111 ~ {1,2,3,5}
  27:  11011 ~ {1,2,4,5}
  29:  11101 ~ {1,3,4,5}
  30:  11110 ~ {2,3,4,5}
  31:  11111 ~ {1,2,3,4,5}
  47: 101111 ~ {1,2,3,4,6}
  55: 110111 ~ {1,2,3,5,6}
  59: 111011 ~ {1,2,4,5,6}
  61: 111101 ~ {1,3,4,5,6}
  62: 111110 ~ {2,3,4,5,6}
(End)
		

Crossrefs

Cf. A181741 (primes), union of A081118 and A000918, apart from initial -1.
For least binary index (instead of rank) we have A001511.
Applying A019565 (Heinz number of binary indices) gives A077011.
For greatest binary index we have A029837 or A070939, opposite A070940.
Row minima of A118462 (binary ranks of strict partitions).
For sum instead of minimum we have A372888, non-strict A372890.
A000009 counts strict partitions, ranks A005117.
A048675 gives binary rank of prime indices, distinct A087207.
A048793 lists binary indices, product A096111, reverse A272020.
A277905 groups all positive integers by binary rank of prime indices.

Programs

  • Haskell
    a089633 n = a089633_list !! (n-1)
    a089633_list = [2 ^ t - 2 ^ k - 1 | t <- [1..], k <- [t-1,t-2..0]]
    -- Reinhard Zumkeller, Feb 23 2012
    
  • Maple
    seq(seq(2^a-1-2^b,b=a-1..0,-1),a=1..11); # Robert Israel, Dec 14 2018
  • Mathematica
    fQ[n_] := DigitCount[n, 2, 0] < 2; Select[ Range[0, 2^10], fQ] (* Robert G. Wilson v, Aug 02 2012 *)
  • PARI
    {insq(n) = local(dd, hf, v); v=binary(n);hf=length(v);dd=sum(i=1,hf,v[i]);if(dd<=hf-2,-1,1)}
    {for(w=0,1536,if(insq(w)>=0,print1(w,", ")))}
    \\ Douglas Latimer, May 07 2013
    
  • PARI
    isoka(n) = #select(x->(x==0), binary(n)) <= 1; \\ Michel Marcus, Dec 14 2018
    
  • Python
    from itertools import count, islice
    def A089633_gen(): # generator of terms
        return ((1<A089633_list = list(islice(A089633_gen(),30)) # Chai Wah Wu, Feb 10 2023
    
  • Python
    from math import isqrt, comb
    def A089633(n): return (1<<(a:=(isqrt((n<<3)+1)-1>>1)+1))-(1<Chai Wah Wu, Dec 19 2024

Formula

A023416(a(n)) <= 1; A023416(a(n)) = A023532(n-2) for n>1;
A000120(a(u)) <= A000120(a(v)) for uA000120(a(n)) = A003056(n).
a(0)=0, n>0: a(n+1) = Min{m>n: BinOnes(a(n))<=BinOnes(m)} with BinOnes=A000120.
If m = floor((sqrt(8*n+1) - 1) / 2), then a(n) = 2^(m+1) - 2^(m*(m+3)/2 - n) - 1. - Carl R. White, Feb 10 2009
A029931(a(n)) = n and A029931(m) != n for m < a(n). - Reinhard Zumkeller, Feb 28 2014
A265705(a(n),k) = A265705(a(n),a(n)-k), k = 0 .. a(n). - Reinhard Zumkeller, Dec 15 2015
a(A014132(n)-1) = 2*a(n-1)+1 for n >= 1. - Robert Israel, Dec 14 2018
Sum_{n>=1} 1/a(n) = A065442 + A160502 = 3.069285887459... . - Amiram Eldar, Jan 09 2024
A019565(a(n)) = A077011(n). - Gus Wiseman, May 24 2024

A246867 Triangle T(n,k) in which n-th row lists in increasing order all partitions lambda of n into distinct parts encoded as Product_{i:lambda} prime(i); n>=0, 1<=k<=A000009(n).

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 14, 15, 13, 21, 22, 30, 17, 26, 33, 35, 42, 19, 34, 39, 55, 66, 70, 23, 38, 51, 65, 77, 78, 105, 110, 29, 46, 57, 85, 91, 102, 130, 154, 165, 210, 31, 58, 69, 95, 114, 119, 143, 170, 182, 195, 231, 330, 37, 62, 87, 115, 133, 138, 187
Offset: 0

Views

Author

Alois P. Heinz, Sep 05 2014

Keywords

Comments

The concatenation of all rows (with offset 1) gives a permutation of the squarefree numbers A005117. The missing positive numbers are in A013929.

Examples

			The partitions of n=5 into distinct parts are {[5], [4,1], [3,2]}, encodings give {prime(5), prime(4)*prime(1), prime(3)*prime(2)} = {11, 7*2, 5*3} => row 5 = [11, 14, 15].
For n=0 the empty partition [] gives the empty product 1.
Triangle T(n,k) begins:
   1;
   2;
   3;
   5,  6;
   7, 10;
  11, 14, 15;
  13, 21, 22, 30;
  17, 26, 33, 35, 42;
  19, 34, 39, 55, 66,  70;
  23, 38, 51, 65, 77,  78, 105, 110;
  29, 46, 57, 85, 91, 102, 130, 154, 165, 210;
  ...
Corresponding triangle of strict integer partitions begins:
                  0
                 (1)
                 (2)
               (3) (21)
               (4) (31)
             (5) (41) (32)
          (6) (42) (51) (321)
        (7) (61) (52) (43) (421)
     (8) (71) (62) (53) (521) (431)
(9) (81) (72) (63) (54) (621) (432) (531). - _Gus Wiseman_, Feb 23 2018
		

Crossrefs

Column k=1 gives: A008578(n+1).
Last elements of rows give: A246868.
Row sums give A147655.
Row lengths are: A000009.
Cf. A005117, A118462, A215366 (the same for all partitions), A258323, A299755, A299757, A299759.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, [1], `if`(i<1, [], [seq(
          map(p->p*ithprime(i)^j, b(n-i*j, i-1))[], j=0..min(1, n/i))]))
        end:
    T:= n-> sort(b(n$2))[]:
    seq(T(n), n=0..14);
  • Mathematica
    b[n_, i_] := b[n, i] = If[n==0, {1}, If[i<1, {}, Flatten[Table[Map[ #*Prime[i]^j&, b[n-i*j, i-1]], {j, 0, Min[1, n/i]}]]]]; T[n_] := Sort[b[n, n]]; Table[T[n], {n, 0, 14}] // Flatten (* Jean-François Alcover, Dec 18 2016, after Alois P. Heinz *)

A066189 Sum of all partitions of n into distinct parts.

Original entry on oeis.org

0, 1, 2, 6, 8, 15, 24, 35, 48, 72, 100, 132, 180, 234, 308, 405, 512, 646, 828, 1026, 1280, 1596, 1958, 2392, 2928, 3550, 4290, 5184, 6216, 7424, 8880, 10540, 12480, 14784, 17408, 20475, 24048, 28120, 32832, 38298, 44520, 51660, 59892, 69230, 79904
Offset: 0

Views

Author

Wouter Meeussen, Dec 15 2001

Keywords

Examples

			The strict integer partitions of 6 are {(6), (5,1), (4,2), (3,2,1)} with sum 6+5+1+4+2+3+2+1 = 24. - _Gus Wiseman_, May 09 2019
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, [1, 0], `if`(i>n, [0$2],
          b(n, i+1)+(p-> p+[0, i*p[1]])(b(n-i, i+1))))
        end:
    a:= n-> b(n, 1)[2]:
    seq(a(n), n=0..80);  # Alois P. Heinz, Sep 01 2014
  • Mathematica
    PartitionsQ[ Range[ 60 ] ]Range[ 60 ]
    nmax=60; CoefficientList[Series[x*D[Product[1+x^k, {k, 1, nmax}], x], {x, 0, nmax}], x] (* Vaclav Kotesovec, Nov 21 2016 *)

Formula

G.f.: sum(n>=1, n*q^(n-1)/(1+q^n) ) * prod(n>=1, 1+q^n ). - Joerg Arndt, Aug 03 2011
a(n) = n * A000009(n). - Vaclav Kotesovec, Sep 25 2016
G.f.: x*f'(x), where f(x) = Product_{k>=1} (1 + x^k). - Vaclav Kotesovec, Nov 21 2016
a(n) = A056239(A325506(n)). - Gus Wiseman, May 09 2019

A209862 Permutation of nonnegative integers which maps A209642 into ascending order (A209641).

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 11, 13, 14, 15, 16, 17, 18, 20, 24, 19, 21, 25, 22, 26, 28, 23, 27, 29, 30, 31, 32, 33, 34, 36, 40, 48, 35, 37, 41, 49, 38, 42, 50, 44, 52, 56, 39, 43, 51, 45, 53, 57, 46, 54, 58, 60, 47, 55, 59, 61, 62, 63, 64, 65, 66, 68, 72, 80, 96, 67, 69, 73, 81, 97, 70, 74, 82, 98, 76, 84, 100, 88, 104, 112, 71, 75, 83
Offset: 0

Views

Author

Antti Karttunen, Mar 24 2012

Keywords

Comments

Conjecture: For all n, a(A054429(n)) = A054429(a(n)), i.e. A054429 acts as a homomorphism (automorphism) of the cyclic group generated by this permutation. This implies also a weaker conjecture given in A209860.
From Gus Wiseman, Aug 24 2021: (Start)
As a triangle with row lengths 2^n, T(n,k) for n > 0 appears (verified up to n = 2^15) to be the unique nonnegative integer whose binary indices are the k-th subset of {1..n} containing n. Here, a binary index of n (row n of A048793) is any position of a 1 in its reversed binary expansion, and sets are sorted first by length, then lexicographically. For example, the triangle begins:
1
2 3
4 5 6 7
8 9 10 12 11 13 14 15
16 17 18 20 24 19 21 25 22 26 28 23 27 29 30 31
Mathematica: Table[Total[2^(Append[#,n]-1)]&/@Subsets[Range[n-1]],{n,5}]
Row lengths are A000079 (shifted right). Also Column k = 1.
Row sums are A010036.
Using reverse-lexicographic order gives A059893.
Using lexicographic order gives A059894.
Taking binary indices to prime indices gives A339195 (or A019565).
The ordering of sets is A344084.
A version using Heinz numbers is A344085.
(End)

Examples

			From _Gus Wiseman_, Aug 24 2021: (Start)
The terms, their binary expansions, and their binary indices begin:
   0:      ~ {}
   1:    1 ~ {1}
   2:   10 ~ {2}
   3:   11 ~ {1,2}
   4:  100 ~ {3}
   5:  101 ~ {1,3}
   6:  110 ~ {2,3}
   7:  111 ~ {1,2,3}
   8: 1000 ~ {4}
   9: 1001 ~ {1,4}
  10: 1010 ~ {2,4}
  12: 1100 ~ {3,4}
  11: 1011 ~ {1,2,4}
  13: 1101 ~ {1,3,4}
  14: 1110 ~ {2,3,4}
  15: 1111 ~ {1,2,3,4}
(End)
		

Crossrefs

Formula

A258323 Sum T(n,k) over all partitions lambda of n into k distinct parts of Product_{i:lambda} prime(i); triangle T(n,k), n>=0, 0<=k<=A003056(n), read by rows.

Original entry on oeis.org

1, 0, 2, 0, 3, 0, 5, 6, 0, 7, 10, 0, 11, 29, 0, 13, 43, 30, 0, 17, 94, 42, 0, 19, 128, 136, 0, 23, 231, 293, 0, 29, 279, 551, 210, 0, 31, 484, 892, 330, 0, 37, 584, 1765, 852, 0, 41, 903, 2570, 1826, 0, 43, 1051, 4273, 4207, 0, 47, 1552, 6747, 6595, 2310
Offset: 0

Views

Author

Alois P. Heinz, May 26 2015

Keywords

Examples

			T(6,2) = 43 because the partitions of 6 into 2 distinct parts are {[5,1], [4,2]} and prime(5)*prime(1) + prime(4)*prime(2) = 11*2 + 7*3 = 22 + 21 = 43.
Triangle T(n,k) begins:
  1
  0,  2;
  0,  3;
  0,  5,   6;
  0,  7,  10;
  0, 11,  29;
  0, 13,  43,  30;
  0, 17,  94,  42;
  0, 19, 128, 136;
  0, 23, 231, 293;
  0, 29, 279, 551, 210;
		

Crossrefs

Row sums give A147655.
T(n*(n+1)/2,n) = A002110(n).
T(n^2,n) = A321267(n).

Programs

  • Maple
    g:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0, expand(
          add(g(n-i*j, i-1)*(ithprime(i)*x)^j, j=0..min(1, n/i)))))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..degree(p)))(g(n$2)):
    seq(T(n), n=0..20);
  • Mathematica
    g[n_, i_] := g[n, i] = If[n==0, 1, If[i<1, 0, Expand[Sum[g[n-i*j, i-1] * (Prime[i]*x)^j, {j, 0, Min[1, n/i]}]]]]; T[n_] := Function[p, Table[ Coefficient[p, x, i], {i, 0, Exponent[p, x]}]][g[n, n]]; Table[T[n], {n, 0, 20}] // Flatten (* Jean-François Alcover, Jan 06 2017, after Alois P. Heinz *)

A339195 Triangle of squarefree numbers grouped by greatest prime factor, read by rows.

Original entry on oeis.org

1, 2, 3, 6, 5, 10, 15, 30, 7, 14, 21, 35, 42, 70, 105, 210, 11, 22, 33, 55, 66, 77, 110, 154, 165, 231, 330, 385, 462, 770, 1155, 2310, 13, 26, 39, 65, 78, 91, 130, 143, 182, 195, 273, 286, 390, 429, 455, 546, 715, 858, 910, 1001, 1365, 1430, 2002, 2145, 2730, 3003, 4290, 5005, 6006, 10010, 15015, 30030
Offset: 0

Views

Author

Gus Wiseman, Dec 02 2020

Keywords

Comments

Also Heinz numbers of subsets of {1..n} that contain n if n>0, where the Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
A019565 in its triangle form, with each row's terms in increasing order. - Peter Munn, Feb 26 2021
From David James Sycamore, Jan 09 2025: (Start)
Alternative definition, with offset = 1: a(1) = 1. For n>1 if a(n-1) = A002110(k), a(n) = prime(k+1). Otherwise a(n) is the smallest novel squarefree number whose prime factors have already occurred as previous terms.
Permutation of A005117, Squarefree version A379746. (End)

Examples

			Triangle begins:
   1
   2
   3   6
   5  10  15  30
   7  14  21  35  42  70  105  210
		

Crossrefs

A011782 gives row lengths.
A339360 gives row sums.
A008578 (shifted) is column k = 1.
A100484 is column k = 2.
A001748 is column k = 3.
A002110 is column k = 2^(n-1).
A070826 is column k = 2^(n-1) - 1.
A209862 takes prime indices to binary indices in these terms.
A246867 groups squarefree numbers by Heinz weight, with row sums A147655.
A261144 divides the n-th row by prime(n), with row sums A054640.
A339116 is the restriction to semiprimes, with row sums A339194.
A005117 lists squarefree numbers, ordered lexicographically by prime factors: A019565.
A006881 lists squarefree semiprimes.
A072047 counts prime factors of squarefree numbers.
A319246 is the sum of prime indices of the n-th squarefree number.
A329631 lists prime indices of squarefree numbers, reversed: A319247.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes.
Cf. A379746.

Programs

  • Maple
    T:= proc(n) option remember; `if`(n=0, 1, (p-> map(
          x-> x*p, {seq(T(i), i=0..n-1)})[])(ithprime(n)))
        end:
    seq(T(n), n=0..6);  # Alois P. Heinz, Jan 08 2025
  • Mathematica
    Table[Prime[n]*Sort[Times@@Prime/@#&/@Subsets[Range[n-1]]],{n,5}]

Formula

For n > 1, T(n,k) = prime(n) * A261144(n-1,k).
a(n) = A019565(A379770(n)). - Michael De Vlieger, Jan 08 2025

Extensions

Row n=0 (term 1) prepended by Alois P. Heinz, Jan 08 2025

A261144 Irregular triangle of numbers that are squarefree and smooth (row n contains squarefree p-smooth numbers, where p is the n-th prime).

Original entry on oeis.org

1, 2, 1, 2, 3, 6, 1, 2, 3, 5, 6, 10, 15, 30, 1, 2, 3, 5, 6, 7, 10, 14, 15, 21, 30, 35, 42, 70, 105, 210, 1, 2, 3, 5, 6, 7, 10, 11, 14, 15, 21, 22, 30, 33, 35, 42, 55, 66, 70, 77, 105, 110, 154, 165, 210, 231, 330, 385, 462, 770, 1155, 2310, 1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 21, 22, 26, 30, 33, 35, 39, 42
Offset: 1

Views

Author

Jean-François Alcover, Nov 26 2015

Keywords

Comments

If we define a triangle whose n-th row consists of all squarefree numbers whose prime factors are all less than prime(k), we get this same triangle except starting with a row {1}, with offset 1. - Gus Wiseman, Aug 24 2021

Examples

			Triangle begins:
1, 2;                        squarefree and 2-smooth
1, 2, 3, 6;                  squarefree and 3-smooth
1, 2, 3, 5, 6, 10, 15, 30;
1, 2, 3, 5, 6,  7, 10, 14, 15, 21, 30, 35, 42, 70, 105, 210;
...
		

Crossrefs

Cf. A000079 (2-smooth), A003586 (3-smooth), A051037 (5-smooth), A002473 (7-smooth), A018336 (7-smooth & squarefree), A051038 (11-smooth), A087005 (11-smooth & squarefree), A080197 (13-smooth), A087006 (13-smooth & squarefree), A087007 (17-smooth & squarefree), A087008 (19-smooth & squarefree).
Row lengths are A000079.
Rightmost terms (or column k = 2^n) are A002110.
Rows are partial unions of rows of A019565.
Row n is A027750(A002110(n)), i.e., divisors of primorials.
Row sums are A054640.
Column k = 2^n-1 is A070826.
Multiplying row n by prime(n+1) gives A339195, row sums A339360.
A005117 lists squarefree numbers.
A056239 adds up prime indices, row sums of A112798.
A072047 counts prime factors of squarefree numbers.
A246867 groups squarefree numbers by Heinz weight, row sums A147655.
A329631 lists prime indices of squarefree numbers, sums A319246.
A339116 groups squarefree semiprimes by greater factor, sums A339194.

Programs

  • Maple
    b:= proc(n) option remember; `if`(n=0, [1],
          sort(map(x-> [x, x*ithprime(n)][], b(n-1))))
        end:
    T:= n-> b(n)[]:
    seq(T(n), n=1..7);  # Alois P. Heinz, Nov 28 2015
  • Mathematica
    primorial[n_] := Times @@ Prime[Range[n]]; row[n_] := Select[ Divisors[ primorial[n]], SquareFreeQ]; Table[row[n], {n, 1, 10}] // Flatten

Formula

T(n-1,k) = A339195(n,k)/prime(n). - Gus Wiseman, Aug 24 2021

A111059 a(n) = Product_{k=1..n} A005117(k), the product of the first n squarefree positive integers.

Original entry on oeis.org

1, 2, 6, 30, 180, 1260, 12600, 138600, 1801800, 25225200, 378378000, 6432426000, 122216094000, 2566537974000, 56463835428000, 1298668214844000, 33765373585944000, 979195833992376000, 29375875019771280000
Offset: 1

Views

Author

Leroy Quet, Oct 07 2005

Keywords

Comments

Do all terms belong to A242031 (weakly decreasing prime signature)? - Gus Wiseman, May 14 2021

Examples

			Since the first 6 squarefree positive integers are 1, 2, 3, 5, 6, 7, the 6th term of the sequence is 1*2*3*5*6*7 = 1260.
From _Gus Wiseman_, May 14 2021: (Start)
The sequence of terms together with their prime signatures begins:
             1: ()
             2: (1)
             6: (1,1)
            30: (1,1,1)
           180: (2,2,1)
          1260: (2,2,1,1)
         12600: (3,2,2,1)
        138600: (3,2,2,1,1)
       1801800: (3,2,2,1,1,1)
      25225200: (4,2,2,2,1,1)
     378378000: (4,3,3,2,1,1)
    6432426000: (4,3,3,2,1,1,1)
  122216094000: (4,3,3,2,1,1,1,1)
(End)
		

Crossrefs

A005117 lists squarefree numbers.
A006881 lists squarefree semiprimes.
A072047 applies Omega to each squarefree number.
A246867 groups squarefree numbers by Heinz weight (row sums: A147655).
A261144 groups squarefree numbers by smoothness (row sums: A054640).
A319246 gives the sum of prime indices of each squarefree number.
A329631 lists prime indices of squarefree numbers (reversed: A319247).

Programs

  • Mathematica
    Rest[FoldList[Times,1,Select[Range[40],SquareFreeQ]]] (* Harvey P. Dale, Jun 14 2011 *)
  • PARI
    m=30;k=1;for(n=1,m,if(issquarefree(n),print1(k=k*n,",")))

Extensions

More terms from Klaus Brockhaus, Oct 08 2005

A325506 Product of Heinz numbers over all strict integer partitions of n.

Original entry on oeis.org

1, 2, 3, 30, 70, 2310, 180180, 21441420, 6401795400, 200984366583000, 41615822944675980000, 10515527757483671302380000, 4919824049783476260137727416400000, 5158181210492841550866520676965246284000000, 29776760895364738730693151196801613158042403043600000000
Offset: 0

Views

Author

Gus Wiseman, May 07 2019

Keywords

Comments

a(n) is the product of row n of A246867 (squarefree numbers arranged by sum of prime indices).
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			The strict integer partitions of 6 are {(6), (5,1), (4,2), (3,2,1)}, with Heinz numbers {13,22,21,30}, with product 13*22*21*30 = 180180, so a(6) = 180180.
The sequence of terms together with their prime indices begins:
                     1: {}
                     2: {1}
                     3: {2}
                    30: {1,2,3}
                    70: {1,3,4}
                  2310: {1,2,3,4,5}
                180180: {1,1,2,2,3,4,5,6}
              21441420: {1,1,2,2,3,4,4,5,6,7}
            6401795400: {1,1,1,2,2,3,3,4,5,5,6,7,8}
       200984366583000: {1,1,1,2,2,2,3,3,3,4,4,5,5,6,6,7,8,9}
  41615822944675980000: {1,1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,5,5,6,6,7,7,8,9,10}
		

Crossrefs

Programs

  • Mathematica
    Table[Times@@Prime/@(Join@@Select[IntegerPartitions[n],UnsameQ@@#&]),{n,0,15}]

Formula

a(n) = Product_{i = 1..A000009(n)} A246867(n,i).
A001222(a(n)) = A015723(n).
A056239(a(n)) = A066189(n).
A003963(a(n)) = A325504(n).
a(n) = A003963(A325505(n)).

A339360 Sum of all squarefree numbers with greatest prime factor prime(n).

Original entry on oeis.org

1, 2, 9, 60, 504, 6336, 89856, 1645056, 33094656, 801239040, 24246190080, 777550233600, 29697402470400, 1250501433753600, 55083063155097600, 2649111037319577600, 143390180403000115200, 8619643674791667302400, 534710099148093259776000, 36412881178052121329664000
Offset: 0

Views

Author

Gus Wiseman, Dec 04 2020

Keywords

Examples

			The initial terms are:
   1 = 1,
   2 = 2,
   9 = 3 + 6,
  60 = 5 + 10 + 15 + 30.
		

Crossrefs

A010036 takes prime indices here to binary indices, row sums of A209862.
A048672 takes prime indices to binary indices in squarefree numbers.
A054640 divides the n-th term by prime(n), row sums of A261144.
A072047 counts prime factors of squarefree numbers.
A339194 is the restriction to semiprimes, row sums of A339116.
A339195 has this as row sums.
A002110 lists primorials.
A005117 lists squarefree numbers.
A006881 lists squarefree semiprimes.
A056239 is the sum of prime indices of n (Heinz weight).
A246867 groups squarefree numbers by weight, with row sums A147655.
A319246 is the sum of prime indices of the n-th squarefree number.
A319247 lists reversed prime indices of squarefree numbers.
A329631 lists prime indices of squarefree numbers.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes.

Programs

  • Maple
    f:= proc(n) local i;
      `if`(n=0, 1, ithprime(n)) *mul(1+ithprime(i),i=1..n-1)
    end proc:
    map(f, [$0..20]); # Robert Israel, Dec 08 2020
  • Mathematica
    Table[Sum[Times@@Prime/@stn,{stn,Select[Subsets[Range[n]],MemberQ[#,n]&]}],{n,10}]

Formula

For n >= 1, a(n) = A054640(n-1) * prime(n).

Extensions

a(0)=1 prepended by Alois P. Heinz, Jan 08 2025
Showing 1-10 of 24 results. Next