Публикации с меткой «icfpc»

Андрей Власовских

Команда spb-archlinux на ICFPC 2009


Никак не возьму в привычку чаще писать в блог. Но на этот раз повод действительно есть. Наша команда spb-archlinux принимала участие в программерском конкурсе ICFP Contest 2009 :)

Ниже рассказ об этом. Если вкратце, было весело! :)

Upd: В финальной таблице результатов мы заняли 95 место с 1826 очками.

Один из членов команды Jkff уже описал наше участие в своём отчёте, так что здесь я остановлюсь на частностях и личном впечатлении, а за общим описанием отсылаю вас в его блог.

Итак, описание Jkff вы уже прочитали. Я бы хотел рассказать про отдельные моменты конкурса. Мне наиболее интересно то, почему мы заняли относительно невысокое место (за 4 часа до окончания мы были 84-ые из 317). Поэтому ниже проводится «разбор полётов». Но я хочу сразу сказать, что мне лично было приятно участвовать с нашей командой! Мы хорошо сработались и я с удовольствием поучаствую в следующем году. Ниже я расскажу также о весёлых моментах в конкурсе :)

Задание

Итак, во время прочтения задания мы сразу и обрадовались, и обеспокоились :) Я, в общем-то, правильно предсказал то, что нас ждёт. До конкурса на каналах Jabber и IRC много говорилось об аэрокосмической теме и о времени 13:00:16. Окончательно раскусить оргов так и не смогли. Мне было ясно следующее: форматы времени и вообще точное время будут играть большую роль (про важность дискретного времени мы вспомнили, к сожалению, слишком поздно). Другой моей догадкой была очередная виртуальная машина (VM), которая обеспечит независимость от платформы: орги слишком уклончиво отвечали на вопрос о том, как надо было отправлять результаты.

В самом начале конкурса выяснилось, что 1995-06-29T13:00:16 (часовой пояс?) — это время первой стыковки шаттла и станции «Мир». Управлять спутниками — это звучало весело :) Но мы сразу стали опасаться сложностей с VM: Jkff вычитал что-то про машинные инструкции, записывающие результат в память по адресу инструкции, и мы стали думать про самомодифицирующийся код и пр. В реальности всё оказалось очень просто. VM была с гарвардской архитектурой: память команд и данных у неё была раздельная. К тому же позднее нам стало понятно, что платформонезависимость была единственной целью включения VM в задание. Способ и эффективность реализации VM на решение почти что не влияли. Мы же напрасно уделили этому внимание, в тактике нашей команды постепенно начала накапливаться ошибка подхода к заданию.

Интерфейс ctl-sim

Почти с самого начала мы стали решать задачу в духе Unix. Вырисовывалась пара взаимодействующих программ: sim — симулятор, выполняемый на VM, и ctl — контроллер, ведущий симулируемый спутник. В ходе небольшого обсуждения я изложил своё видение взаимодействия этой пары. ctl и sim будут обмениваться командами чтения/записи в порты IO через неименованные каналы ОС. Формат обмена FMT было решено сделать текстовым, что давало нам читаемый и отлаживаемый интерфейс в стиле Unix. Например, имея готовую трассу команд на вход симулятора, вывод можно было, в принципе, увидеть картинку полёта спутников следующим образом:

$ cat trace.fmt | sim | vis

В таком же стиле можно было вручную вести спутник (номеров портов было мало и мы помнили их наизусть):

$ sim
0 0x3e80=1002
0 0x0=0 0x1=10000 0x2=-6352278.930 0x3=6361717.566 0x4=21082000
1
1 0x0=0 0x1=10000 0x2=-6347554.360 0x3=6366431.626 0x4=21082000
2 0x2=5000
2 0x0=0 0x1=5000 0x2=-6345326.291 0x3=6371142.177 0x4=21082000

По общему мнению, интерфейс взаимодействия контроллеров и симуляторов получился хороший, он позволил нам писать их на разных языках: Python, Haskell, C++ и Java (на ней был написан визуализатор). Однако сразу же надо высказаться критически об этом интерфейсе. Он способствовал написанию контроллеров на разных языках, что уменьшало шансы совместной работы над кодом и повторного использования кода в следующих задачах. Мы обезопасили себя от проблем с отладкой VM, нашего физического эмулятора и разных контроллеров. Однако это дало начало ненужному в таком коротком конкурсе языковому разброду. К тому же отладка интерфейса обошлась нам в 4 человеко-часа. Но в целом, всё же, интерфейс был положительным моментом.

