A316855 -id:A316855 - cp's OEIS frontend

A051908 Number of ways to express 1 as the sum of unit fractions such that the sum of the denominators is n.

1, 0, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 3, 0, 1, 1, 1, 1, 2, 3, 2, 2, 1, 2, 2, 2, 4, 5, 5, 2, 4, 5, 5, 9, 4, 4, 6, 4, 4, 7, 8, 4, 10, 9, 9, 11, 8, 13, 13, 15, 16, 21, 18, 16, 22, 19, 18, 30, 24, 19, 26, 28, 26, 29, 35, 29, 44, 28, 47, 48

Offset: 1

Jud McCranie, Dec 16 1999

nonn

Also the number of partitions of n whose reciprocal sums to 1; "exact partitions". - Robert G. Wilson v, Sep 30 2009

			1 = 1/2 + 1/2, the sum of denominators is 4, and this is the only expression of 1 as unit fractions with denominator sum 4, so a(4)=1.
The a(22) = 3 partitions whose reciprocal sum is 1 are (12,4,3,3), (10,5,5,2), (8,8,4,2). - _Gus Wiseman_, Jul 16 2018

Derrick Niederman, "Number Freak, From 1 to 200 The Hidden Language of Numbers Revealed", a Perigee Book, Penguin Group, NY, 2009, pp. 82-83. [From Robert G. Wilson v, Sep 30 2009]

David A. Corneth, Table of n, a(n) for n = 1..200 (terms a(1)-a(86) from Jud McCranie, a(87)-a(88) from Robert G. Wilson v, a(89)-a(100) from Seiichi Manyama)
David A. Corneth, Tuples up to n = 170
Gus Wiseman, Sequences counting and ranking integer partitions by their reciprocal sums
Index entries for sequences related to Egyptian fractions

A028229 lists n such that a(n)=0.

Cf. A002966, A058360, A270599, A316854, A316855, A316888-A316904.

(* first do *) << "Combinatorica`"; (* then *) f[n_] := Block[{c = i = 0, k = PartitionsP@n, p = {n}}, While[i < k, If[1 == Plus @@ (1/p), c++ ]; i++; p = NextPartition@p]; c]; Array[f, 88] (* Robert G. Wilson v, Sep 30 2009 *)
Table[Length[Select[IntegerPartitions[n],Sum[1/m,{m,#}]==1&]],{n,30}] (* Gus Wiseman, Jul 16 2018 *)

def partition(n, min, max)
  return [[]] if n == 0
  [max, n].min.downto(min).flat_map{|i| partition(n - i, min, i).map{|rest| [i, *rest]}}
end
def A051908(n)
  ary = [1]
  (2..n).each{|m|
    cnt = 0
    partition(m, 2, m).each{|ary|
      cnt += 1 if ary.inject(0){|s, i| s + 1 / i.to_r} == 1
    }
    ary << cnt
  }
  ary
end
p A051908(100) # Seiichi Manyama, May 31 2016

a(n) > 0 for n > 23.

A316854 Number of integer partitions of n whose reciprocal sum is the reciprocal of an integer.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 2, 3, 2, 2, 2, 1, 1, 1, 4, 2, 4, 1, 5, 1, 5, 1, 3, 4, 2, 5, 6, 5, 5, 4, 5, 5, 4, 8, 10, 9, 7, 5, 9, 10, 6, 12, 10, 8, 7, 6, 9, 13, 15, 8, 19, 13, 19, 19, 19, 18, 22, 26, 28, 28, 29, 22, 33, 29, 28, 38, 34, 26, 40, 32, 43, 39, 51, 38, 62, 46

Offset: 1

Gus Wiseman, Jul 14 2018

nonn

The reciprocal sum of (y_1, ..., y_k) is 1/y_1 + ... + 1/y_k.

			The a(36) = 10 partitions:
  (36),
  (30,6), (24,12), (18,18),
  (12,12,12),
  (12,12,6,6),
  (15,10,4,4,3), (12,12,6,3,3), (12,8,8,6,2),
  (6,6,6,6,6,6).

Giovanni Resta, Table of n, a(n) for n = 1..200 (first 100 terms from Robert G. Wilson v)

Cf. A000041, A051908, A058360, A316855, A316856, A316857.

Table[Length[Select[IntegerPartitions[n],IntegerQ[1/Sum[1/m,{m,#}]]&]],{n,30}]
ric[n_, p_, s_] := If[n == 0, If[IntegerQ[1/s], c++], Do[If[s + 1/i <= 1, ric[n - i, Append[p, i], s + 1/i]], {i, Min[p[[-1]], n], 1, -1}]]; a[n_] := (c = 0; Do[ric[n - j, {j}, 1/j], {j, n}]; c); Array[a, 80] (* after Giovanni Resta in A316898, Robert G. Wilson v, Jul 23 2018 *)

a(n)={my(s=0); forpart(p=n, if(frac(1/sum(i=1, #p, 1/p[i]))==0, s++)); s} \\ Andrew Howroyd, Jul 15 2018

a(51)-a(77) from Giovanni Resta, Jul 15 2018

A316857 Heinz numbers of integer partitions whose reciprocal sum is the reciprocal of an integer.

Original entry on oeis.org

2, 3, 5, 7, 9, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 49, 53, 59, 61, 65, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 125, 127, 131, 137, 139, 147, 149, 151, 157, 163, 167, 169, 173, 179, 181, 191, 193, 195, 197, 199, 211, 223, 227, 229, 233, 239

Offset: 1

Gus Wiseman, Jul 14 2018

nonn

The reciprocal sum of (y_1, ..., y_k) is 1/y_1 + ... + 1/y_k.

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

Cf. A000041, A051908, A056239, A058360, A072411, A296150, A316854, A316855, A316856.

Select[Range[2,100],IntegerQ[1/Sum[m[[2]]/PrimePi[m[[1]]],{m,FactorInteger[#]}]]&]

A316856 Heinz numbers of integer partitions whose reciprocal sum is an integer.

Original entry on oeis.org

1, 2, 4, 8, 9, 16, 18, 32, 36, 64, 72, 81, 125, 128, 144, 147, 162, 195, 250, 256, 288, 294, 324, 390, 500, 512, 576, 588, 648, 729, 780, 1000, 1024, 1125, 1152, 1176, 1296, 1323, 1458, 1560, 1755, 2000, 2048, 2250, 2304, 2352, 2401, 2592, 2646, 2916, 3120

Offset: 1

Gus Wiseman, Jul 14 2018

nonn

The reciprocal sum of (y_1, ..., y_k) is 1/y_1 + ... + 1/y_k.

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

			195 is the Heinz number of (6,3,2), which has reciprocal sum 1/6 + 1/3 + 1/2 = 1, which is an integer, so 195 belongs to the sequence.
The sequence of all integer partitions whose reciprocal sum is an integer begins: (), (1), (11), (111), (22), (1111), (221), (11111), (2211), (111111), (22111), (2222).

Cf. A000041, A051908, A056239, A058360, A072411, A296150, A316854, A316855, A316857.

Select[Range[1000],IntegerQ[Sum[m[[2]]/PrimePi[m[[1]]],{m,If[#==1,{},FactorInteger[#]]}]]&]

A316888 Heinz numbers of aperiodic integer partitions into relatively prime parts whose reciprocal sum is 1.

Original entry on oeis.org

2, 195, 3185, 6475, 10527, 16401, 20445, 20535, 21045, 25365, 46155, 164255, 171941, 218855, 228085, 267883, 312785, 333925, 333935, 335405, 343735, 355355, 414295, 442975, 474513, 527425, 549575, 607475, 633777, 691041, 711321, 722425, 753865, 804837, 822783

Offset: 1

Gus Wiseman, Jul 16 2018

nonn

The reciprocal sum of (y_1, ..., y_k) is 1/y_1 + ... + 1/y_k.
The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
A partition is aperiodic if its multiplicities are relatively prime.
Does not contain 29888089, which belongs to A316890 and is the Heinz number of a periodic partition.

			The partition (6,4,4,3) with Heinz number 3185 is aperiodic, has relatively prime parts, and 1/6 + 1/4 + 1/4 + 1/3 = 1, so 3185 belongs to the sequence.
The sequence of partitions whose Heinz numbers belong to the sequence begins: (1), (6,3,2), (6,4,4,3), (12,4,3,3), (10,5,5,2), (20,5,4,2), (15,10,3,2), (12,12,3,2), (18,9,3,2), (24,8,3,2), (42,7,3,2).

Gus Wiseman, Sequences counting and ranking integer partitions by their reciprocal sums

Cf. A000837, A002966, A007916, A051908, A100953, A289509, A296150, A316855, A316856, A316857, A316888-A316904.

Select[Range[2,100000],And[GCD@@FactorInteger[#][[All,2]]==1,GCD@@PrimePi/@FactorInteger[#][[All,1]]==1,Sum[m[[2]]/PrimePi[m[[1]]],{m,FactorInteger[#]}]==1]&]

A316904 Heinz numbers of aperiodic integer partitions into relatively prime parts whose reciprocal sum is an integer.

Original entry on oeis.org

2, 18, 72, 162, 195, 250, 288, 294, 390, 500, 588, 648, 780, 1125, 1152, 1176, 1458, 1560, 1755, 2000, 2250, 2352, 2592, 2646, 3120, 3185, 3510, 4000, 4500, 4608, 4704, 4802, 5292, 6240, 6370, 6475, 7020, 8450, 9000, 9408, 10125, 10368, 10527, 10584, 12480

Offset: 1

Gus Wiseman, Jul 16 2018

nonn

The reciprocal sum of (y_1, ..., y_k) is 1/y_1 + ... + 1/y_k.
The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
A partition is aperiodic if its multiplicities are relatively prime.

			The sequence of partitions whose Heinz numbers belong to this sequence begins: (1), (221), (22111), (22221), (632), (3331), (2211111), (4421), (6321), (33311), (44211), (2222111).

Gus Wiseman, Sequences counting and ranking integer partitions by their reciprocal sums

Cf. A000837, A002966, A051908, A058360, A100953, A296150, A316854, A316855, A316856, A316857, A316888-A316904.

Select[Range[2,20000],And[GCD@@FactorInteger[#][[All,2]]==1,GCD@@PrimePi/@FactorInteger[#][[All,1]]==1,IntegerQ[Sum[m[[2]]/PrimePi[m[[1]]],{m,FactorInteger[#]}]]]&]

A325618 Numbers m such that there exists an integer partition of m whose reciprocal factorial sum is 1.

Original entry on oeis.org

1, 4, 11, 18, 24, 31, 37, 44, 45, 50, 52, 57, 58, 65, 66, 70, 71, 73, 76, 78, 79, 83, 86, 87, 89, 91, 92, 94, 96, 97, 99, 100, 102, 104, 107, 108, 109, 110, 112, 113, 114, 115, 117, 118, 119, 120, 121, 122, 123, 125, 126, 127, 128, 130, 131

Offset: 1

Gus Wiseman, May 13 2019

nonn

The reciprocal factorial sum of an integer partition (y_1,...,y_k) is 1/y_1! + ... + 1/y_k!.

Conjecture: 137 is the greatest integer not in this sequence. - Charlie Neder, May 14 2019

			The sequence of terms together with an integer partition of each whose reciprocal factorial sum is 1 begins:
   1: (1)
   4: (2,2)
  11: (3,3,3,2)
  18: (3,3,3,3,3,3)
  24: (4,4,4,4,3,3,2)
  31: (4,4,4,4,3,3,3,3,3)
  37: (4,4,4,4,4,4,4,4,3,2)
  44: (4,4,4,4,4,4,4,4,3,3,3,3)
  45: (5,5,5,5,5,4,4,4,3,3,2)
  50: (4,4,4,4,4,4,4,4,4,4,4,4,2)

Charlie Neder, Table of n, a(n) for n = 1..1000

Factorial numbers: A000142, A002982, A007489, A011371, A022559, A064986, A115627, A284605, A325616.

Reciprocal factorial sum: A002966, A051908, A058360, A316854, A316855, A325619, A325620, A325621, A325622, A325623, A325624.

a(11)-a(55) from Charlie Neder, May 14 2019

A325619 Heinz numbers of integer partitions whose reciprocal factorial sum is 1.

Original entry on oeis.org

2, 9, 375, 15625

Offset: 1

Gus Wiseman, May 13 2019

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

The reciprocal factorial sum of an integer partition (y_1,...,y_k) is 1/y_1! + ... + 1/y_k!.

			The sequence of terms together with their prime indices begins:
      1: {}
      2: {1}
      9: {2,2}
    375: {2,3,3,3}
  15625: {3,3,3,3,3,3}

Factorial numbers: A000142, A007489, A022559, A064986, A108731, A115944, A284605, A325508, A325616.

Reciprocal factorial sum: A002966, A051908, A316855, A325618, A325624.

Select[Range[100000],Total[Cases[FactorInteger[#],{p_,k_}:>k/PrimePi[p]!]]==1&]

Contains prime(n)^(n!) for all n > 0, including 191581231380566414401 for n = 4.

A318573 Numerator of the reciprocal sum of the integer partition with Heinz number n.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Aug 29 2018

The reciprocal sum of (y_1, ..., y_k) is 1/y_1 + ... + 1/y_k. The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Antti Karttunen, Table of n, a(n) for n = 1..10000
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537
Gus Wiseman, Sequences counting and ranking integer partitions by their reciprocal sums
Index entries for sequences computed from indices in prime factorization
Index entries for sequences related to Heinz numbers

Positions of 1's are A316857.

Cf. A051908, A056239, A058360, A112798, A289506, A289507, A296150, A316854, A316855, A316856, A318574, A325704.

Table[Sum[pr[[2]]/PrimePi[pr[[1]]],{pr,If[n==1,{},FactorInteger[n]]}],{n,100}]//Numerator

A318573(n) = { my(f=factor(n)); numerator(sum(i=1,#f~,f[i, 2]/primepi(f[i, 1]))); }; \\ Antti Karttunen, Nov 17 2019

If n = Product prime(x_i)^y_i is the prime factorization of n, then a(n) is the numerator of Sum y_i/x_i.

More terms from Antti Karttunen, Nov 17 2019

A318584 Number of integer partitions of n whose sum of reciprocals squared is 1.

Original entry on oeis.org

0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 1, 0, 2, 2, 0, 0, 0, 2, 0, 2, 1, 2, 2, 2, 1, 1, 2, 3, 0, 1, 1, 6, 2, 3, 2, 6, 2, 2, 3, 2, 6, 7, 2, 4, 3, 9, 4, 7, 5, 8, 8, 7, 9, 9, 11, 12, 7, 9, 11, 17, 9, 13, 12, 17, 16, 13, 15, 20, 26, 27, 18, 23

Offset: 0

Gus Wiseman, Aug 29 2018

nonn

The a(16) = 1 integer partition:
  (6,3,3,2,2,2)
The a(48) = 2 integer partitions:
  (18,9,9,3,3,2,2,2)
  (6,6,6,6,3,3,3,3,3,3,3,3)
The a(56) = 3 integer partitions:
  (12,6,6,4,4,4,4,4,4,4,2,2)
  (10,6,5,5,5,5,5,5,3,3,2,2)
  (6,6,4,4,4,4,4,4,4,4,3,3,3,3)
The a(60) = 6 integer partitions:
  (12,12,12,12,3,3,2,2,2)
  (8,8,8,8,6,4,4,4,3,3,2,2)
  (6,6,6,6,6,6,6,6,6,2,2,2)
  (12,12,12,4,3,3,3,3,3,3,2)
  (10,5,5,5,5,5,5,4,4,4,4,2,2)
  (6,4,4,4,4,4,4,4,4,4,4,4,4,3,3)

Cf. A000041, A051908, A058360, A289506, A289507, A316854, A316855.

Cf. A318585, A318586, A318587, A318588, A318589.

Table[Length[Select[IntegerPartitions[n],Total[#^(-2)]==1&]],{n,30}]

a(61)-a(100) from Alois P. Heinz, Aug 30 2018

cp's OEIS Frontend

A051908 Number of ways to express 1 as the sum of unit fractions such that the sum of the denominators is n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Ruby

Formula

A316854 Number of integer partitions of n whose reciprocal sum is the reciprocal of an integer.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A316857 Heinz numbers of integer partitions whose reciprocal sum is the reciprocal of an integer.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

A316856 Heinz numbers of integer partitions whose reciprocal sum is an integer.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A316888 Heinz numbers of aperiodic integer partitions into relatively prime parts whose reciprocal sum is 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A316904 Heinz numbers of aperiodic integer partitions into relatively prime parts whose reciprocal sum is an integer.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A325618 Numbers m such that there exists an integer partition of m whose reciprocal factorial sum is 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Extensions

A325619 Heinz numbers of integer partitions whose reciprocal factorial sum is 1.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs