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

A358793 Lexicographically earliest sequence of positive and unique integers such that 2*Sum_{k = 1..n} a(k) = Sum_{k = 1..n} a(a(k)) for n > 1 and a(1) = 1.

Original entry on oeis.org

1, 3, 7, 5, 10, 8, 14, 16, 11, 20, 22, 13, 26, 28, 17, 32, 34, 19, 38, 40, 23, 44, 46, 25, 50, 52, 29, 56, 58, 31, 62, 64, 35, 68, 70, 37, 74, 76, 41, 80, 82, 43, 86, 88, 47, 92, 94, 49, 98, 100, 53, 104, 106, 55, 110, 112, 59, 116, 118, 61, 122, 124, 65, 128
Offset: 1

Views

Author

Thomas Scheuerle, Dec 01 2022

Keywords

Comments

There is a second version of this sequence possible if we change the definition to a(1) = 2 and a(n) > 1, then the sequence will start 2, 4, 5, 8, 10, 7, 14, ... . It will after this continue in the same way as our actual sequence does (and would also extend the valid range of the recurrence formulas).
Start a(1) = 2 and value 1 allowed is A257794.

Crossrefs

Programs

  • Mathematica
    LinearRecurrence[{0, 0, 1, 0, 0, 1, 0, 0, -1}, {1, 3, 7, 5, 10, 8, 14, 16, 11, 20, 22, 13, 26, 28, 17}, 100] (* Paolo Xausa, Jan 23 2025 *)
  • PARI
    a(n) = {my(v = [1, 3, 7, 5, 10, 8]);if(n < 7, v[n], n*(1+min(1, n%3))+(n%3 == 0)+(n%6 == 3))}

Formula

G.f.: x*(1 + 3*x + 7*x^2 + 4*x^3 + 7*x^4 + x^5 + 8*x^6 + 3*x^7 - 4*x^8 + 2*x^9 - x^10 + x^11 - 3*x^12 + x^14)/(1 - x^3 - x^6 + x^9).
a(n) = a(n-3) + a(n-6) - a(n-9) for n >= 16.
a((3*(2*n-1) - (-1)^n)/4) = (3*(2*n-1) - (-1)^n)/2, for n > 3.
a(6*n) = 6*n+1, for n > 1.
a(6*n+3) = 6*n+5, for n > 0.
a(n) = 30*n - 2*a(n-1) - 3*a(n-2) - 3*a(n-3) - 3*a(n-4) - 3*a(n-5) - 2*a(n-6) - a(n-7) - 96, for n > 13.

A381464 Lexicographically earliest positive integer sequence satisfying a(n) = a(a(n))/n.

Original entry on oeis.org

1, 3, 6, 5, 20, 18, 8, 56, 10, 90, 12, 132, 14, 182, 16, 240, 19, 108, 323, 100, 22, 462, 24, 552, 26, 650, 28, 756, 30, 870, 32, 992, 34, 1122, 36, 1260, 38, 1406, 40, 1560, 42, 1722, 44, 1892, 46, 2070, 48, 2256, 50, 2450, 52, 2652, 54, 2862, 57, 448, 3135, 59, 3422, 61, 3660, 63, 3906, 65, 4160
Offset: 1

Views

Author

Thomas Scheuerle, Feb 24 2025

Keywords

Comments

While extending the sequence at a(k) we will check if k equals a previous term in the sequence. If such a term a(m) = k is found a(k) is determined as a(k) = a(m)*m. If no previous term matches k we may choose a(k) = k+c with the least c such that c > 0 and k+c does not equal any previous term in the sequence. It is conjectured that this check is sufficient. Reasoning behind this conjecture:
The greatest common divisor of two consecutive Fibonacci numbers is 1, thus we know that (k-1)^F(m)*k^F(m+1) and (t-1)^F(n)*t^F(n+1) are all different for some m,n > 1 if k and t are chosen such that for m or n < 2 no solution for (k-1)^F(m)*k^F(m+1) = (t-1)^F(n)*t^F(n+1) exist, because this cannot be equal if (k-1)*k and (t-1)*t have different prime numbers as divisors and if the only difference is the exponent of the prime factors, then the distribution of these between (t-1) and t and thus their progression F(n) or F(n+1) is individually distinct. In this sequence we need also to consider the more general case (k-c)^F(m)*k^F(m+1) = (t-1)^F(n)*t^F(n+1) because sometimes we need to set a(k) = k+c. It is conjectured that in this case c is bounded to be < 3.

Crossrefs

Cf. A099267 ( a(n) = a(a(n))-n ).

Programs

  • PARI
    listA(max_n) = {my(v=[1, 0], t=1); for(k=2, max_n, if(v[k]==0, t=1; if(k+t<#v, while(v[k+t]>0, t++)); v[k]=k+t); v=concat(v, vector(max(0, v[k]+1-#v))); if(v[v[k]]>0, print("The conjecture that a single forward check is sufficient failed at:", k)); v[v[k]]=k*v[k]); v[1..max_n]}

Formula

Let b(n, m) be m times recursion into a(n), for example b(3, 2) = a(a(3)).
b(3, m) = A000304(m+1) for m > 0.
b(n, m+2) = b(n, m)*b(n, m+1).
b(5, m) = 4^F(m)*5^F(m+1), where F(m) = A000045(m).
b(k, m) = (k-1)^F(m)*k^F(m+1), for all k where k+1 = a(k).
Showing 1-2 of 2 results.