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

A109298 Primal codes of finite idempotent functions on positive integers.

Original entry on oeis.org

1, 2, 9, 18, 125, 250, 1125, 2250, 2401, 4802, 21609, 43218, 161051, 300125, 322102, 600250, 1449459, 2701125, 2898918, 4826809, 5402250, 9653618, 20131375, 40262750, 43441281, 86882562, 181182375, 362364750, 386683451, 410338673, 603351125, 773366902, 820677346
Offset: 1

Views

Author

Jon Awbrey, Jul 06 2005

Keywords

Comments

Finite idempotent functions are identity maps on finite subsets, counting the empty function as the idempotent on the empty set.
From Gus Wiseman, Mar 09 2019: (Start)
Also numbers whose ordered prime signature is equal to the distinct prime indices in increasing order. A prime index of n is a number m such that prime(m) divides n. The ordered prime signature (A124010) is the sequence of multiplicities (or exponents) in a number's prime factorization, taken in order of the prime base. The case where the prime indices are taken in decreasing order is A324571.
Also numbers divisible by prime(k) exactly k times for each prime index k. These are a kind of self-describing numbers (cf. A001462, A304679).
Also Heinz numbers of integer partitions where the multiplicity of m is m for all m in the support (counted by A033461). The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
Also products of distinct elements of A062457. For example, 43218 = prime(1)^1 * prime(2)^2 * prime(4)^4.
(End)

Examples

			Writing (prime(i))^j as i:j, we have the following table of examples:
Primal Codes of Finite Idempotent Functions on Positive Integers
` ` ` 1 = { }
` ` ` 2 = 1:1
` ` ` 9 = ` ` 2:2
` ` `18 = 1:1 2:2
` ` 125 = ` ` ` ` 3:3
` ` 250 = 1:1 ` ` 3:3
` `1125 = ` ` 2:2 3:3
` `2250 = 1:1 2:2 3:3
` `2401 = ` ` ` ` ` ` 4:4
` `4802 = 1:1 ` ` ` ` 4:4
` 21609 = ` ` 2:2 ` ` 4:4
` 43218 = 1:1 2:2 ` ` 4:4
`161051 = ` ` ` ` ` ` ` ` 5:5
`300125 = ` ` ` ` 3:3 4:4
`322102 = 1:1 ` ` ` ` ` ` 5:5
`600250 = 1:1 ` ` 3:3 4:4
From _Gus Wiseman_, Mar 09 2019: (Start)
The sequence of terms together with their prime indices begins as follows. For example, we have 18: {1,2,2} because 18 = prime(1) * prime(2) * prime(2) has prime signature {1,2} and the distinct prime indices are also {1,2}.
       1: {}
       2: {1}
       9: {2,2}
      18: {1,2,2}
     125: {3,3,3}
     250: {1,3,3,3}
    1125: {2,2,3,3,3}
    2250: {1,2,2,3,3,3}
    2401: {4,4,4,4}
    4802: {1,4,4,4,4}
   21609: {2,2,4,4,4,4}
   43218: {1,2,2,4,4,4,4}
  161051: {5,5,5,5,5}
  300125: {3,3,3,4,4,4,4}
  322102: {1,5,5,5,5,5}
  600250: {1,3,3,3,4,4,4,4}
(End)
		

Crossrefs

Programs

  • Mathematica
    Select[Range[10000],And@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>PrimePi[p]==k]&]
  • PARI
    is(n) = my(f = factor(n)); for(i = 1, #f~, if(prime(f[i, 2]) != f[i, 1], return(0))); 1 \\ David A. Corneth, Mar 09 2019

Formula

Sum_{n>=1} 1/a(n) = Product_{n>=1} (1 + 1/prime(n)^n) = 1.6807104966... - Amiram Eldar, Jan 03 2021

Extensions

Offset set to 1, missing terms inserted and more terms added by Alois P. Heinz, Mar 08 2019

A304678 Numbers with weakly increasing prime multiplicities.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 18, 19, 21, 22, 23, 25, 26, 27, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 41, 42, 43, 46, 47, 49, 50, 51, 53, 54, 55, 57, 58, 59, 61, 62, 64, 65, 66, 67, 69, 70, 71, 73, 74, 75, 77, 78, 79, 81, 82, 83
Offset: 1

Views

Author

Gus Wiseman, May 16 2018

Keywords

Comments

Complement of A112769.

Examples

			12 = 2*2*3 has prime multiplicities (2,1) so is not in the sequence.
36 = 2*2*3*3 has prime multiplicities (2,2) so is in the sequence.
150 = 2*3*5*5 has prime multiplicities (1,1,2) so is in the sequence.
		

Crossrefs

Programs

  • Maple
    q:= n-> (l-> (t-> andmap(i-> l[i, 2]<=l[i+1, 2],
            [$1..t-1]))(nops(l)))(sort(ifactors(n)[2])):
    select(q, [$1..120])[];  # Alois P. Heinz, Nov 11 2019
  • Mathematica
    Select[Range[200],OrderedQ[FactorInteger[#][[All,2]]]&]
    Select[Range[90],Min[Differences[FactorInteger[#][[;;,2]]]]>=0&] (* Harvey P. Dale, Jan 28 2024 *)
  • PARI
    isok(n) = my(vm = factor(n)[,2]); vm == vecsort(vm); \\ Michel Marcus, May 17 2018

A324572 Number of integer partitions of n whose multiplicities (where if x < y the multiplicity of x is counted prior to the multiplicity of y) are equal to the distinct parts in decreasing order.

Original entry on oeis.org

1, 1, 0, 0, 2, 0, 1, 0, 1, 1, 2, 0, 3, 0, 2, 0, 4, 1, 2, 1, 4, 1, 3, 1, 5, 3, 5, 1, 6, 2, 6, 1, 7, 2, 7, 2, 11, 4, 8, 3, 11, 5, 10, 4, 13, 5, 11, 5, 16, 8, 14, 5, 19, 8, 18, 6, 22, 8, 22, 7, 26, 10, 25, 8, 33, 12, 29, 11, 36, 13, 34, 12, 40, 16, 41, 14, 47, 17, 45, 16, 55
Offset: 0

Views

Author

Gus Wiseman, Mar 08 2019

Keywords

Comments

These are a kind of self-describing partitions (cf. A001462, A304679).
The Heinz numbers of these partitions are given by A324571.
The case where the distinct parts are taken in increasing order is counted by A033461, with Heinz numbers given by A109298.

Examples

			The first 19 terms count the following integer partitions:
   1: (1)
   4: (22)
   4: (211)
   6: (3111)
   8: (41111)
   9: (333)
  10: (511111)
  10: (322111)
  12: (6111111)
  12: (4221111)
  12: (33222)
  14: (71111111)
  14: (52211111)
  16: (811111111)
  16: (622111111)
  16: (4444)
  16: (442222)
  17: (43331111)
  18: (9111111111)
  18: (7221111111)
  19: (533311111)
		

Crossrefs

Sequences related to self-description: A000002, A001462, A079000, A079254, A276625, A304360.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Union[#]==Length/@Split[#]&]],{n,0,30}]

Extensions

More terms from Alois P. Heinz, Mar 08 2019

A304686 Numbers with strictly decreasing prime multiplicities.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 12, 13, 16, 17, 19, 20, 23, 24, 25, 27, 28, 29, 31, 32, 37, 40, 41, 43, 44, 45, 47, 48, 49, 52, 53, 56, 59, 61, 63, 64, 67, 68, 71, 72, 73, 76, 79, 80, 81, 83, 88, 89, 92, 96, 97, 99, 101, 103, 104, 107, 109, 112, 113, 116, 117, 121
Offset: 1

Views

Author

Gus Wiseman, May 16 2018

Keywords

Examples

			10 = 2*5 has prime multiplicities (1,1) so is not in the sequence.
20 = 2*2*5 has prime multiplicities (2,1) so is in the sequence
90 = 2*3*3*5 has prime multiplicities (1,2,1) so is not in the sequence.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[200],Greater@@FactorInteger[#][[All,2]]&]
  • PARI
    isok(n) = my(vm = factor(n)[,2]); vm == vecsort(vm,,4) && (#vm == #Set(vm)); \\ Michel Marcus, May 17 2018
    
  • PARI
    list(lim)=my(v=List()); forfactored(n=1,lim\1, if(n[2][,2]==vecsort(n[2][,2],,8), listput(v,n[1]))); Vec(v) \\ Charles R Greathouse IV, Oct 28 2021

Formula

a(n) ~ n log n. - Charles R Greathouse IV, Oct 28 2021

A324525 Numbers divisible by prime(k)^k for each prime index k.

Original entry on oeis.org

1, 2, 4, 8, 9, 16, 18, 27, 32, 36, 54, 64, 72, 81, 108, 125, 128, 144, 162, 216, 243, 250, 256, 288, 324, 432, 486, 500, 512, 576, 625, 648, 729, 864, 972, 1000, 1024, 1125, 1152, 1250, 1296, 1458, 1728, 1944, 2000, 2048, 2187, 2250, 2304, 2401, 2500, 2592
Offset: 1

Views

Author

Gus Wiseman, Mar 08 2019

Keywords

Comments

These are a kind of self-describing numbers (cf. A001462, A304679).
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. The prime signature of a number is the multiset of multiplicities (or exponents) in its prime factorization.
Also Heinz numbers of integer partitions where the multiplicity of k is at least k (A117144). The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Examples

			The sequence of terms together with their prime indices begins as follows. For example, 36 = prime(1) * prime(1) * prime(2) * prime(2) is a term because the prime multiplicities are {2,2}, which are greater than or equal to the prime indices {1,2}.
    1: {}
    2: {1}
    4: {1,1}
    8: {1,1,1}
    9: {2,2}
   16: {1,1,1,1}
   18: {1,2,2}
   27: {2,2,2}
   32: {1,1,1,1,1}
   36: {1,1,2,2}
   54: {1,2,2,2}
   64: {1,1,1,1,1,1}
   72: {1,1,1,2,2}
   81: {2,2,2,2}
  108: {1,1,2,2,2}
  125: {3,3,3}
  128: {1,1,1,1,1,1,1}
		

Crossrefs

Sequences related to self-description: A000002, A001462, A079000, A079254, A276625, A304360.

Programs

  • Maple
    q:= n-> andmap(i-> i[2]>=numtheory[pi](i[1]), ifactors(n)[2]):
    select(q, [$1..3000])[];  # Alois P. Heinz, Mar 08 2019
  • Mathematica
    Select[Range[1000],And@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>k>=PrimePi[p]]&]
    seq[max_] := Module[{ps = {2}, p, s = {1}, s1, s2, emax}, While[ps[[-1]]^Length[ps] < max, AppendTo[ps, NextPrime[ps[[-1]]]]]; Do[p = ps[[k]]; emax = Floor[Log[p, max]]; s1 = Join[{1}, p^Range[k, emax]]; s2 = Select[Union[Flatten[Outer[Times, s, s1]]], # <= max &]; s = Union[s, s2], {k, 1, Length[ps]}]; s]; seq[3000] (* Amiram Eldar, Nov 23 2020 *)

Formula

Closed under multiplication.
Sum_{n>=1} 1/a(n) = Product_{k>=1} 1 + 1/(prime(k)^(k-1) * (prime(k)-1)) = 2.35782843100111139159... - Amiram Eldar, Nov 23 2020

A324524 Numbers where every prime index divides its multiplicity in the prime factorization. Numbers divisible by a power of prime(k)^k for each prime index k.

Original entry on oeis.org

1, 2, 4, 8, 9, 16, 18, 32, 36, 64, 72, 81, 125, 128, 144, 162, 250, 256, 288, 324, 500, 512, 576, 648, 729, 1000, 1024, 1125, 1152, 1296, 1458, 2000, 2048, 2250, 2304, 2401, 2592, 2916, 4000, 4096, 4500, 4608, 4802, 5184, 5832, 6561, 8000, 8192, 9000, 9216
Offset: 1

Views

Author

Gus Wiseman, Mar 07 2019

Keywords

Comments

These are a kind of self-describing numbers (cf. A001462, A304679).
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. The prime signature of a number is the multiset of multiplicities (or exponents) in its prime factorization.
Also Heinz numbers of integer partitions in which every part divides its multiplicity (counted by A001156). The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
Also products of elements of A062457.

Examples

			The sequence of terms together with their prime indices begins as follows. For example, we have 18: {1,2,2} because 18 = prime(1) * prime(2) * prime(2).
    1: {}
    2: {1}
    4: {1,1}
    8: {1,1,1}
    9: {2,2}
   16: {1,1,1,1}
   18: {1,2,2}
   32: {1,1,1,1,1}
   36: {1,1,2,2}
   64: {1,1,1,1,1,1}
   72: {1,1,1,2,2}
   81: {2,2,2,2}
  125: {3,3,3}
  128: {1,1,1,1,1,1,1}
  144: {1,1,1,1,2,2}
  162: {1,2,2,2,2}
  250: {1,3,3,3}
  256: {1,1,1,1,1,1,1,1}
		

Crossrefs

Range of values of A090884.
Sequences related to self-description: A000002, A001462, A079000, A079254, A276625, A304360.

Programs

  • Maple
    q:= n-> andmap(i-> irem(i[2], numtheory[pi](i[1]))=0, ifactors(n)[2]):
    select(q, [$1..10000])[];  # Alois P. Heinz, Mar 08 2019
  • Mathematica
    Select[Range[1000],And@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>Divisible[k,PrimePi[p]]]&]
    v = Join[{1}, Prime[(r = Range[10])]^r]; n = Length[v]; vmax = 10^4; s = {1}; Do[v1 = v[[k]]; rmax = Floor[Log[v1, vmax]]; s1 = v1^Range[0, rmax]; s2 = Select[Union[Flatten[Outer[Times, s, s1]]], # <= vmax &]; s = Union[s, s2], {k, 2, n}]; Length[s] (* Amiram Eldar, Sep 30 2020 *)

Formula

Closed under multiplication.
Sum_{n>=1} 1/a(n) = Product_{k>=1} 1/(1-prime(k)^(-k)) = 2.26910478689594012492... - Amiram Eldar, Sep 30 2020

A324571 Numbers whose ordered prime signature is equal to the set of distinct prime indices in decreasing order.

Original entry on oeis.org

1, 2, 9, 12, 40, 112, 125, 352, 360, 675, 832, 1008, 2176, 2401, 3168, 3969, 4864, 7488, 11776, 14000, 19584, 29403, 29696, 43776, 44000, 63488, 75600, 104000, 105984, 123201, 151552, 161051, 214375, 237600, 267264, 272000, 335872, 496125, 561600, 571392, 608000
Offset: 1

Views

Author

Gus Wiseman, Mar 08 2019

Keywords

Comments

These are a kind of self-describing numbers (cf. A001462, A304679). The increasing case is A109298.
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. The ordered prime signature (A124010) is the sequence of multiplicities (or exponents) in a number's prime factorization, taken in order of the prime base.
Also Heinz numbers of the integer partitions counted by A324572. The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
Each finite set of positive integers determines a unique term with those prime indices. For example, corresponding to {1,2,4,5} is 1397088 = prime(1)^5 * prime(2)^4 * prime(4)^2 * prime(5)^1.

Examples

			The sequence of terms together with their prime indices begins as follows. For example, we have 40: {1,1,1,3} because 40 = prime(1) * prime(1) * prime(1) * prime(3).
      1: {}
      2: {1}
      9: {2,2}
     12: {1,1,2}
     40: {1,1,1,3}
    112: {1,1,1,1,4}
    125: {3,3,3}
    352: {1,1,1,1,1,5}
    360: {1,1,1,2,2,3}
    675: {2,2,2,3,3}
    832: {1,1,1,1,1,1,6}
   1008: {1,1,1,1,2,2,4}
   2176: {1,1,1,1,1,1,1,7}
   2401: {4,4,4,4}
   3168: {1,1,1,1,1,2,2,5}
   3969: {2,2,2,2,4,4}
   4864: {1,1,1,1,1,1,1,1,8}
   7488: {1,1,1,1,1,1,2,2,6}
  11776: {1,1,1,1,1,1,1,1,1,9}
  14000: {1,1,1,1,3,3,3,4}
  19584: {1,1,1,1,1,1,1,2,2,7}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[1000],Reverse[PrimePi/@First/@If[#==1,{},FactorInteger[#]]]==Last/@If[#==1,{},FactorInteger[#]]&]
Showing 1-7 of 7 results.