С упомянутой отладкой интерфейса отдельная история. Несмотря на то, что программы взаимодействовали через каналы (pipes), они были крайне сильно связаны. Контроллер ожидал ответа от симулятора после каждой своей команды. Фактически, это был протокол с обменом сообщениями «запрос-ответ». А для этого только каналов недостаточно, нужно организовать подчиенение процессов. Эта типичная схема описана у Реймонда в “The Art of Unix Programming”. Мы решили оставить всё на каналах и просто сделать программу-драйвер fmtdump, которая связывала каналы двух программ крест-на-крест и заодно писала логи их взаимодействия в общем формате FMT. Всё бы хорошо, да только мы не задумывались над обязательностью flush(2) после write(2) и над тем, что строки текствого протокола мы не имеем права буферизовать при чтении: строку неизвестной длины без блокировки можно читать только по одному байту за раз.

После того, как минут за 30 всё это было осознано, началась борьба по отключению буферизации ввода во всех языках, на которых писались контроллеры. В Python это оказалось сделать на удивление сложно. Способ, описанный в документации, не работал и одно время мне пришлось использовать свою функцию чтения строки по байту. Позже правильный способ всё-таки был найден: открыть файл в бинарном режиме, указать размер буфера 0 и читать строку через readline, а не итерироваться по строкам файла. В других вариантах буферизация всё равно происходит.

Вот с этими техническими деталями VM и интерфейса симулятора к контроллеру мы провозились до середины субботы. Сама VM получилась тоже интересно, см. блог Jkff. Напомню, что для получения симулятора использовалось взаимодействие языков байткода VM, Python, Haskell и C++.

Все пишут контроллеры

Закопавшись с VM, мы сели за саму задачу, т. е. за контроллеры, только к полудню субботы. До этого я писал транслятор байткода в Haskell-ассемблер VM, а также драйвер пары контроллер-симулятор fmtdump, Kerzum писал компилятор sim с Haskell-ассемблера в C++, Jkff писал визуализацию выходного лога в формате FMT, а Cfr и Incubos писали альтернативный sim, моделирующий мир по формулам из задания.

Cfr уже без Incubos (который пошёл спать) начал писать свой контроллер на Haskell, однако скоро от его кода ответвился Kerzum, который решил попробовать свою реализацию. Параллельно Jkff начал ещё один контроллер на Python. Я занимался скриптами сборки всех проектов и небольшими исправлениями к fmtdump. Не очень понятно, что мешало всем вместе вначале решить задачу математически и алгоритмически. Наверное, всем просто хотелось уже наконец попробовать подойти к задачам практически: симулятор и драйвер были как раз готовы. Начавшееся разделение хаскеловских и питоновских контроллеров в результате затормозит нас на довольно большое время, хотя и обезопасит от отдельных неудач с конкретными контроллерами. Как и в случае с интерфейсом контроллер-симулятор, мы опять подошли к проблеме слишком защищёно.

Вечером в субботу свой контроллер начал и Incubos, взяв за основу мои наброски (об этом речь ниже). Так что в воскресенье практически все работали над контроллерами, но своими, а не общим. Знания тоже локализовались в коде и обмен ими шёл не особо эффективно.

Эволюция контроллеров на Python

После того, как к полудню субботы я решил вопросы с драйвером, я взялся за контроллеры. Однако вместо того, чтобы познакомиться наконец детально с физикой задачи (а не только с VM и форматами), я решил придумать хороший интерфейс контроллеров на Python, чтобы их было легко писать и комбинировать.

Утром я написал для ребят простой парсер формата FMT на Python и дёрнул меня тогда чёрт сделать его интерфейс потоковым! Возможно, это была главная моя ошибка за весь конкурс. Контроллер представлялся в этом первом интерфейсе функцией типа Iterable(Ports) -> Iterable(Ports). Чуть позже, думая об универсальном интерфейсе, я написал вспомогательную функцию buffered :: ([Ports] -> Ports), int -> (Iterable(Ports) -> Iterable(Ports)), которая позволяла работать с одним шагом симуляции, но иметь буфер портов за предыдущие ходы. Тогда мне показалось, что такая потоковая обработка наряду с буферизацией будет удобной для решения всех задач. Но я почему-то забыл с самого начала подумать про суть обмена данными с симулятором и про композицию алгоритмов управления в контроллерах. С таким интерфейсом это получалось крайне неестественно.

