Блог пользователя Burunduk1

Автор Burunduk1, 13 лет назад, По-русски

** UPD : Считается, что все формулы уже исправлены **

163A - Подстрока и подпоследовательность

Кратко о решении: динамика.

Одно из авторских решений: 1415300 (автор = levlam)

Задача решалась следующей динамикой. Пусть f[i, j] — количество равных пар вида "подстрока начинающаяся в позиции i" и "подпоследовательность подстроки t[j..|t|]".

Тогда:

f[i, j] = f[i, j + 1];
if (s[i] == t[j])
  add(f[i, j], f[i + 1, j + 1] + 1)

Answer = f[i,0]

163B - Лемминги

Кратко о решении: сортировка + бинпоиск по ответу.

Одно из авторских решений: 1415306 (автор: Burunduk1)

Нам нужно найти минимальное время T. Давайте сделаем бинарный поиск по ответу.

После этого можно расставлять леммингов жадно или снизу вверх, или сверху вниз. В этом разборе я воспользуюсь способом <снизу вверх>. Из тех леммингов, которые могут забраться на 1-й уступ, выберем лемминга с минимальной массой, а из таких — с минимальной скоростью. Почему так можно сделать? Пусть мы поставили на первый уступ лемминга с большой массой, значит, всех леммингов с меньшей массой мы использовать уже не можем, значит, можно было поставить на первый уступ любого лемминга с меньшей массой. А из леммингов с одинаковой массой всегда выгодно оставить тех, чья скорость больше, они могут залезть на более высокие уступы.

Чтобы быстро расставлять леммингов в таком порядке, отсортируем их заранее, сравнивая сперва по массе, затем по скорости. После этого мы рассматриваем всех леммингов в порядке возрастания массы ровно один раз и каждого или берем, или не берем. Второй раз рассматривать лемминга не имеет смысла, т.к. или мы его уже взяли на i-й уступ (он занят) или мы его на i-й уступ взять не смогли, тогда на i + 1-й и т.д. тоже не сможем.

Итого решение = сортировка + бинпоиск по ответу с линейным проходом внутри. Время работы = O(Nlog?)

Но написать бинпоиск и верную жадность в этой задаче было не достаточно для получения Accepted. Нужно было еще не попасться на проблемы с точностью. Давайте точно оценим число итераций бинпоиска.

0 ≤ T ≤ H * K (максимальное время не больше 109). Нужно понять, насколько могут отличаться два времени. В случае "N = K = 105, H = 104, V порядка 109" времена, за которые некоторые два лемминга поднимаются на какие-то уступы, могут отличаться на 10 - 18. Это следует из того, что времена = рациональные дроби , где X и Y порядка 109.

Т.е. нам нужно сделать число итераций = log21027 = 90. На практике моему решению 70 итераций не хватало, а вот 75 было уже достаточно. Рассуждения выше показывают, что 90 точно хватит.

163C - Конвейер

Кратко о решении: события на окружности или бинпоиск.

Одно из авторских решений: 1415310 (автор: Burunduk1)

Одно из авторских решений: 1415316 (автор: arseny30)

Антон, откуда бы он ни начал, пробежит по конвейеру отрезок длины D = (v2 * l) / (v1 + v2). Рассмотрим одну конфету. Чтобы ее скушать, нужно попасть на конвейер в любой момент из отрезка [ ai — D .. ai ]. Рассмотрим все точки вида a[i] — D и a[i] , если какие-то из них отрицательны, прибавим 2 * l. Также не забудем добавить в массив точек числа 0 и 2 * l. Отсортируем эти точки. Рассмотрим теперь две соседние — x[i] и x[i+1]. В какой бы момент времени между этими двумя точками Антон ни начал, он съест одно и то же количество конфет.

С этого момента можно получить одно из двух решений:

1) Запуститься от каждой середины отрезка M = (x[i] + x[i+1]) / 2 и бинпоиском по отсортированному удвоенному массиву a найти количество конфет на отрезке [M..M + D].

2) Сказать, что у нас есть события на окружности вида (a[i]-D,+1) и (a[i],-1). Нужно пройти по кругу два раза, и на втором проходе говорить, что мы знаем, сколько раз покрыт текущая точка конфетами-отрезками.

Оба решения работают за время O(NlogN).

163D - Большой холодильник

Кратко о решении: перебор с отсечением по ответу

Одно из авторских решений: 1415320 (автор: Burunduk1)

Решение состоит из трех главных идей:

  1. V = ABC, A ≤ B ≤ C тогда A ≤ N1 / 3, B ≤ (N / A)1 / 2.

  2. A и B — делители V, а поскольку нам уже дано разложение числа V на простые множители, можно перебирать только делители

  3. Если зафиксировано A, то оптимальные вещественные B и C = (N / A)1 / 2 (обозначим эту величину за X). Т.е. площадь поверхности получится в любом случае  ≥ 2(2AX + X2). Таким образом можно использовать отсечение по ответу.

Если применение в лоб этих идей давало TL, то можно было воспользоваться следующими оптимизациями:

  1. Отсечение по ответу лучше работает, если перебирать A в порядке убывания

  2. A и B — 32-битные целые.

  3. Если мы для числа V уже считали ответ, то нужно запомнить ответ (таким образом мы чуть уменьшили max тест).

Если кому-то интересно теоретическое обоснование "Почему это все работает за нужное время?", приведу некоторую статистику:

  1. Максимальное количество делителей у чисел от 1 до 1018 = 103860 (у числа 897612484786617600 = 2834527211 ...)

  2. Если использовать только 2 первые оптимизации, количество чисел A, которое переберет наше решение = 10471 (на числе с max количеством делителей)

  3. Если использовать только 2 первые оптимизации, количество пар A и B, которое переберет наше решение = 128264 (на числе с max количеством делителей).

163E - Электронное правительство

Кратко о решении: Ахо-Корасик + сумма на пути в дереве суффиксных ссылок

Одно из авторских решений: 1415345 (автор: arseny30)

Предполагается, что читатель знаком с алгоритмом Ахо-Корасик (прочитать о нем можно здесь http://e-maxx.ru/algo/aho_corasick).

Пусть у нас есть бор имен и суффиксные ссылки на этом боре. Также для каждой вершины v посчитано количество имен, кончающихся в этой вершине (end[v]).

Тогда добавить имя i: end[v[i]] += 1, где v[i] — вершина-конец i-го имени

Удалить имя i: end[v[i]] -= 1.

Чтобы ответить на запрос "политизированность текста", нужно увидеть, что суффиксные ссылки образуют дерево. Когда мы проходимся по тексту и параллельно ходим по бору с суффиксными ссылками, если мы находимся в вершине бора v, нужно прибавить к политизированности сумму end[v[i]] на пути по суффиксным ссылкам от v до корня дерева. Теперь нужно научиться быстро решать задачу <сумма весов вершин на пути до корня, веса меняются>.

Это можно сделать например так:

  1. Говорим, что вес не в вершине, а на ребре, соединяющем вершину с ее отцом.

  2. Храним Эйлеров обход дерева, в котором при спуске вниз записан плюс вес ребра, а при подъеме наверх минус вес ребра.

  3. Сумма на пути до корня = сумма на отрезке в Эйлеровом обходе (концы отрезка — произвольные вхождения двух вершин в Эйлеров обход).

Быстрый и удобный способ считать сумму на отрезке — дерево Фенвика.

Мы получили решение работающее за время O(|SumLen| * log). Памяти используется |SumLen| * 26 (бор мы храним не сжатым).

Разбор задач VK Cup 2012 Раунд 2
  • Проголосовать: нравится
  • +42
  • Проголосовать: не нравится

»
13 лет назад, # |
  Проголосовать: нравится +2 Проголосовать: не нравится

Есть такой спецэффект, что решения авторов контеста нельзя посмотреть. Вот совсем. И в списке посылок их тоже нет.

»
13 лет назад, # |
  Проголосовать: нравится 0 Проголосовать: не нравится

Thanks for the editorial. Please provide an English edition. Meanwhile I noticed that the author's submissions aren't accessible at this moment.

»
13 лет назад, # |
  Проголосовать: нравится -8 Проголосовать: не нравится

Я в D придумал решение за (количество делителей * количество различных простых) на тест. Правда набажил и еще не заработало) Как вы думаете, пройдет?

  • »
    »
    13 лет назад, # ^ |
      Проголосовать: нравится 0 Проголосовать: не нравится

    Зная то, как ваша команда умеет запихивать все, конечно:)

  • »
    »
    13 лет назад, # ^ |
      Проголосовать: нравится 0 Проголосовать: не нравится

    Количество делителей — это уже очень много. Примерно 10^5. Мое решение делает 3*10^4 итераций на тест.

»
13 лет назад, # |
  Проголосовать: нравится 0 Проголосовать: не нравится

А почему решение №2 задачи C — это ? Вроде же O(N), у нас 2N событий и 2 прохода.

  • »
    »
    13 лет назад, # ^ |
      Проголосовать: нравится 0 Проголосовать: не нравится

    А отсортировать?

    • »
      »
      »
      13 лет назад, # ^ |
        Проголосовать: нравится 0 Проголосовать: не нравится

      А, ну да.

      • »
        »
        »
        »
        13 лет назад, # ^ |
          Проголосовать: нравится 0 Проголосовать: не нравится

        данные на входе уже отсортированы. и если делать события не в лоб, то сложность будет O(N)

    • »
      »
      »
      13 лет назад, # ^ |
        Проголосовать: нравится +3 Проголосовать: не нравится

      Что сортировать? Данные уже на вход отсортированные даются. У меня чистый O(N). Линейно двигаю отрезок по циклу при помощи метода 2 указателей.

  • »
    »
    13 лет назад, # ^ |
      Проголосовать: нравится 0 Проголосовать: не нравится

    а сортировка?

  • »
    »
    13 лет назад, # ^ |
      Проголосовать: нравится 0 Проголосовать: не нравится

    ну их еще упорядочить надо. события то. хотя, можно их сливанием упорядоченных массивов посортить... тогда будет O(N). но зачем усложнять себе жизнь?

    • »
      »
      »
      13 лет назад, # ^ |
        Проголосовать: нравится 0 Проголосовать: не нравится

      Ну я, например, решал с помощью двух указателей. Хотя алгоритмически это то же самое, что и слияние, только результат нигде не хранится. Преимущество в том, что не нужно создавать отдельные объекты для событий.

      • »
        »
        »
        »
        13 лет назад, # ^ |
          Проголосовать: нравится 0 Проголосовать: не нравится

        да, указателями тоже можно. когда я писал комментарий, идеи с указателями мне в голову не приходило (в том числе и из других комментариев).

        я лично решал сортировкой событий за NlogN.

  • »
    »
    13 лет назад, # ^ |
      Проголосовать: нравится +2 Проголосовать: не нравится

    К сожалению, нет время описывать подробно. Мое решение вроде как O(N): http://mirror.codeforces.com/contest/163/submission/1414357 Мне казалось, что все так решали и не надо никаких O(NlogN)... :)

»
13 лет назад, # |
Rev. 2   Проголосовать: нравится 0 Проголосовать: не нравится

In problem B " So we need to make log (10^27) = 90 iterations " where from 10^27 came ?.....i just want to know how you concluded that 75 iterations will be enough.

  • »
    »
    13 лет назад, # ^ |
      Проголосовать: нравится 0 Проголосовать: не нравится

    You have a range of doubles [0, 1e9] and two doubles here are considered "equal" if their substraction is less than 1e-18. Then you actually have 1e27 "different" double values, don't you?

»
13 лет назад, # |
Rev. 2   Проголосовать: нравится +16 Проголосовать: не нравится

На мой взгляд, разбор задачи B (про леммингов) не полон.

Максимальный ответ равен 10^9, а расстояние между оптимальным ответом и следующим за ним не меньше 10^-18. Таким образом, если у вы используете числа с как минимум 27 десятичными знаками относительной точности, можно делать 90 итераций и быть спокойным...

На контесте я сдал (на ACC) решение с double (16 десятичных знаков относительной точности) и очень боялся, что оно упадет по точности. Очень странно, что рассуждения на тему ТОЧНОСТИ (а не количества итераций) отсутствуют в разборе жюри.

Как я понимаю, нужно доказать, что если оптимальный ответ равен A, а какой-то другой ответ равен B, то |B-A|/A >= 10^-16. Если это неверно, то многие засчитанные решения (в том числе мое) могут быть убиты хаком.

  • »
    »
    13 лет назад, # ^ |
    Rev. 4   Проголосовать: нравится +11 Проголосовать: не нравится

    Вообще, подобное рассуждение в контест напугало меня настолько, что я придумал неправильное "точное" решение и отправил его.

    А вот как на самом деле. Ответом может быть некоторая дробь j/v[i], 1 <= j <= k. Параметр h вообще в задаче лишний, примем его за единицу. Пусть у нас есть i1/v1 и i2/v2 как возможные ответы задачи. Относительное отличие 1 - (v1/v2)*(i2/i1). Тогда, если взять i1 = A, i2 = A+D, v1 = B-C, v2 = B, получим 1 - ((1-C/B)*(1+D/A)) = C/B-D/A+(CD)/(AB), что, вроде, положительным меньше 10^(-9)*10^(-5)=10^(-14) не бывает.

  • »
    »
    13 лет назад, # ^ |
    Rev. 2   Проголосовать: нравится 0 Проголосовать: не нравится

    Таким образом, если у вы используете числа с как минимум 27 десятичными знаками относительной точности, можно делать 90 итераций и быть спокойным…

    А я правильно понимаю, что среди стандартных в С такого типа не предусмотрено? В доках GCC __float128 обзывают гнутым расширением.

    • »
      »
      »
      13 лет назад, # ^ |
        Проголосовать: нравится 0 Проголосовать: не нравится

      Да, такого типа в большинстве языков нет.

      Но это теперь неважно, поскольку I_love_natalia уже показал, что достаточно 14 десятичных знаков относительной точности, то есть обычного double хватает=)

»
13 лет назад, # |
  Проголосовать: нравится 0 Проголосовать: не нравится

А можно более подробно как-нибудь описать решение первой задачи. А то, что-то не совсем доходит какая там динамика. f[i][j] — количество пар подстрок что строка S начинается с позиции i а строка T, не сильно понимаю. (Заранее благодарен).

  • »
    »
    13 лет назад, # ^ |
      Проголосовать: нравится +2 Проголосовать: не нравится

    Попробую(правда, у меня слегка другое решение). Итак, пусть f[i][j] = число равных пар "подстрока-подпоследовательность" таких, что подстрока закончилась ровно в символе i, подпоследовательность закончилась не правее символа j. (Всё инициализируется нулями). Тогда:
    1)s[i] != t[j]
    f[i][j] = f[i][j — 1] (поскольку ясно, что сделать иначе мы не можем)
    2)s[i] == t[j]
    f[i][j] = (f[i][j — 1] + 1 + f[i — 1][j — 1]) % MOD
    (первое слагаемое — см. случай 1, второе слагаемое — новые подстрока и подпоследовательность длины 1, третье слагаемое — возьмём все "старые" пары и припишем новый символ).
    Как "собрать" ответ? Очевидно, что искомые подстроки заканчиваются на любом символе, а подпоследовательности — гарантированно не правее, чем на самом правом сиволе. Таким образом,
    for(int i = 0; i < n; i++){
    ans += f[i][m — 1];
    ans %= MOD; }

»
13 лет назад, # |
  Проголосовать: нравится 0 Проголосовать: не нравится

А будет разбор задач для Div.2?

  • »
    »
    13 лет назад, # ^ |
      Проголосовать: нравится 0 Проголосовать: не нравится

    Врядли

    • »
      »
      »
      13 лет назад, # ^ |
      Rev. 2   Проголосовать: нравится 0 Проголосовать: не нравится

      Прикольно. Обычный рейтинговый контест — а разбора не будет? Не в духе codeforces.
      Они даже для Testing Round делают разборы)

      • »
        »
        »
        »
        13 лет назад, # ^ |
        Rev. 2   Проголосовать: нравится 0 Проголосовать: не нравится

        думаю для А и В (Div2) вполне достаточно посмотреть решения других участников.

  • »
    »
    13 лет назад, # ^ |
      Проголосовать: нравится 0 Проголосовать: не нравится

    А что тут разбирать та? A) отсортировать массив и вывести h[n — a] + h[n — a — 1]; B) пройтись по цифрам числа а и пытаться жадно увеличить каждую цифру.

    ПС: Получил на бинпоске по С WA, когда сделал 100 итерации получил АС, вместо того чтобы стать фиолетовым, стал зеленым Т_Т

  • »
    »
    13 лет назад, # ^ |
      Проголосовать: нравится +5 Проголосовать: не нравится

    Может быть, но точно не от меня :-) Я даже условий не читал.

»
13 лет назад, # |
Rev. 3   Проголосовать: нравится 0 Проголосовать: не нравится

please help!! i tried to submit my solution to the problem E many times, but my solution couldn't be received. the screen says "Ooops! Something has broken down in Codeforces. Do not panic, you can try to reload the page or return Home. Anyway we will carefully read megabytes of logs, analyze stacktraces and fix the problem.". i tried submitting a empty-code to it,and the code was received(of course, wrong answer). and i tried log-out and log-in,but the situation didn't improve. ??why?? my code is no compile-error, or memory limit exeed. also, i tried the custom test,and the custom test was no problem.

thank you for fixing the problem ^_^ my cause of trouble is using the "%lld". now,the program doesn't be received,too,however,the screen told me a proper warning.

»
12 лет назад, # |
  Проголосовать: нравится 0 Проголосовать: не нравится

Мой ответ на тест задачи В циферка в циферку похож на ответы тех, у кого полное решение, но почему-то мой не проходит. Как так может быть?

  • »
    »
    12 лет назад, # ^ |
      Проголосовать: нравится 0 Проголосовать: не нравится

    В чем-то разница точно есть ;-)

    Возможно, ответы даны с символом "..." посередине. Значит, часть пропущена. Возможно, ошибка именно там.

»
10 лет назад, # |
  Проголосовать: нравится -8 Проголосовать: не нравится

Could someone explain me, why mathematicaly 90 iterations is enough for 163B — Lemmings? I have done a solution with 40

»
7 лет назад, # |
  Проголосовать: нравится 0 Проголосовать: не нравится

can someone please elaborate the dp states and transition of A.

»
6 лет назад, # |
Rev. 8   Проголосовать: нравится 0 Проголосовать: не нравится

In problem 163A - Substring and Subsequence, should the formula for Answer be the summation over i of f[i,|t|] instead of f[i,0]? 39291217

Update: The answer is NO! The editorial code reverses the direction of the j loop, and the original formula is correct.

»
2 месяца назад, # |
  Проголосовать: нравится 0 Проголосовать: не нравится

Why is the time limit so tight for E problem, I tried a suffix array approach and a lazy segtree to maintain the contributions, but got TLE at testcase 4.