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-6 of 6 results.

A054973 Number of numbers whose divisors sum to n.

Original entry on oeis.org

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

Views

Author

Henry Bottomley, May 16 2000

Keywords

Comments

a(n) = frequency of values n in A000203(m), where A000203(m) = sum of divisors of m. a(n) >= 1 for such n that A175192(n) = 1, a(n) >= 1 if A000203(m) = n for any m. a(n) = 0 for such n that A175192(n) = 0, a(n) = 0 if A000203(m) = n has no solution. - Jaroslav Krizek, Mar 01 2010
First occurrence of k: 2, 1, 12, 24, 96, 72, ..., = A007368. - Robert G. Wilson v, May 14 2014
a(n) is also the number of positive terms in the n-th row of triangle A299762. - Omar E. Pol, Mar 14 2018
Also the number of integer partitions of n whose parts form the set of divisors of some number (necessarily the greatest part). The Heinz numbers of these partitions are given by A371283. For example, the a(24) = 3 partitions are: (23,1), (15,5,3,1), (14,7,2,1). - Gus Wiseman, Mar 22 2024

Examples

			a(12) = 2 since 11 has factors 1 and 11 with 1 + 11 = 12 and 6 has factors 1, 2, 3 and 6 with 1 + 2 + 3 + 6 = 12.
		

Crossrefs

Cf. A000203 (sum-of-divisors function).
For partial sums see A074753.
The non-strict version is A371284, ranks A371288.
These partitions have ranks A371283, unsorted version A275700.
A000005 counts divisors, row-lengths of A027750.
A000041 counts integer partitions, strict A000009.

