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

A026477 a(1) = 1, a(2) = 2, a(3) = 3; and for n > 3, a(n) = smallest number > a(n-1) and not of the form a(i)*a(j)*a(k) for 1 <= i < j < k < n.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 9, 11, 13, 16, 17, 19, 23, 25, 29, 31, 37, 41, 43, 47, 49, 53, 59, 61, 67, 71, 73, 79, 81, 83, 89, 97, 101, 103, 107, 109, 113, 120, 121, 127, 131, 137, 139, 149, 151, 157, 163, 167, 168, 169, 173, 179, 181, 191, 193, 197, 199, 210, 211, 216, 223
Offset: 1

Views

Author

Keywords

Comments

From Bob Selcoe, Aug 25 2016: (Start)
Once a term with a given prime signature S (i.e., multiset of prime exponents) appears, then all numbers with the same prime signature follow. So either all or no terms with the same prime signature appear (first conjectured by Charles R Greathouse IV).
Proof:
i. Let S = {i,j,k..,y}, i>=j>=k.. be prime signatures, i.e., numbers of the form (p^i*q^j*r^k..*z^y) where {p,q,r,..,z} are distinct primes; denote S = {} as the signature for p^0 = 1.
ii. By definition, 1 and all primes p appear in the sequence; so terms where S = {k} (i.e., one-digit signatures denoting prime powers p^k) are k=0 and k = A026474(n) = {1,2,4,8,15,22,29..}, because k cannot be the sum of any combination of 3 smaller k. So the only prime power terms are p^0 = 1 and p^A026474(n), or S = {}, {1}, {2}, {4}, {8}, {15}, {22}, {29}...
iii. By induction, once S = {k} is determined, all other terms with the same prime signature appear (or do not) depending on the various combinations of signature exponents of previous terms with smaller signature sums. So for example, terms with the same two-digit signature S = {k,k} (i.e. (pq)^k) are constrained by the following: since p = {1} and q = {1} appear, then (pq) = {1,1} does not because in this case {1}*{1}*{} = {1,1}; since p^2 = {2} and q^2 = {2} appear, then (pq)^2 = {2,2} does not; since {4} appears but {3} does not, then {3,3} appears but {4,4} and {5,5} do not (the latter because {3,3}*{2}*{2} = {5,5} when p^2,q^2 = {2}). Note that {3,3} would not appear if {3,2} appeared because {3,2}*{1}*{} = {3,3} when q = {1}; but because p = {1} and p^2,q^2 = {2} appear, then {2}*{2}*{1} = {3,2} does not. {6,6} does not appear because {6,5} appears (by virtue of other constraints) and {6,5}*{1}*{} = {6,6} when q = {1}. Determining which signatures appear and which do not becomes increasingly complicated as the sequence increases.
Don Reble offered a proof on the Sequence Fans Mailing List which seems to be different (and certainly more formal) than mine. Perhaps mine is more of an "explanation" than a "proof"? (End)

Crossrefs

There are six related sequences: A026477: 1 <= i < j < k < n starting 1,2,3; A026478: 1 <= i <= j <= k < n starting 1,2,3; A026479: 1 <= i < j < k < n starting 1,2,4; A026480: 1 <= i <= j <= k < n starting 1,2,4; A026481: 1 <= i < j < k < n starting 1,3,4; A026482: 1 <= i <= j <= k < n starting 1,3,4.