Что же было естественным способом управления симулятором? Старое доброе ООП. Наша задача сводится к синхронному обмену сообщениями (message passing) между процессами с локальным изменяемым состоянием.

К сожалению, Jkff и Incubos уже вовсю писали свой код на потоковом интерфейсе (код Incubos особенно критично зависел от этого), а ОО-интерфейс я придумал только днём в воскресенье, после того, как, наконец, выспался после 36-часового бодрствования.

Новый интерфейс заключался вот в чём. Есть интерфейс ООП Sim с единственными методом: interact :: Ports -> Ports, который в ответ на управляющие сигналы в порты ввода выдаёт сигналы портов вывода. Всё очень банально! Зато сразу стало возможно делать следующие вещи: объединять контроллеры (процедуры аргумента Sim) последовательной композицией, организовывать внутри них конечные автоматы, вызывать внутри контроллеров другие контроллеры (вложенные конечные автоматы), делать классы-прокси для дампа логов взаимодействия, запускать альтернативный симулятор внутри процесса для предсказания физики и пр.

Посетила бы нас эта простая мысль раньше — наш код был бы повторно используемым, более выразительным и наверняка содержал бы меньше ошибок.

Последнее улучшение интерфейса контроллера состояло в том, что специальный прокси TraceSim добавлял поле trace :: [(Ports, Ports)] с историей команд ввода-вывода, что нужно почти всем алгоритмам.

Физика и геометрия

Отдельной темой стоят наши обсуждения физики и геометрии задачи. К сожалению, в субботу их было очень мало, но зато они были интересные :) и в воскресенье они позволили нам прийти, фактически, к решению задачи. Особенно эффективными мне показались обсуждения физики у доски, в которых участвовали я, Jkff и пару раз Cfr с Kerzum. На мой взгляд, примерно так и стоило бы решать задачу: обсуждать и вместе писать алгоритмы, а не распределять усилия на несвязанный код. Хоть в итоге именно распределение позволило набрать нам наши очки, но оно и помешало набрать нам больше.

Надо заметить, что ещё в субботу вечером у нас был код, выводящий спутники первых задач на нужную орбиту, но нам упорно не давали очков. Оказалось, мы невнимательно читали задание, допустили несколько опечаток в расчётах, выводили данные со значительной потерей точности, также совсем забыли, что мы летаем в дискретном времени с шагом 1 сек и скоростью порядка 10 км/сек. Это значит, что без учёта дискретизации (например, параллельным расчётом по непрерывной формуле и дискретной симуляцией с коррекцией после этого) мы не могли достичь нужной точности в орбитальных переходах.

В воскресенье часть этих проблем была решена и при помощи многих подхаков мы решили все задачи, кроме последней, главной. До неё мы, к сожалению, не успели дойти из-за нехватки времени и неизвестного бага в симуляторе, проявлявшегося только на этой задаче.

Причины неуспеха

По моему мнению, нам помешало следующее:

  • Не пытались писать алгоритмы контроллеров вместе
  • Не провели своевременный поиск доков по предметной области
  • Долго не вспоминали о дискретном времени
  • Писали код на разных языках, не думали о композиции алгоритмов
  • Потеряли полдня на баг с выводом неточных чисел с плавающей точкой
  • Плохо организовали математический код, допускали много опечаток

Что мне понравилось

Ну и после такого анализа хочу поделиться приятными впечатлениями :) Вот что было интересно:

  • Мы запрограммили классный многоязыковой процесс получения симулятора
  • Картинка concept map на доске в самом начале позволила быстро найти связи, которые ещё никто не понимал
  • Unix-интерфейсы рулят! Мы игрались с интерфейсом ctl-sim, рисуя бешеные трассы, отлавливая баги во взаимодействии: всё было текстовым и простым
  • Мы выпили в ходе обсуждения в первую ночь кучу кофе и продержались где-то до 20 часов субботы
  • У нас были интересные обсуждения у доски. Было здорово подсказывать, ловить баги друг друга, придумывать новые подходы :)
  • Система управления версиями Mercurial всё время рулила
  • Мы хотели и продумывали (но так и не сделали) интерактивный контроллер в визуализаторе, чтобы играться в спутник с клавиатуры

Цифры

Ну и совсем под конец несколько цифр:

  • 219 ревизий в Mercurial
  • 11 исполняемых программ
  • 1701 строка на Python, 630 на C++, 475 на Java, 460 на Haskell, 171 на C, 111 на make, 67 на sh
  • На момент заморозки таблицы 84 место, 1654 очка (до отправки последних результатов)

Спасибо членам команды за классно проведённое время! :)

Tagged: 2009, c++, haskell, icfpc, java, python, spb-archlinux

Murkt’s codehole

ICFPC'09: Космоопера

  • Имя команды: Concrete mixers
  • Итоговый балл: Weighted Total Score 2852.2285 (14 problems solved)

Начало

Узнал я про проведение ICFPC в течении 27-30 июня внезапно и всего за пару дней. Отправился в конференцию pythonua@cjr, и быстренько нашёл там среди знакомых имён тех, кто собирался участвовать. Собирались мы в секретной комнате на том же conference.jabber.ru :) Это был наш первый опыт участия в ICFPC, а для меня — и вообще в подобных соревнованиях (xa4a и tilarids уже участвовали в Sapka contest). Забегая наперёд, я думаю, что справились мы очень хорошо, даже отлично и должны быть где-то в первых 10% итоговой таблицы. Конечно, теоретически можно было и лучше выступить, но теоретически все могли лучше выступить :)

День 0, вечер

xa4a заранее создал приватный репозитарий на битбакете. Контест начинался в девять часов по киевскому времени, и, конечно же, ровно в девять его сайт практически перестал работать и или отваливался по таймауту, или выдавал страничку в течении пары-тройки минут. Кто-то из первых скачавших спецификацию сразу выложил её на зеркале и показал ссылку в irc-канале #icfp-contest@freenode, поэтому особых проблем с получением задания не было. На сам сайт первым пробился tilarids через несколько минут после начала, и зарегистрировал нашу команду — мы были уже 129-ми.

Суть задания — нужно управлять спутником на околоземных орбитах. Чуть больше часа мы читали задание, вникали в суть и обсуждали как это нужно делать. Если вкратце — у нас есть двумерное пространство (x-y), в котором у нас есть Земля и спутники. Всего четыре типа заданий, и для каждого есть четыре разных набора данных:

  • Перевести спутник с одной круговой орбиты на другую. Задание засчитывается после 900 секунд полёта в пределах одного километра от заданного радиуса.
  • Состыковаться со спутником на другой круговой орбите. Задание засчитывается после 900 секунд полёта в радиусе одного километра от цели.
  • Состыковаться со спутником на эллиптической орбите, причём стартовая орбита тоже может быть эллиптической. Засчитывается так же, как второе.
  • Пролететь в радиусе километра от одиннадцати спутников. Держаться 900 секунд возле каждого не нужно.

Организаторы выдали три бинарника — соответственно пёрвым трём типам задач, которые нужно выполнять на виртуальной машине, и собственно спецификацию к виртуальной машине, форматам файлов, и всему прочему. VM элементарна — четыре арифметические операции, операции сравнения с нулём, копирование данных и ввод-вывод. Четвёртый бинарник пообещали выдать спустя сутки — после окончания lightning round.

Где-то в 22:20 я начал писать парсер бинарников, через двадцать минут он был готов, а ещё минут через пятнадцать в нём были выловлены все допущенные ошибки :) Ещё через час (в полночь) у меня уже был готов базовый вариант виртуальной машины — хотя я писал её первый раз в жизни :) VM запустилась, и смогла выполнить каждый из бинарников не упавши. Ура!

Спецификация написана в некоторых местах очень мутно, и спустя час была уже третья версия — в соответствии с которой я исправил и свою VM. Всё это время мы параллельно обсуждали как управлять спутником, какой должен быть интерфейс общения решателей проблем (надо было их называть мистерами Вулфами ;) ), виртуальной машины и чего-нибудь ещё буде оно появится.

Получилось так — есть порты для входящих данных, и есть для исходящих. Даём виртуальной машине входящие данные, исполняем полностью весь код задачи и получаем исходящие данные. Каждый такой проход — секунда внутреннего времени. В портах ввода мы задаём номер задания (только для первой итерации), а также задающие воздействия на актуаторы — двигатели нашего спутника. Ускорение, которое произойдёт в течении последующей итерации — то есть, разница скоростей между временами t и t+1 задаётся в м/с и указывается вектором — отдельно для координат x и y. Например, мы можем сказать спутнику — «по оси x изменить скорость на 40 м/с, по оси y — на 452,4 м/с».

В выходных портах мы можем прочесть заработанные очки, количество оставшегося топлива, свои координаты и радиус целевой орбиты. Оставшееся топливо меряется просто напросто в м/с — удобно. Заработанные очки равны нулю, пока ничего интересного не произошло, -1 если у нас закончилось топливо или спутник погиб при ударе о Землю, и какому-то определённому числу, если мы выполнили задание. Они зависят от количества использованного топлива и от времени выполнения задания.

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

В общем, ещё спустя минут сорок-пятьдесят усилиями xa4a мы увидели картинку в гнуплоте, а усилиями tilarids у нас был визуализатор на pygame, с анимацией и наблюдением в реальном времени. Рисовало как наш спутник летает вокруг Земли, и толстой зелёной линией целевую орбиту. Красота! :)

В то самое время у некоторых людей из конференции icfpc@cjr VM была в зачаточном состоянии и они ещё даже не пытались её запускать. В этот момент я порадовался что мы продвигаемся достаточно неплохо и отправился спать.

День 1

Я проснулся в половину восьмого утра и пошёл смотреть что же успели накрутить товарищи за то время, пока я спал. Оказалось, что у нас за ночь появился нерабочий логгер трейсов и код, который пытался перевести спутник между орбитами по Гоману… Но спутничек почему-то не долетал пару тысяч километров до нужного радиуса. Тем временем спецификации выросли уже до седьмой версии, да так, что опять пришлось исправлять VM :)

Ещё за час я в несколько раз ускорил визуализатор tilarids’а и поменял константы (гравитационная постоянная, масса Земли) в считалке орибт — там использовались данные из реального мира, а нужны были такие, как в спецификации. Таким образом мы получили решения для трёх задач — 1001, 1002 и 1003. В четвёртой задаче спутник изначально вращался в другую сторону, что создавало некоторые проблемы для того решения, что было, и я решил пока переключиться на второй тип задач.

К полудню начали активизироваться другие участники, и у нас появилось решение для 1004-й, потом появился рабочий логгер — и мы заслали на сайт решения для всех задач первого типа. А к трём часам дня подоспел подарок от tilarids’а — он переписал за полчаса-час виртуальную машину с обычного интерпретирования на генерацию питоновского кода, чем ускорил её раз в десять-двадцать :)

Во второй задаче была такая проблема — если мы сразу переходим на целевую орбиту, то до целевого спутника может быть от шестой части оборота до половины. Выхода нашлось два:

  • Применять специальные методы-фазеры. Если мы опережаем цель, то нужно выходить на внешнюю орбиту, если отстаём от неё — на внутреннюю. Чем ниже орбита, тем быстрее по ней движется спутник и тем меньше длина орбиты, ну и наоборот :)
  • Ждать некоторое время на своей орбите, и в чётко выверенный момент начать переходить на целевую орбиту.

Из-за того, что управление у нас дискретное, а скорости значительные (единицы километров в секунду), то в некоторых случаях погрешности на выходе составляли не сильно удобоваримые единицы километров — а нам нужно находиться ближе километра! При помощи таких-то слов и ручной подстройки времени отлёта (плюс-минус 1-3 секунды) в двух задачах из четырёх я смог подобраться на расстояние меньше километра и продержаться так в течении 1200 секунд… И ничего! Никаких очков нам это не давало! В #icfp-contest и icfpc@cjr наблюдалось несколько человек с подобными проблемами, причём в icfpc@cjr был даже человек, который подбирался к цели буквально на расстояние в миллиметр — и крутился так в течении многих тысяч секунд, и никакого результата это не давало.

Кстати, в некоторых случаях метод задержки давал ошибку в 350 километров, что значительно — ошибался со временем старта почти на 300 секунд. В чём может быть проблема мы так до самого конца соревнования и не поняли, а жаль.

Потратив несколько часов чуть ли не на биение головой об стенку, около семи часов вечера я нашёл корень проблемы. И он оказался в мутной спецификации — оказалось, что на самом деле нам указывается сдвиг Земли относительно спутника, а не сдвиг спутника относительно Земли. То есть, наши координаты нужно инвертировать по знаку. После того, как я это сделал пришлось ещё починить определялку направления вращения. Ну и, естественно, нам стали давать очки за вторые задания, которые мы тут же отправили на сайт.

Приблизительно в это же время выяснилось, что организаторами в первой задаче была допущена ошибка в подсчёте очков — там давали больше очков за сжигание топлива, а не за его сохранение. Вместо того чтоб подправить задание, орги подправили спецификацию — оригинальный путь :) Некоторые люди обсуждали стратегии выжигания топлива аж до конца третьих суток. Предлагали выходить на далёкую орбиту и потом возвращаться на целевую, предлагали дёргать орбиту туда-сюда внизу под тем предлогом, что внизу скорости больше и затраты топлива тоже будут больше… Я сделал самое очевидное — вышел на нужную траекторию, разделил всё оставшееся топливо на две части, и сделал два импульса — половиной толкнул вперёд, половиной толкнул назад, за пределы допустимого километра такой манёвр спутник не выбросил, и отлично. Так мы получили за первые задания 95-96 баллов вместо 60-67 начальных.

Мы начали экспериментировать с разными трансферами и третьим заданием, где-то в этот момент закончился lightning round — к сожалению, в его топ20 мы не вошли, и на тот момент таблица результатов была только для топ20, полную вывесили гораздо позже. Мы набрали где-то 945 баллов, а для попадания в топ20 нужно было около 1100. Минут через двадцать выложили бинарник для четвёртого типа задач, и в очередной раз обновили спецификацию — оказывается, там ещё и Луна есть в четвёртом. xa4a начал рефакторинг наших тупых решателей в более-менее удобоваримые state-machine — до этого у нас был просто метод run() с циклом и переменной state, и длинной простынёй условий, потом начал делать негомановский трансфер.

Я всё это время ещё медитировал над хоть каким-то стабилизатором орбиты или чем-нибудь подобным. К полуночи я придумал и написал совершенно убогую хрень, которая, тем не менее, смогла как-то догонять цели, которые летели в нескольких километрах от нашего спутника и у нас были решения для всех четырёх задач второго типа — это выразилось в увеличении счётчика баллов до 1127 очков. Я ещё минут двадцать пообщался в разных конференциях и около часа ночи опять отправился поспать.

День 2

После того, как я ушёл спать, tilarids ночью сделал канонически правильный фазер, который без всяких проблем достаточно точно подстраивался даже при различиях в треть оборота и работал очень симпатично. В двенадцать дня обнаружилось, что уже доступна вся таблица результатов, и мы в ней были на 36-м месте из более чем двухсот участников, которые сделали хотя бы одно задание.

Я скачал и начал изучать замечательную книгу Orbital Mechanics, которая к тому моменту была уже у всех, кто хотел :) Потом мы решили провести ещё один рефакторинг и собрать весь код в одном месте, плюс к тому я написал удобную запускалку, которая вызывалась просто как ./run.py <task number>. Я ещё написал контроллер-визуализатор для четвёртой задачи, который пока ничего не делал. Но зато он мог показывать все одиннадцать спутников и их траектории движения, собранные спутники выделял другим цветом и даже рисовал где находится заправочная станция.

Да-да, в четвёртой задаче была заправочная станция — если пролететь мимо неё, то она восстанавливала запас топлива в баках нашего космического коня :) В этом задании топлива у нас было не в пример меньше. Луну рисовать я обломался, потому что в большинстве случаев она была далеко за пределами изображаемой области. Такой мы выбирали масштаб, чтоб было видно более подробно что же происходит в околоземном пространстве.

Информацию о всех искусственных и натуральных спутниках Земли, в том числе и АЗС мы доставали из всё тех же портов вывода. Там указывались все координаты относительно нашего управляемого спутника, “поприветствован” спутник или нет, запас топлива на орбитальной АЗС и координаты Луны тоже. На всякий случай уточню — все спутники кроме управляемого летают по одним и тем же траекториям и никак не пытаются упасть на Землю, сгореть в Солнце или отправиться колонизировать Марс или Альфу Центавру :)

Практически весь оставшийся день мы все втроём усиленно читали Orbital Mechanics, пытались как-то определять параметры орбиты, пытались как-то переходить между эллиптическими орбитами и ещё всякая такая муть. У нас ничего не получалось, там постоянно вылазили проблемы с углами, куда что направлено, что где считать и всё такое. Где-то к половине второго меня осенило, и я вспомнил что в питоне есть встроенный тип данных — комплексные числа. С их помощью очень легко производить все векторные вычисления на плоскости, определять углы, радиусы и т.п. На их основе я за полчаса сделал отличную догонялку под названием naive_chase(), которая презрев все законы астрофизики просто напросто летела в сторону цели. Работает догонялка так:

  • определяет векторные скорости цели и управляемого спутника;
  • вычисляет разницу скоростей, которую нужно скомпенсировать для того чтоб выравнять их скорости, то есть результирующие векторы скоростей будут направлены в одну сторону
  • добавляет к этой разнице небольшой вектор в сторону цели (3..50 м/с в зависимости от расстояния);
  • ждать пока не догоним, если расстояние начало увеличиваться — пересчитать и скорректировать траекторию.

Результаты применения оказались замечательными, хотя и топливо она могла жрать немеряно. Пример неумеренного употребления можно посмотреть здесь — красный спутник по Гоману переходит на внешнюю орбиту — а зелёный (цель) к этому моменту ещё даже до нижней точки не долетел :) Поэтому догонялка начинает выделывать такие кренделя.

В результате её применения к 2002-й у нас выросло количество очков с 1127 до 1133, потом ей же решил 3001 и 3004 — получилось 1897 баллов, и ещё 3002 и 3003, и улучшенные 3001 и 3004 — 2364 балла! За двадцать минут мы резко рванули вверх в турнирной таблице :) И опять на такой радостной ноте я пошёл отдыхать :)

А если не делать никаких задержек, и включить сразу Гомановский переход, а потом догонялку, то получилась бы такая картинка.

День 3

Проснувшись с утра, я увидел, что наши результаты за ночь ещё чуть-чуть улучшились, и я ещё чуть-чуть покрутил решения третьих задач, получив 2588.7917 очка и 21 место в таблице — из трёхсот.

Опять почти весь день мы пытались что-то сделать с математикой определения эллиптических орбит и перехода между ними, я пытался развернуть спутник на 4004 задаче… Кстати о развороте спутников, есть несколько способов (для круговой орбиты):

  • На месте где находится спутник просто развернуть тангенциальную скорость (касательную). Для этого нужен импульс, равный двойной скорости. В 4004-й спутник на старте летел по низкой орбите со скоростью около 7 км/с — импульс нужен величиной в 14 км/с, а у нас топлива есть только на 10 км/с. При попытке развернуться - фейл.
  • перейти на вытянутую эллиптическую орбиту и где-нибудь в средине обратного движения переломить траекторию на другую эллиптическую орбиту, и обойти Землю с другой стороны. Таким образом основную работу по развороту мы перекладываем на Землю — чем-то похоже на гравитационную пращу. Метод применить не удалось — опять запутался в математике.
  • Выйти на вытянутую эллиптическую орбиту и в её апогее опять же развернуть тангенциальную скорость. Можно легко совместить с переходом на другую орбиту — очень похоже на би-эллиптический переход. Топлива требуется значительно меньше, чем в первом варианте — раза в два-три (чем более вытянутая орбита, тем меньше топлива)

Я смог без особых проблем реализовать третий метод, но ничего на его основе я уже сделать не успел. За четыре часа до окончания контеста таблицу с достижениями заморозили, и уже никаких движений видно не было. Последние часа два-три мы пытались хоть что-то подкрутить и заработать хоть какие-то очки в дополнение к тем, что были, и нам кое-что удалось — мы поймали три спутника в 4001 и два спутника в 4003 задании, что выразилось в итоговых 2852.2285 баллах.

Да, в четвёртом задании не обязательно было ловить всё — нужно было поймать хоть какой-нибудь спутник. Последний спутник (третий в 4001) поймал tilarids и закинул на сайт буквально за несколько минут до окончания контеста. Всё, хана! :-D

Итоги подведём

Море фана, удовольствие когда что-то начинает получаться и летать, и всё такое. Что можно было сделать лучше:

  • раньше вспомнить про комплексные числа;
  • попытаться во все формулы вручную подставлять числа и смотреть, что же там не так;
  • было бы неплохо всем вместе находиться в одном месте — если бы можно было обсуждать математику с маркером у доски, а не переписываясь по джабберу, то было бы гораздо проще.

Я считаю что мы выступили отлично, и особенно — как для первого раза. Если кроме нас особо никто за последние четыре часа ничего не придумал, то мы заняли 21 место, но это вряд ли. Так что, учитывая некоторое творчество других команд, я думаю что мы занимаем где-то 25-30 место (из 320). Ура! В следующем году тоже постараюсь поучаствовать.

P.S.Отчёт tilarids, xa4a.

Метки

