urn:lj:livejournal.com:atom1:janatem Артем Артем 2009-12-06T21:28:56Z urn:lj:livejournal.com:atom1:janatem:63576 2009-12-06T21:25:35Z 2009-12-06T21:28:56Z <span class='ljuser ljuser-name_plakhov' lj:user='plakhov' style='white-space: nowrap;'><a href='http://plakhov.livejournal.com/profile'><img src='http://l-stat.livejournal.com/img/userinfo.gif' alt='[info]' width='17' height='17' style='vertical-align: bottom; border: 0; padding-right: 1px;' /></a><a href='http://plakhov.livejournal.com/'><b>plakhov</b></a></span> подкинул <a href="http://plakhov.livejournal.com/116046.html">задачку</a> размять мозг.<br /><br /><div style="margin-left: 40px;">Есть входящий поток объектов. Сколько их, заранее неизвестно. Для каждого объекта известен его вес (float). Нужно придумать алгоритм, выбирающий N из этих объектов (число N задано свыше) таким образом, чтобы вероятность оказаться выбранным была для каждого объекта пропорциональна его весу, и не зависела от его места в очереди. Память <s>константная</s> O(N) (при этом число объектов K &gt;&gt; N). За один проход. Естественно, алгоритм имеет право использовать random() - иначе непонятно, о каких &quot;вероятностях&quot; идет речь, или уж выбрал, или уж нет.</div><br />После некоторых препирательств по поводу формулировки задачку я решил, но в ней обнаружилось двойное дно, ради чего, собственно, пишу отдельным постом.<br /><br />Как обычно, буду решать математическую, а не программистскую задачу. Поэтому сразу выкидываю расплывчатое и ненужное условие K&gt;&gt;N, ограничиваясь необходимым для осмысленности K&gt;=N. Также забиваю на требование однопроходности и ограничения по памяти. Это потом само получится. ;)<br /><br />На самом деле всё же однопроходность учитывается некоторым образом. А именно, я считаю, что при поступлении k-го объекта уже требуется иметь заготовленный ответ на случай, если этот объект окажется последним. Поэтому в обозначениях ниже понятие <img src="http://www.codecogs.com/eq.latex?A_i^k" alt="" /> является осмысленным для N&lt;=k&lt;=K.<br /><br /><strong>Обозначения.</strong> Пусть <img src="http://www.codecogs.com/eq.latex?\{w_i\}" alt="" /> -- веса объектов; <img src="http://www.codecogs.com/eq.latex?W_k = \sum_{i=1}^k w_i" alt="" />; <img src="http://www.codecogs.com/eq.latex?A_i^k" alt="" /> для k&gt;=N обозначает вероятностное событие &quot;i-й объект был выбран среди k первых&quot;.<br /><br /><a name="cutid1"></a><strong>Теорема.</strong> Для k&gt;N верно <img src="http://www.codecogs.com/eq.latex?p(A_i^k)=\frac{Nw_i}{W_k}" alt="" />. (А для k=N эта вероятность тривиальным образом равна единице.)<br /><br />Из этой теоремы, в частности, выводятся ограничения на веса, при которых задача имеет решение. А именно, выражение справа всегда должно иметь смысл вероятности, т.е. не превосходить единицу.<br /><br />Идея доказательства в следующем. Разобьем вероятностное пространство на элементарные события -- выбраны конкретные N элементов из k возможных -- коих будет <img src="http://www.codecogs.com/eq.latex?C_k^N" alt="" /> штук. Обозначим их вероятности как <img src="http://www.codecogs.com/eq.latex?p_{i_1\dots i_N}" alt="" />, где индекс <img src="http://www.codecogs.com/eq.latex?i_j" alt="" /> соответствует наличию <img src="http://www.codecogs.com/eq.latex?i_j" alt="" />-го объекта в выборке. Тогда полная вероятность выбора i-го объекта равна сумме элементарных вероятностей, у которых присутствует индекс i. По условию задачи эта вероятность равна <img src="http://www.codecogs.com/eq.latex?\lambda w_i" alt="" /> для некоторого значения лябмда. Если просуммировать всё подряд, то получится <img src="http://www.codecogs.com/eq.latex?\sum_{i=1}^k p(A_i^k) = \sum_{\dots} p_{i_1\dots i_N}" alt="" />.<br /><br />Можно заметить, что в правой части каждое элементарное событие будет представлено ровно N раз. Поэтому, учитывая, что сумма вероятностей элементарных событий равна 1, сумма справа равна N. Отсюда находим лябмда и получаем утверждение теоремы.<br /><br /><a name="cutid2"></a><strong>Теперь самое интересное</strong> -- про множественность решений. Хотя мы определили вероятность выбора каждого объекта, однозначно приписать вероятности элементарным событиям в общем случае не удается, потому что уравнений получается порядка N, а элементарных событий <img src="http://www.codecogs.com/eq.latex?C_k^N" alt="" />.<br /><br />Я не поленился расписать случай N=2, K=4 и полностью решить систему уравнений на <img src="http://www.codecogs.com/eq.latex?p_{i_1\dots i_N}" alt="" /> -- действительно, пространство решений двухпараметрично (6 неизвестных, 5 уравнений, общая система имеет ранг 4).<br /><br />Это говорит о том, что есть множество решений, которые дают правильные вероятности, но разное совместное распределение для выбранных объектов. В условии задачи ничего не говорится о совместном распределении, но, поскольку задача взята из какой-то прикладной области (так ведь?), это может иметь значение.<br /><br />Из условия алгоритмической эффективности легко строится одно из решений (привожу его ниже), но я не выяснил, существуют ли другие решения, удовлетворяющие алгоритмическим требованиям. Наверно, достаточно о-большое сделать побольше (в смысле, взять побольше памяти, не нарушив при этом асимптотику), чтобы построить решение с другим совместным распределением.<br /><br /><a name="cutid3"></a><strong>Алгоритм.</strong> Заводим пул на N объектов. Заводим регистр, который хранит W_k на k-м шаге (т.е. при поступлении k-го объекта). Первые N объектов закидываем в пул. Далее, беря очередной k+1-й объект и имея, таким образом, N+1 объект, вспоминаем, что Боливар не вынесет (N+1)-рых -- решаем, кого выбросить: новый объект или что-нибудь из пула. Делается это следующим образом. Пусть <img src="http://www.codecogs.com/eq.latex?B_i^k" alt="" /> -- вероятность выбросить i-й объект на k-м шаге. Тогда по формуле Байеса <br /><img src="http://www.codecogs.com/eq.latex?p(B_i^{k+1}|A_i^k) = \frac{p(B_i^{k+1}\cap A_i^k)}{p(A_i^k)} = \frac{p(A_i^{k+1})}{p(A_i^k)}" alt="" />.<br /><br />Смысл условной вероятности слева таков: &quot;на k+1-м шаге i-й объект будет оставлен при условии, что он дожил до этого момента&quot;. То есть в точности то, что нам нужно при принятии решения на k+1-м шаге! Надо только заметить, что N+1-й шаг особый, потому что вероятность события <img src="http://www.codecogs.com/eq.latex?A_i^N" alt="" /> равна единице, в то время, как на следующих шагах она нетривиальна.<br /><br />Чтобы наш алгоритм заработал, надо сконструировать &quot;бросание жребия&quot; таким образом, чтобы из N+1 объекта (тех, что в пуле, и одного нового) был выброшен ровно один. При этом мы уже нашли вероятности выбрасывания объектов из пула (для i-го объекта это будет <img src="http://www.codecogs.com/eq.latex?1-\frac{p(A_i^{k+1})}{p(A_i^k)}=1-\frac{W_k}{W_{k+1}}" alt="" />; NB: она не зависит от объекта, только от номера шага!) и знаем вероятность выбрасывания свежепоступившего объекта -- она должна быть равна <img src="http://www.codecogs.com/eq.latex?1-p(A_{k+1}^{k+1}) = 1-\frac{Nw_{k+1}}{W_{k+1}}" alt="" />. (Учет-контроль: сумма всех должна быть равна единице. Я даже не поленился проверить апостериори вероятность выживания объекта в конце работы алгоритма.)<br /><br />Как-то так. urn:lj:livejournal.com:atom1:janatem:63326 2009-11-21T15:04:41Z 2009-11-21T15:04:41Z <div style="text-align: right;">Детям нельзя в интернет. От детей интернет тупеет. (<a href="http://bash.org.ru/quote/40084">bash</a>)</div><br />Не знаю, тупеет ли интернет и от чего именно, но, по-моему, это явный признак деградации.<br /><br />Мне понадобилось решить нехитрый дифур, а я позабыл, как это делается. Полез гуглить, а там по запросу типа &quot;методы решения ОДУ&quot; и т.п. все ссылки на первой странице -- про численные методы, маткады и т.п. Как будто люди забыли математику (точнее, ее раздел -- дифуры). (Может быть кто-нибудь покажет мне правильный запрос?)<br /><br />После некоторого напряжения мозга я увидел, как сделать замену, и решил уравнение сам, но осадок остался. Мне неудивительно, что я сам слегка отупел, всё-таки давно не брал в руки шашек. Но куда пропали ссылки на учебники и лекции по дифурам? Может их там и не было? urn:lj:livejournal.com:atom1:janatem:63142 2009-11-21T10:54:28Z 2009-11-21T15:07:17Z Был такой выпуск Ералаша, где чел на уроке труда сделал табуретку трехметровой высоты, &laquo;...сказали &laquo;семь раз отмерь, один раз отрежь&raquo;, ну я так и сделал -- отмерил семь раз&raquo;.<br /><br />В книжке по оригами дети нашли модель лебедя и сделали:<br /><br /><img src="http://xoft.ru/~art/lj-share/lebed.jpg" alt="" /><br /><br />Там было написано, что понадобится 20 листов А4 (это довольно много, большинство моделей требуют 1-2 листа). Но дальше было написано то, что дети поленились читать. А именно: нарезать листы А4 на куски вроде 6х8 см. Ну вот и получилось такое чудище. ;)<br /><br /><strong>Upd.</strong> Он еще и крылья умеет расправлять:<a name="cutid1"></a><br /><br /><img src="http://xoft.ru/%7Eart/lj-share/lebed2.jpg" alt="" /> urn:lj:livejournal.com:atom1:janatem:62734 2009-11-12T21:20:05Z 2009-11-12T21:20:05Z Странное дело, как нелогично мыслят дизайнеры софта. В данном случае речь пойдет о почтовых клиентах.<br /><br />Недавно решил разгрести авгиевы конюшни и попытаться привести в порядок свою почту -- создать правильные фильтры, метки и т.д. Одно из вопиющих неудобств -- это то, что мои (т.е. отправленные мною) письма сохраняются в отдельное место (которое называется &quot;Отправленные&quot;), вместо того, чтобы сортироваться по общему принципу. Например, если пишу в группу рассылки, то должно попадать туда же, куда сваливается от других членов группы.<br /><br /> Хотел было по-быстрому исправить это, но оказалось, что задача нетривиальна. Я пользуюсь evolution и браузером для gmail-а. Ни там, ни там я не нашел естественных способов сделать, что хочу. Раньше я пользовался емаксовым gnus-ом, и там была возможность автоматически добавлять что-нибудь к заголовкам, поэтому там работает добавление bcc себе. Правда, это решение, хоть и работает, попахивает извращением. Как сделать это в гугле или в evolution, я не так и не придумал.<br /><br />Разве я хочу странного? По-моему, естественно, чтобы фильтры работали единообразно как для входящих, так и для исходящих, и чтобы то, что гмейл называет цепочками, никогда не ломалось. А функция &quot;найти все мои письма&quot; и так тривиально выражается через общие механизмы; ну ладно, если кому-то очень хочется, путь на нее будет специальная кнопка или виртуальная папка. urn:lj:livejournal.com:atom1:janatem:62468 2009-10-25T17:07:53Z 2009-10-25T17:07:53Z Есть идея о том, как контролировать точность компьютерных вычислений. Представление с плавающей точкой не дает никакой гарантии точности,&nbsp;поэтому для каждого алгоритма нужно отдельно доказывать корректность. Интервальную арифметику придумали для того, чтобы доказывалось само,&nbsp;но у нее есть проблема с тем, что длины интервалов слишком быстро растут. Если бы имелся механизм, учитывающий функциональную зависимость между аргументами, то, возможно, всё было бы не так плохо.<br /><br /> Я предлагаю другой, возможно, новый подход, который, по идее, дает стохастическое доказательство корректности и учитывает зависимости. Я потратил некоторое время, чтобы проверить, нет ли уже чего-нибудь в этом роде. Вроде бы прослеживается некоторая связь с нечеткой логикой (точнее, с нечеткой интервальной арифметикой -- обобщением интервальной арифметики). Может быть, стоит копать куда-то в строну так называемой вероятностной логики (product logic), чтобы увидеть, что предлагаемое мной уже изобретено. Тем не менее попытаюсь изложить концепцию.<br /><br />Вместо того, чтобы представлять вещественное число в виде интервала или в формате IEEE 754 (что в общем-то тоже в некотором смысле есть интервал), предлагается задавать число в виде случайной величины с заданным распределением. Для этого надо сделать две вещи:<ol><li>Зафиксировать некоторое параметризованное семейство распределений, так чтобы набор параметров можно было представить в виде битовой последовательности (скорее всего ее длина должна быть фиксированной и равно какому-то круглому количеству битов).</li><li>Определить частичный порядок на нашем семестве распределений, так чтобы можно было сказать, что одно распределение &quot;точнее&quot; другого.</li></ol>Имея частичный порядок, всегда можно можно перейти от &quot;реального&quot; представления <img src="http://www.codecogs.com/eq.latex?a(p_1, \dots, p_n)" alt="" /> с вещественными параметрами к &quot;дискретному&quot; <img src="http://www.codecogs.com/eq.latex?a(\tilde p_1, \dots, \tilde p_n)" alt="" />, где параметры <img src="http://www.codecogs.com/eq.latex?\tilde p_i" alt="" /> конечно-преставимы в машине.<br /><br />Как определить порядок &quot;точнее&quot;? Мне не удалось найти универсального определения, но можно ввести такое понятие, назовем его уровнем значимости, -- положительное число (близкое к нулю), которое выбирается для каждой задачи заранее из каких-то физических или &quot;вычислительных&quot; соображений. Тогда для заданного уровня значимости <img src="http://www.codecogs.com/eq.latex?\varepsilon" alt="" /> два распределения a и b связаны отношением <img src="http://www.codecogs.com/eq.latex?a\prec b" alt="" /> (a &quot;точнее&quot; b), если для любого интервала <img src="http://www.codecogs.com/eq.latex?B\subset R,\quad \mu_b(B) \ge 1-\varepsilon" alt="" /> существует интервал <img src="http://www.codecogs.com/eq.latex?A\subset B,\quad \mu_a(A) \ge 1-\varepsilon" alt="" />.<br /><br />Для примера можно попробовать двупараметрическое семейство нормальных распределений N(c, s), пополненное дельта-функциями (формально записываемыми как N(c, 0)). Здесь при любом уровне значимости N(c, 0) всегда точнее N(c, s), и вообще <img src="http://www.codecogs.com/eq.latex?N(c, s_1)\prec N(c, s_2)" alt="" /> при s1 &lt; s2. Если взять разные c, то сравнимость будет зависеть от уровня значимости.<br /><br />Используя порядок, можно сформулировать задачу выполнения арифметической операции как нахождение какого-нибудь представимого распределения, так чтобы истинное распределение было &quot;точнее&quot; найденного. При этом желательно подыскивать не какое попало, а &quot;поточнее&quot;, т.е. &quot;наименьшее&quot; из тех, которые &quot;больше&quot; истинного. (Пока не обращаем внимание на то, что существование &quot;наименьшего&quot;, вообще говоря, под вопросом.) Это задача кажется трудоемкой, поэтому, во-первых, надо позаботиться о рациональном выборе представляемого семейства и, во-вторых, о том, чтобы эффективно реализовать арифметику, впоть до привлечения специальных аппаратных средств.<br /><br />Теперь о зависимостях. Можно считать, что задача формулируется в следующем виде: дано совместное распределение n чисел, результат надо выразить в виде совместного распределения m чисел. Если подойти наивно, то количество информации, требуемое для представления вектора, зависит экспоненциально от длины вектора. Так, скорее всего, не годится. Можно заметить, что почти все элементарные операции над числами бинарны (ну или унарны). Тогда можно хранить только попарные зависимости, коих O(n^2), что гораздо лучше, чем экспонента. Получается, что исходный вектор (а также результатный и все промежуточные) хранится в виде набора n распределий и некоего обобщения ковариационной матрицы -- набора попарных распределений. (Здесь наверно есть избыточность -- вовсе незачем хранить попарное распределение двух независимых чисел...)<br /><br />Проблем здесь, конечно, куча.&nbsp;Я пока не готов описать арифметику для вышеупомянутого семейства гауссиан даже для простых случаев (ну разве что a+b для независимых). Кроме того, совершенно неясно, что делать с операциями арности больше 2 (хотя бы if, у которого условием стоит a&gt;0).<br /><br /><br /> urn:lj:livejournal.com:atom1:janatem:62266 2009-10-19T13:30:17Z 2009-10-19T13:30:17Z Только что наблюдал прикольный глюк. Когда идешь на gmail.com, то запрос обычно перенаправляется куда-то вроде <a href="https://www.google.com/accounts/ServiceLogin?чего-то_там=mail.google.com">https://www.google.com/accounts/ServiceLogin?чего-то_там=mail.google.com</a> . Смысл глюка был в том, что на этом процесс не останавливался, а в URL-е дальше происходила замена этого mail.google.com . Причем это происходило рекурсивно; до тех пор, пока не упиралось в ограничение браузера на длину URL-а.<br /><br />Минут через 10 всё пофиксилось само, но осадок остался... urn:lj:livejournal.com:atom1:janatem:62125 2009-10-04T11:02:01Z 2009-10-04T17:21:42Z Вот какой забавный подарок крестного достался детям на день рожденья -- <a href="http://xoft.ru/~art/lj-share/incredible_machine.avi">видео</a> (95 МБ). Правда, там написано &quot;от 15 лет&quot;, да еще инструкция не содержит русских букв, так что мое участие в сборке было существенным.<br /><br />Можно было бы мертвую петлю побольше сделать, если бы шарик был не обычный,&nbsp;стальной, а, например, пластмассовый со свинцовым ядром, чтобы момент у него поменьше был. Но всё равно, и так круто.<br /> urn:lj:livejournal.com:atom1:janatem:61749 2009-10-01T20:35:55Z 2009-10-01T20:35:55Z Мне тут напомнили, что у меня скоро юбилей,&nbsp;очень круглый. 12 октября исполняется миллиард секунд! Наверно надо что-то замутить, ибо столь круглых дат почти не бывает. ;)<br /><br />Много это или мало? Если пренебречь тем, что гига- отличается от гиби-,&nbsp;то это половина 32-битной эпохи. (Не четверть, а половина,&nbsp;потому что один бит зарезервирован под знак -- можно отмерять от начала эпохи не только вперед,&nbsp;но и назад.) На самом деле 10^9 заметно отличается от 2^30, так что бинарный юбилей будет на два года позже... urn:lj:livejournal.com:atom1:janatem:61559 2009-09-29T20:41:05Z 2009-10-02T08:06:12Z Давненько я здесь не оставлял записей.&nbsp;И&nbsp;не потому что иссякла мысля и нечего выразить, и даже не из-за лени. Просто образовалось перманентное ощущение цейтнота...<br /><br />Для начала -- просто лытдыбр.<br /><br />Я недавно вернулся с конференции в Абрау-Дюрсо &quot;Научный сервис в сети Интернет&quot;.&nbsp;Несмотря на название,&nbsp;инета там не было, а мой телефон почему-то не давал gprs (и я пока не осилил разобраться, видимо придется пешком идти в офис мтс). В результате образовалось немного скучных, без компромата, <a href="http://xoft.ru/~art/pics/2009-09-xx.Abrau.byArt/">фоток</a>. Также в процессе сбора воедино фотографии, сделанные участниками (пока в доступности вижу только <a href="http://picasaweb.google.com/katrin.elinor/200903?feat=embedwebsite#">эти</a> от <span class='ljuser ljuser-name_katrin_elinor' lj:user='katrin_elinor' style='white-space: nowrap;'><a href='http://katrin-elinor.livejournal.com/profile'><img src='http://l-stat.livejournal.com/img/userinfo.gif' alt='[info]' width='17' height='17' style='vertical-align: bottom; border: 0; padding-right: 1px;' /></a><a href='http://katrin-elinor.livejournal.com/'><b>katrin_elinor</b></a></span> ).<br /><br />Еще я занимался нехарактерной для моей специальности деятельностью, а именно -- художественной. От чего получил моральное удовлетворение по полной программе, пытаясь раскрыть тему сисек. (Вроде бы можно выложить результат в открытый доступ, но я пока не получил четких санкций.) Также открыл для себя занятный инструмент для самовыражения в векторной графике -- <a href="http://www.inkscape.org/">inkscape</a>. (Конечно же довольно быстро обнаружился ряд недостатков и недоделок, но всё равно штучка очень хорошая; можно давать детям.)<br /><br /><strong>Upd</strong>. Кроме фотографий зацените модное доменное имя: <a href="http://exascale.ru/abrau/09/">http://exascale.ru/abrau/09/</a> .<br /> urn:lj:livejournal.com:atom1:janatem:61240 2009-09-09T08:44:21Z 2009-09-09T08:44:21Z На xkcd есть картинка, от которой меня просто прет.<br /><br /><a href="http://xkcd.com/55/"><img src="http://imgs.xkcd.com/comics/useless.jpg" alt="" /></a><br /><br />Правда там нелогично написано преобразование Фурье (непонятная штука была использована как связанная переменная), к тому же с ошибкой (минус потерялся). Но всё равно -- для меня это больше чем картинка-анекдот; она обозначает тип мышления.<br /><br />А вот есть <a href="http://xkcd.com/150/">другая картинка</a>, которая тоже в топе. Но я не понимаю ее смысла. Объясните мне кто-нибудь, а то я вижу там только сцену инцеста.<br /><br /><br /> urn:lj:livejournal.com:atom1:janatem:60972 2009-08-20T19:50:48Z 2009-08-20T19:50:48Z Я, кажется,&nbsp;изобрел полезную штуку. Настолько полезную, что не понимаю, почему это не внедрили. (Ну уж совсем невероятно, что это еще не придумали.)<br /><br />Проблема такая. Кто ездит на автомобиле ночью, знает,&nbsp;что есть товарищи, которые нам совсем не товарищи,&nbsp;и не удосуживаются нормально настроить фары. Также, бывает, забывают вовремя выключить дальний свет. Ну и вообще,&nbsp;даже если все всё делают как надо, свет чужих фар доставляет неудобства.<br /><br /><a name="cutid1"></a>Идея в следующем: пусть свет фар всех автомобилей будет поляризованным, а все зеркала и лобовое стекло отражают/пропускают свет с перпендикулярным направлением поляризации. (Может быть со стеклом возникнут технические сложности,&nbsp;тогда достаточно поляризационных очков для водителя.)<br /><br />Конечно, направление поляризации света фар должно быть у всех одинаковым, а то выйдет фигня, как с право- и левостороннестью движения. Возможно, есть некоторые соображения,&nbsp;из которых следует, что вертикальная (или горизонтальная) поляризация по каким-то соображениям лучше (например, из-за естественной поляризации солнечного света атмосферой -- получаем халявный бонус к защите от солнца у горизонта). Ну тем проще выбрать.<br /><br />Теперь о том,&nbsp;как внедрить. Ясно, что покупать очки (или менять зеркала/стекла) -- дело выгодное исключительно владельцу. Заставить же произвести замену фар сложнее, придется простимулировать автопроизводителей. Однако мы видим, что различные системы безопасности, контроля выхлопа и т.п., напрямую не нужные покупателю, более-менее успешно внедряются,&nbsp;поэтому можно поверить, что необходимые механизмы воздействия существуют. В результате поляризующие фары должны постепенно вытеснить обычные.<br /><br />Я нашел пока только одну техническую сложность: если тупо поставить поляризационные фильтры к обычным лампам, то получим 50% потери света. Поэтому хорошо бы как-то усовершенствовать конструкцию фары в целом. Также, возможно, придется считаться с тем, что отраженный от дороги свет частично сохранит начальную поляризацию, и из-за этого эффективность освещения немного упадет.<br /> urn:lj:livejournal.com:atom1:janatem:60812 2009-08-16T09:59:42Z 2009-08-16T09:59:42Z Снилось, будто сторговал с Маратом большой шлем. Раздача была такая: Т??, ТДхх, Т??, x (распределение фосок в трефе и черве не помню). Торговля шла очень глупо -- я обсчитался в количестве тузов и спасовал на открытии. Партнер открылся в мажоре, потом убежал в 3 бубей. Тут я обнаружил потерявшегос я туза и, кроме того, заметил, что есть все задержки и почти все они первого класса. Вполне шлемовая карта. Учитывая, что партнер еще не знает про бубновый фит, но ленясь считать силу, говорю сразу 6 бубей. Он поднимает до 7, и я просыпаюсь. Зная, что партнер склонен к импульсивности, ощущаю беспокойство -- что ж там дальше? Долго пытаюсь снова заснуть, чтобы посмотреть продолжение. Наконец, мне удается вернуться в сон, и я смотрю расклад. Вижу у партнера Вхх в пике при почти равномерном раскладе и от расстройства сразу просыпаюсь,&nbsp;даже не проверив, берется ли малый шлем. :( urn:lj:livejournal.com:atom1:janatem:60475 2009-08-13T14:38:29Z 2009-08-13T14:38:29Z Сегодня продолжаю писать всякие гадости. На этот раз, запоздало выполняя просьбу <span class='ljuser ljuser-name_yurikl' lj:user='yurikl' style='white-space: nowrap;'><a href='http://yurikl.livejournal.com/profile'><img src='http://l-stat.livejournal.com/img/userinfo.gif' alt='[info]' width='17' height='17' style='vertical-align: bottom; border: 0; padding-right: 1px;' /></a><a href='http://yurikl.livejournal.com/'><b>yurikl</b></a></span>&nbsp;, расскажу про работу наших интернет-провайдеров.<br /><br />Было так: наш домашний интернет был подключен через ADSL&nbsp;от domolink'а. В июле там задумали менять какое-то оборудование. Из-за этого сеть пропала не на день, как можно было ожидать,&nbsp;а более чем на две недели. Даже не горела лампочка DSL. Собственно,&nbsp;из-за этого я <a href="http://janatem.livejournal.com/60067.html">возился со скайлинком</a>.<br /><br /><a name="cutid1"></a>В августе я не выдержал и переподключился к другой конторе -- billing-lan.net . Раньше мы имели с ними дело,&nbsp;но потом нашлось несколько причин (в том числе высокие тарифы) уйти к домолинку. У биллинга вроде бы всё сделано красиво -- в дом заходит прямо-таки ethernet. Но ихний VPN -- это гемор еще тот. Мне кое-как удалось настроить соединение с компьютера,&nbsp;но ведь не хочется же держать его,&nbsp;большой и гудящий, включенным 24/7 для того, чтобы жила вся домашняя сеть. Надо устанавливать VPN-соединение с рутера. Однако рутер -- устройство довольно тупое,&nbsp;поэтому заставить его делать слегка нестандартные действия непросто.<br /><br /> Проблема в том, что VPN-сервер находится не в локальной сети, и нужно ходить к нему через шлюз. При этом, если установлен дефолтный маршрут, то поднявшийся PPP перезатирает его своей обычной опцией replacedefaultroute, в результате срубая сук, на котором сидит. Поэтому приходится держать специальный route для доступа подсети, где живет vpn-сервер. Моя коробочка не имеет в конфиге отдельного поля для локального шлюза, она думает, что vpn-сервер всегда живет в локальной подсети. Поэтому пришлось брать в руки бубен и настраивать маршруты вручную. При этом из-за ошибки в прошивке мои настройки часто не могут пережить ресета или даже временного обрыва VPN-соединения -- почему-то адрес интерфейся WAN самопроизвольно заменяется на адрес ppp-интерфейса. (Железка DLink Dir-320 с прошивками dd-wrt, я пробовал три версии прошивки,&nbsp;но все они с этой ошибкой.) Я пока не осилил создание хака для этой прошивки,&nbsp;чтобы запускать в нужный момент загрузки мой скрипт из nvram. Т.е. я не хочу осваивать весь toolchain для создания новой прошивки, а пытаюсь перебиться тем, что можно записывать в nvram.<br /><br /> В общем этот виндовозно заточенный протокол PPTP меня бесит. Неужели нельзя сделать чего-нибудь более цивильного? Хотя бы тот же openvpn, который работает понятно и надежно. Вообще ведь ничто не мешает держать набор сервисов, чтобы клиент мог выбрать наиболее ему подходящий.<br /><br />А еще биллинговский саппорт отжигает. У них довольно часто отваливается то vpn-сервер, то шлюз. В очередной раз звонил им по этому поводу, говорю, что шлюз не работает. Оператор-девушка начинает что-то впаривать про то, устанавливается ли соединение и т.п. Я пытаюсь донести, что если шлюз не работает, то и VPN-соединение в принципе не может быть установлена. Она вообще создает впечатление достаточно грамотного для такой должности айтишника,&nbsp;но здесь как будто специально начинает прикидываться дурой. &quot;Мы не отвечаем за настройки маршрутизатора, попробуйте с компьютера&quot; и т.д. Я думаю, хрен с вами, цепляю нетбук,&nbsp;меняю мак-адрес -- вижу шлюз есть,&nbsp;и даже vpn поднимается;&nbsp;переподключаю опять рутер -- тоже всё работает. Т.е.&nbsp;они в саппорте приучены зачем-то сначала лапшу вешать,&nbsp;а потом,&nbsp;пользуясь временным замешательством клиента,&nbsp;быстренько проверяют свои железки.&nbsp;Часто за это время удается поправить все косяки, и в результате -- клиент дурак, а саппорт весь в белом.<br /><br />Однако, возвращаясь к домолинку. По совету <span class='ljuser ljuser-name_yurikl' lj:user='yurikl' style='white-space: nowrap;'><a href='http://yurikl.livejournal.com/profile'><img src='http://l-stat.livejournal.com/img/userinfo.gif' alt='[info]' width='17' height='17' style='vertical-align: bottom; border: 0; padding-right: 1px;' /></a><a href='http://yurikl.livejournal.com/'><b>yurikl</b></a></span>&nbsp; я попробовал изменить настроки в разделе ADSL: бывают какие-то ADSL2, ADSL2+ и т.д. (точнее, я тупо включил все галки, чтобы он сам выбирал). В результате случилось чудо: и светодиод DSL зажегся,&nbsp;и, более того, оно соединилось и заработало. Правда,&nbsp;как оказалось,&nbsp;работает плохо, довольно часто отваливается. Да и всё равно подло с их стороны после смены железа не сказать о том,&nbsp;что надо какие-то заточки сменить.<br /><br />В результате имеем два интернета, оба плохо работающие,&nbsp;но плохо по-разному. (Хех,&nbsp;что получается! Все хорошо работающие устройства работают одинаково хорошо,&nbsp;каждое плохо работающее устройство плохо работает по-своему.)<br /><br />Интересно, это везде так или где-то всё же могут обеспечить нормальный сервис?&nbsp;И&nbsp;связаны ли текущие проблемы с периодом отпусков?<br /> urn:lj:livejournal.com:atom1:janatem:60161 2009-08-13T11:46:51Z 2009-08-13T11:46:51Z <a href="http://fregimus.livejournal.com/66234.html">Суть</a> в том,&nbsp;что <span class='ljuser ljuser-name_fregimus' lj:user='fregimus' style='white-space: nowrap;'><a href='http://fregimus.livejournal.com/profile'><img src='http://l-stat.livejournal.com/img/userinfo.gif' alt='[info]' width='17' height='17' style='vertical-align: bottom; border: 0; padding-right: 1px;' /></a><a href='http://fregimus.livejournal.com/'><b>fregimus</b></a></span>&nbsp; неточно процитировал некоего журналиста <span class='ljuser ljuser-name_scottishkot' lj:user='scottishkot' style='white-space: nowrap;'><a href='http://scottishkot.livejournal.com/profile'><img src='http://l-stat.livejournal.com/img/userinfo.gif' alt='[info]' width='17' height='17' style='vertical-align: bottom; border: 0; padding-right: 1px;' /></a><a href='http://scottishkot.livejournal.com/'><b>scottishkot</b></a></span>&nbsp;(причем,&nbsp;кажется,&nbsp;более-менее правильно передав дух высказывания, но изменив букву). А тот прямо таки в суд подал за клевету (или пока только пригрозил).<br /><br />На месте fregimus'а мог бы оказаться каждый. Даже я, хотя с меньшей вероятностью, -- я бы, конечно, в такое непонятное месиво не полез, однако иногда позволяю себе и в более извращенной форме перевирать чужие высказывания.<br /><br />Поэтому мне представляется полезным проанализировать причины столь резкой и явно неадекватной реакции. Может быть, журналист преследовал какие-то свои, мне непонятные, пиар-цели. А может он просто в каком-то непонятном порыве написал гневный комментарий,&nbsp;а потом пожалел об этом, да поздно было ,-- получилось, взял сам себя на понт. Кстати, сам <a href="http://bars-of-cage.livejournal.com/562731.html?thread=8269867#t8269867">текст</a>, помимо легко извлекаемого смысла содержит в себе куски чистого словоблудия, столь искусного,&nbsp;что я, имея необходимость на него отвечать,&nbsp;боролся бы с искушением выразить свой ответ на каком-нибудь языке программирования (например, брейнфаке), с тем, чтобы отомстить автору нарочитой непонятностью за нарочитую непонятность.<br /><br />Я, признаюсь, не осилил разобраться в сути флейма -- почитал лишь по диагонали. Ну собачатся люди себе и ладно. Но стоило fregimus'у вынести свое замечаение в другое место... Кстати, может быть в этом дело? Наверно людям не нравится,&nbsp;когда о них говорят в третьем лице. Я сам сталкивался с подобной, столь же импульсивной, реакцией -- после вялотекущей перебранки в чужом журнале оставил о ней несколько язвительное упоминание в своем журнале.<br /> urn:lj:livejournal.com:atom1:janatem:60067 2009-08-05T23:40:45Z 2009-08-05T23:40:45Z Мне попалась в руки такая игрушка -- модем C-motech CNU-680. Я хотел быстренько справить интернет-нужду, но не тут-то было.<br /><br /> Мы уже привыкли, что многие производители железа не заботятся о пользователях, которые живут на линуксе. Чисто формально по задумке разработчиков этот CNU-680 должен работать и с линуксовым компом, даже без специального ядерного модуля. Реально это сделано так: usb-устройство при втыкании обнаруживается как несколько устройств: (1) обычный флеш-накопитель, (2) CD-диск, а также (3) некое дополнительное SCSI generic устройство (я так понял, три штуки /dev/sg* отвечают соответственно за (1) и (2), а также за еще одно дополнительное).<br /><br />Первое на основную функциональность не влияет, просто дает флешку (правда, я не смог ее подмонтирова, там какая-то неизвестная таблица разделов,&nbsp;ну и фиг с ней, не стал разбираться).<br /><br /> На CD-диске находятся драйвера для разных ОС,&nbsp;причем то, что лежит в директории Linux, заслуживает отдельного комментария. Там есть бинарник с именем RDEVCHG, который надо запустить из-под рута и скрипт,&nbsp;который запускает wvdial с якобы правильным конфиг-файлом. Еще (что само по себе замечательно) там есть исходник этого RDEVCHG.<br /><br />Ни бинарник, ни скрипт у меня не заработали, но это само по себе не примечательно. Важно другое -- созерцание исходника и скрипта вызывает ощущение, что разработчики намеревались плюнуть в душу всем пользователям, и это, похоже, им удалось. Я даже не уверен, что автор кода дошел до стадии &quot;works for me&quot; -- там реально написан полный бред. Что я выяснил: требуется всего-то сделать один вызов ioctl() в применении к устройству /dev/sg* (номер 3 в моем перечне). Хитрость лишь в том, чтобы угадать правильный номер устройства (вместо звездочки). Результатом этого ioctl() должно быть появление нового устройства /dev/ttyACM*,&nbsp;которое является одним стандартным устройств соответствующих модему; с ним уже работает wvdial (ну или напрямую pppd).<br /><br />В общем, они пытаются перебором проверить все /dev/sg*, но делают это неправильно -- условие успеха не то. Кроме того, они неправильно интерпретируют возвращаемое значение ioctl(). Не знаю, есть ли дистрибутив линукса, на котором этот код работает. Даже если бинарник где-то сработает,&nbsp;то скрип точно нерабочий. Еще они ввели меня в заблуждение относительно цели выполнения ioctl(),&nbsp;поэтому я не смог правильно подхачить код, чтобы всё заработало.<br /><br />Потом уже я нашел на sourceforge маленький проект <a href="http://cmotech-tools.sourceforge.net/">cmotech-tools</a> -- тот же самый вызов ioctl(), только правильно и аккуратно обернут.&nbsp;Плюс к тому примочка с udev, которая запускает всё автоматически при включении (а также, похоже, само восстанавливает соединение при обрыве (а обрывы, как я успел убедиться, довольно часты)). В общем, это cmotech-tools прекрасно заработал у меня (правда я не стал заморачиваться с udev -- всё вручную). (Только меня подстерегли еще одни грабли -- долго не мог понять, где брать логин/пароль для скайлинка.)<br /><br />PS. А производители этой железки -- козлы. Отобрали у меня кучу времни, причем совершенно не по делу.<br /><br /> urn:lj:livejournal.com:atom1:janatem:59771 2009-08-05T22:45:14Z 2009-08-05T23:46:49Z Во-первых, я теперь считаюсь больным и сижу дома, на работу не хожу. После похода в Хибины у меня заболела спина и я, наконец, добрался до врача. Теперь вот лечусь.<br /><br />Во-вторых дома уже две недели нет инета. Кто-то -- то ли провайдер, то ли АТС -- менял оборудование, и теперь DSL не работает. В Домолинке придумали клевую отмазку: согласно договору они имеют право &quot;починять&quot; (т.е. с чистой совестью вообще ничего не делать) в течение двух рабочих недель. И это вроде как написано в договоре. (Почитать договор что ли?)<br /><br />Еще я убил кучу времени на то, чтобы подключиться через gprs. Всё безрезультатно -- не помогл даже мой старый телефон, который я точно использовал как модем. Да, тут еще проблема курицы и яйца -- чтобы подключиться к инету, надо вкурить мануалы, которые водятся в инете... В результате я на несколько дней оказался почти в полной изоляции.<br /><br />Из более оптимистичного:<ul><li>Я взял погонять скайлинковский модем и кое-как подключился через него&nbsp;(но эти грабли -- отдельная тема).</li><li>Наконец-то подошел к Хаселю на интимное расстояние (благо необходимый софт был заранее установлен, а также распечатана брошюрка). Я восхищен! Если раньше я осознавал, что алгебраические полиморфные типы данных в ML по выразительности сравнимы с классами в ООП, то сейчас вижу, что ООП как парадигма программирования просто сосет у хаскельных классов (особенно, если включить какие-нибудь расширения, например GADT -- &quot;обобщенные&nbsp; алгебраические типы&quot;).</li><li>Я самостоятельно заботал, как делать внутримышечную инъекцию. И, собственно, сделал. Надеюсь Хе, когда узнает, меня не прибьет за то, что я постеснялся обращаться к соседу-врачу, а сделал всё сам. ;)</li></ul> urn:lj:livejournal.com:atom1:janatem:59474 2009-07-11T18:08:08Z 2009-07-11T18:08:08Z ...и переждем там недели полторы.<br /><br />Это примерно здесь:&nbsp;<a href="http://xoft.ru/~art/hibiny/01.jpg">http://xoft.ru/~art/hibiny/01.jpg</a> (осторожно,&nbsp;19 МБ)<br /><br />Стартовые массы: 34.0,&nbsp;22.7, 6.0, 6.0 кг. urn:lj:livejournal.com:atom1:janatem:59311 2009-07-03T19:33:56Z 2009-07-03T19:33:56Z Роясь в блоге замечательной харьковской команды THIRTEEN (я уже упоминал ее), увидел пост с вызывающим текстом: <a href="http://thirteen.kharkov.ru/2009/07/01/icfpc-2009-thirteen-testcase-4001-video/">Здесь базу посетили _после_ Луны</a>.<br /><br />Созерцание ролика вызвало у меня культурный шок.&nbsp;Они реально жгут! Без дозаправки облетают все цели. Главное обидно, что мы, выведя все нужные формулы, не доперли до очень простых вещей (которые становятся наглядными на видео):<ol><li>Для изменения энергии орбиты выгодней всего ускоряться в перигее. Причина довольно очевидна: скорость максимальна =&gt; приращение квадрата скорости максимально. Энергия однозначно задает a -- большую полуось эллипса.</li><li>Для изменения других параметров выгоднее ускоряться в апогее. Это верно по крайней мере для угла наклона оси эллиса относительно координатной оси. (Где лучше менять эксцентриситет, не меняя энергии, так сходу не скажу, но это вроде и не требуется.)</li><li>Почему-то я не верил, что перехват может быть успешным и при большой разности векторов скорости. Из-за этого мы рассматривали только орбиты касательные к целевым. Сейчас я почти уверен, что мой алгоритм локальной доводки будет успешно работать не только на касательных. По крайней мере вдалеке от Земли -- там линеаризация будет достаточно точной.</li></ol> urn:lj:livejournal.com:atom1:janatem:59019 2009-07-01T21:32:03Z 2009-07-01T21:32:03Z Еще раз о пользе и вреде чистой, непрограммистской математики в условиях соревнования.<br /><br />Как восторженно и поэтично <a href="http://thirteen.kharkov.ru/2009/06/30/san-esse-o-polze-matematiki/">пишет</a> об этом участник-нематематик (команда THIRTEEN):<br /><br /><div style="margin-left: 40px;">Понимающие математику! О, они подобны музыканту, берущему нотную запись незнакомого музыкального произведения, и начинающему играть ноты, претворяя их в жизнь. Понимающие математику берут с нужной полочки формулы, и знают, что в них к чему! Для того, чтобы перевести тело из одной любой точки в другую в орбитальном пространстве земли, нет нужды методом приближений искать нужный вектор приложения сил, запуская каждый раз прилагаемую программу, и наблюдая, чем это закончится, а потом чуть-чуть изменять вектор и пробовать еще раз, приближаясь к нужному варианту. Так делаем мы, люди, но не математики. [...]</div><br />А вот слова математика <span class='ljuser ljuser-name_janatem' lj:user='janatem' style='white-space: nowrap;'><a href='http://janatem.livejournal.com/profile'><img src='http://l-stat.livejournal.com/img/userinfo.gif' alt='[info]' width='17' height='17' style='vertical-align: bottom; border: 0; padding-right: 1px;' /></a><a href='http://janatem.livejournal.com/'><b>janatem</b></a></span>'а:<br /><br /><div style="margin-left: 40px;">В очередной раз напоминаю (в первую очередь себе!), что аналитические методы на ICFPC -- это зло.&nbsp;Конечно, приятно аналитически решить задачу и получить удовольствие, но всегда есть риск, что это не удастся сделать за разумное время. [...]</div><br /> urn:lj:livejournal.com:atom1:janatem:58784 2009-07-01T14:11:05Z 2009-07-01T17:02:15Z <strong>Событие</strong>:&nbsp;http://icfpcontest.org/<br /><strong>Команда</strong>: Error 404<br /><strong>Состав</strong>:<a href="http://yurikl.livejournal.com/profile"><img width="17" height="17" style="border: 0pt none ; vertical-align: bottom; padding-right: 1px;" alt="[info]" src="http://l-stat.livejournal.com/img/userinfo.gif" /></a><a href="http://yurikl.livejournal.com/"><b>yurikl</b></a>&nbsp;, <span class='ljuser ljuser-name_janatem' lj:user='janatem' style='white-space: nowrap;'><a href='http://janatem.livejournal.com/profile'><img src='http://l-stat.livejournal.com/img/userinfo.gif' alt='[info]' width='17' height='17' style='vertical-align: bottom; border: 0; padding-right: 1px;' /></a><a href='http://janatem.livejournal.com/'><b>janatem</b></a></span>&nbsp;, Антон и еще Антон&nbsp;(ощутимо не хватало <span class='ljuser ljuser-name_druxa_druxa' lj:user='druxa_druxa' style='white-space: nowrap;'><a href='http://druxa-druxa.livejournal.com/profile'><img src='http://l-stat.livejournal.com/img/userinfo.gif' alt='[info]' width='17' height='17' style='vertical-align: bottom; border: 0; padding-right: 1px;' /></a><a href='http://druxa-druxa.livejournal.com/'><b>druxa_druxa</b></a></span>)<br /><a href="http://icfpcontest.org/scoreboard.php"><strong>Предварительный результат</strong></a>:&nbsp;2735.2781 баллов,&nbsp;14 задач (25-е место).<br /><br /><span style="font-size: large;">Собственно задание</span><br /><br /> <a name="cutid1"></a>Дан космический аппарат в поле тяжести Земли, и надо им управлять так, чтобы наиболее полно выполнить заданную миссию, причем побыстрее и потратив поменьше топлива. Худший результат достигается, если аппарат врезается в Землю или расходует всё топливо -- тогда за раунд дают ноль баллов. Вот перечень заданий-миссий:<br /><ol><li>Перейти с одной круговой орбиты на другую круговую. (Здесь организаторы облажались, и большее количество очков получали те, кто сжигал <u>больше</u> топлива. Исправлять они отказались, ибо &quot;так даже прикольней получилось&quot;.)</li><li>То же самое, но надо еще состыковаться с другим спутником на целевой орбите.</li><li>Обобщение (2): начальная и целевая орбиты эллиптические.</li><li>(Задание было выложено по прошествии суток -- после завершения lightning-раунда.) Нужно было собирать космический мусор -- пачка спутников на различных орбитах. Чтобы жизнь малиной не казалась, добавили заправочную станцию на низкой орбите и Луну. (Опять же без багов со стороны организаторов не обошлось -- Луна не действовала на наш КА, только на мусор. Из-за этого было нереально поймать спутник на окололунной орбите. Но его всё-таки исправили.)</li></ol><span style="font-size: large;">Организация</span><br /><br />Найдите десять отличий от ICFPC-2008 (там надо было управляться с марсоходом). ;)<br /><br />Моя критика в адрес организаторов состоит в основном в том, что была дана известная и скорей всего давно решеная задача -- те, кто в теме и &quot;держит в руках&quot; соответсвующие алгоритмы, могли бы довольно быстро и без особого труда всё решить. Ну разве что оставлось бы место для творчества в четвертом задании -- нужно было придумать стратегию очередности обхода спутников и заправки. В этом свете не слишком впечатляющие результаты команды ИПМ выглядят особенно бледно. В наше оправдание замечу, что некоторая часть аппаратов серии <a href="http://ru.wikipedia.org/wiki/%D0%9B%D1%83%D0%BD%D0%B0_%28%D0%BA%D0%BE%D1%81%D0%BC%D0%B8%D1%87%D0%B5%D1%81%D0%BA%D0%B0%D1%8F_%D0%BF%D1%80%D0%BE%D0%B3%D1%80%D0%B0%D0%BC%D0%BC%D0%B0%29">Луна</a> была утеряна из-за расчетных ошибок.<br /><br />Тем, кто внимательно и до конца прочитал условие, обломился бонус: ссылка на википедийную статью про <a href="http://en.wikipedia.org/wiki/Hohmann">маневр Хохманна</a>. (На самом деле в русской транскипции Hohmann вовсе не Хохманн, а Гоман -- мы потом уже нашли советскую книгу докосмической эры (!), где было написано по теме почти всё.) В этой статье есть все три формулы, чтобы без напряженя мозга решить первое и второе задания. (Ну на самом деле не совсем так: точности не всегда хватало, поэтому нужен был еще алгоритм finetune для локальной доводки.)<br /><br /><span style="font-size: large;">Сам процесс</span><br /><br />В течение всего контеста команда почти полностью была&nbsp; в одном месте в реале. Это хороший бонус. Благодаря удачному распределению сил нашей команде удалось быстро реализовать (1) виртуальную машину, (2) маневр Гомана и (3) наивный finetune, что позволило в первую же ночь закомитить пристойные результаты для первого и второго задания.<br /><br />На второй день я научился аналитически стоить эллипс по трем точкам.Очень удобно уравнение эллипса в полярных координатах:<br /><img src="http://www.codecogs.com/eq.latex?%5Crho=%5Cfrac%7Bp%7D%7B1-e%5Ccos%28%5Cphi-%5Calpha%29%7D" alt="" /><br />Здесь фокус эллипса находится в начале координат, как раз там, где у нас Земля. (По сравнению с <a href="http://ru.wikipedia.org/wiki/%D0%AD%D0%BB%D0%BB%D0%B8%D0%BF%D1%81#.D0.A3.D1.80.D0.B0.D0.B2.D0.BD.D0.B5.D0.BD.D0.B8.D0.B5_.D0.B2_.D0.BF.D0.BE.D0.BB.D1.8F.D1.80.D0.BD.D1.8B.D1.85_.D0.BA.D0.BE.D0.BE.D1.80.D0.B4.D0.B8.D0.BD.D0.B0.D1.82.D0.B0.D1.85">канонической формулой</a> здесь добавлен угол alpha между большой полуосью эллипса и осью координат.) Таким образом, любая эллиптическая орбита задается тремя параметрами.<br /><br />Потом, <span class='ljuser ljuser-name_yurikl' lj:user='yurikl' style='white-space: nowrap;'><a href='http://yurikl.livejournal.com/profile'><img src='http://l-stat.livejournal.com/img/userinfo.gif' alt='[info]' width='17' height='17' style='vertical-align: bottom; border: 0; padding-right: 1px;' /></a><a href='http://yurikl.livejournal.com/'><b>yurikl</b></a></span> научился (тоже аналитически) строить эллипс, касательный двум данным эллипсам (правда, с некоторыми ограничениями), а мне удалось допилить алгоритм локальной доводки до почти идеального.<br /><br />Дальше пошли танцы с бубнами: всю эту машинерию так просто применить не всегда удается из-за потери точности (особенно на высоких орбитах), а также из-за влияния Луны в четвертом задании. Численно решать ограниченную задачу трех тел как-то не хотелось (точнее,&nbsp;обратную к ней), хотя наверно, если покопаться, нужные алгоритмы можно было найти. Поэтому решили считать Луну всего лишь возмущением для эллиптических орбит.<br /><br />Лирическое отсупление -- попутно нашел прикол в википедийной статье про <a href="http://ru.wikipedia.org/wiki/%D0%97%D0%B0%D0%B4%D0%B0%D1%87%D0%B0_%D1%82%D1%80%D0%B5%D1%85_%D1%82%D0%B5%D0%BB">задачу трех тел</a>: <div style="margin-left: 40px;">[...] Это подтолкнуло Пуанкаре и Зундмана искать решение не в виде функций от <span class="texhtml"><i>t</i></span>, а в виде рядов от некоторого параметра. [...] Эти ряды были найдены [...] К несчастью, как показал Д. Белорицкий, по крайней мере в случае Лагранжа для нужд вычислительной астрономии в сходящихся рядах Зудмана <u>нужно брать как минимум</u> <img src="http://www.codecogs.com/eq.latex?10%5E%7B8%5Ccdot%7B10%5E6%7D%7D" alt="" /> членов.</div><div style="margin-left: 40px;">&nbsp;</div><br />В общем, на Луну мы даже и не пытались летать, и занимались всякой эвристической фигней типа подбора констант. И довольно быстро скатились в таблице с третьих мест до двадцатых. Поскольку в последние четыре часа мы комитили мало, ожидается, что окончательный результат будет хуже предварительного (25-е место).<br />Да, можно отметить такой полезный,&nbsp;но простой как валенок,&nbsp;хак: если на очередной маневр не хватает топлива, и неясно, что делать,&nbsp;то надо тупо отключить двигатели и ждать окончания раунда, надеясь не врезаться в Землю. Тогда исключается возможность получить ноль баллов за опустошенный бак.<br /><br /><span style="font-size: large;">Замечания на будущее</span><br /><br /> Язык программирования Java показался вполне удовлетворительным. Правда к концу мероприятия бОльшая часть кода превратилась в плохочитаемую лапшу. В связи с этим надо подумать, как с этим бороться. Наверно в определенный момент надо было потратить время на рефакторинг. Тем более, что идея использовать continuation passing style мне пришла в голову не в последний момент,&nbsp;а когда еще можно было задешево переоформить имеющееся хозяйство. Задним числом видно, что это было бы полезно:&nbsp;появилась бы так необходимая модульность. Конечно, нет повода не искать и другие языки.<br /><br />Полезно заранее подготовить инструментарий. Об этом говорилось уже сто раз, но опять наступаем на грабли: к началу контеста выяснилось,&nbsp;что я потерял свой доступ к репозиторию,&nbsp;и ни у кого нет админского доступа. (К слову сказать, я умудрился забыть дома компьютер.) При этом не было защиты от падения интернета и всё хозяйство хранилось на удаленном сервере. Чудесным образом проблем с инетом не было.<br /><br />Да, еще мы так и не нашли удовлетворительного средства для реализации блога+вики+чего-еще-надо в одном флаконе. Хотя заранее потратили изрядно времени. Это отдельная тема для обсуждения.<br /><br />В очередной раз напоминаю (в первую очередь себе!), что аналитические методы на ICFPC -- это зло.&nbsp;Конечно, приятно аналитически решить задачу и получить удовольствие, но всегда есть риск, что это не удастся сделать за разумное время.&nbsp;Таким образом,&nbsp;правильная стратегия -- иметь в запасе набор алгоритмов разной степени универсальности для решения следующих типовых задач: решение системы алгебраических уравнений, поиск экстремумов,&nbsp;поиск путей в графах и т.д. Если задача не берется аналитически за небольшое время (полчаса-час), то сразу переходить на численные методы. Уж потом можно распараллелиться и пытаться заменить плохое численное решение на хорошее аналитическое.<br /><br />Хорошо бы заранее озаботиться готовкой правильной еды. Для этого не жалко завести специального члена команды.&nbsp;;)&nbsp;А то всякие пиццы, пусть даже и горячие, меня совсем не радуют... Ну и надо подумать о режиме сна. То ли спать понемногу, но часто,&nbsp;то ли наоборот, не знаю. А то я часто ловил себя на том,&nbsp;что мозг уже не соображает, и только зря время тратится.<br /><br /><strong>Upd</strong>. Да, для подобных задач обязательно иметь под руками библу для работы с комплексными числами. А то по какому-то недоразумению в стандартной яве ее нет. В результате и в этом контесте, и в прошлогоднем дело было так: вначале всё писалось на отдельных координатах, потом в какой-то момент меня заколебывало вспоминать, как сторить перпендикуляр или уравнение касательной,&nbsp;и сначала вся аналитика, а потом и программа переделывалась под использование комплексных чисел. Вероятно, из-за этого немного падает эффективность, но ведь продуктивность гораздо важнее!<br /><br />PS. Вот еще один <a href="http://users.livejournal.com/_adept_/94839.html">отчет</a> для сравнения.<br /> urn:lj:livejournal.com:atom1:janatem:58518 2009-06-25T07:20:08Z 2009-06-25T07:20:08Z Насчет того, как <a href="http://janatem.livejournal.com/55881.html">сделать проигрыватель на рутере</a>:<br /><div style="margin-left: 40px;"><em>К сожалению,&nbsp;я обнаружил, что выбор usb-звуковушек совсем не велик;&nbsp;кроме того,&nbsp;всё что сейчас есть в магазинах -- очень наворочено и поэтому слишком дорого.</em></div><br />На самом деле правильные устройства конечно же есть, например некая <a href="http://mightyohm.com/blog/2008/10/inside-the-syba-sd-cm-uaud-usb-stereo-audio-adapter/">SYBA SD-CM-UAUD</a>. И стоит она по 50 копеек пучок (точнее $8). Только это за морем телушка -- полушка... А в российских магазинах я не нашел даже упоминания о ней.<br /><h2>&nbsp;</h2> urn:lj:livejournal.com:atom1:janatem:58287 2009-06-24T13:59:41Z 2009-06-24T13:59:41Z Почитывал описание правил фишеровских шахмат в русской википедии и заметил, что приведенная там <a href="http://ru.wikipedia.org/w/index.php?title=Случайные_шахматы&amp;oldid=15406093#.D0.9D.D0.B0.D1.87.D0.B0.D0.BB.D1.8C.D0.BD.D0.B0.D1.8F_.D0.BF.D0.BE.D0.B7.D0.B8.D1.86.D0.B8.D1.8F">процедура выбора начальной расстановки</a> явно кривая.<br /><br /><div style="margin-left: 40px;">...<br />Случайным образом выбрать вертикаль для белого короля из диапазона b1-g1.<br />Случайным образом выбрать позиции для ладей: одну из диапазона a1-позиция короля, другую из диапазона от позиции короля до h1.<br />Из оставшихся белых клеток первой горизонтали выбрать случайную позицию для белопольного слона.<br />Из оставшихся чёрных клеток первой горизонтали выбрать случайную позицию для чёрнопольного слона.<br />...</div><div><br />Невооруженным взглядом видно, что, хотя все 960 позиций достижимы, они выпадают не равновероятно. Если король попадет на b1 или g1, то вариантов расположения ладей будет 6*1=6, а если попадет в центр -- то 3*4=12. (Здесь не учитывается ограничение на разнопольность слонов, но всё равно понятно, что будет криво.)<br />&nbsp;</div>Полез разбираться, что пишут в английской версии,&nbsp;и нашел целый список процедур установки. Там были и правильные (например, <a href="http://en.wikipedia.org/wiki/Chess960_starting_position#Bodlaender.27s_dice-rolling_method">бодландеровская</a>), и кривые, в том числе попавшая в русскую (<a href="http://en.wikipedia.org/wiki/Chess960_starting_position#Rovida.27s_dice-rolling_method">by Ravida</a>).<br /><br /> В связи с этим возникает вопрос: зачем существуют кривые способы, если они ничем не лучше правильных (ну разве что вместо 6-7 бросков кубика достаточно 5-6)? Причем кривых способов много разных. Неужели люди,&nbsp;пользуясь ими, не чувствуют что что-то не так?!<br /> urn:lj:livejournal.com:atom1:janatem:58021 2009-06-17T11:55:04Z 2009-06-17T11:58:34Z Мы с Хеленкой отметились посещением Dream Theater в Лужниках. Кратко о впечатлениях. Никакой особой программы, тем более концепта, у них не было,&nbsp;просто подборка из разных альбомов, ну и, само собой, новенькое из не вышедшего пока альбома. Да, конечно же были и <strike>пьяные вариации</strike> импровизации, что характерно для DT. :)<br /><br />Вот только организаторы подкачали. Звук бы отвратен. Еще для такого огромного зала просто необходимы были экраны. Ну и цены на билеты -- в прошлом году у меня был шанс за ту же сумму попасть на _трехдневный_ фестиваль Graspop, где DT был лишь одним из многих пунктов программы. Хотя цены всё же можно назвать рыночно адекватными, поскольку зал был полон несмотря ни на что.<br /><br /> Я уже давно понял, что не стоит ходить на монстров, собирающих большие залы (прямо-таки стадионы на десятки килозрителей). Эти стадионы и спортивные залы не заточены под концерты, в них неизбежно поганый звук, и даже на хороших местах нет ощущения присутствия. (Однако ж опять наступаю на старые грабли.) Хочется чего-нибудь более камерного.<br /><br />В связи с этим я придумал новый жанр развлечений, дешевый и бесхитростный. Вот есть театр -- живая игра,&nbsp;а есть кино -- мертвая. Предлагаю аналогично жанр мертвой музыки: пусть будет маленький зальчик на 20-100 чел с хорошей аудиоаппаратурой и акустически правльной геометрией. Важно, чтобы затея не выродилась в обычный бар с фоновой музыкой;&nbsp;т.е. должно быть расписание,&nbsp;что и в какое время дают. Такие зальчики могут находиться в любых местах города,&nbsp;так чтобы людям было удобно, например, после работы потусоваться,&nbsp;выпить пивка и заценить новый (или наоборот)&nbsp;альбомчик любимых авторов.<br /><br /> urn:lj:livejournal.com:atom1:janatem:57708 2009-06-02T11:29:38Z 2009-06-02T11:29:38Z Настройка женского линукса:<br /><br /><div style="margin-left: 40px;">$ echo 'alias lady=man' &gt;&gt; ~/.bashrc</div> urn:lj:livejournal.com:atom1:janatem:57527 2009-05-24T21:02:57Z 2009-05-24T21:04:21Z Кажется, так линукс еще никто не опускал. ;)<br /><br />Мне звонят из магазина:<br />-- Вы заказывали ноутбук такой-то?<br />-- Да.<br />-- Доставка туда-то в такое-то время. Так?<br />-- Да, всё правильно.<br />-- Только вы знаете, что он <em>без операционной системы</em>?<br />-- Как?! На сайте было написано, что там линукс установлен.<br />-- Да, но там нет графики. И вообще...<br />-- Я знаю, что такое линукс. [Только не втыкаю, почему &quot;нет графики&quot;, и почему бы ясно не сказать, какой там дистрибутив.]<br />-- Ну ладно, я вас предупредил.<br /><br />Видимо следует считать, что в данном контексте linux = free dos = &quot;просто бесплатная ОС&quot;. Просто для того, чтобы ноутбук мог при включении &quot;помигать лампочками&quot;. Может быть по закону запрещено продавать компьютеры без предустановленной ОС?