A002033 -id:A002033 - cp's OEIS frontend

A325781 Heinz numbers of complete integer partitions.

1, 2, 4, 6, 8, 12, 16, 18, 20, 24, 30, 32, 36, 40, 42, 48, 54, 56, 60, 64, 72, 80, 84, 90, 96, 100, 108, 112, 120, 126, 128, 132, 140, 144, 150, 160, 162, 168, 176, 180, 192, 198, 200, 210, 216, 220, 224, 234, 240, 252, 256, 260, 264, 270, 280, 288, 294, 300

Offset: 1

Gus Wiseman, May 21 2019

nonn

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

The sum of prime indices of n is A056239(n). A number is in this sequence iff its divisors have sums of prime indices covering an initial interval of nonnegative integers. For example, the divisors of 60 are {1, 2, 3, 4, 5, 6, 10, 12, 15, 20, 30, 60}, with respective sums of prime indices {0, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 7}, so 60 is in the sequence.

			The sequence of terms together with their prime indices begins:
     1: {}
     2: {1}
     4: {1,1}
     6: {1,2}
     8: {1,1,1}
    12: {1,1,2}
    16: {1,1,1,1}
    18: {1,2,2}
    20: {1,1,3}
    24: {1,1,1,2}
    30: {1,2,3}
    32: {1,1,1,1,1}
    36: {1,1,2,2}
    40: {1,1,1,3}
    42: {1,2,4}
    48: {1,1,1,1,2}
    54: {1,2,2,2}
    56: {1,1,1,4}
    60: {1,1,2,3}
    64: {1,1,1,1,1,1}

Cf. A002033, A056239, A103295, A112798, A126796, A188431, A299702, A325684, A325770, A325780, A325782.