Programs

  • Mathematica
    nn = 105; t = Table[0, {nn}]; k = 1; While[k < 6 nn^(3/2)/Pi^2, d = DivisorSigma[1, k]; If[d < nn + 1, t[[d]]++]; k++]; t (* Robert G. Wilson v, May 14 2014 *)
    Table[Length[Select[IntegerPartitions[n],#==Reverse[Divisors[Max@@#]]&]],{n,30}] (* Gus Wiseman, Mar 22 2024 *)
  • PARI
    a(n)=v = vector(0); for (i = 1, n, if (sigma(i) == n, v = concat(v, i));); #v; \\ Michel Marcus, Oct 22 2013
    
  • PARI
    a(n)=sum(k=1,n,sigma(k)==n) \\ Charles R Greathouse IV, Nov 12 2013
    
  • PARI
    first(n)=my(v=vector(n),t); for(k=1,n, t=sigma(n); if(t<=n, v[t]++)); v \\ Charles R Greathouse IV, Mar 08 2017
    
  • PARI
    A054973(n)=#invsigma(n) \\ See Alekseyev link for invsigma(). - M. F. Hasler, Nov 21 2019

Formula

Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = A308039. - Amiram Eldar, Dec 23 2024

Extensions

Incorrect comment deleted by M. F. Hasler, Nov 21 2019

A275700 a(n) = Product_{d|n} prime(d).

Original entry on oeis.org

2, 6, 10, 42, 22, 390, 34, 798, 230, 1914, 62, 101010, 82, 4386, 5170, 42294, 118, 547170, 134, 951258, 12410, 14694, 166, 170807910, 2134, 24846, 23690, 3285114, 218, 660741510, 254, 5540514, 42470, 49206, 55726, 21399271530, 314, 65526, 68470, 3126785046, 358
Offset: 1

Views

Author

Jaroslav Krizek, Aug 05 2016

Keywords

Comments

a(n) mod n = 0 for n: 1, 2, 6, 30, 78, 330, 390, 870, 1410, 3198, ...

Examples

			a(4) = 42 because the divisors of 4 are: 1, 2 and 4; and prime(1) * prime(2) * prime(4) = 2 * 3 * 7 = 42.
		

Crossrefs

Cf. A007445 (Sum_{d|n} prime(d)).
A version for binary indices is A034729.
Partitions of this type are counted by A054973, strict case of A371284.
The sorted version is A371283, squarefree case of A371288.
These numbers have products A371286, unsorted version A371285.
A000005 counts divisors, row-lengths of A027750.
A027746 lists prime factors, indices A112798, length A001222.

Programs

  • Magma
    [(&*[NthPrime(d): d in Divisors(n)]): n in [1..100]]
    
  • Mathematica
    Table[Times@@(Prime[#]&/@Divisors[n]),{n,50}] (* Harvey P. Dale, Jun 16 2017 *)
  • PARI
    a(n) = my(d=divisors(n)); prod(i=1, #d, prime(d[i])) \\ Felix Fröhlich, Aug 05 2016
    
  • Perl
    use ntheory ":all"; sub a275700 { vecprod(map { nth_prime($) } divisors($[0])); } # Dana Jacobsen, Aug 09 2016

A371288 Numbers whose distinct prime indices form the set of divisors of some positive integer.

Original entry on oeis.org

2, 4, 6, 8, 10, 12, 16, 18, 20, 22, 24, 32, 34, 36, 40, 42, 44, 48, 50, 54, 62, 64, 68, 72, 80, 82, 84, 88, 96, 100, 108, 118, 124, 126, 128, 134, 136, 144, 160, 162, 164, 166, 168, 176, 192, 200, 216, 218, 230, 236, 242, 248, 250, 252, 254, 256, 268, 272, 288
Offset: 1

Views

Author

Gus Wiseman, Mar 22 2024

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Examples

			The prime indices of 694782 are {1,2,2,5,5,5,10} with distinct elements {1,2,5,10}, which form the set of divisors of 10, so 694782 is in the sequence.
The terms together with their prime indices begin:
    2: {1}
    4: {1,1}
    6: {1,2}
    8: {1,1,1}
   10: {1,3}
   12: {1,1,2}
   16: {1,1,1,1}
   18: {1,2,2}
   20: {1,1,3}
   22: {1,5}
   24: {1,1,1,2}
   32: {1,1,1,1,1}
   34: {1,7}
   36: {1,1,2,2}
   40: {1,1,1,3}
   42: {1,2,4}
   44: {1,1,5}
   48: {1,1,1,1,2}
		

Crossrefs

The squarefree case is A371283, unsorted version A275700.
Partitions of this type are counted by A371284, strict A054973.
Products of squarefree terms are A371286, unsorted version A371285.
A000005 counts divisors.
A001221 counts distinct prime factors.
A027746 lists prime factors, indices A112798, length A001222.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],Union[prix[#]]==Divisors[Max@@prix[#]]&]

A371283 Heinz numbers of sets of divisors of positive integers. Numbers whose prime indices form the set of divisors of some positive integer.

Original entry on oeis.org

2, 6, 10, 22, 34, 42, 62, 82, 118, 134, 166, 218, 230, 254, 314, 358, 382, 390, 422, 482, 554, 566, 662, 706, 734, 798, 802, 862, 922, 1018, 1094, 1126, 1174, 1198, 1234, 1418, 1478, 1546, 1594, 1718, 1754, 1838, 1914, 1934, 1982, 2062, 2126, 2134, 2174, 2306
Offset: 1

Views

Author

Gus Wiseman, Mar 21 2024

Keywords

Comments

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Examples

			The terms together with their prime indices begin:
     2: {1}
     6: {1,2}
    10: {1,3}
    22: {1,5}
    34: {1,7}
    42: {1,2,4}
    62: {1,11}
    82: {1,13}
   118: {1,17}
   134: {1,19}
   166: {1,23}
   218: {1,29}
   230: {1,3,9}
   254: {1,31}
   314: {1,37}
   358: {1,41}
   382: {1,43}
   390: {1,2,3,6}
		

Crossrefs

Partitions of this type are counted by A054973.
The unsorted version is A275700.
These numbers have products A371286, unsorted version A371285.
Squarefree case of A371288, counted by A371284.
A000005 counts divisors.
A001221 counts distinct prime factors.
A027746 lists prime factors, A112798 indices, length A001222.
A355731 counts choices of a divisor of each prime index, firsts A355732.
A355741 counts choices of a prime factor of each prime index.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n], {p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[2,100],SameQ[prix[#],Divisors[Last[prix[#]]]]&]

A371285 Heinz number of the multiset union of the divisor sets of each prime index of n.

Original entry on oeis.org

1, 2, 6, 4, 10, 12, 42, 8, 36, 20, 22, 24, 390, 84, 60, 16, 34, 72, 798, 40, 252, 44, 230, 48, 100, 780, 216, 168, 1914, 120, 62, 32, 132, 68, 420, 144, 101010, 1596, 2340, 80, 82, 504, 4386, 88, 360, 460, 5170, 96, 1764, 200, 204, 1560, 42294, 432, 220, 336
Offset: 1

Views

Author

Gus Wiseman, Mar 21 2024

Keywords

Comments

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Examples

			The prime indices of 105 are {2,3,4}, with divisor sets {{1,2},{1,3},{1,2,4}}, with multiset union {1,1,1,2,2,3,4}, with Heinz number 2520, so a(105) = 2520.
The terms together with their prime indices begin:
          1: {}
          2: {1}
          6: {1,2}
          4: {1,1}
         10: {1,3}
         12: {1,1,2}
         42: {1,2,4}
          8: {1,1,1}
         36: {1,1,2,2}
         20: {1,1,3}
         22: {1,5}
         24: {1,1,1,2}
        390: {1,2,3,6}
         84: {1,1,2,4}
         60: {1,1,2,3}
         16: {1,1,1,1}
         34: {1,7}
         72: {1,1,1,2,2}
		

Crossrefs

Product of A275700 applied to each prime index.
The squarefree case is also A275700.
The sorted version is A371286.
A000005 counts divisors.
A001221 counts distinct prime factors.
A027746 lists prime factors, A112798 indices, length A001222.
A355731 counts choices of a divisor of each prime index, firsts A355732.
A355741 counts choices of a prime factor of each prime index.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n], {p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Times@@Prime/@Join@@Divisors/@prix[n],{n,100}]

Formula

If n = prime(x_1)*...*prime(x_k) then a(n) = A275700(x_1)*...*A275700(x_k).

A371286 Products of elements of A275700 (Heinz numbers of divisor sets). Numbers with a (necessarily unique) factorization into elements of A275700.

Original entry on oeis.org

1, 2, 4, 6, 8, 10, 12, 16, 20, 22, 24, 32, 34, 36, 40, 42, 44, 48, 60, 62, 64, 68, 72, 80, 82, 84, 88, 96, 100, 118, 120, 124, 128, 132, 134, 136, 144, 160, 164, 166, 168, 176, 192, 200, 204, 216, 218, 220, 230, 236, 240, 248, 252, 254, 256, 264, 268, 272, 288
Offset: 1

Views

Author

Gus Wiseman, Mar 22 2024

Keywords

Comments

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Examples

			The terms together with their prime factorizations and unique factorizations into terms of A275700 begin:
   1 =             = ()
   2 = 2           = (2)
   4 = 2*2         = (2*2)
   6 = 2*3         = (6)
   8 = 2*2*2       = (2*2*2)
  10 = 2*5         = (10)
  12 = 2*2*3       = (2*6)
  16 = 2*2*2*2     = (2*2*2*2)
  20 = 2*2*5       = (2*10)
  22 = 2*11        = (22)
  24 = 2*2*2*3     = (2*2*6)
  32 = 2*2*2*2*2   = (2*2*2*2*2)
  34 = 2*17        = (34)
  36 = 2*2*3*3     = (6*6)
  40 = 2*2*2*5     = (2*2*10)
  42 = 2*3*7       = (42)
  44 = 2*2*11      = (2*22)
  48 = 2*2*2*2*3   = (2*2*2*6)
  60 = 2*2*3*5     = (6*10)
  62 = 2*31        = (62)
  64 = 2*2*2*2*2*2 = (2*2*2*2*2*2)
  68 = 2*2*17      = (2*34)
  72 = 2*2*2*3*3   = (2*6*6)
  80 = 2*2*2*2*5   = (2*2*2*10)
  82 = 2*41        = (82)
  84 = 2*2*3*7     = (2*42)
  88 = 2*2*2*11    = (2*2*22)
  96 = 2*2*2*2*2*3 = (2*2*2*2*6)
		

Crossrefs

Products of elements of A275700.
The squarefree case is A371283.
The unsorted version is A371285.
A000005 counts divisors.
A001221 counts distinct prime factors.
A027746 lists prime factors, A112798 indices, length A001222.
A355731 counts choices of a divisor of each prime index, firsts A355732.
A355741 counts choices of a prime factor of each prime index.

Programs

  • Mathematica
    nn=100;
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n], {p_,k_}:>Table[PrimePi[p],{k}]]]];
    facs[n_]:=If[n<=1, {{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]], {d,Rest[Divisors[n]]}]];
    s=Table[Times@@Prime/@Divisors[n],{n,nn}];
    m=Max@@Table[Select[Range[2,k],prix[#] == Divisors[Last[prix[#]]]&],{k,nn}];
    Join@@Position[Table[Length[Select[facs[n], SubsetQ[s,Union[#]]&]],{n,m}],1]
Showing 1-6 of 6 results.