Programs

  • Mathematica
    a = {1, 2, 3}; no = {1 2 3};
    Do[x = SelectFirst[Range[Last[a] + 1, 1000], ! MemberQ[no, #] &]; AppendTo[a, x]; no = Union[Times @@@ Subsets[a, {3}]], 200]; a (* Robert Price, May 26 2019 *)
  • PARI
    list(lim)=my(v=List(),n,d,k); while(n++<=lim, d=divisors(n); for(i=1,#d-2, if(!setsearch(v,d[i]), next); for(j=i+1,#d-1, if(!setsearch(v,d[j]), next); k=n/(d[i]*d[j]); if(d[j]>=k, break); if(denominator(k)==1 && setsearch(v,k), next(3)))); listput(v,n)); Vec(v) \\ Charles R Greathouse IV, Sep 16 2015

Extensions

More terms from Christian G. Bower, Nov 15 1999

A026478 a(n) = least positive integer > a(n-1) and not of form a(i)*a(j)*a(k) for 1<=i<=j<=k

Original entry on oeis.org

1, 2, 3, 5, 7, 11, 13, 16, 17, 19, 23, 24, 29, 31, 36, 37, 40, 41, 43, 47, 53, 54, 56, 59, 60, 61, 67, 71, 73, 79, 81, 83, 84, 88, 89, 90, 97, 100, 101, 103, 104, 107, 109, 113, 126, 127, 128, 131, 132, 135, 136, 137, 139, 140, 149, 150, 151, 152, 156, 157, 163, 167
Offset: 1

Views

Author

Keywords

Comments

Number of repeated prime divisors of a(n) is a square for a(n) < 128. - Jason Earls, Jul 01 2001
Numbers a(n) such that A001222(a(n)) is not a square are 128, 192, 288, 320, 432, 448, 480, 648, 672, 704, 720, 800, 832, 972, ... - Altug Alkan, Sep 26 2016

Examples

			13 is an obvious term because it is a prime.
15 is not a term because it is a semiprime; 15 = a(1)*a(3)*a(4) = 1*3*5.
		

Crossrefs

There are six related sequences: A026477: 1 <= i < j < k < n starting 1,2,3; A026478: 1 <= i <= j <= k < n starting 1,2,3; A026479: 1 <= i < j < k < n starting 1,2,4; A026480: 1 <= i <= j <= k < n starting 1,2,4; A026481: 1 <= i < j < k < n starting 1,3,4; A026482: 1 <= i <= j <= k < n starting 1,3,4.

Programs

  • PARI
    print1(1, ", "); for(n=1, 1e3, if(bigomega(n) % 3 == 1, print1(n, ", "))); \\ Altug Alkan, Sep 26 2016

Formula

1 together with numbers with 3m+1 prime factors (for some m >= 0).

Extensions

Formula from Henry Bottomley, Feb 09 2000
Definition corrected by and more terms from Larry Reeves (larryr(AT)acm.org), Mar 24 2000
Terms corrected (128 removed) by Charles R Greathouse IV, Aug 25 2016
As pointed out by Don Reble, 128 IS a member of this sequence. - N. J. A. Sloane, Sep 23 2016

A026479 a(n) = least positive integer > a(n-1) and not a(i)*a(j)*a(k) for 1<=i

Original entry on oeis.org

1, 2, 4, 5, 6, 7, 9, 11, 13, 15, 16, 17, 19, 21, 23, 25, 27, 29, 31, 33, 37, 39, 41, 43, 47, 49, 51, 53, 57, 59, 61, 67, 69, 71, 73, 79, 81, 83, 87, 89, 93, 97, 101, 103, 107, 109, 111, 113, 121, 123, 127, 129, 131, 137, 139, 141, 149, 151, 157, 159, 163
Offset: 1

Views

Author

Keywords

Crossrefs

There are six related sequences: A026477: 1 <= i < j < k < n starting 1,2,3; A026478: 1 <= i <= j <= k < n starting 1,2,3; A026479: 1 <= i < j < k < n starting 1,2,4; A026480: 1 <= i <= j <= k < n starting 1,2,4; A026481: 1 <= i < j < k < n starting 1,3,4; A026482: 1 <= i <= j <= k < n starting 1,3,4.

Programs

  • Mathematica
    a = {1, 2, 4}; no = {1 2 4};
    Do[x = SelectFirst[Range[Last[a] + 1, 1000], ! MemberQ[no, #] &]; AppendTo[a, x]; no = Union[Times @@@ Subsets[a, {3}]], 200]; a (* Robert Price, May 26 2019 *)

Extensions

a(31)-a(61) from Robert Price, May 26 2019

A026480 a(1) = 1, a(2) = 2, a(3) = 4, and for n > 3, a(n) is the least positive integer > a(n-1) not of form a(i)*a(j)*a(k) for 1 <= i <= j <= k < n.

Original entry on oeis.org

1, 2, 4, 5, 6, 7, 9, 11, 13, 15, 17, 19, 21, 23, 27, 29, 31, 33, 37, 39, 41, 43, 47, 51, 53, 57, 59, 61, 67, 69, 71, 73, 79, 83, 87, 89, 93, 97, 101, 103, 107, 109, 111, 113, 123, 127, 128, 129, 131, 137, 139, 141, 149, 151, 157, 159, 160, 163, 167, 173, 177, 179, 181
Offset: 1

Views

Author

Keywords

Comments

From David A. Corneth, Sep 24 2016 (Start):
Let t be a term of this sequence. Let v = (i, j, m) be a vector of three elements where i is A007814(t), j = A007949(t) and m = A001222(t) - i - j. Then possible vectors v for t up to 100000 are {[0, 0, 0], [0, 0, 1], [0, 0, 4], [0, 0, 7], [0, 1, 1], [0, 1, 4], [0, 2, 0], [0, 3, 0], [0, 4, 3], [0, 5, 3], [0, 6, 2], [0, 7, 2], [0, 8, 1], [0, 9, 1], [0, 10, 0], [1, 0, 0], [1, 0, 3], [1, 0, 6], [1, 1, 0], [1, 1, 3], [1, 1, 6], [1, 4, 2], [1, 5, 2], [1, 6, 1], [1, 7, 1], [1, 8, 0], [1, 9, 0], [2, 0, 0], [2, 3, 2], [2, 5, 1], [2, 7, 0], [3, 0, 2], [3, 0, 5], [3, 1, 2], [3, 1, 5], [3, 3, 1], [3, 4, 1], [3, 5, 0], [3, 6, 0], [4, 2, 1], [4, 2, 4], [4, 4, 0], [5, 0, 1], [5, 0, 4], [5, 1, 1], [5, 1, 4], [5, 2, 0], [5, 3, 0], [6, 0, 3], [6, 1, 0], [6, 1, 3], [7, 0, 0], [7, 3, 2], [8, 0, 2], [8, 1, 2], [8, 3, 1], [8, 5, 0], [9, 2, 1], [9, 4, 0], [10, 0, 1], [10, 1, 1], [10, 2, 0], [10, 3, 0], [11, 1, 0], [12, 0, 0]}.
(End)

Crossrefs

There are six related sequences: A026477: 1 <= i < j < k < n starting 1,2,3; A026478: 1 <= i <= j <= k < n starting 1,2,3; A026479: 1 <= i < j < k < n starting 1,2,4; A026480: 1 <= i <= j <= k < n starting 1,2,4; A026481: 1 <= i < j < k < n starting 1,3,4; A026482: 1 <= i <= j <= k < n starting 1,3,4.

Programs

  • Mathematica
    a = {1, 2, 4}; no = {1 2 4};
    Do[x = SelectFirst[Range[Last[a] + 1, 1000], ! MemberQ[no, #] &]; AppendTo[a, x]; no = Union[Times @@@ Tuples[a, {3}]], 60]; a (* Robert Price, May 26 2019 *)

Extensions

Name corrected by Charles R Greathouse IV, Sep 23 2016

A026481 a(n) = least positive integer > a(n-1) and not a(i)*a(j)*a(k) for 1<=i

Original entry on oeis.org

1, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 16, 17, 19, 22, 23, 25, 26, 29, 31, 34, 37, 38, 41, 43, 46, 47, 49, 53, 58, 59, 61, 62, 67, 71, 73, 74, 79, 81, 82, 83, 86, 89, 94, 97, 101, 103, 106, 107, 109, 113, 118, 121, 122, 127, 131, 134, 137, 139, 142, 146, 149, 151, 157
Offset: 1

Views

Author

Keywords

Crossrefs

There are six related sequences: A026477: 1 <= i < j < k < n starting 1,2,3; A026478: 1 <= i <= j <= k < n starting 1,2,3; A026479: 1 <= i < j < k < n starting 1,2,4; A026480: 1 <= i <= j <= k < n starting 1,2,4; A026481: 1 <= i < j < k < n starting 1,3,4; A026482: 1 <= i <= j <= k < n starting 1,3,4.

Programs

  • Mathematica
    a = {1, 3, 4}; no = {1 3 4};
    Do[x = SelectFirst[Range[Last[a] + 1, 1000], ! MemberQ[no, #] &]; AppendTo[a, x]; no = Union[Times @@@ Subsets[a, {3}]], 62]; a (* Robert Price, May 26 2019 *)

Extensions

More terms from Todd Will (towill(AT)davidson.edu).
Showing 1-5 of 5 results.