normQ[m_]:=Or[m=={},Union[m]==Range[Max[m]]];
hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
Select[Range[1000],normQ[hwt/@Rest[Divisors[#]]]&]

A065036 Product of the cube of a prime (A030078) and a different prime.

Original entry on oeis.org

24, 40, 54, 56, 88, 104, 135, 136, 152, 184, 189, 232, 248, 250, 296, 297, 328, 344, 351, 375, 376, 424, 459, 472, 488, 513, 536, 568, 584, 621, 632, 664, 686, 712, 776, 783, 808, 824, 837, 856, 872, 875, 904, 999, 1016, 1029, 1048, 1096, 1107, 1112

Offset: 1

Alford Arnold, Nov 04 2001

This sequence appears on row 8 of the list illustrated in A064839 and is similar to A054753 which appears on row 6. Previous rows are generated by A000007, A000040, A001248, A006881, A030078 respectively.

Or, the numbers n such that 20=number of perfect partitions of n. - Juri-Stepan Gerasimov, Sep 26 2009

			a(4)= 56 since 56 = 2*2*2*7.

Cf. A000007, A000040, A001248, A002033, A006881, A030078, A054753, A064839, A089233.

Select[ Range[1500], Sort[ Transpose[ FactorInteger[ # ]] [[2]]] == {1, 3} & ]
Module[{upto=1200},Select[(Union[Flatten[{#[[1]]^3 #[[2]],#[[1]]#[[2]]^3}&/@Subsets[Prime[Range[upto/8]],{2}]]]),#<=upto&]] (* Harvey P. Dale, May 23 2015 *)

list(lim)=my(v=List(),t);forprime(p=2,(lim\2)^(1/3),t=p^3; forprime(q=2,lim\t,if(p==q,next);listput(v,t*q)));vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 20 2011

is(n)=my(f=factor(n)[,2]); f==[3,1]~||f==[1,3]~ \\ Charles R Greathouse IV, Oct 15 2015

from sympy import primepi, primerange, integer_nthroot
def A065036(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x-sum(primepi(x//p**3) for p in primerange(integer_nthroot(x,3)[0]+1))+primepi(integer_nthroot(x,4)[0])
    return bisection(f,n,n) # Chai Wah Wu, Feb 21 2025

A002033(a(n)) = 20. - Juri-Stepan Gerasimov, Sep 26 2009
A089233(a(n)) = 3. - Reinhard Zumkeller, Sep 04 2013
A000005(a(n)) = 8. - Altug Alkan, Nov 11 2015

A089723 a(1)=1; for n>1, a(n) gives number of ways to write n as n = x^y, 2 <= x, 1 <= y.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1

Offset: 1

Naohiro Nomoto, Jan 07 2004

This function depends only on the prime signature of n. - Franklin T. Adams-Watters, Mar 10 2006
a(n) is the number of perfect divisors of n. Perfect divisor of n is divisor d such that d^k = n for some k >= 1. a(n) > 1 for perfect powers n = A001597(m) for m > 2. - Jaroslav Krizek, Jan 23 2010
Also the number of uniform perfect integer partitions of n - 1. An integer partition of n is uniform if all parts appear with the same multiplicity, and perfect if every nonnegative integer up to n is the sum of a unique submultiset. The Heinz numbers of these partitions are given by A326037. The a(16) = 3 partitions are: (8,4,2,1), (4,4,4,1,1,1), (1,1,1,1,1,1,1,1,1,1,1,1,1,1,1). - Gus Wiseman, Jun 07 2019
The record values occur at 1 and at 2^A002182(n) for n > 1. - Amiram Eldar, Nov 06 2020

			144 = 2^4 * 3^2, gcd(4,2) = 2, d(2) = 2, so a(144) = 2. The representations are 144^1 and 12^2.
From _Friedjof Tellkamp_, Jun 14 2025: (Start)
n:          1, 2, 3, 4, 5, 6, 7, 8, 9, ...
----------------------------------------------------
1st powers: 1, 1, 1, 1, 1, 1, 1, 1, 1, ... (A000012)
Squares:    1, 0, 0, 1, 0, 0, 0, 0, 1, ... (A010052)
Cubes:      1, 0, 0, 0, 0, 0, 0, 1, 0, ... (A010057)
Quartics:   1, 0, 0, 0, 0, 0, 0, 0, 0, ... (A374016)
...
Sum:       oo, 1, 1, 2, 1, 1, 1, 2, 2, ...
a(1)=1:     1, 1, 1, 2, 1, 1, 1, 2, 2, ... (= this sequence). (End)

N. J. A. Sloane, Table of n, a(n) for n = 1..20000
Solomon W. Golomb, A new arithmetic function of combinatorial significance, J. Number Theory, Vol. 5, No. 3 (1973), pp. 218-223. 1973JNT.....5..218G
Jan Mycielski, Sur les représentations des nombres naturels par des puissances à base et exposant naturels, Colloquium Mathematicum, Vol. 2 (1951), pp. 254-260. See gamma(n).
N. J. A. Sloane, Maple programs for A007916, A278028, A278029, A052409, A089723, A277564

Cf. A000005, A277564, A278028.
Cf. A000961, A002033, A002182, A007916, A047966, A070941, A108917, A126796, A276024, A326035, A326036, A326037.
Cf. A010052, A010057, A374016, A259362.

with(numtheory):
A089723 := proc(n) local t1,t2,g,j;
if n=1 then 1 else
t1:=ifactors(n)[2]; t2:=nops(t1); g := t1[1][2];
for j from 2 to t2 do g:=gcd(g,t1[j][2]); od:
tau(g); fi; end;
[seq(A089723(n),n=1..100)]; # N. J. A. Sloane, Nov 10 2016

Table[DivisorSigma[0, GCD @@ FactorInteger[n][[All, 2]]], {n, 100}] (* Gus Wiseman, Jun 12 2017 *)

a(n) = if (n==1, 1, numdiv(gcd(factor(n)[,2]))); \\ Michel Marcus, Jun 13 2017

from math import gcd
from sympy import factorint, divisor_sigma
def a(n):
    if n == 1: return 1
    e = list(factorint(n).values())
    g = e[0]
    for ei in e[1:]: g = gcd(g, ei)
    return divisor_sigma(g, 0)
print([a(n) for n in range(1, 105)]) # Michael S. Branicky, Jul 15 2021

If n = Product p_i^e_i, a(n) = d(gcd()). - Franklin T. Adams-Watters, Mar 10 2006
Sum_{n=1..m} a(n) = A255165(m) + 1. - Richard R. Forberg, Feb 16 2015
Sum_{n>=2} a(n)/n^s = Sum_{n>=2} 1/(n^s-1) = Sum_{k>=1} (zeta(s*k)-1) for all real s with Re(s) > 1 (Golomb, 1973). - Amiram Eldar, Nov 06 2020
For n > 1, a(n) = Sum_{i=1..floor(n/2)} floor(n^(1/i))-floor((n-1)^(1/i)). - Wesley Ivan Hurt, Dec 08 2020
Sum_{n>=1} (a(n)-1)/n = 1 (Mycielski, 1951). - Amiram Eldar, Jul 15 2021
From Friedjof Tellkamp, Jun 14 2025: (Start)
a(n) = 1 + A259362(n) = 1 + A010052(n) + A010057(n) + A374016(n) + (...), for n > 1.
G.f.: x + Sum_{j>=2, k>=1} x^(j^k). (End)

A034729 a(n) = Sum_{ k, k|n } 2^(k-1).

Original entry on oeis.org

1, 3, 5, 11, 17, 39, 65, 139, 261, 531, 1025, 2095, 4097, 8259, 16405, 32907, 65537, 131367, 262145, 524827, 1048645, 2098179, 4194305, 8390831, 16777233, 33558531, 67109125, 134225995, 268435457, 536887863, 1073741825, 2147516555, 4294968325, 8590000131

Offset: 1

Erich Friedman

nonn

Dirichlet convolution of b_n=1 with c_n = 2^(n-1).
Equals row sums of triangle A143425, & inverse Möbius transform (A051731) of [1, 2, 4, 8, ...]. - Gary W. Adamson, Aug 14 2008
Number of constant multiset partitions of normal multisets of size n, where a multiset is normal if it spans an initial interval of positive integers. - Gus Wiseman, Sep 16 2018

			From _Gus Wiseman_, Sep 16 2018: (Start)
The a(4) = 11 constant multiset partitions:
  (1)(1)(1)(1)
    (11)(11)
    (12)(12)
     (1111)
     (1222)
     (1122)
     (1112)
     (1233)
     (1223)
     (1123)
     (1234)
(End)

T. D. Noe, Table of n, a(n) for n = 1..1000

Cf. A000005, A002033, A003238, A008480, A018783, A047968, A051731.
Cf. A052409, A055895, A078392, A143425, A215366, A245282, A248906.
Cf. A289508.
Sums of the form Sum_{d|n} q^(d-1): this sequence (q=2), A034730 (q=3), A113999 (q=10), A339684 (q=4), A339685 (q=5), A339686 (q=6), A339687 (q=7), A339688 (q=8), A339689 (q=9).

A034729:= func< n | (&+[2^(d-1): d in Divisors(n)]) >;
[A034729(n): n in [1..40]]; // G. C. Greubel, Jun 26 2024

seq(add(2^(k-1),k=numtheory:-divisors(n)), n = 1 .. 100); # Robert Israel, Aug 22 2014

Rest[CoefficientList[Series[Sum[x^k/(1-2*x^k),{k,1,30}],{x,0,30}],x]] (* Vaclav Kotesovec, Sep 08 2014 *)

A034729(n) = sumdiv(n,k,2^(k-1)) \\ Michael B. Porter, Mar 11 2010

{a(n)=polcoeff(sum(m=1,n,2^(m-1)*x^m/(1-x^m +x*O(x^n))),n)}
for(n=1,40,print1(a(n),", ")) \\ Paul D. Hanna, Aug 21 2014

{a(n)=local(A=x+x^2);A=sum(m=1,n,x^m*sumdiv(m,d,1/(1 - x^(m/d) +x*O(x^n))^d) );polcoeff(A,n)}
for(n=1,40,print1(a(n),", ")) \\ Paul D. Hanna, Aug 21 2014

from sympy import divisors
def A034729(n): return sum(1<<(d-1) for d in divisors(n,generator=True)) # Chai Wah Wu, Jul 15 2022

def A034729(n): return sum(2^(k-1) for k in (1..n) if (k).divides(n))
[A034729(n) for n in range(1,41)] # G. C. Greubel, Jun 26 2024

G.f.: Sum_{n>0} x^n/(1-2*x^n). - Vladeta Jovovic, Nov 14 2002
a(n) = 1/2 * A055895(n). - Joerg Arndt, Aug 14 2012
G.f.: Sum_{n>=1} 2^(n-1) * x^n / (1 - x^n). - Paul D. Hanna, Aug 21 2014
G.f.: Sum_{n>=1} x^n * Sum_{d|n} 1/(1 - x^d)^(n/d). - Paul D. Hanna, Aug 21 2014
a(n) ~ 2^(n-1). - Vaclav Kotesovec, Sep 09 2014
a(n) = Sum_{k in row n of A215366} A008480(k) * A000005(A289508(k)). - Gus Wiseman, Sep 16 2018
a(n) = Sum_{c is a composition of n} A000005(gcd(c)). - Gus Wiseman, Sep 16 2018

A188431 The number of n-full sets, F(n).

Original entry on oeis.org

1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 2, 2, 1, 2, 1, 2, 3, 4, 5, 7, 7, 8, 9, 11, 10, 13, 14, 17, 20, 25, 28, 34, 40, 46, 54, 62, 69, 80, 90, 102, 115, 131, 144, 167, 186, 213, 239, 273, 304, 349, 388, 441, 495, 563, 625, 710, 790, 890, 990, 1114, 1232, 1387, 1530, 1713, 1894, 2119, 2330, 2605, 2866, 3192, 3512, 3910, 4289, 4774, 5237, 5809, 6377, 7068, 7739

Offset: 0

Madjid Mirzavaziri, Mar 31 2011

nonn

Let A be a set of positive integers. We say that A is n-full if (sum A)=[n] for a positive integer n, where (sum A) is the set of all positive integers which are a sum of distinct elements of A and [n]={1,2,...,n}. Then F(n) denotes the number of n-full sets.
Also the number of distinct and complete partitions of n, by definition, which are counted by A000009 and A126796. - George Beck, Nov 06 2017
An integer partition of n is complete (see also A325781) if every number from 0 to n is the sum of some submultiset of the parts. The Heinz numbers of these partitions are given by A325986. - Gus Wiseman, May 31 2019

			a(26) = 10, because there are 10 26-full sets: {1,2,4,5,6,8}, {1,2,3,5,7,8}, {1,2,3,5,6,9}, {1,2,3,4,7,9}, {1,2,3,4,6,10}, {1,2,3,4,5,11}, {1,2,4,8,11}, {1,2,4,7,12}, {1,2,4,6,13}, {1,2,3,7,13}.
G.f.: 1 = 1/(1+x) + 1*x/((1+x)*(1+x^2)) + 0*x^2/((1+x)*(1+x^2)*(1+x^3)) + 1*x^3/((1+x)*(1+x^2)*(1+x^3)*(1+x^4)) +...+ a(n)*x^n / Product_{k=1..n+1} (1+x^k) +...

Alois P. Heinz, Table of n, a(n) for n = 0..10000 (terms n=1..1000 from Reinhard Zumkeller)
Mohammad Saleh Dinparvar, Python program
L. Naranjani and M. Mirzavaziri, Full Subsets of N, Journal of Integer Sequences, 14 (2011), Article 11.5.3.
Arash Sal Moslehian, JavaScript p5 program for three different algorithms

Cf. A188429, A188430, A126796.

Cf. A002033, A103295, A108917, A276024, A325763, A325765, A325781, A325782, A325788, A325986, A325790, A325791.

import Data.MemoCombinators (memo2, integral, Memo)
a188431 n = a188431_list !! (n-1)
a188431_list = map
   (\x -> sum [fMemo x i | i <- [a188429 x .. a188430 x]]) [1..] where
   fMemo = memo2 integral integral f
   f _ 1 = 1
   f m i = sum [fMemo (m - i) j |
                j <- [a188429 (m - i) .. min (a188430 (m - i)) (i - 1)]]
-- Reinhard Zumkeller, Aug 06 2015

sums:= proc(s) local i, m;
          m:= max(s[]);
         `if`(m<1, {}, {m, seq([i, i+m][], i=sums(s minus {m}))})
       end:
a:= proc(n) local b;
      b:= proc(i,s) local si;
            if i=1 then `if`(sums(s)={$1..n}, 1, 0)
          else si:= s union {i};
               b(i-1, s)+ `if`(max(sums(si)[])>n, 0, b(i-1, si))
            fi
          end; b(n, {1})
    end:
seq(a(n), n=1..40);  # Alois P. Heinz, Apr 03 2011
# second Maple program:
b:= proc(n, i) option remember; `if`(i*(i+1)/2n or i>n-i+1, 0, b(n-i, i-1))))
    end:
a:= n-> b(n$2):
seq(a(n), n=0..80);  # Alois P. Heinz, May 20 2017

Sums[s_] := Sums[s] = With[{m = Max[s]}, If[m < 1, {}, Union @ Flatten @ Join[{m}, Table[{i, i + m}, {i, Sums[s ~Complement~ {m}]}]]]];
a[n_] := Module[{b}, b[i_, s_] := b[i, s] = Module[{si}, If[i == 1, If[Sums[s] == Range[n], 1, 0], si = s ~Union~ {i}; b[i-1, s] + If[Max[ Sums[si]] > n, 0, b[i - 1, si]]]]; b[n, {1}]];
Table[an = a[n]; Print["a(", n, ") = ", an]; an, {n, 1, 80}] (* Jean-François Alcover, Apr 12 2017, after Alois P. Heinz *)
Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&Union[Total/@Union[Subsets[#]]]==Range[0,n]&]],{n,30}] (* Gus Wiseman, May 31 2019 *)

/* As coefficients in g.f. */
{a(n)=local(A=[1]); for(i=1, n+1, A=concat(A,0); A[#A]=polcoeff(1 - sum(m=1,#A,A[m]*x^m/prod(k=1, m, 1+x^k +x*O(x^#A) )), #A) ); A[n+1]}
for(n=0, 50, print1(a(n),", ")) /* Paul D. Hanna, Mar 06 2012 */

F(n) = Sum_(i=L(n) .. U(n), F(n,i)), where F(n,i) = Sum_(j=L(n-i) .. min(U(n-i),i-1), F(n-i,j) ) and L(n), U(n) are defined in A188429 and A188430, respectively.
G.f.: 1 = Sum_{n>=0} a(n)*x^n / Product_{k=1..n+1} (1+x^k), with a(0)=1. - Paul D. Hanna, Mar 08 2012
a(n) ~ c * exp(Pi*sqrt(n/3)) / n^(3/4), where c = 0.03316508... - Vaclav Kotesovec, Oct 21 2019

More terms from Alois P. Heinz, Apr 03 2011

a(0)=1 prepended by Alois P. Heinz, May 20 2017

A326020 Number of complete subsets of {1..n}.

Original entry on oeis.org

1, 2, 3, 4, 6, 9, 15, 27, 50, 95, 185, 365, 724, 1441, 2873, 5735, 11458, 22902, 45789, 91561, 183102, 366180, 732331, 1464626, 2929209, 5858367, 11716674, 23433277, 46866473, 93732852, 187465596, 374931067, 749861989, 1499723808, 2999447418

Offset: 0

Gus Wiseman, Jun 04 2019

nonn

A set of positive integers summing to n is complete if every nonnegative integer up to n is the sum of some subset.

			The a(0) = 1 through a(6) = 15 subsets:
  {}  {}   {}     {}       {}         {}           {}
      {1}  {1}    {1}      {1}        {1}          {1}
           {1,2}  {1,2}    {1,2}      {1,2}        {1,2}
                  {1,2,3}  {1,2,3}    {1,2,3}      {1,2,3}
                           {1,2,4}    {1,2,4}      {1,2,4}
                           {1,2,3,4}  {1,2,3,4}    {1,2,3,4}
                                      {1,2,3,5}    {1,2,3,5}
                                      {1,2,4,5}    {1,2,3,6}
                                      {1,2,3,4,5}  {1,2,4,5}
                                                   {1,2,4,6}
                                                   {1,2,3,4,5}
                                                   {1,2,3,4,6}
                                                   {1,2,3,5,6}
                                                   {1,2,4,5,6}
                                                   {1,2,3,4,5,6}

Charlie Neder, Table of n, a(n) for n = 0..300
Andrzej Kukla and Piotr Miska, On practical sets and A-practical numbers, arXiv:2405.18225 [math.NT], 2024.

Cf. A002033, A103295, A108917, A126796, A188431, A276024.

Cf. A325684, A325781, A325790, A325791, A325986, A325988, A326016, A326021, A326022, A326036.

Table[Length[Select[Subsets[Range[n]],Union[Plus@@@Subsets[#]]==Range[0,Total[#]]&]],{n,0,10}]

a(17)-a(34) from Charlie Neder, Jun 05 2019

A330972 Sorted list containing the least number with each possible nonzero number of factorizations into factors > 1.

Original entry on oeis.org

1, 4, 8, 12, 16, 24, 36, 48, 60, 72, 96, 120, 128, 144, 180, 192, 216, 240, 256, 288, 360, 384, 420, 432, 480, 576, 720, 768, 840, 864, 900, 960, 1024, 1080, 1152, 1260, 1440, 1680, 1728, 1800, 1920, 2048, 2160, 2304, 2520, 2592, 2880, 3072, 3360, 3456, 3600

Offset: 1

Gus Wiseman, Jan 06 2020

nonn

This is the sorted list of positions of first appearances in A001055 of each element of the range (A045782).

			Factorizations of n for n = 4, 8, 12, 16, 24, 36, 48, 60:
  4    8      12     16       24       36       48         60
  2*2  2*4    2*6    2*8      3*8      4*9      6*8        2*30
       2*2*2  3*4    4*4      4*6      6*6      2*24       3*20
              2*2*3  2*2*4    2*12     2*18     3*16       4*15
                     2*2*2*2  2*2*6    3*12     4*12       5*12
                              2*3*4    2*2*9    2*3*8      6*10
                              2*2*2*3  2*3*6    2*4*6      2*5*6
                                       3*3*4    3*4*4      3*4*5
                                       2*2*3*3  2*2*12     2*2*15
                                                2*2*2*6    2*3*10
                                                2*2*3*4    2*2*3*5
                                                2*2*2*2*3

R. E. Canfield, P. Erdős and C. Pomerance, On a Problem of Oppenheim concerning "Factorisatio Numerorum", J. Number Theory 17 (1983), 1-28.

All terms belong to A025487
Includes all highly factorable numbers A033833.
Factorizations are A001055, with image A045782.
The least number with A045782(n) factorizations is A045783(n).
The least number with n factorizations is A330973(n).
The strict version is A330997.
Cf. A001222, A002033, A007716, A045778, A318284, A325238, A330935, A330936, A330976, A330977, A330989, A330991, A330992.

nn=1000;
facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
nds=Length/@Array[facs,nn];
Table[Position[nds,i][[1,1]],{i,First/@Gather[nds]}]

A045782 Number of factorizations of n for some n (image of A001055).

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 9, 11, 12, 15, 16, 19, 21, 22, 26, 29, 30, 31, 36, 38, 42, 45, 47, 52, 56, 57, 64, 66, 67, 74, 77, 92, 97, 98, 101, 105, 109, 118, 135, 137, 139, 141, 162, 165, 171, 176, 181, 189, 195, 198, 203, 212, 231, 249, 250, 254, 257, 267, 269, 272, 289

Offset: 1

David W. Wilson

nonn

Also the image of A318284. - Gus Wiseman, Jan 11 2020

Florian Luca, Anirban Mukhopadhyay and Kotyada Srinivas, On the Oppenheim's "factorisatio numerorum" function, arXiv:0807.0986 [math.NT], 2008.

Factorizations are A001055 with image this sequence and complement A330976.
Strict factorizations are A045778 with image A045779 and complement A330975.
The least number with exactly a(n) factorizations is A045783(n).
The least number with exactly n factorizations is A330973(n).
Cf. A002033, A007716, A033833, A318284, A325238, A330935, A330936, A330977, A330989, A330991, A330992, A330997.

terms = 61; m0 = 10^5; dm = 10^4;
f[1, ] = 1; f[n, k_] := f[n, k] = Sum[f[n/d, d], {d, Select[Divisors[n], 1 < # <= k &]}];
Clear[seq]; seq[m_] := seq[m] = Sort[Tally[Table[f[n, n], {n, 1, m}]][[All, 1]]][[1 ;; terms]]; seq[m = m0]; seq[m += dm]; While[Print[m]; seq[m] != seq[m - dm], m += dm];
seq[m] (* Jean-François Alcover, Oct 04 2018 *)

The Luca et al. paper shows that the number of terms with a(n) <= x is x^{ O( log log log x / log log x )}. - N. J. A. Sloane, Jun 12 2009

Name edited by Gus Wiseman, Jan 11 2020

A325862 Number of integer partitions of n such that every set of distinct parts has a different sum.

Original entry on oeis.org

1, 1, 2, 3, 5, 7, 10, 14, 19, 26, 34, 46, 58, 77, 93, 122, 146, 188, 217, 282, 327, 410, 470, 596, 673, 848, 947, 1178, 1325, 1629, 1798, 2213, 2444, 2962, 3247, 3935, 4292, 5149, 5579, 6674, 7247, 8590, 9221, 10964, 11804, 13870, 14843, 17480, 18675, 21866

Offset: 0

Gus Wiseman, May 31 2019

nonn

A knapsack partition (A108917, A299702) is an integer partition such that every submultiset has a different sum. The one non-knapsack partition counted under a(4) is (2,1,1).

			The a(1) = 1 through a(7) = 14 partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)
       (11)  (21)   (22)    (32)     (33)      (43)
             (111)  (31)    (41)     (42)      (52)
                    (211)   (221)    (51)      (61)
                    (1111)  (311)    (222)     (322)
                            (2111)   (411)     (331)
                            (11111)  (2211)    (421)
                                     (3111)    (511)
                                     (21111)   (2221)
                                     (111111)  (4111)
                                               (22111)
                                               (31111)
                                               (211111)
                                               (1111111)
The three non-knapsack partitions counted under a(6) are:
  (2,2,1,1)
  (3,1,1,1)
  (2,1,1,1,1)

Dominates A108917.

Cf. A002033, A034444, A196723, A275972, A276024, A299702, A325592, A325856, A325863, A325864, A325865, A325877.

Table[Length[Select[IntegerPartitions[n],UnsameQ@@Plus@@@Subsets[Union[#]]&]],{n,0,20}]

A076274 2p-1 where p is 1 or a prime.

Original entry on oeis.org

1, 3, 5, 9, 13, 21, 25, 33, 37, 45, 57, 61, 73, 81, 85, 93, 105, 117, 121, 133, 141, 145, 157, 165, 177, 193, 201, 205, 213, 217, 225, 253, 261, 273, 277, 297, 301, 313, 325, 333, 345, 357, 361, 381, 385, 393, 397, 421, 445, 453, 457, 465, 477, 481, 501, 513, 525, 537, 541

Offset: 1

Reinhard Zumkeller, Oct 04 2002

Vincenzo Librandi, Table of n, a(n) for n = 1..1000

Positions of squares in A076271.
Subsequence of A186621 and of A005408.
Cf. A000010, A000720, A002033, A008578.

[2*NthPrime(n-1)-1: n in [1..60]]; // Vincenzo Librandi, Jul 31 2017

Join[{1},2*Prime[Range[100]]-1] (* Harvey P. Dale, Nov 29 2015 *)

a(n)=if(n>1,2*prime(n-1)-1,1) \\ Charles R Greathouse IV, Apr 29 2015

A076271(a(n)) = A008578(n)^2 for all n.
a(n+1) - a(n) = 2*A075526(n) for all n.
G.f.: 1 + Sum_{ j>=1 } Sum_{ k>=1 } x^(pi(k/2)+j), where pi(k) is the prime counting function. - Benedict W. J. Irwin, Jun 13 2016
a(n) = 2*A008578(n) - 1. - Alois P. Heinz, Nov 28 2023

New definition from Thomas Baruchel, Oct 21 2003

Offset changed from Vincenzo Librandi, Jul 31 2017

cp's OEIS Frontend

A325781 Heinz numbers of complete integer partitions.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A065036 Product of the cube of a prime (A030078) and a different prime.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Python

Formula

A089723 a(1)=1; for n>1, a(n) gives number of ways to write n as n = x^y, 2 <= x, 1 <= y.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

A034729 a(n) = Sum_{ k, k|n } 2^(k-1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

PARI

Python

SageMath

Formula

A188431 The number of n-full sets, F(n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

Extensions

A326020 Number of complete subsets of {1..n}.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs