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

A304360 Lexicographically earliest infinite sequence of numbers m > 1 with the property that none of the prime indices of m are in the sequence.

Original entry on oeis.org

2, 4, 5, 8, 10, 13, 16, 17, 20, 23, 25, 26, 31, 32, 34, 37, 40, 43, 46, 47, 50, 52, 61, 62, 64, 65, 67, 68, 73, 74, 79, 80, 85, 86, 89, 92, 94, 100, 103, 104, 107, 109, 113, 115, 122, 124, 125, 128, 130, 134, 136, 137, 146, 148, 149, 151, 155, 158, 160, 163
Offset: 1

Views

Author

Gus Wiseman, Aug 16 2018

Keywords

Comments

A self-describing sequence.
The prime indices of m are the numbers k such that prime(k) divides m.
The sequence is monotonically increasing, since once a number is rejected it stays rejected. Sequence is closed under multiplication for a similar reason. - N. J. A. Sloane, Aug 26 2018

Examples

			After the initial term 2, the next term cannot be 3 because 3 has prime index 2, and 2 is already in the sequence. The next term could be 10, which has prime indices 1 and 3, but 4 (with prime index 1) is smaller. So a(2) = 4.
		

Crossrefs

For first differences see A317963, for primes see A317964.

Programs

  • Maple
    A:= NULL:
    P:= {}:
    for n  from 2 to 1000 do
      pn:= numtheory:-factorset(n);
      if pn intersect P = {} then
        A:= A, n;
        P:= P union {ithprime(n)};
      fi
    od:
    A; # Robert Israel, Aug 26 2018
  • Mathematica
    gaQ[n_]:=Or[n==0,And@@Cases[FactorInteger[n],{p_,k_}:>!gaQ[PrimePi[p]]]];
    Select[Range[100],gaQ]

Extensions

Added "infinite" to definition. - N. J. A. Sloane, Sep 28 2019

A324695 Lexicographically earliest sequence of positive integers whose prime indices are not already in the sequence.

Original entry on oeis.org

1, 3, 7, 9, 11, 13, 19, 21, 27, 29, 33, 37, 39, 43, 47, 49, 53, 57, 59, 61, 63, 71, 77, 79, 81, 83, 87, 89, 91, 97, 99, 101, 107, 111, 113, 117, 121, 127, 129, 131, 133, 139, 141, 143, 147, 149, 151, 159, 163, 169, 171, 173, 177, 179, 181, 183, 189, 193, 197
Offset: 1

Views

Author

Gus Wiseman, Mar 10 2019

Keywords

Comments

A self-describing sequence, similar to A304360.
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 sequence of terms together with their prime indices begins:
   1: {}
   3: {2}
   7: {4}
   9: {2,2}
  11: {5}
  13: {6}
  19: {8}
  21: {2,4}
  27: {2,2,2}
  29: {10}
  33: {2,5}
  37: {12}
  39: {2,6}
  43: {14}
  47: {15}
  49: {4,4}
  53: {16}
  57: {2,8}
  59: {17}
  61: {18}
  63: {2,2,4}
		

Crossrefs

Programs

  • Mathematica
    aQ[n_]:=And@@Cases[If[n==1,{},FactorInteger[n]],{p_,k_}:>!aQ[PrimePi[p]]];
    Select[Range[100],aQ]

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

A277098 Finitary primes. Primes of finitary index.

Original entry on oeis.org

2, 3, 5, 11, 13, 29, 31, 41, 47, 79, 101, 109, 113, 127, 137, 167, 179, 211, 257, 271, 293, 313, 317, 397, 401, 421, 449, 487, 491, 547, 599, 601, 617, 677, 709, 733, 773, 811, 823, 829, 907, 929, 977, 991, 1033, 1063, 1109, 1187, 1231, 1259, 1297, 1361, 1429, 1483, 1489, 1559, 1609, 1621, 1741, 1759, 1831, 1871, 1889
Offset: 1

Views

Author

Gus Wiseman, Sep 29 2016

Keywords

Comments

Definition: prime(n) is a finitary prime iff n is a product of distinct finitary primes, where prime = A000040. This sequence could be described as a "multiplicative Aronson transform" of A005117. (Aronson transforms such as A079000 satisfy "n is in a if and only if a(n) is in b".
The composite bijection (finitary primes -> finitary numbers -> finite sets of finitary primes) can be used to construct a natural linear extension SET : N -> F where F is the partially ordered inverse limit of all finite Boolean algebras of finite sets of positive integers. Then a(n) = prime(Prod_p a(p)) where the product is over SET(n).

Examples

			The sequence of all nonempty finite sets of positive integers (a=1 b=2.. *=27) begins:
0,a,b,c,ab,ac,d,e,bc,ad,ae,f,abc,
g,bd,be,h,i,cd,af,ag,ce,abd,abe,
j,ah,bf,bg,ai,k,l,acd,m,bh,n,ace,
o,bi,de,cf,cg,aj,bcd,p,abf,q,abg,
bce,ak,ch,r,al,am,ci,bj,abh,an,s,
t,ao,abi,ade,acf,u,bk,acg,v,w,df,
bl,abcd,ap,bm,dg,aq,ef,bn,abce,
cj,x,y,eg,ach,bo,z,ar,bde,bcf,*
		

Crossrefs

Subsequence of A302491.

Programs

  • PARI
    has(p)=if(p<7, 1, my(f=factor(primepi(p))); if(vecmax(f[,2])>1, return(0)); for(i=1,#f~, if(!has(f[i,1]), return(0))); 1)
    is(n)=isprime(n) && has(n) \\ Charles R Greathouse IV, Aug 03 2023

Formula

a(n) = A000040(A276625(n)).

A324694 Lexicographically earliest sequence of positive integers divisible by prime(m) for some m not already in the sequence.

Original entry on oeis.org

2, 4, 5, 6, 8, 10, 12, 14, 15, 16, 17, 18, 20, 22, 23, 24, 25, 26, 28, 30, 31, 32, 34, 35, 36, 38, 40, 41, 42, 44, 45, 46, 48, 50, 51, 52, 54, 55, 56, 58, 60, 62, 64, 65, 66, 67, 68, 69, 70, 72, 73, 74, 75, 76, 78, 80, 82, 84, 85, 86, 88, 90, 92, 93, 94, 95
Offset: 1

Views

Author

Gus Wiseman, Mar 10 2019

Keywords

Comments

A self-describing sequence, similar to A304360.
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 sequence of terms together with their prime indices begins:
   2: {1}
   4: {1,1}
   5: {3}
   6: {1,2}
   8: {1,1,1}
  10: {1,3}
  12: {1,1,2}
  14: {1,4}
  15: {2,3}
  16: {1,1,1,1}
  17: {7}
  18: {1,2,2}
  20: {1,1,3}
  22: {1,5}
  23: {9}
  24: {1,1,1,2}
  25: {3,3}
  26: {1,6}
  28: {1,1,4}
  30: {1,2,3}
		

Crossrefs

Programs

  • Mathematica
    aQ[n_]:=!And@@Cases[If[n==1,{},FactorInteger[n]],{p_,k_}:>aQ[PrimePi[p]]];
    Select[Range[100],aQ]

A324696 Lexicographically earliest sequence containing 1 and all numbers divisible by prime(m) for some m not already in the sequence.

Original entry on oeis.org

1, 3, 6, 7, 9, 11, 12, 14, 15, 18, 19, 21, 22, 24, 27, 28, 29, 30, 33, 35, 36, 38, 39, 41, 42, 44, 45, 48, 49, 51, 53, 54, 55, 56, 57, 58, 59, 60, 63, 66, 69, 70, 71, 72, 75, 76, 77, 78, 81, 82, 83, 84, 87, 88, 90, 91, 93, 95, 96, 97, 98, 99, 101, 102, 105
Offset: 1

Views

Author

Gus Wiseman, Mar 10 2019

Keywords

Comments

A self-describing sequence, similar to A304360.
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 sequence of terms together with their prime indices begins:
   1: {}
   3: {2}
   6: {1,2}
   7: {4}
   9: {2,2}
  11: {5}
  12: {1,1,2}
  14: {1,4}
  15: {2,3}
  18: {1,2,2}
  19: {8}
  21: {2,4}
  22: {1,5}
  24: {1,1,1,2}
  27: {2,2,2}
  28: {1,1,4}
  29: {10}
  30: {1,2,3}
  33: {2,5}
  35: {3,4}
  36: {1,1,2,2}
		

Crossrefs

Programs

  • Mathematica
    aQ[n_]:=n==1||Or@@Cases[FactorInteger[n],{p_,k_}:>!aQ[PrimePi[p]]];
    Select[Range[100],aQ]

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

A324704 Lexicographically earliest sequence containing 1 and all numbers > 2 divisible by prime(m) for some m already in the sequence.

Original entry on oeis.org

1, 4, 6, 7, 8, 10, 12, 13, 14, 16, 17, 18, 19, 20, 21, 22, 24, 26, 28, 29, 30, 32, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 46, 48, 49, 50, 51, 52, 53, 54, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 70, 71, 72, 73, 74, 76, 77, 78, 79, 80, 82, 84
Offset: 1

Views

Author

Gus Wiseman, Mar 11 2019

Keywords

Comments

A self-describing sequence, similar to A304360.
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 sequence of terms together with their prime indices begins:
   1: {}
   4: {1,1}
   6: {1,2}
   7: {4}
   8: {1,1,1}
  10: {1,3}
  12: {1,1,2}
  13: {6}
  14: {1,4}
  16: {1,1,1,1}
  17: {7}
  18: {1,2,2}
  19: {8}
  20: {1,1,3}
  21: {2,4}
  22: {1,5}
  24: {1,1,1,2}
  26: {1,6}
  28: {1,1,4}
		

Crossrefs

Programs

A324698 Lexicographically earliest sequence containing 2 and all numbers > 1 whose prime indices already belong to the sequence.

Original entry on oeis.org

2, 3, 5, 9, 11, 15, 23, 25, 27, 31, 33, 45, 47, 55, 69, 75, 81, 83, 93, 97, 99, 103, 115, 121, 125, 127, 135, 137, 141, 155, 165, 197, 207, 211, 225, 235, 243, 249, 253, 257, 275, 279, 291, 297, 309, 341, 345, 347, 363, 375, 379, 381, 405, 411, 415, 419, 423
Offset: 1

Views

Author

Gus Wiseman, Mar 10 2019

Keywords

Comments

A self-describing sequence, similar to A304360.
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 sequence of terms together with their prime indices begins:
   2: {1}
   3: {2}
   5: {3}
   9: {2,2}
  11: {5}
  15: {2,3}
  23: {9}
  25: {3,3}
  27: {2,2,2}
  31: {11}
  33: {2,5}
  45: {2,2,3}
  47: {15}
  55: {3,5}
  69: {2,9}
  75: {2,3,3}
  81: {2,2,2,2}
  83: {23}
  93: {2,11}
  97: {25}
  99: {2,2,5}
		

Crossrefs

Programs

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
Showing 1-10 of 36 results. Next