Долгие сообщения в RabbitMQ
Почему стандартная Python библиотека для RabbitMQ не работает с длинными сообщениями и как это исправить.
Блог 529
Project Euler и остальное
Долгие сообщения в RabbitMQ
Почему стандартная Python библиотека для RabbitMQ не работает с длинными сообщениями и как это исправить.
Экспорт партий с Lichess
Сохраняем себе шахматные партии с Lichess.org. Рабочий скрипт на Python прилагается.
Моё решение задачи 134
Краткое условие: назовём порождающим для двух последовательных простых \(p_1 < p_2\) наименьшее натуральное число, что оно закачивается на \(p_1\) и при этом делится на \(p_2\). Необходимо найти сумму порождающих для всех \(p_1 \in \left[ 5; 10^6 \right]\)
-
Моё решение задачи 146
Краткое условие: необходимо найти сумму всех натуральных \(n\), что \(n^2+1\), \(n^2+3\), \(n^2+7\), \(n^2+9\), \(n^2+13\), и \(n^2+27\) будут последовательными простыми числами.
Моё решение задачи 146
Краткое условие: необходимо найти сумму всех натуральных \(n\), что \(n^2+1\), \(n^2+3\), \(n^2+7\), \(n^2+9\), \(n^2+13\), и \(n^2+27\) будут последовательными простыми числами.
-
Нахождение суммы k-ых степеней
Как придумать формулу для суммы \(1^5 + 2^5 + 3^5 + \ldots + n^5\) и есть ли она вообще?
Нахождение суммы k-ых степеней
Как придумать формулу для суммы \(1^5 + 2^5 + 3^5 + \ldots + n^5\) и есть ли она вообще?
Wallabag и реальная жизнь
Как я устанавливал известнейшее свободное решение для отложенного чтения, и что из этого получилось.
Как я шахматного бота писал
Как я проверял шахматное приложение на “ботоустойчивость”.
CrossGen v1.0
Читая хабр, случайно натолкнулся на идею сделать программу, которая по заданной кроссвордной сетке находит способ её заполнить. В этом посте вкратце напишу про моё решение и первую версию приложения.
Моё решение задачи 60
Краткое условие: необходимо найти множество из пяти простых чисел с минимальной суммой такое, что после “склеивания” в любом порядке любых двух чисел из него тоже будет простое число.
Ещё одно вычисление выражений
На хабре когда-то увидел статью про то, что в Яндексе двум сотрудникам дали задачу на написание приложения, для вычисления выражений. Менеджер справился за 4 часа, а программист за два. Я решил попробовать свои силы.
Wallabag и реальная жизнь
Как я устанавливал известнейшее свободное решение для отложенного чтения, и что из этого получилось.
Как я шахматного бота писал
Как я проверял шахматное приложение на “ботоустойчивость”.
CrossGen v1.0
Читая хабр, случайно натолкнулся на идею сделать программу, которая по заданной кроссвордной сетке находит способ её заполнить. В этом посте вкратце напишу про моё решение и первую версию приложения.
Моё решение задачи 60
Краткое условие: необходимо найти множество из пяти простых чисел с минимальной суммой такое, что после “склеивания” в любом порядке любых двух чисел из него тоже будет простое число.
\ No newline at end of file + \ No newline at end of file diff --git a/index.html.gz b/index.html.gz index d1577fc..ff2d00f 100644 Binary files a/index.html.gz and b/index.html.gz differ diff --git a/index2.html b/index2.html index 2a11213..4b34753 100644 --- a/index2.html +++ b/index2.html @@ -1,6 +1,6 @@ -
Блог 529
Project Euler и остальное
Мой первый пост или зачем этот блог?
Кратко о том, почему я сделал блог.
Здравствуйте, я Алексей Лобанов
Сейчас я студент ВМК МГУ. Немного занимаюсь фрилансом и машинным обучением.
Умею писать:
- Бекенды на
Python(+flask) - Фронтенды на
jQuery - Парсеры на
scrapy - Кроссплатформенные приложения на
C++/Pythonс использованием wxWidgets - И интегрировать во всё это машинное обучение на
scikit-learn,xgboostиTensorFlow.
Если вы хотите связаться со мной и обсудить что-нибудь, то можете выбрать удобный способ связи:
Блог 529
Project Euler и остальное
Мои проекты
Значительная часть моих проектов есть на GitHub или BitBucket вместе с открытым исходным кодом.
CrossGen
Что использовалось
- C++
- wxWidgets
Что реализовано
- Графический интерфейс
- Автоматическая локализация всего интерфейса на русский и английский языки
- Быстрая генерация кроссворда по заданной сетке с использованием эвристики
- Автоматический бенчмарк для измеренеия производительноти генерации Подробнее я писал тут
За подробностями пишите мне на почту
Блог 529
Project Euler и остальное
CrossGen v1.0
Начать, наверное, нужно с того, что реальную практическую значимость я осознал после того, как реализовал 90% того, что есть сейчас. Сейчас мне кажется, что единственное применение данного приложения лишь в том, чтобы создавать кроссворды очень сложной или необычной формы. Зачем это надо обычному человеку я вообще не знаю.
Собственно, алгоритм генерации изначально был примитивным: простой рекурсивный поиск с отсечением. Скорость генерации более-менее сложных сеток была ужасной (для перцентиля 30% это примерно 40 минут на этой сетке, в общем случае, время генерации непредсказуемо.
В дальнейшем были выполнены некоторые оптимизации. Первой более-менее значимой стала замена передачи сетки в юникоде (во внутреннем цикле) на передачу сетки в однобайтовой кодировке, таким образом, языки с алфавитом больше ~200 букв пролетают. Впрочем, мне кажется, что им не слишком сильно требуются кроссворды. Такая оптимизация дала ~35% прироста при значимом времени перебора (больше секунды).
В какой-то момент, мне показалось, что оптимизация структуры данных для хранения сетки тоже могла сильно увеличить производительность, но это будет заметно только на разряженных сетках, которые и так достаточно быстро генерируются. На сложных же плотных сетках, прирост скорости может быть минимальным, вплоть до отрицательного. Таким образом, используются просто двумерный массив.
Наибольший прирост, как и ожидалось, дало упорядочивание словаря по некоторому критерию. Таким образом, при переборе, “плохие” варианты будут попадаться редко. На это ушло немногим больше дня. В конечном итоге, целевой функцией, которая стала критерием сортировки, стала такая \(\sqrt[n]{\prod_{i=1}^{n}{f\left( a_i \right)}}\) (иначе говоря, среднее геометрическое отлично подошло), где \(n\) — количество букв в слове, \(a_i\) — буквы слова, а \(f \left( x \right)\) есть вероятность появления буквы \(x\). Я отсекаю примерно половину самых “плохих” слов, а из того, что осталось я составляю кроссворд. Производительность, в среднем, увеличилась в ~1000 раз. На тестовых сетках для перцентиля 90% длительность работы около секунды.
Что реализовано
- Прозрачная интернационализация (пока только русский и английский языки)
- Быстрая генерация
- Экспорт в текстовом формате
- Простая смена словаря
- Относительно простое редактирование сеток
Что может быть реализовано
- Визуальный редактор сетки
- Экспорт в HTML, PDF, Markdown
- Экспорт сетки в файл картинки
Скриншоты
- Версия для Xubuntu:
- Версия для Windows 7:
\ No newline at end of file + \ No newline at end of file diff --git a/posts/crossgen-v10/index.html.gz b/posts/crossgen-v10/index.html.gz index 66f53cd..86cd1fb 100644 Binary files a/posts/crossgen-v10/index.html.gz and b/posts/crossgen-v10/index.html.gz differ diff --git a/posts/dolgie-soobshcheniia-v-rabbitmq/index.html b/posts/dolgie-soobshcheniia-v-rabbitmq/index.html index 2d98fb4..783faae 100644 --- a/posts/dolgie-soobshcheniia-v-rabbitmq/index.html +++ b/posts/dolgie-soobshcheniia-v-rabbitmq/index.html @@ -1,138 +1,138 @@ -
Блог 529
Project Euler и остальное
Долгие сообщения в RabbitMQ
Предположим, что у вас появилось желание перекодировать фильмы на вашем медиасервере, и вы решили использовать production-ready решение для хранения заданий. Вы взяли RabbitMQ для управления очередями сообщений и Python для их обработки. Но почему-то сообщения обрабатываются нестабильно, клиент падает без всяких видимых причин. Попробуем понять почему такое может быть.
Возьмём готовый код из официального туториала RabbitMQ и немного его модифицируем, чтобы он обрабатывал сообщения за различное время.
Producer
#!/usr/bin/python -# -*- coding: utf-8 -*- - -__author__ = "Aleksey Lobanov" -__license__ = "MIT" - -import sys - -import pika - -if __name__ == "__main__": - connection = pika.BlockingConnection( - pika.ConnectionParameters(host='localhost') - ) - channel = connection.channel() - - - channel.queue_declare(queue='demo.hello') - - # чтение первого аргумента командной строки, если он есть - if len(sys.argv) > 1: - delay_to_send = int(sys.argv[1]) - else: - delay_to_send = 3 - channel.basic_publish( - exchange='', - routing_key='demo.hello', - body=str(delay_to_send) - ) - print(" [x] Sent " + str(delay_to_send)) - connection.close() -
Consumer
#!/usr/bin/python -# -*- coding: utf-8 -*- - -__author__ = "Aleksey Lobanov" -__license__ = "MIT" - -import sys -import time -import math - -import pika - - -def callback(ch, method, properties, body): - print(" [x] Received %r" % body) - delay = int(body) - begin_at = time.time() - time.sleep(delay) - print(" [x] Finished {}".format(body)) - - -if __name__ == "__main__": - connection = pika.BlockingConnection(pika.ConnectionParameters(host='localhost')) - channel = connection.channel() - - channel.queue_declare(queue='demo.hello') - - channel.basic_consume( - callback, - queue='demo.hello', - no_ack=True - ) - - print(' [*] Waiting for messages. To exit press CTRL+C') - channel.start_consuming() -
Пока мы посылаем маленькие числа, всё будет хорошо работать. Но если послать какое-то большое (в моём случае достаточно 200), то consumer потеряет соединение с сервером, будет ошибка. Скорее всего это будет
pika.exceptions.ConnectionClosed: (-1, 'EOF')илиSocket Error 104(тут есть обсуждение на GitHub библиотеки, но установкаprefetch_count=1тоже не поможет). Эта проблема актуальна для обоих веток Python.Реальная причина в том, что при обработке сообщения не происходит необходимого взаимодействия с RabbitMQ, не отправляются hearbeats, а без них сервер считает, что этот клиент погиб окончательно. Нужно отметить, что переход на другой тип соединения не помогает. Например использование примера на Twisted из официальной документации ничего не изменит.
У этой проблемы есть много решений. Будем считать, что сообщение, которое мы обрабатываем не разделяется на подзадачи и рассмотрим некоторые из них:
- Отключить hearbeats/увеличить их интервал так, чтобы самprefetch_countая долгая обработка сообщения вела к потере не более, чем одного. Это самое простое решение, но в таком случае мы теряем в надёжности. Внешнем сервисам мониторинга будет сложнее понять, обрабатываются ли сейчас сообщения или уже нет. И чем больше интервал, тем серьёзнее проблема. Не подходит
- Разбить обработку сообщения на несколько этапов. Часто это хорошее решение, но в предположении, что данная обработка не разбивается на более мелкие тоже не подходит.
- Использовать
connection.sleepвместоtime.sleep, а также регулярный вызовBlockingConnection.process_data_events. Оба эти решения помогают, но их использование плохо по многим причинам. Во-первых это явное протекании абстракции, когда код для обработки сообщения вынужден работать с очередью. Во-вторых не всегда можно гарантировать, что эти функции вызываются достаточно часто, а это главная проблема. Не подходит - Выделить отдельный процесс под
pika. Вероятно, единственно универсальное решение. Если исходный код был правильно написан, то адаптация будет простой. Но у этого решения есть минусы, обязательные при использовании нескольких потоков/процессов. Также сама библиотекаpikaне является потокобезопасной.
Возможное решение будет заключаться в вынесение обработки данных в отдельный метод:
def real_work(body): - delay = int(body) - begin_at = time.time() - time.sleep(delay) -
и небольшой доработке callback-метода:
# Конструктор, при запуске процесс запустит -# функцию с указанными параметрами -work_process = multiprocessing.Process( - target=real_work, - args=(body, ) -) -work_process.start() # процесс нужно явно запустить -while True: - ch.connection.sleep(1) - # пока поток работает, используем функцию pika - # для обработки необходимых сообщений - if not work_process.is_alive(): - break -
Тогда готовый код будет выглядеть так:
#!/usr/bin/python -# -*- coding: utf-8 -*- - -__author__ = "Aleksey Lobanov" -__license__ = "MIT" - -import sys -import time -import math -import multiprocessing - -import pika - - -def real_work(body): - delay = int(body) - begin_at = time.time() - time.sleep(delay) - - -def callback(ch, method, properties, body): - print(" [x] Received %r" % body) - work_process = multiprocessing.Process( - target=real_work, - args=(body, ) - ) - work_process.start() - while True: - ch.connection.sleep(1) - if not work_process.is_alive(): - break - - print(" [x] Finished " + body) - - -if __name__ == "__main__": - connection = pika.BlockingConnection(pika.ConnectionParameters(host='localhost')) - channel = connection.channel() - - channel.queue_declare(queue='demo.hello') - - channel.basic_consume( - callback, - queue='demo.hello', - no_ack=True - ) - - print(' [*] Waiting for messages. To exit press CTRL+C') - channel.start_consuming() -
Блог 529
Project Euler и остальное
Экспорт партий с Lichess
Одну партию с Lichess скачать довольно просто, это можно сделать на странице с самой игрой. Скачать все игры тоже несложно, для этого есть специальный раздел.
Если же хочется скачать свои партии, то я нашёл лишь одно решение. Последний коммит был сделан 2 года назад, также оно не учитывает ограничения API, и больше ~2000 партий, за раз скачать не получится.
Моей целью было написание альтернативы, которая бы работала быстрее и работала корректно.
Проблема формата решена довольно просто: API отдаёт PGN, значит его поддержка уже есть. Но это плохой формат для хранения и обработки. Поэтому нужно было выбрать что-то ещё.
Существует, практически, только Scid, если вы хотите хранить свои шахматные партии, используя открытый софт. В его дистрибутиве уже есть необходимые утилиты для конвертации pgn в свой формат.
Репозиторий здесь, а здесь можно скачать файлы сразу.
Возможности
- Многопоточность
- Сохранение как в PGN, так и в формат Scid
- Для работы необходим только Python
- Работает при любом числе партий
Установка
pip install grequests -Запуск
python lichess.py -n hippo23 -t pgn -o hippo23.pgn --threads 6 -Здесь
hippo23мой ник, аpgnформат сохранения (доступен такжеscid).
Блог 529
Project Euler и остальное
Ещё одно вычисление выражений
Задачка кажется не очень сложной, даже, если не знать как её делать (я не знал). Целью является быстрое вычисление чего-то типа
4 * ( 5 + 7 ^ 4). Для это я парсил исходную строку в список токенов, а затем непосредственно вычислял, что получится.Я решил, что проще всего будет реализовать (а мне потом и понять) алгоритм, когда после каждого действия будет выполняться некий “хороший” инвариант. Первое что приходит в голову — это то, что истинность выражение после выполнения операции не меняется (TITO соблюдается). То есть выражение
3 + 5можно заменить на8или хотя бы на4 * 2.Непосредственно сама обработка является несколькими проходами, так что в каждом проходе мы избавляемся от операций одного приоритета.
4 + 5 * 3заменяется на4 + 15, 7 - 5 * 2^3заменяется на7 - 5*8. Таким образом, каждый цикл тривиален, и легко задавать приоритеты операций.Если использовать один список как контейнер для токенов и при работе изменять непосредственно его, сохраняя указанные инварианты, то сложность получается \(O\left( N \right)\), где \(N\) — число токенов.
Времени на непосредственно кодирование ушло часа три-четыре, но в это время не входит продумывание мелких деталей.
Всё написано на C++11. Исходники лежат на GitHub и BitBucket.
\ No newline at end of file + \ No newline at end of file diff --git a/posts/eshchio-odno-vychislenie-vyrazhenii/index.html.gz b/posts/eshchio-odno-vychislenie-vyrazhenii/index.html.gz index a7374af..fd378be 100644 Binary files a/posts/eshchio-odno-vychislenie-vyrazhenii/index.html.gz and b/posts/eshchio-odno-vychislenie-vyrazhenii/index.html.gz differ diff --git a/posts/kak-ia-delaiu-bekapy-pochty/index.html b/posts/kak-ia-delaiu-bekapy-pochty/index.html new file mode 100644 index 0000000..ac08934 --- /dev/null +++ b/posts/kak-ia-delaiu-bekapy-pochty/index.html @@ -0,0 +1,14 @@ +
Блог 529
Project Euler и остальное
-
Как я делаю бэкапы почты
“Люди делятся на две категории: кто еще не делает бэкапы, и кто их уже делает” (с)
С бэкапами файлов всё, в целом, довольно понятно — есть множество решений разной сложности, как открытые, так и коммерческие. Если взять что-нибудь более специфическое, даже обычную почту, то открытые решения найти уже сложнее.
Мне кажется, что хорошее решение точно должно удовлетворять следующим требованиям:
- Дружественность file-based системам синхронизации: Dropbox, Яндекс.Диск и другие.
- Простота потенциального переноса данных на другое решение бэкапов. Не хочется сталкиваться с vendor lock-in в случае проблем.
- Бэкапы должны быть зашифрованы.
- Естественная поддержка нескольких почтовых ящиков.
- Хочется легко извлекать результаты простых запросов, например, “все письма с таким автором”, без сложных процедур импорта.
До этого у меня на домашнем сервере был настроен бэкап по IMAP по каждому из аккаунтов. Хотелось получить какой-то более простой интерфейс к поиску, чем grep и не хранить лишнее.
В итоге я пришёл к самописному решению, которое удовлетворяло всем требованиям, оставаясь компактным при этом. Схема работы такая:
- Есть два места для сохранения результатов: небольшая sqlite база с метаданными и некоторое файловое хранилище (далее просто хранилище), которое может быть S3 корзиной, локальной директорией или чем-нибудь ещё. Хранилище используется только на запись, поэтому поверх него можно добавить GPG шифрование.
- Сначала используется mail-parser для извлечения простых метаданных из данных письма.
- По этим метаданным добавляем запись в sqlite базу с адресами участников, темой письма и временем. Этого достаточно для большинства простых запросов.
- Извлекаем вложения, метаданные по ним сохраняем в sqlite, а сами файлы — в хранилище.
- Сохраняем полные данные письма, кроме уже сохранённых вложений, в хранилище.
Проще всего поддерживать сохранение объектов в локальную директорию, поэтому сейчас поддерживается только оно. Слой шифрования достаточно естественно добавляется ещё одним решением, gocryptfs reverse mode. Мой полный бэкап почты остаётся достаточно компактным, чтобы можно было не разделять метаданные и объекты. В таком случае процесс бэкапа целиком выглядит так:
- Получить mbox по всем аккаунтам
- Обработать mbox и сохранить новые объекты, добавить новые метаданные.
- Добавить слой шифрования с помощью reverse mode в gocryptfs.
- Сохранить на бэкап хранилища.
Актуальная версия моего решения вместе с подробной информацией по использованию лежит на GitHub.
Блог 529
Project Euler и остальное
Как я шахматного бота писал
Лет 5 назад я достаточно активно играл в живые шахматы. Потом времени на это стало не хватать и постепенно перешёл на редкие партии в онлайне. Сейчас для игры я использую одно из самых популярных приложений вк. Это проще, чем использовать, к примеру, FICS. Предмет обсуждения появился из-за того, что я как-то раз встретился с соперником, который на все ходы потратил порядка 10 секунд, при этом не допустив значимых ошибок. Тогда я решил написать своего бота, чтобы узнать что с ним будет и столкнусь ли я с какими-нибудь подводными камнями.
Целью было максимально быстрое написание максимально простого решения. Поэтому от разбора протокола я сразу отказался, тем более у меня не было подобного опыта ранее. Была мысль работать с FICS (у меня есть библиотека для работы с их протоколом), но поскольку я там не играю, то и результаты были бы не так интересны, во всяком случае, для меня. Таким образом, я писал простого кликера для приложения вк.
Изначально хотелось найти доску и определить положения всех фигур, это было бы достаточно универсально, хотя и привязало бы меня к OpenCV. Тем не менее, решил не усложнять: можно определять только последний ход, а это можно сделать, проверяя цвет только одного пикселя.
В самом скрипте около 200 строк на python. Очень сильно помогла библиотека chess, которая взяла на себя общение с движком (я использовал stockfish), проверку на допустимые ходы и определение мата. Некоторое время я уделил тому, чтобы сделать бота максимально похожим на человека, чтобы было невозможно выявить, что это бот полностью автоматическими средствами. По пунктам:
- Клик по полю в случайном месте, с распределением по Гауссу, центр которого не совпадает с центром клетки
- Случайное время хода, длительность которого распределена по Гауссу, причём средняя длительность хода изменяется, в зависимости от номера текущего хода.
- Прокладываются дополнительные точки, с распределением по Гауссу, при перемещении курсора от точки к точке.
По факту, всё это было лишним, бан получить не удалось даже при простом клике из начальной точки в конечную.
Примеры работы скрипта можно посмотреть тут и тут, анализы двух сыгранных игр лежат тут и тут (оппонент имеет рейтинг около 2100).
На момент публикации аккаунт вполне жив. Рейтинг достиг некоторого потолка (около 2200), после которого найти игроков примерно равного рейтинга, не являющихся ботами, очень сложно. Сражаться же с ботами сильно сложнее, такую цель я не ставил. Интересно, хоть и ожидаемо, что при наборе рейтинга было достаточно личностей, для которых возможность того, что их нагло обманули, и они играли с ботом была столь неприятна, что они не могли сдерживаться. Например:
p.s. Уже после создания рабочей версии от одного из оппонентов узнал про lichess.org. Это отличный ресурс на котором кроме, собственно, платформы для игры в шахматы (поддерживается большое количество их вариантов), есть тренировки по дебютам, анализ игр. Самое интересное — ресурс полностью открытый, все исходники есть на github.
В процессе подготовки данного материала узнал про InternetChessKiller, который делает, фактически, тоже самое, что и мой скрипт, но без привязки к какой-то одной игровой площадке. Исходники старых версий можно найти, например, в этом репозитории.
p.p.s. Боты на серверах, предназначенных для людей, играющие в игры, которые предназначены для людей очень сильно мешают людям. Тем не менее, проверять, насколько сильно они мешают не нужно!
Блог 529
Project Euler и остальное
Мой первый пост или зачем этот блог?
Есть несколько причин появления этого блога. Самая главная — хочется как-то использовать домен, который я купил изначально для красивой почты. Но она не единственная. Мне также хочется делиться частью того, что я делаю на платформе, которую я могу контролировать. Ну и контактные данные. Хорошо, когда все в одном месте.
Пока я планирую публиковать свои, возможно не лучшие, но рабочие решения для задач из Project Euler (projecteuler.net). На самом деле, решения значительного числа задач уже есть в Сети, но на английском. Хотя, это и не есть большая проблема.
Блог 529
Project Euler и остальное
Моё решение задачи 134
Назовём порождающим для двух последовательных простых \(p_1 < p_2\) наименьшее натуральное число, что оно закачивается на \(p_1\) и при этом делится на \(p_2\). Необходимо найти сумму порождающих для всех \(p_1 \in \left[ 5; 10^6 \right]\)
Например, если \(p_1 = 19\), то следующее простое \(p_2 = 23\). Тогда порождающим будет число \(1219\), при этом \(1219 \: \vdots \: 23\).
Полное условие можно найти тут
Несмотря на то, что сложность задачи 45%, для её решения достаточно выписать условие.
Пусть \(p_1\) содержит в себе \(k\) цифр, т.е. \(n = r \cdot 10^k + p_1\), где \(r\) — какое-то натуральное число с отрезка \(\left[ 1; p_2-1 \right]\)
Давайте посчитаем остатки по модулю \(p_2\): \(n \equiv r \cdot 10^k + p_1 \equiv 0\). Отсюда получим явную формулу для \(r\):
$$ r \equiv -p_1 \cdot 10^{-k} \equiv -p_1 \cdot 10^{p_2 -1-k} $$Комментарии:
- Так как \(a^p \equiv a \mod p\), то верно что \(a^{-k} \equiv a^{p -1-k} \mod p\)
- Это всё бессмысленно, если не знать про алгоритм быстрого возведения в степень, который делает асимптотическую сложность возведения в степень логарифмической.
У нас есть явная формула для порождающего, и мы знаем как её быстро посчитать. Ниже приведён код на Python с использованием sympy.
from sympy import primerange # для получения простых чисел - -# быстрое возведение в степень по модулю -def fast_pow(x, y, modulo): - if y == 0: - return 1 - p = fast_pow(x, y // 2, modulo) - p = (p * p) % modulo - if y % 2: - p = (p * x) % modulo - return p - -# нам нужно первое простое, которое больше 10^6 -- 10^6+3 -primes = list(primerange(5,10**6+4)) - -sm = 0 - -for i in range(len(primes) - 1): - digs = len(str(primes[i])) # количество цифр - r = (primes[i+1]**2 - primes[i] * fast_pow(10, primes[i+1] - 1 - digs, primes[i+1])) % primes[i+1] - sm += r * 10**digs + primes[i] - -print('Result is {}'.format(sm)) -
Ответ: 18613426663617118
\ No newline at end of file + \ No newline at end of file diff --git a/posts/moio-reshenie-zadachi-134/index.html.gz b/posts/moio-reshenie-zadachi-134/index.html.gz index a8150b0..909e82e 100644 Binary files a/posts/moio-reshenie-zadachi-134/index.html.gz and b/posts/moio-reshenie-zadachi-134/index.html.gz differ diff --git a/posts/moio-reshenie-zadachi-146/index.html b/posts/moio-reshenie-zadachi-146/index.html index d50b4f7..0872d66 100644 --- a/posts/moio-reshenie-zadachi-146/index.html +++ b/posts/moio-reshenie-zadachi-146/index.html @@ -1,112 +1,112 @@ -
Блог 529
Project Euler и остальное
Моё решение задачи 146
Необходимо найти сумму всех натуральных \(n\), что \(n^2+1\), \(n^2+3\), \(n^2+7\), \(n^2+9\), \(n^2+13\), и \(n^2+27\) будут последовательными простыми числами.
Полное условие можно найти тут
Хочется отметить, что сложность у задачи 50%, а на текущий момент её решило меньше 4000 человек. Тем не менее, мне она показалось простой. Простейшее решение отработало очень быстро.
Для начала, можно отметить, что в лоб проверять условие очень долго. Проверять на простоту числа порядка \(10^{15}\) достаточно сложно, поэтому их нужно как-то отсеять.
Самое простое — не рассматривать те \(n\), что хотя бы одно из \(n^2+1\), \(n^2+3\), \(n^2+7\), \(n^2+9\), \(n^2+13\), и \(n^2+27\) будет заведомо делиться на какое-то маленькое простое число. Это даёт достаточно хорошие результаты: из 150 миллионов чисел, после отсеивания по простым числам \(< 3000\) (этот параметр я подбирал уже после решения задач: если он слишком маленький, то будет слишком много проверок на простоту, если же слишком большой, то мы делаем слишком много работы, чтобы отсеять несколько чисел), останется меньше \(2000\) чисел. Их уже можно проверить непосредственно.
Тогда алгоритм может быть таким:
- Находим простые числа меньше \(3000\).
- Для каждого из них находим допустимые остатки.
- Для каждого из чисел от \(1\) до \(n\) проверяем, что остатки по всем простым хорошие.
- Непосредственно проверяем условие. Важно не забыть проверить непростоту оставшихся нечётных чисел из диапазона \(n^2 + 1 \ldots n^2 + 27\) там могут быть (и будут!) другие простые числа.
Непосредственно сам поиск такой клики можно реализовать тривиально. Ниже мой код на C++11 с использованием библиотек Flint и primesieve. Распараллеливание хоть и просится, но смысла не имеет, т.к. я получил ответ менее, чем за 5 секунд.
/* - * Problem 146 on Project Euler - * Aleksey Lobanov (c) 2016 - */ - -#include <iostream> -#include <vector> -#include <cstdint> -#include <set> -#include <iomanip> -#include <algorithm> - -#include "fmpzxx.h" -#include "arithxx.h" - -#include "primesieve.hpp" - -using namespace std; -using namespace flint; - -bool is_prime(int64_t num) { - fmpz_factorxx fact; - fact.set_factor(num); - if ( fact.size() != 1 ) - return false; - if ( fact.exp(0) != 1) - return false; - return true; -} - -bool is_possible(int64_t num, const vector<int64_t> &to_add) { - for (auto &&add: to_add) - if ( !is_prime(num*num + add) ) - return false; - - // primes must be consecutive - // so we need check, that other numbers like n^2 + i is not primes - vector<int64_t> other_adds; - for (size_t i = to_add[0] + 2; i < to_add[to_add.size() - 1]; i += 2) - if ( !binary_search(to_add.begin(), to_add.end(), i) ) - other_adds.push_back(i); - - for (auto &&add: other_adds) - if ( is_prime(num*num + add) ) - return false; - - return true; -} - -int main() { - const vector<int64_t> to_add = {1, 3, 7, 9, 13, 27}; - - const int64_t MAX_N = 1l*150*1000*1000; - const int64_t MAX_PRIME = 3000; - - vector<int64_t> sieve_primes; - primesieve::generate_primes(MAX_PRIME, &sieve_primes); - - vector< vector <int64_t> > good_remainders; - for (auto &&prime: sieve_primes) { - set<int64_t> remainders; - for(int64_t rem = 0; rem < prime; ++rem) - remainders.insert(rem); - - set<int64_t> base_remainders; - for (auto &&base: to_add) - base_remainders.insert(base % prime); - - for (int64_t rem = 0; rem < prime; ++rem) - for (auto &&base_rem: base_remainders) - if ( (rem*rem + base_rem) % prime == 0 ) - remainders.erase(rem); - - good_remainders.push_back(vector<int64_t>(remainders.begin(), remainders.end())); - } - - size_t cnt = 0; - size_t sum = 0; - - /* WARNING - * for small n can be that - * n^2 + [1 or 3 or .. ] is prime in sieve primes - * but there is only one n < 315410 is 10, - * so we need add 10 to n - */ - sum += 10; - - for (int64_t i = 1; i < MAX_N; ++i) { - bool is_good = true; - for (size_t prime_ind = 0; prime_ind < sieve_primes.size(); ++prime_ind) { - if ( !binary_search(good_remainders[prime_ind].begin(), good_remainders[prime_ind].end(), i % sieve_primes[prime_ind]) ) { - is_good = false; - break; - } - } - if ( is_good ) { - cnt++; - if ( is_possible(i, to_add) ) { - sum += i; - } - } - } - cout << "count = " << cnt << endl; - cout << "Result is: " << sum << endl; - return 0; -} -
Ответ: 676333270
\ No newline at end of file + \ No newline at end of file diff --git a/posts/moio-reshenie-zadachi-146/index.html.gz b/posts/moio-reshenie-zadachi-146/index.html.gz index 405fb3a..2db1e65 100644 Binary files a/posts/moio-reshenie-zadachi-146/index.html.gz and b/posts/moio-reshenie-zadachi-146/index.html.gz differ diff --git a/posts/moio-reshenie-zadachi-60/index.html b/posts/moio-reshenie-zadachi-60/index.html index d497ca8..32e81d4 100644 --- a/posts/moio-reshenie-zadachi-60/index.html +++ b/posts/moio-reshenie-zadachi-60/index.html @@ -1,171 +1,171 @@ -
Блог 529
Project Euler и остальное
Моё решение задачи 60
Необходимо найти множество из пяти простых чисел с минимальной суммой такое, что после “склеивания” в любом порядке любых двух чисел из него тоже будет простое число. Здесь под процедурой “склеивания” чисел \(a\) и \(b\) подразумевается получения из \(a = \overline{a_1 a_2 \ldots a_n}\) и \(b = \overline{b_1 b_2 \ldots b_m}\) некоторого \(c\) так, что \(c = \overline{a_1 a_2 \ldots a_n b_1 b_2 \ldots b_m}\).
Полное условие можно найти тут
Для начала, можно понять, что непосредственный перебор “в лоб” слишком медленный и нужного результата не даст. Поэтому хочется уйти от, как мне кажется, не самого формализуемого условия к чему-то более простого, с чем проще работать. Давайте сначала поймём, какие вообще числа могут быть в одном множестве. Для этого достаточно перебрать все разбиения на два подчисла всех простых чисел. Это достаточно быстро, порядка \(O\left( N \right)\) операций. Важно не забыть, что мы можем разбивать число \(p\) только на \(\overline{p_1 p_2}\), между числами не может быть нулей! То есть если число 37 разбивается на 3 и 7, то 307 нет.
Пусть мы получили набор таких разбиений, то есть набор пар вида \(\left( p, q \right)\). Давайте составим из них граф, где вершинки это простые числа, а ориентированное ребро из \(p\) в \(q\) означает, что есть пара \(\left( p, q \right)\). Из того, что порядок склеивания чисел произвольный сразу следует, что рассматриваемый граф должен быть неориентированным. Таким образом, все пары \(\left( p, q \right)\) для которых нет пары \(\left( q, p \right)\) необходимо выкинуть, а из оставшихся построить граф.
Теперь задача стало гораздо понятнее: достаточно выбрать клику размера 5, что сумма значений её вершин минимальна. В общем случае, это достаточно ресурсоёмкая(как мне кажется) задача, но в реальном графе количество рёбер не слишком большое. В худшем случае, по теореме Турана, количество рёбер в графе лишь с одной такой кликой примерно на 10% меньше числа рёбер в полном графе.
Непосредственно сам поиск такой клики можно реализовать тривиально. Ниже мой код на C++11 с использованием библиотеки Boost Graph Library (BGL).
#include <iostream> -#include <algorithm> -#include <vector> -#include <set> -#include <map> -#include <cstdint> - -#include <boost/graph/adjacency_list.hpp> -#include "primesieve.hpp" - -using namespace std; -using namespace boost; - -typedef adjacency_list<listS, vecS, undirectedS> Graph; - -template <class T> -vector< T > getDigs(T num) { - vector< T > res; - while ( num != 0 ) { - res.push_back(num % 10); - num = num / 10; - } - reverse(res.begin(),res.end()); - return res; -} - -template <class T> -T getNum(const vector<T> &digs){ - T res(0); - T cur_mult(1); - for (auto it = digs.rbegin(); it != digs.rend(); ++it){ - res += *it * cur_mult; - cur_mult *= T(10); - } - return res; -} - -template <class T> -vector< pair<T,T> > getPrimePairs(T prime, set<T> &primes) { - vector< pair<T,T> > res; - auto digs = getDigs<T>(prime); - for (size_t i = 1; i <= digs.size() - 1; ++i){ - T q1 = getNum(vector<T>(digs.begin(), digs.begin() + i)); - T q2 = getNum(vector<T>(digs.begin() + i, digs.end())); - - if ( (primes.find(q1) != primes.end()) - && (primes.find(q2) != primes.end()) && (digs.at(i) != 0) ){ - res.push_back(make_pair(q1,q2)); - } - } - return res; -} - -// removes all (p,q) that (p,q) in edges but (q,p) is not -template <class T> -void delBadEdges(vector< pair<T,T> > &edges){ - map< int64_t, int > mults; // contains counts of p * q - - for (auto p: edges) { - int64_t key = p.first * p.second; - if ( mults.find(key) != mults.end() ) - mults[key] += 1; - else - mults[key] = 1; - } - vector< pair<T,T> > res; - for (auto p: edges) { - int64_t key = p.first * p.second; - if ( mults[key] == 2 ){ - res.push_back(p); - // edge (p,q) is already exists - mults[key] = 0; - } - } - swap(res, edges); -} - -template <class T> -set<T> getUniqNums(vector< pair<T,T> > &edges){ - set<T> res; - for (auto p: edges) { - res.insert(p.first); - res.insert(p.second); - } - return res; -} - -template <class T> -void getNeedSet(Graph &g, Graph::vertex_descriptor st,int depth, - set<Graph::vertex_descriptor> cur_set, const map<Graph::vertex_descriptor, T> &vert_prime, - vector< set<Graph::vertex_descriptor> > &v_out){ - auto edges = out_edges(st, g); - vector< Graph::vertex_descriptor > adj_verts; - for (auto e_it = edges.first; e_it != edges.second; ++e_it) - adj_verts.push_back(target(*e_it,g)); - - for (auto prev_v: cur_set) - if ( find(adj_verts.begin(), adj_verts.end(), prev_v) == adj_verts.end() ) - return; - - cur_set.insert(st); - - if ( depth == 1 ){ - v_out.push_back(cur_set); - return; - } - - for (auto v: adj_verts){ - if ( (vert_prime.at(v) < vert_prime.at(st)) ) - continue; - getNeedSet(g,v,depth - 1, cur_set, vert_prime, v_out); - } -} - -int main(int argc, char **argv) { - const int MAX_VAL = 100000000; - vector< long > v_primes; - primesieve::generate_primes(MAX_VAL, &v_primes); - set< long > s_primes(v_primes.begin(), v_primes.end()); - - cout << "Number of primes is " << v_primes.size() << endl; - v_primes.clear(); - - vector< pair< long,long > > all_edges; - for (auto pr: s_primes){ - auto t_v = getPrimePairs(pr, s_primes); - copy(t_v.begin(),t_v.end(),back_inserter(all_edges)); - } - cout << "There are " << all_edges.size() << " edges" << endl; - delBadEdges(all_edges); - cout << "After deleting bad edges: " << all_edges.size() << " edges" << endl; - // We only need this primes - s_primes = getUniqNums(all_edges); - map< Graph::vertex_descriptor, long > vert_prime; - map< long, Graph::vertex_descriptor > prime_vert; - Graph g; - for (auto pr: s_primes){ - vert_prime[add_vertex(g)] = pr; - } - for (auto v: vert_prime){ - prime_vert[v.second] = v.first; - } - for (auto p: all_edges) - add_edge(prime_vert.at(p.first),prime_vert.at(p.second),g); - auto all_vertices = vertices(g); - vector< set<Graph::vertex_descriptor> > goods; - for (auto v_it = all_vertices.first; v_it != all_vertices.second; ++v_it){ - getNeedSet(g, *v_it,5, set<Graph::vertex_descriptor>(), vert_prime, goods); - } - vector< pair<size_t, size_t> > good_sums; - for (size_t i = 0; i < goods.size(); ++i){ - cout << "Good set #" << i+1 << ": "; - size_t cur_sum = 0; - for (auto el: goods[i]){ - el = vert_prime.at(el); - cur_sum += el; - cout << el << " "; - } - cout << endl; - good_sums.push_back(make_pair(cur_sum, i)); - } - if ( !good_sums.empty() ) - cout << "Result is " << min_element(good_sums.begin(),good_sums.end(), [](pair<size_t, size_t> a, pair<size_t, size_t> b){return a.first < b.first;})->first << endl; - return 0; -} -
Ответ: 26033
UPD. Кажется,
getNeedSetможно было реализовать гораздо проще, использую BGL посетителей (Visitors). Но, к сожалению, на этапе написания кода я про это забыл.
\ No newline at end of file + \ No newline at end of file diff --git a/posts/moio-reshenie-zadachi-60/index.html.gz b/posts/moio-reshenie-zadachi-60/index.html.gz index 2087f3c..6c58ec1 100644 Binary files a/posts/moio-reshenie-zadachi-60/index.html.gz and b/posts/moio-reshenie-zadachi-60/index.html.gz differ diff --git a/posts/nakhozhdenie-summy-k-ykh-stepenei/index.html b/posts/nakhozhdenie-summy-k-ykh-stepenei/index.html index 60cee34..dcb6d5d 100644 --- a/posts/nakhozhdenie-summy-k-ykh-stepenei/index.html +++ b/posts/nakhozhdenie-summy-k-ykh-stepenei/index.html @@ -1,6 +1,6 @@ -
Блог 529
Project Euler и остальное
Нахождение суммы k-ых степеней
Давайте сразу обобщим задачу до нахождения \(f_k\left( n \right)\), где
$$f_k\left( n \right) = 1^k + 2^k + \ldots + n^k$$Для \(k=1\) формула известна всем школьникам: \(f_1\left( n \right) = \frac{n\left(n+1 \right)}{2}\). Формулу для \(k=2\) знают уже не все, но всё же в школе её найти можно (я видел на обложке учебника по алгебре): \(f_2\left( n \right) = \frac{n\left(n+1 \right) \left( 2n + 1 \right)}{6}\)
Интуиция может подсказать, что \(f_k \left( n \right)\) есть некий полином со степенью \(k+1\). Если это так, то его нахождение тривиально. Например, можно посчитать его в явном виде, используя полином Лагранжа. Осталось показать, что наша функция представима в таком виде.
Для начала введём обозначение. “Нижней степенью”, \(x^{\underline{k}}\), будем обозначать такое выражение:
$$x^{\underline{k}} = x(x-1)\cdot \ldots \cdot (x-k+1)$$.
Далее, заметим следующее, если \(a_i = A_{i+1} - A_i\), где \(\lbrace a_i \rbrace\) и \(\lbrace A_i \rbrace\) — некие последовательности, то \(\sum_{i=1}^n = A_{n+1}-A_1\) (телескопирование, можно посмотреть тут, с. 6).
Теперь посчитаем сумму \(\sum_{i=1}^n i^{\underline{k}}\). Для этого достаточно понять, что \(\left( x+1 \right)^{\underline{k+1}} - \left(x \right)^{\underline{k+1}} = \left( k+ 1 \right)x^{\underline{k}}\). Отсюда сразу получаем, что
$$\sum_{i=1}^n i^{\underline{k}} = \frac{\left( n+1 \right)^{\underline{k+1}}}{k+1}$$Осталось показать, что “нормальные” степени выражаются через нижние. Начнём со степени \(k=1\), тут всё просто:
$$x = x^{\underline{1}}$$С бОльшими степенями сделаем следующее: считая, что все степени, меньше, чем \(k\) мы выражать умеем, раскроем скобки в определении нижней степени. Теперь поймём, что старший коэффициент \(1\): \(x^{\underline{k}} = x^k + \sum_{i=1}^k a_ix^i\) или \(x^k = \sum_{i=1}^k a_ix^i - x^{\underline{k}}\). Осталось понять, что каждое из слагаемых вида \(a_ix^i\) мы умеем выражать через нижние степени. Таким образом, можно получить следующее:
$$ \sum_{i=1}^{n} i^k = \sum_{i=1}^{n} \sum_{j=1}^{k} a_j i^{\underline{k}} = \sum_{j=1}^{k} \sum_{i=1}^{n} a_j i^{\underline{k}} = \sum_{j=1}^{k} \frac{a_j \left(n+1 \right)^{\underline{k+1}}}{k+1}$$Кстати, формула для суммы в самом начале такая:
$$ \sum_{i=1}^n i^5 = \frac{1}{12} n^2 \left(n+1 \right)^2 \left(2n^2 + 2n-1 \right) $$
\ No newline at end of file + \ No newline at end of file diff --git a/posts/nakhozhdenie-summy-k-ykh-stepenei/index.html.gz b/posts/nakhozhdenie-summy-k-ykh-stepenei/index.html.gz index ae77b1d..d3fe57f 100644 Binary files a/posts/nakhozhdenie-summy-k-ykh-stepenei/index.html.gz and b/posts/nakhozhdenie-summy-k-ykh-stepenei/index.html.gz differ diff --git a/posts/wallabag-i-realnaia-zhizn/index.html b/posts/wallabag-i-realnaia-zhizn/index.html index 47e705f..c048b1a 100644 --- a/posts/wallabag-i-realnaia-zhizn/index.html +++ b/posts/wallabag-i-realnaia-zhizn/index.html @@ -1,6 +1,6 @@ -
Блог 529
Project Euler и остальное
Wallabag и реальная жизнь
Начать следует с того, что Wallabag действительно является самым популярным среди открытых приложений для отложенного чтения. Можно взять, например, alternativeto:
Первый релиз вышел почти два года назад. Тем не менее, мне сложно назвать продукт зрелым. Последняя, на момент написания, версия
2.0.0-beta.2не может похвастаться простым процессом установки. Вариант просто выполнить команды из мануалов по очереди у меня не получился. В этом соперничать с тем же Pocket, очевидно, бессмысленно.Стандартная тема, material, ужасно выглядит на моём ноутбуке с разрешением 1366x768, элементы явно рассчитаны на большую диагональ. Ещё часть места отъедает неубирающаяся плашка внизу страницы, предупреждающая о том, что баги в бета версии не есть что-то плохое. Официальное Android приложение упорно не может найти сервер.
После волевого решения перейти на стабильную версию (то есть откатиться в равзвитии на полгода назад), дела улучшились, но не сильно. Количество настроек минимально, если не сказать, что их вообще нет. Но оно работает, вроде.
Сложности, впрочем, только начались. Обещанная синхронизация с Pocket работает, мягко говоря, неоптимально. После загрузки экспортированного html файла со ссылками, wallabag почти час выкачивал мои 750 статей. Где-то на 500-й статье он выделил слишком много памяти и упал с ошибкой. После этого пришлось руками искать потенциально проблемные статьи и удалять их по одной — удалить сразу несколько элементов невозможно. Затем надо было выкачивать ещё 250 статей, периодически посматривая, чтобы ничего снова не упало.
Кажется, что все проблемы закончились, но нет. Дальше синхронизация с телефоном. Она заняла почти столько же времени, при этом, начальное очевидное предположение о том, что токен безопасности вбивать не нужно (он сам заполняется в приложении) стоило где-то 15 минут поиска. Но и это не всё. После того, как база загрузилась, небольшие изменения на телефоне (такие как удаление статьи) синхронизировались больше минуты!
Резюме: пользоваться можно, но советовать кому-либо это использовать я точно не стану.
Этого бы поста не было бы, если бы я сказал, что Wallabag плохой, а OTHER_PRODUCT хороший и можно пользоваться им. Но я так написать не могу. Ни я, ни alternativeto других решений не знают. Значит нужно их создать. На моём слабеньком VPS уже почти год трудится Syncthing, управляя значительным количеством файлов с минимальной нагрузкой на ЦП. Поэтому мне кажется, что Go подойдёт идеально.
Если написать подобный продукт на Go, то многие проблемы даже не появятся:
- не нужно разрешать большое количество зависимостей — достаточно одного бинарника
- скорость генерации контента (например, создание pdf, epub) будет значительно выше
- небольшой оверхед позволит всё хранить в памяти
- работа почти на чём угодно
В совокупности, будет достигнуто, как мне кажется, самое главное — удобство для обычного пользователя, который сможет без проблем развернуть это на любом VPS.
Блог 529
Project Euler и остальное
-
Как я делаю бэкапы почты
Решение, которое я сейчас исользую для бэкапа своей почты
Page 1 / 1
Блог 529
Project Euler и остальное
Моё решение задачи 60
Краткое условие: необходимо найти множество из пяти простых чисел с минимальной суммой такое, что после “склеивания” в любом порядке любых двух чисел из него тоже будет простое число.
Page 1 / 1
Блог 529
Project Euler и остальное
Мой первый пост или зачем этот блог?
Кратко о том, почему я сделал блог.
Page 1 / 1
Блог 529
Project Euler и остальное
Как я шахматного бота писал
Как я проверял шахматное приложение на “ботоустойчивость”.
Page 1 / 1
Блог 529
Project Euler и остальное
Моё решение задачи 146
Краткое условие: необходимо найти сумму всех натуральных \(n\), что \(n^2+1\), \(n^2+3\), \(n^2+7\), \(n^2+9\), \(n^2+13\), и \(n^2+27\) будут последовательными простыми числами.
CrossGen v1.0
Читая хабр, случайно натолкнулся на идею сделать программу, которая по заданной кроссвордной сетке находит способ её заполнить. В этом посте вкратце напишу про моё решение и первую версию приложения.
Моё решение задачи 60
Краткое условие: необходимо найти множество из пяти простых чисел с минимальной суммой такое, что после “склеивания” в любом порядке любых двух чисел из него тоже будет простое число.
Page 1 / 1
CrossGen v1.0
Читая хабр, случайно натолкнулся на идею сделать программу, которая по заданной кроссвордной сетке находит способ её заполнить. В этом посте вкратце напишу про моё решение и первую версию приложения.
Моё решение задачи 60
Краткое условие: необходимо найти множество из пяти простых чисел с минимальной суммой такое, что после “склеивания” в любом порядке любых двух чисел из него тоже будет простое число.
Page 1 / 1
\ No newline at end of file
+
\ No newline at end of file
diff --git a/tag/c.html.gz b/tag/c.html.gz
index 0712fac..b9db8f2 100644
Binary files a/tag/c.html.gz and b/tag/c.html.gz differ
diff --git a/tag/flint.html b/tag/flint.html
index 9807932..a249079 100644
--- a/tag/flint.html
+++ b/tag/flint.html
@@ -1,6 +1,6 @@
- Блог 529
Project Euler и остальное
Моё решение задачи 146
Краткое условие: необходимо найти сумму всех натуральных \(n\), что \(n^2+1\), \(n^2+3\), \(n^2+7\), \(n^2+9\), \(n^2+13\), и \(n^2+27\) будут последовательными простыми числами.
Page 1 / 1
Page 1 / 1
\ No newline at end of file
+
\ No newline at end of file
diff --git a/tag/flint.html.gz b/tag/flint.html.gz
index aa915db..395eeb7 100644
Binary files a/tag/flint.html.gz and b/tag/flint.html.gz differ
diff --git a/tag/go.html b/tag/go.html
index d9d30be..ad171a4 100644
--- a/tag/go.html
+++ b/tag/go.html
@@ -1,6 +1,6 @@
- Блог 529
Project Euler и остальное
Wallabag и реальная жизнь
Как я устанавливал известнейшее свободное решение для отложенного чтения, и что из этого получилось.
Page 1 / 1
Блог 529
Project Euler и остальное
Экспорт партий с Lichess
Сохраняем себе шахматные партии с Lichess.org. Рабочий скрипт на Python прилагается.
Page 1 / 1
Блог 529
Project Euler и остальное
Нахождение суммы k-ых степеней
Как придумать формулу для суммы \(1^5 + 2^5 + 3^5 + \ldots + n^5\) и есть ли она вообще?
Page 1 / 1
Page 1 / 1
\ No newline at end of file
+
\ No newline at end of file
diff --git a/tag/matematika.html.gz b/tag/matematika.html.gz
index a00e53d..3e66477 100644
Binary files a/tag/matematika.html.gz and b/tag/matematika.html.gz differ
diff --git a/tag/open-source.html b/tag/open-source.html
index 67f60a4..d09b244 100644
--- a/tag/open-source.html
+++ b/tag/open-source.html
@@ -1,6 +1,6 @@
- Блог 529
Project Euler и остальное
Wallabag и реальная жизнь
Как я устанавливал известнейшее свободное решение для отложенного чтения, и что из этого получилось.
Page 1 / 1
Блог 529
Project Euler и остальное
Экспорт партий с Lichess
Сохраняем себе шахматные партии с Lichess.org. Рабочий скрипт на Python прилагается.
Как я шахматного бота писал
Как я проверял шахматное приложение на “ботоустойчивость”.
CrossGen v1.0
Читая хабр, случайно натолкнулся на идею сделать программу, которая по заданной кроссвордной сетке находит способ её заполнить. В этом посте вкратце напишу про моё решение и первую версию приложения.
Ещё одно вычисление выражений
На хабре когда-то увидел статью про то, что в Яндексе двум сотрудникам дали задачу на написание приложения, для вычисления выражений. Менеджер справился за 4 часа, а программист за два. Я решил попробовать свои силы.
Page 1 / 1
Блог 529
Project Euler и остальное
Моё решение задачи 134
Краткое условие: назовём порождающим для двух последовательных простых \(p_1 < p_2\) наименьшее натуральное число, что оно закачивается на \(p_1\) и при этом делится на \(p_2\). Необходимо найти сумму порождающих для всех \(p_1 \in \left[ 5; 10^6 \right]\)
-
Моё решение задачи 146
Краткое условие: необходимо найти сумму всех натуральных \(n\), что \(n^2+1\), \(n^2+3\), \(n^2+7\), \(n^2+9\), \(n^2+13\), и \(n^2+27\) будут последовательными простыми числами.
Моё решение задачи 146
Краткое условие: необходимо найти сумму всех натуральных \(n\), что \(n^2+1\), \(n^2+3\), \(n^2+7\), \(n^2+9\), \(n^2+13\), и \(n^2+27\) будут последовательными простыми числами.
Моё решение задачи 60
Краткое условие: необходимо найти множество из пяти простых чисел с минимальной суммой такое, что после “склеивания” в любом порядке любых двух чисел из него тоже будет простое число.
Page 1 / 1
Моё решение задачи 60
Краткое условие: необходимо найти множество из пяти простых чисел с минимальной суммой такое, что после “склеивания” в любом порядке любых двух чисел из него тоже будет простое число.
Page 1 / 1
\ No newline at end of file
+
\ No newline at end of file
diff --git a/tag/project-euler.html.gz b/tag/project-euler.html.gz
index ddb807e..8be5024 100644
Binary files a/tag/project-euler.html.gz and b/tag/project-euler.html.gz differ
diff --git a/tag/python.html b/tag/python.html
index edcf1df..cc05f21 100644
--- a/tag/python.html
+++ b/tag/python.html
@@ -1,6 +1,6 @@
- Блог 529
Project Euler и остальное
Долгие сообщения в RabbitMQ
Почему стандартная Python библиотека для RabbitMQ не работает с длинными сообщениями и как это исправить.
Моё решение задачи 134
Краткое условие: назовём порождающим для двух последовательных простых \(p_1 < p_2\) наименьшее натуральное число, что оно закачивается на \(p_1\) и при этом делится на \(p_2\). Необходимо найти сумму порождающих для всех \(p_1 \in \left[ 5; 10^6 \right]\)
Page 1 / 1
Page 1 / 1
\ No newline at end of file
+
\ No newline at end of file
diff --git a/tag/python.html.gz b/tag/python.html.gz
index 832eb0f..4dd14b8 100644
Binary files a/tag/python.html.gz and b/tag/python.html.gz differ
diff --git a/tag/rabbitmq.html b/tag/rabbitmq.html
index dedc719..5d01560 100644
--- a/tag/rabbitmq.html
+++ b/tag/rabbitmq.html
@@ -1,6 +1,6 @@
- Блог 529
Project Euler и остальное
Долгие сообщения в RabbitMQ
Почему стандартная Python библиотека для RabbitMQ не работает с длинными сообщениями и как это исправить.
Page 1 / 1
Блог 529
Project Euler и остальное
Экспорт партий с Lichess
Сохраняем себе шахматные партии с Lichess.org. Рабочий скрипт на Python прилагается.
Как я шахматного бота писал
Как я проверял шахматное приложение на “ботоустойчивость”.
Page 1 / 1
Блог 529
Project Euler и остальное
Моё решение задачи 134
Краткое условие: назовём порождающим для двух последовательных простых \(p_1 < p_2\) наименьшее натуральное число, что оно закачивается на \(p_1\) и при этом делится на \(p_2\). Необходимо найти сумму порождающих для всех \(p_1 \in \left[ 5; 10^6 \right]\)
Page 1 / 1
Page 1 / 1
\ No newline at end of file
+
\ No newline at end of file
diff --git a/tag/sympy.html.gz b/tag/sympy.html.gz
index 370782d..473bd25 100644
Binary files a/tag/sympy.html.gz and b/tag/sympy.html.gz differ
diff --git a/tag/wallabag.html b/tag/wallabag.html
index f268508..a2bc2df 100644
--- a/tag/wallabag.html
+++ b/tag/wallabag.html
@@ -1,6 +1,6 @@
- Блог 529
Project Euler и остальное
Wallabag и реальная жизнь
Как я устанавливал известнейшее свободное решение для отложенного чтения, и что из этого получилось.
Page 1 / 1
Блог 529
Project Euler и остальное
CrossGen v1.0
Читая хабр, случайно натолкнулся на идею сделать программу, которая по заданной кроссвордной сетке находит способ её заполнить. В этом посте вкратце напишу про моё решение и первую версию приложения.
Page 1 / 1