dimview: (Default)
[personal profile] dimview
1, 11, 21, 1211, 111221, 312211, 13112221, ...

Следующий элемент должен быть...



...1113213211.

Читаем 13112221 как «одна единица, одна тройка, две единицы, три двойки, одна единица». Эта последовательность называется look-and-say sequence, A005150 в OEIS.

Это была приcказка. Сказка в том, что каждое следующее число примено на 30% длиннее. Если быть точным, то если в десятичном представлении i-го элемента последовательности Li цифр, то Li+1/Li стремится к λ≈1.303577 при i, стремящемся к бесконечности.

Доказательство тривиально — достаточно перечислить 92 возможных подпоследовательности, нарисовать матрицу переходов, рассмотреть её как рекуррентное соотношение для длины элементов последовательности, использовать тот факт, что предел отношения элементов последовательности равен спектральному радиусу матрицы переходов, выписать соответствующий многочлен

x71
-x69-2x68-x67+2x66+2x65+x64-x63-x62-x61-x60
-x59+2x58+5x57+3x56-2x55-10x54-3x53-2x52+6x51+6x50
+x49+9x48-3x47-7x46-8x45-8x44+10x43+6x42+8x41-5x40
-12x39+7x38-7x37+7x36+x35-3x34+10x33+x32-6x31-2x30
-10x29-3x28+2x27+9x26-3x25+14x24-8x23-7x21+9x20
+3x19-4x18-10x17-7x16+12x15+7x14+2x13-12x12-4x11-2x10
+5x9+x7-7x6+7x5-4x4+12x3-6x2+3x-6


и найти его единственный положительный вещественный корень λ, например, с помощью wxMaxima:

    float(realroots(x^71
    -x^69-2*x^68-x^67+2*x^66+2*x^65+x^64-x^63-x^62-x^61-x^60
    -x^59+2*x^58+5*x^57+3*x^56-2*x^55-10*x^54-3*x^53-2*x^52+6*x^51+6*x^50
    +x^49+9*x^48-3*x^47-7*x^46-8*x^45-8*x^44+10*x^43+6*x^42+8*x^41-5*x^40
    -12*x^39+7*x^38-7*x^37+7*x^36+x^35-3*x^34+10*x^33+x^32-6*x^31-2*x^30
    -10*x^29-3*x^28+2*x^27+9*x^26-3*x^25+14*x^24-8*x^23-7*x^21+9*x^20
    +3*x^19-4*x^18-10*x^17-7*x^16+12*x^15+7*x^14+2*x^13-12*x^12-4*x^11-2*x^10
    +5*x^9+x^7-7*x^6+7*x^5-4*x^4+12*x^3-6*x^2+3*x-6, 1e-8));


Источник: Nathaniel Johnston
This account has disabled anonymous posting.
If you don't have an account you can create one now.
HTML doesn't work in the subject.
More info about formatting

Profile

dimview: (Default)
dimview

November 2017

S M T W T F S
   12 34
56789 1011
1213141516 1718
192021 22232425
26 27 282930  

Style Credit

Expand Cut Tags

No cut tags
Page generated Jun. 15th, 2025 09:06 pm
Powered by Dreamwidth Studios