.net .NET C# .sort 1.2 2009 2010 404 error admin ajax amazon analytics and apache api archlinux asp.net async asynchronous autocomplete bash blender blog blogengine blogs book bootstrap bot bpython buildout byteflow bzr C c plus plus C++ cache cbv Chaco checkio chrome ci ckeditor class based views clojure closure cms cms с удобной админкой code coding style collectd COM comet competition conference ConfigParser contest Context continuous integration CouchDB coverage CppCMS cpyext cpython crud csrf CSS ctypes curl custom model fields cx_freeze cython database db dbm dbqueries debian debug debugging decorator decorators deploy deployment descriptor design dev devconf developers development diveintopython Django django 1.2 django 1.3 django advent django framework django template django trunk django weblog django-admin-tools django-cms django-compressor django-hosts django-piston django-registration django-sphinx django.admin djangoadvent djangocms djangodash doc documentation drupal e-legion eclipse EGit emacs encoding Enthought epoll erlang event exception ExtJS fabric facebook fastcgi finaloption fixtures fonts forms formset fp framework freebsd freeswitch fs2web ftp fun funcparserlib functional gae gamin gandi generic views gettext gevent gil git github gitosis Google Google App Engine google picasa Google Translate google wave Google Web Toolkit grab grablab greenlet gtd gui haskell hg hgshelve highlighter host hosting how-to howto html html5lib Hudson humor i18n icfpc ide idiomatic image-scripting improvements Internet interpreter ipython ironpython izmenimsya.ru jabber java javascript jenkins jetbrains JIT job jquery json jstree jython kde kiev kiyv kyivpy l10n ldap library libs Life Links linux Linux & Unix LLVM logging logs lxml Mac OS X magic mail markdown Matplotlib Mayavi maybe mediavirus meetup memcache Memcached memory messages metaclass middleware migration mikrotik mkd model models mod_python mod_wsgi mongodb monitoring mptt musicmans.ru musicx mvc my-projects mysql netCDF networkx newforms newforms-admin news nginx Nhibernate nix nose NoSQL numpy oop open source OpenID openoffice opster optimization oracle orm os pagination parsing path patterns pdf PDF-принтер PEP PEP8 performance performance optimization perl personality photo php picture-driven computing PIL pinax pingback pip plasma plone plugin plugins postgresql programming progress bar psycopg2 py2exe pybb pybbm pycamp pycharm pycon pycow pycurl pydev pygtk pylons PyNGL pypy pyqt PyQt4 pyrad pyramid PySide Python Python 2.5 python 2.7 python 3 python c api python speed python-mssql python3 pywinauto Qt Qt4 queue rabbitmq radius raw sql re redis redsolution redsolution cms regexp regular expressions release repoze.bfg RequestContext reusable apps robokassa rss ru ruby ruby-on-rails sample satchmo scalability SciPy scraping screencast search selenium self.error seo server setattr settings setuptools shell sikuli sms snippet socket.io software sorting south sphinx spider sql sqlalchemy sqlite ssh startup step-by-step subdomain subversion svn SyntaxHighlighter system tags tdd tddspry teh drama template templates templatetags test testing thinkpad threading threads tips tips and tricks tools tornadio tornado tornado server tricks tutorial tweepy twisted twitter typography uapycon Ubuntu ucsvlog uml Uncategorized unicode unit test unit testing UnitTest Unladen Swallow upload urllib urls utf-8 uwsgi validation vcs versioning video vim virtualenv Visual Studio vkontakte voip wave web web-devel web-services web-разработка webdev webfaction webkit webpy websockets webtest widget widgets Win API windows Wirbel work wrapper wsgi wxPython wxWidgets wysiwyg xapian xml xmonad xmpp xpath yandex youtube zip zomg zope [cdata[cbv]] [cdata[ci]] [cdata[class based views]] [cdata[continuous integration]] [cdata[django framework]] [cdata[django-sphinx]] [cdata[django]] [cdata[nginx]] [cdata[python]] [cdata[virtualenv]] [cdata[программирование]] автоматизация администрирование администрирование django админка алгоритмы архитектура атрибуты базы данных Без рубрики безопасность библиотеки блоге бот веб-разработка видео Визуализация данных вконтакте Все записи гвидо ван россум граббер графика графы декоратор декораторы дескриптор дескрипторы документация заметки игра жизнь идея интересное киев Клиентам книги конференция личное математика метаклассы модели модули монады морфология мысли невозможное новости о облачные вычисления обо мне Обработка данных оптимизация оптимизация кода Основная лента основы парсинг парсинг сайтов перевод песочница Питон поебень поиск правила кодирования программирование Проектирование производительность работа рабочее размышлизмы Разное разработка разработка приложений разработки регулярные выражения сайт событие события ссылки статьи тестирование тесты Тюмень убунтариум фигня философия формы форум Хабрахабр хакинг хостинг шаблоны шаблоны проектирования эксперимент Эксперименты юмор я пиарюсь Яндекс