Веб-версияОткрыть в Telegram
ООбсуждения компьютерной математики

Обсуждения компьютерной математики

@compmathchat · группа · Образование · в индексе с 2026-07-13
6участников
828сообщений в индексе
К
Ссылка
нажмите — покажем
рассмотрим последовательность a(1) = 7 a(n) = a(n-1)+НОД(n, a(n-1)) Rowland доказал, что каждый раз число увеличивается либо на 1, либо на простое число (но появляются ли так все нечетные простые, неизвестно) ( и такая задача предлагалась, как научили в комментариях, на Турнире городов — problems.ru/view_problem_details_new.php?id=64532 ) коллега Медведь поделился забавным родственником этой последовательности, который для разных начальных условий (гипотетически) генерирует всё большие простые числа-близнецы: https://math.stackexchange.com/q/5142627/152 (по ссылке как раз эксперименты и обсуждения)
11 · 1.4K ·
  1. М
    Эффект Баадера — Майнхоф: вчера перед сном как раз смотрел видео Rowland об этой задаче.)
  2. D
    Еще это задача старого тургора и там a1 = 6
      1. G
        ну да, там та же последовательность, просто индексы сдвинули типа. забавно что в решениях и тп вроде не упоминается Rowland (думал даже что мб на тургоре появилось раньше — но вроде нет)
        1. N
          Фотография
          нажмите — покажем
          зато, наоборот, у Роуланда упоминается Франк
      2. D
        Ну это я по памяти. Как обозначение начального элемента. Так то там вообще в терминах шагов :)
Вся ветка · 12 ответов →
К
Ссылка
нажмите — покажем
упомянутый в прошлом посте Rowland (относительно) недавно рассказывал, оказывается, на семинаре по экспериментальной математике вот про что возмьем квадратную матрицу неотрицательных чисел. и будем нормировать строки-столбцы: разделим каждую строку на сумму чисел в ней, потом каждый столбец, потом снова каждую строку… к чему это сойдется («Sinkhorn limit»)? вот, например, для матрицы [4 1\\ 2 1] можете сообразить, что это за числа получаются? можно посмотреть на приближенные значения: import numpy as np A = np.array([[4,1],[2,1]]) for _ in range(5): A = A / A.sum(axis=1, keepdims=True) A = A / A.sum(axis=0, keepdims=True) print(A) (для матриц произвольного размера жизнь быстро усложняется — и доклад был как раз про разные экспериментальные гипотезы по этому поводу… мб напишу позже какие-то подробности)
7 · 993 ·
    1. G
      именно так ну тут ясно, что в пределе должно быть два числа с суммой 1; указание: придумать еще инвариант, который сохраняется при операциях
  1. A
    Попробовал с другими числами, там так хорошо не укладывается 🤔🤔
Вся ветка · 3 ответа →
К
Ссылка
нажмите — покажем
во время ЛШСМ на компьютерные развлечения не хватает энергии, так что вот пока вместо моего поста — пост Тао: terrytao.wordpress.com/2026/07/14/visualizing-the-gilbreath-expectation-sequence/
5 · 470 ·
К
Видео
buffon-heptagon.mp4 · 90 КБ · нажмите — покажем
будем переходить от многоугольника к новому многоугольнику с вершинами в серединах сторон исходного («отображение Бюффона») если итерировать это отображение, то все вершины стремятся к одной точке… но что при этом происходит с формой N-угольника? оказывается, он становится аффинно правильным (это, кстати, достаточно просто доказать); в приземленных терминах: после нескольких итераций он практически становится «параллельником» — (практически) параллельны те же стороны и диагонали, что у правильного N-угольника (для 4-угольников уже после первой итерации получается параллелограмм… ну вот Бюффона можно считать таким обобщением Вариньона) можно посмотреть пример на гифке или попробовать разные многоугольника по ссылке dev.mccme.ru/~merzon/compmath/midpoints.html *** услышал про это в начале лекции А.П.Веселова на ЛШСМ-2026 — дальше обсуждался трехмерный случай, где все сложнее upd: появилось, кстати, видео — mathnet.ru/rus/present50983
15 · 1.1K ·
  1. Т
    А есть какое-то понимание станет многоугольник после конечного числа итераций точным «параллельником» или нет?
    1. G
      Для четного N есть еще «полупараллельники», которые за один ход становятся параллельниками (при N=4 каждый является полупараллельником, а дальше это очень рядкая ситуация — подпространство размерности 3 в пространстве размерности N). А для нечетного N так вообще не бывает. Ну там вот на картинке видны фурье-компоненты многоугольника (что бы это ни означало). Каждая домножается на свое собственное значение¹. Поэтому ненулевая координата редко может стать нулевой (компонента с нулевым собственным зачения бывает только одна и только при четном N). ¹ На картинке еще все отнормированно, поэтому длины столбиков меняются не буквально так.
  2. D
    А для большого количества точек он превращается в конце концов в эллипс
Вся ветка · 3 ответа →
К
Ссылка
нажмите — покажем
история про Rowland'а и Sinkhorn limit немного повисла в воздухе — вернемся ненадолго матрица 2×2, у которой суммы по строкам и столбцам единицы, имеет вид [x 1-x\\ 1-x x] чтобы понять, к какой из матриц такого вида мы сойдемся, нужно еще найти инвариант процесса — ну вот он такой соверешнно в духе мат. кружка: произведение чисел на черных клетках делить на произведение чисел на белых, ad/bc — получается уравнение на x, из которого x = √ad/(√ad+√bc) — вот эти квадратные корни и были видны в эксперименте конечно аналогичные инварианты можно писать для 2×2 подматриц матрицы 3×3 — и в итоге на элементы предельной матрицы получаются алгебраические уравнения степени 6 ну на этом история не заканчивается, можно почитать дальше доклад
5 · 563 ·
  1. D
    Если я правильно понимаю, не так давно они таки написали правильный многочлен с коэффициентами выражающимися через миноры изначальной матрицы, у которого x^{Sinkhorn}_{1, 1} должно быть корнем, но почему-то не получается это доказать. (знаки в суммах подсказал ChatGPT o3-mini-high...) Инвариантность относительно операций достаточно тривиальна, так что вопрос сводится к уравнению на элементы doubly stochastic matrix. Так как многочлен определяется довольно сложно через комбинаторику, то вопрос не самый простой, но мне кажется, что должно быть можно разобраться Я спросил LLM, что в этой задаче сейчас происходит, и она мне очень nonchalantly сообщила, что в статье https://arxiv.org/pdf/2506.06338, которая доказывает, что алгебраическая степень x^{Sinkhorn}_{1, 1} (над \Q с изначальными элементами или просто над \Q, если изначально матрица рациональная) не больше binom(2n-2, n-1), есть ошибка. Вроде и правда ошибка: в середине страницы 8 там пишется равенство вида dim(C\cap L)=dim(C)+dim(L)-N для проективных многообразий, хотя без пояснений верна только оценка снизу, и надо отедльно доказывать, что пересечение какое-то generic. Но я в деталях не разбирался, может LLM-ка чего-то не учла. Но ничего, говорит LLM, не кручинься, в статье 2017 года без каких либо упоминаний Sinkhorn limit уже доказали эту оценку! https://arxiv.org/pdf/1703.02251 Там явно этой теоремы нет, но вроде результат сразу следует из этой статьи. Думаю, сейчас AI может разобраться во всей этой комбинаторике и дорешать задачу. Какой в этом будет толк сказать сложно...
Вся ветка · 1 ответ →

Открытая публичная лента из поискового индекса ChatCrawler — «Google по публичному Telegram»; обновляется по мере обхода площадки. Время — UTC.

Только публичный контент, официальный API Telegram. О проекте · Вопросы · Чего мы не делаем · Убрать страницу из выдачи · Каталог · Поиск · Как мы считаем