Решение цветных японских кроссвордов со скоростью света

    Японские кроссворды (также нонограммы) — логические головоломки, в которых зашифровано пиксельное изображение. Разгадывать кроссворд нужно с помощью чисел, расположенных слева от строк и сверху от столбцов.


    Размер кроссвордов может доходить до 150x150. Игрок с помощью специальных логических приемов вычисляет цвет каждой клетки. Решение может занять как пару минут на кроссвордах для начинающих, так и десятки часов на сложных головоломках.


    Хороший алгоритм может решить задачу намного быстрее. В тексте описано, как с помощью наиболее подходящих алгоритмов (которые вообще приводят к решению), а также их оптимизаций и использования особенностей C++ (которые уменьшают время работы в несколько десятков раз) написать решение, работающее почти мгновенно.



    В мобильной версии Хабра формулы в тексте могут не показываться (не во всех мобильных браузерах) — пользуйтесь полной версией или другим мобильным браузером, если заметили проблемы с формулами


    Правила игры


    Изначально холст кроссворда белый. Для ванильных черно-белых кроссвордов нужно восстановить местоположения черных клеток.


    В черно-белых кроссвордах количество чисел для каждой строки (слева от холста) или для каждого столбца (сверху от холста) показывает, сколько групп черных клеток находятся в соответствующих строке или столбце, а сами числа — сколько слитных клеток содержит каждая из этих групп. Набор чисел $[3, 2, 5]$ значит, что в определенном ряду есть три последовательные группы из $3$, $2$ и $5$ черных клеток подряд. Группы могут быть расположены в ряду как попало, не нарушая относительный порядок (цифры задают их длину, а позицию надо угадать), но они обязательно должны разделяться хотя бы одной белой клеткой.



    Корректный пример


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


    Что не является японским кроссвордом


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


    Такие кроссворды являются не кроссвордами, а графоманией. Считается, что корректный кроссворд — такой, что логическим путем можно прийти к единственному решению.


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



    Неправильная нонограмма — решение единственное, но решить нормально нельзя. Оранжевым помечены "нерешаемые" клетки. Взято из Wikipedia.


    Такое ограничение объясняется так — в самом общем случае японский кроссворд это NP-полная задача. Однако, NP-полной задачей разгадывание не становится, если есть алгоритм, в каждый момент времени однозначно указывающий, какие клетки открыть далее. Все методы разгадывания кроссвордов, применяемые человеком (за исключением метода Монте-Карло проб и ошибок), основываются именно на этом.


    У наиболее православных кроссвордов ширина и высота делится на 5, нет рядов, которых можно посчитать мгновенно (такие, где либо цветные клетки забивают все места, либо их нет совсем), и ограничено количество дополнительных цветов. Эти требования не обязательные.


    Наипростейший неправильный кроссворд:


      |1 1|
    --+---+
     1|   |
     1|   |
    --+---+

    Часто не решаются взад закодированные пиксельные арты, в которых используется "шахматный порядок" для имитации градиента. Лучше понять будет, если вы наберете в поиске "pixelart gradient". Градиент как раз похож на этот фейловый кроссворд 2x2.



    Возможные варианты решений


    У нас есть размер кроссворда, описание цветов и всех строк и столбцов. Как можно собрать из этого картинку:


    Полный перебор


    Для каждой клетки перебираем все возможные состояния и проверяем на удовлетворительность. Сложность такого решения для черно-белого кроссворда $O(N*M*{2}^{N*M})$, для цветного $O(N*M*{colors}^{N*M})$. По аналогии с клоунской сортировкой это решение можно тоже назвать клоунским. Оно годится для размера 5x5.


    Backtracking


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


    Отдельно для каждого ряда


    Это решение намного лучше и оно истинно верное. Мы можем проанализировать каждую строку и каждый столбец по отдельности. У каждого ряда попытаемся раскрыть максимум информации.


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


    Истинно верное решение


    Одна строка, два цвета


    Эффективное отгадывание черно-белого "однострочника", для которого некоторые клетки уже отгаданы — весьма жесткая задача. Она встречалась в таких местах, как:


    • Четвертьфинал ACM ICPC 2006задачу можно попробовать решить самому. Тайм-лимит 1 секунда, ограничение количества групп 400, длины ряда тоже 400. Имеет сильно высокий уровень сложности по сравнению с другими задачами.
    • International Olympiad in Informatics 2016условие, сдать задачу. Тайм-лимит 2 секунды, ограничение кол-ва групп 100, длины ряда 200'000. Такие ограничения оверкилл, но решение задачи с ограничениями ACM набирает 80/100 баллов в этой задаче. Тут тоже уровень не подкачал, школьники со всего мира с жестоким уровнем IQ тренируются несколько лет решать разную жесть, потом проходят на эту олимпиаду (пройти должны только 4 человека от страны), решают 2 тура по 5 часов и в случае epic win (бронза в разные годы за 138-240 баллов из 600) поступление в Оксфорд, потом офферы от понятных компаний в отдел поиска, долгая и счастливая жизнь в обнимку с дакимакурой.

    Монохромный однострочник тоже можно решать по-разному, и за $O(N*2^N)$ (перебор всех вариантов, проверка на корректность, выделение клеток, которые имеют один и тот же цвет во всех вариантах), и еще как-то менее тупо.


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


    Псевдокод
    def EpicWin(group, cell):                                                          
        if cell >= last_cell: # Условие выигрыша                                       
            return group == group_size                                              
    
        win = False                                                                 
    
        # Можем ли вставить группу в этой позиции                                   
        if group < group_size  # Группы еще есть                                    
                and CanInsertBlack(cell, len[group])  # Вставка черной группы       
                and CanInsertWhite(cell + len[group], 1)  # Вставка белой клетки    
                and EpicWin(group + 1, cell + len[group] + 1):  # Можно заполнить далее
            win = True                                                              
            can_place_black[cell .. (cell + len[group] - 1)] = True                 
            can_place_white[cell + len[group]] = True;                              
    
        # Можем ли вставить белую клетку в этой позиции                                
        if CanInsertWhite(cell, 1)                                                     
                and EpicWin(group, cell + 1):                                          
            win = True                                                                 
            can_place_white[cell] = True                                               
    
        return win
    
    EpicWin(0, 0)        

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


    Функции CanInsertBlack/CanInsertWhite нужны, чтобы проверить, можно ли теоретически поставить группу нужного размера в нужное место. Все что им надо сделать — проверить, что в указанном интервале клеток нет "100% белой" клетки (для первой функции) или "100% черной" (для второй). Значит, они могут работать за $O(1)$, это можно сделать с помощью частичных сумм:


    CanInsertBlack
    white_counter = [ 0, 0, ..., 0 ]  # Массив длины n                              
    
    def PrecalcWhite():                                                             
        for i in range(0, (n-1)):                                                        
            if is_anyway_white[i]:  # 1, если i-ая клетка 100% белая                
                white_counter[i]++                                                  
            if i > 0:                                                               
                white_counter[i] += white_counter[i - 1]                            
    
    def CanInsertBlack(cell, len):                                                  
        # Опускаем проверку корректности границ [cell, cell + len - 1]              
        ans = white_counter[cell + len - 1]                                         
        if cell > 0:                                                                
            ans -= white_counter[cell - 1]                                          
        # В ans - количество белых клеток от cell до (cell + len - 1)               
        return ans == 0        

    Такое же колдунство с частичными суммами ждет строки вида


    can_place_black[cell .. (cell + len[group] - 1)] = True

    Тут можно вместо = True увеличивать число на 1. А если нам надо произвести много прибавлений на отрезке в неком массиве $array$, и притом этот массив мы никак не используем перед разными прибавлениями (про такое говорят, что эта задача "решается оффлайн"), то вместо цикла:


    # Много раз такой код
    for i in range(L, R + 1):
        array[i]++

    Можно сделать так:


    # Много раз такой код
    array[L]++
    array[R + 1]--
    # После всех прибавлений
    for i in range(1, n):
        array[i] += array[i - 1]

    Таким образом, работает весь алгоритм за $O(k*n)$, где $k$ — количество групп черных клеток, $n$ — длина строки. И мы проходим на полуфинал ACM ICPC или получаем бронзу межнара. Решение ACM (Java). Решение IOI (C++).


    Одна строка, много цветов


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



    Источник — Методы решения японских кроссвордов


    В то время как двухцветные нонограммы нормально обходились без этого, им не надо было оглядываться на ортогональных друзей.
    На картинке видно, что у левого примера три крайние правые клетки пустые, потому что поломается конфигурация, если окрасить эти клетки в те цвета, в которых окрасить себя не-белый цвет.


    Но этот прикол математически разрешим — надо каждой клетке выдать число, где $i$-й бит будет означать, можно ли дать этой клетке $i$-й цвет. Изначально у всех клеток значение $2^{colors} - 1$. Решение динамики поменяется не очень сильно.


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


    Если считать нулевой бит "белым", первый "синим", второй "зеленый", то строка вычислила для последней клетки состояние $011$, а столбец $101$. Значит, реально у этой клетки состояние $011\&101 = 001$


    Псевдокод
    source = [...]  # Начальные состояния
    result = [0, 0, ..., 0]  # Итоговые состояния
    len = [...]  # Длины групп клеток
    clrs = [...]  # Цвета групп клеток
    
    def CanInsertColor(color, cell, len):                                                  
        for i in range(cell, cell + len):
            if (source[i] & (1 << color)) == 0:  # В некоторой клетке цвет color поставить не можем
                return False;
        return True
    
    def PlaceColor(color, cell, len):                                                                                                                      
        for i in range(cell, cell + len):                                           
            result[i] |= (1 << color)  # Добавляем цвет color операцией OR
    
    def EpicWinExtended(group, cell):                                               
        if cell >= last_cell: # Условие выигрыша                                    
            return group == group_size                                              
    
        win = False                                                                 
    
        # Можем ли вставить группу в этой позиции                                   
        if group < group_size  # Группы еще есть                                    
                and CanInsertColor(clrs[group], cell, len[group])  # Вставка черной группы
                and SequenceCheck()  # Если следующая группа имеет тот же цвет, проверяем
                                     # что можем после этой группы поставить белую клетку
                and EpicWin(group + 1, cell + len[group]):  # Можно заполнить далее 
            win = True                                                              
            PlaceColor(clrs[group], cell, len[group])                               
            PlaceColor(0, cell + len[group], 1)  # Белая клетка - только если нужно 
    
        # Можем ли вставить белую клетку в этой позиции                             
        # Белая клетка - бит 0                                                      
        if CanInsertWhite(0, cell, 1)                                               
                and EpicWinExtended(group, cell + 1):                                                                                                      
            win = True                                                              
            PlaceColor(0, cell, 1)                                                  
    
        return win                      

    Много строк, много цветов


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


    Первые результаты


    Допустим, что код мы писали не как дятлы, то есть никуда по ошибке объекты не копируем вместо передачи по ссылке, нигде в алгоритме не косячим, велосипедов не изобретаем, для битовых операций используем __builtin_popcount и __builtin_ctz (особенности gcc), все аккуратно и чисто.


    Посмотрим на время работы программы, которая считывает из файла головоломку и решает ее полностью. Стоит оценить достоинства машинки, на которой все это добро писалось и тестировалась:


    Спецификация движка моего 1932 Harley Davidson Model B Classic Motorcycle
    RAM - 4GB
    Процессор - AMD E1-2500, частота 1400MHz
    Кэш L1 - 128KiB, 1GHz
    Кэш L2 - 1MiB, 1GHz
    Ядер - 2, потоков - 2
    Представлен как «малопроизводительная SoC для компактных ноутбуков» в середине 2013 года
    Стоимость процессора 28 долларов

    Понятно, что такой суперкомпьютер был выбран, потому что оптимизации на нем имеют больший эффект, чем на летающей шайтан-машине.


    Итак, после прогона нашего алгоритма на самом сложном кроссворде (по версии nonograms.ru), получаем не очень хороший результат — от 47 до 60 секунд (в это входит считывание из файла и решение). Надо заметить, что "сложность" на сайте посчитана хорошо — этот же кроссворд во всех версиях программы так же был самым тяжелым, другие наиболее сложные кроссворды по мнению архива держались в топе по времени.


    Цветной кроссворд №9596, наибольшая сложность

    Для быстрого тестирования была сделана функциональность для бенчмарка. Для получения данных для него я специальным скриптом спарсил 4683 цветных кроссворда (из 10906) и 1406 черно-белых (из 8293) с nonograms.ru (это один из крупнейших архивов нонограмм в интернете) и сохранил их в формате, понятном программе. Можно считать, что эти кроссворды являются случайной выборкой, поэтому бенчмарк показал бы адекватные средние значения. Также номера пары дюжин самых "сложных" кроссвордов (также самых больших по размеру и количеству цветов) записал в скрипт для загрузки ручками.



    Оптимизация


    Здесь показаны возможные приемы для оптимизации, которые были сделаны (1)во время написания всего алгоритма, (2)для ужимания времени работы с полминуты до долей секунды, (3)те оптимизации, которые могут быть полезны далее.


    Во время написания алгоритма


    • Специальные функции для битовых операций, в нашем случае __builtin_popcount для подсчета единиц в двоичной записи числа, и __builtin_ctz для количества нулей перед первой самой младшей единицей. Таких функций может не оказаться в некоторых компиляторах. Для Windows подойдут такие аналоги:

    Windows popcount
    #ifdef _MSC_VER                                                                                                
    #  include <intrin.h>                                                                                          
    #  define __builtin_popcount __popcnt                                                                          
    #endif   

    • Организация массивов — меньший размер стоит вначале. Например, лучше использовать массив [2][3][500][1024], чем [1024][500][3][2].
    • Самое главное — общая адекватность кода и избегание лишних забот для вычислений.

    Что уменьшает время работы


    • Флаг -O2 при компиляции.
    • Чтобы не бросать в алгоритм полностью решенную строку/столбец заново, можно в отдельном std::vector/массиве завести флаги для каждого ряда, помечать их при полном решении, и не давать идти дальше, если решать уже нечего.
    • Специфика многоповторного решения задачи на динамику предполагает, что специальный массив, который содержит флаги, помечающие уже "вычисленные" куски задачи, следует обнулять каждый раз при новом решении. В нашем случае это двумерный массив/вектор, где первое измерение — количество групп, второе — текущая клетка (см. псевдокод EpicWin сверху, где этого массива нет, но идея ясна). Вместо обнуления можно сделать так — пусть у нас будет переменная-"таймер", а массив состоит из чисел, показывающих последнее "время", когда этот кусок пересчитывался в последний раз. При запуске новой задачи "таймер" увеличивается на 1. Вместо проверки булевого флага следует проверять равенство элемента массива и "таймера". Это эффективно особенно в тех случаях, когда далеко не все возможные состояния покрываются алгоритмом (а значит, обнуление массива "считали ли мы уже это" занимает здоровый кусок времени).
    • Замена несложных однотипных операций (циклы с присваиванием и т.д.) на аналоги в STL или более адекватные вещи. Иногда работает быстрее велосипеда.
    • std::vector<bool> в С++ сжимает все булевые значения до отдельных битов в числах, что работает при доступе чуть медленнее, чем если бы это было обычное значение по адресу. Если программа ну очень-очень часто обращается к таким значениям, то замена bool на целочисленный тип может хорошо повлиять на производительность.
    • Остальные слабые места можно искать через профайлеры и править их. Я сам использую Valgrind, его анализ производительности удобно смотреть через kCachegrind. Профайлеры встроены во многие IDE.

    Этих правок оказалось достаточно, чтобы получить такие данные на бенчмарке:


    Цветные нонограммы
    izaron@izaron:~/nonograms/build$ ./nonograms_solver -x ../../nonogram/source2/ -e
    [2018-08-04 22:57:40.432] [Nonograms] [info] Starting a benchmark...
    [2018-08-04 22:58:03.820] [Nonograms] [info] Average time: 0.00497556, Median time: 0.00302781, Max time: 0.215925
    [2018-08-04 22:58:03.820] [Nonograms] [info] Top 10 heaviest nonograms:
    [2018-08-04 22:58:03.820] [Nonograms] [info] 0.215925 seconds, file 9596
    [2018-08-04 22:58:03.820] [Nonograms] [info] 0.164579 seconds, file 4462
    [2018-08-04 22:58:03.820] [Nonograms] [info] 0.105828 seconds, file 15831
    [2018-08-04 22:58:03.820] [Nonograms] [info] 0.08827 seconds, file 18353
    [2018-08-04 22:58:03.820] [Nonograms] [info] 0.0874717 seconds, file 10590
    [2018-08-04 22:58:03.820] [Nonograms] [info] 0.0831248 seconds, file 4649
    [2018-08-04 22:58:03.820] [Nonograms] [info] 0.0782811 seconds, file 11922
    [2018-08-04 22:58:03.820] [Nonograms] [info] 0.0734325 seconds, file 4655
    [2018-08-04 22:58:03.820] [Nonograms] [info] 0.071352 seconds, file 10828
    [2018-08-04 22:58:03.820] [Nonograms] [info] 0.0708263 seconds, file 11824

    Черно-белые нонограммы
    izaron@izaron:~/nonograms/build$ ./nonograms_solver -x ../../nonogram/source/ -e -b
    [2018-08-04 22:59:56.308] [Nonograms] [info] Starting a benchmark...
    [2018-08-04 23:00:09.781] [Nonograms] [info] Average time: 0.0095679, Median time: 0.00239274, Max time: 0.607341
    [2018-08-04 23:00:09.781] [Nonograms] [info] Top 10 heaviest nonograms:
    [2018-08-04 23:00:09.781] [Nonograms] [info] 0.607341 seconds, file 5126
    [2018-08-04 23:00:09.781] [Nonograms] [info] 0.458932 seconds, file 19510
    [2018-08-04 23:00:09.781] [Nonograms] [info] 0.452245 seconds, file 5114
    [2018-08-04 23:00:09.781] [Nonograms] [info] 0.19975 seconds, file 18627
    [2018-08-04 23:00:09.781] [Nonograms] [info] 0.163028 seconds, file 5876
    [2018-08-04 23:00:09.781] [Nonograms] [info] 0.161675 seconds, file 17403
    [2018-08-04 23:00:09.781] [Nonograms] [info] 0.153693 seconds, file 12771
    [2018-08-04 23:00:09.781] [Nonograms] [info] 0.146731 seconds, file 5178
    [2018-08-04 23:00:09.781] [Nonograms] [info] 0.142732 seconds, file 18244
    [2018-08-04 23:00:09.781] [Nonograms] [info] 0.136131 seconds, file 19385

    Можно заметить, что в среднем черно-белые кроссворды "сложнее" цветных. Это подтверждает наблюдения любителей игры, которые также считают, что "цветные" решаются в среднем легче.


    Таким образом, без радикальных правок (таких как переписывание всего кода на на С или ассемблерных вставок с fastcall и опусканием указателя фрейма) можно достичь высокой производительностью, заметим, на весьма скромном компьютере. К оптимизациям может быть применим принцип Парето — окажется, что мелкая оптимизация влияет сильно, потому что этот кусок кода критичен и вызывается очень часто.


    Дальнейшая оптимизация


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


    • Переписывание кода на C-style и прочий 1972 год. Заменяем ifstream на сишный аналог, векторы на массивы, учим все опции компилятора и боремся за каждый такт процессора.
    • Распараллеливание. Например, в коде есть кусок, где последовательно обновляются строки, потом столбцы:

    bool Puzzle::UpdateState(...)
        ...
        if (!UpdateGroupsState(solver, dead_rows, row_groups, row_masks)) {                                                                               
            return false;                                                              
        }                                                                              
        if (!UpdateGroupsState(solver, dead_cols, col_groups, col_masks)) {           
            return false;                                                              
        }                                                                             
        return true;

    Эти функции независимы друг от друга и у них нет общей памяти, кроме переменной solver (тип OneLineSolver), так что можно создать два объекта "решателя" (здесь очевидно только один — solver) и запустить два потока в этой же функции. Эффект такой — в цветных кроссвордах "самый тяжелый" решается в два раза быстрее, в черно-белых такой же на треть быстрее, а среднее время увеличилось, за счет сравнительно больших затрат на создание потоков.


    Но вообще я бы не советовал делать прямо так в текущей программе и спокойно использовать — во-первых, создание потоков затратная операция, не стоит постоянно создавать потоки для микросекундных задач, во-вторых, при некоторой комбинации аргументов программы, потоки могут обращаться одновременно к какой-то внешней памяти, например при создании картинок решения — это надо бы учесть и обезопасить.


    Если бы задача была серьезной и у меня было бы много данных и многоядерные машины, я бы пошел еще дальше — можно завести несколько постоянных потоков, у каждого будет свой объект OneLineSolver, и еще одна thread-safe структура, которая рулит распределением работы и по запросу к ней выдает референс на очередную строку/столбец для решения. Потоки после решения очередной задачи обращаются к структуре заново, чтобы решать что-то еще. Какую-то задачу-нонограмму в принципе можно начать решать, не закончив предыдущей, например когда эта структура занимается взаимным AND всех клеток, и тогда какие-то потоки свободны и ничего не делают. Еще распараллеливание можно провести на графическом процессоре через CUDA — вариантов много.


    • Использование классических структур данных. Обратите внимание — когда я показывал псевдокод решения для цветных нонограмм, функции CanInsertColor и PlaceColor работают вовсе не за $O(1)$, в отличие от черно-белых нонограмм. Выглядит в программе это так:

    CanPlaceColor и SetPlaceColor
    bool OneLineSolver::CanPlaceColor(const vector<int>& cells, int color,          
            int lbound, int rbound) {                                               
        // Went out of the border                                                   
        if (rbound >= cells.size()) {                                               
            return false;                                                           
        }                                                                           
    
        // We can paint a block of cells with a certain color if and only if it is  
        // possible for all cells to have this color (that means, if every cell     
        // from the block has color-th bit set to 1)                                                                                                       
        int mask = 1 << color;                                                                                                                             
        for (int i = lbound; i <= rbound; ++i) {                                    
            if (!(cells[i] & mask)) {                                               
                return false;                                                       
            }                                                                       
        }                                                                           
        return true;
    }
    
    void OneLineSolver::SetPlaceColor(int color, int lbound, int rbound) {                                                                                 
        // Every cell from the block now can have this color                        
        for (int i = lbound; i <= rbound; ++i) {                                    
            result_cells_[i] |= (1 << color);                                       
        }                                                                           
    } 

    То есть работает за линию, за $O(N)$ (Позже будет объяснение смысла именно такого кода).


    Посмотрим, как можно получить лучшую сложность. Возьмем CanPlaceColor. Эта функция проверяет, что среди всех чисел вектора в отрезке $[lbound, rbound]$ бит номер $color$ установлен в 1. Эквивалентно этому можно взять $AND$ всех чисел этого отрезка и проверить бит номер $color$. Используя тот факт, что операция $AND$ коммутативная, также как сумма, минимум/максимум, произведение, или операция $XOR$, для быстрого подсчета $AND$ всего отрезка можно использовать почти любую структуру данных, работающую с коммутативными операциями на отрезке. Это:


    1. SQRT-декомпозиция. Предподсчет $O(\sqrt{N})$, запрос $O(\sqrt{N})$. Статья на Хабре.
    2. Дерево отрезков. Предподсчет $O(N\log{N})$, запрос $O(\log{N})$. Сотни статей в интернете.
    3. Разреженная таблица (Sparse Table). Предподсчет $O(N\log{N})$, запрос $O(1)$. Статья.

    К сожалению, особо сильные колдунства как алгоритм Фарака-Колтона и Бендера (предподсчет $O(N)$, запрос $O(1)$) использовать нельзя, так как вкурив статьи, можно понять, они созданы только для таких операций $\varphi$, что $\varphi(\alpha, \beta) \in \{\alpha, \beta\}$, то есть результат коммутативной операции — один из аргументов (жаль, а так хотелось...)


    Теперь возьмем SetPlaceColor. Тут надо произвести операцию на отрезке массива. Это тоже можно делать с SQRT-декомпозицией или деревом отрезков с ленивым обновлением (оно же "с проталкиванием"). А для обеих функций одновременно можно еще использовать убер-структуру декартово дерево по неявному ключу с обновлением и запросом за логарифм.


    Еще можно использовать расширение алгоритма для черного-белого кроссворда — частичные суммы для всех цветов.


    Итак, возникает вопрос — почему мы не используем все это богатство комплюктерн саенс, а делаем за линию? На это есть несколько ответов:


    1. Меньшая сложность вычислений не значит меньшее время работы. Алгоритм за $\log$ может потребовать различных преподсчетов, выделений памяти, другой тряски ресурсов — у этого алгоритма может быть довольно высокая константа (не в смысле magic number в нейроночках, а аффект по времени работы). Очевидно, что если $N=10^{5}$, то условный алгоритм за $O(N^2)$ будет работать условные 10 секунд, а условный алгоритм $O(N\log{N})$ за условные 0.150 секунд, но все может поменяться при достаточно маленьких $N$, особенно когда таких вычислений много. Еще непонятнее, когда сложности очень похожие и одной сложности сложно перебить другую сложность (сложные приколы): $O(N\sqrt{N})$ versus $O(N\log^2{N})$. В нашей задаче $N$ (длина ряда) очень маленькое — в среднем около 15-30.
    2. Запросов может оказаться достаточно мало, чтобы предпосчеты были бесполезными и жрали ресурсы просто так.

    То есть объяснение простое — оба этих пункта выполняются и вставка этих чудес программерской мысли вместо тупого алгоритма либо очень слабо оптимизирует программу, либо только увеличивает время работы из-за специфики вычислений — очень маленького $N$ и не такого большого количества запросов, чтобы они грузили процессор. Факт про запросы доказывает то, что по мнению профайлера те функции занимают ~25% и ~11% времени соответственно, то есть довольно мало для такого потенциально слабого места программы. Даже если у нас из-за этого возникает большая оценка сложности алгоритма, стоит понимать, что в таких типах задач это оценка сверху, а реальное время работы на рандомном кроссворде всегда премного ниже.


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


    • Правки алгоритма. Может оказаться так, что в среднем алгоритм на этих входных данных очень неплохо начинает работать, если поменять что-то неочевидное. В нашем случае этим может быть такое: логично же предположить, что если мы успешно обновили данные в строке, то обновленные в нем клетки скорее всего "триггерят" соответствующие столбцы? Значит, лучше эти столбцы обновить быстрее всех, сразу после этой строки! То есть образуется очередь из таких подзадач. Я не пробовал именно такой алгоритм, может это реально быстрее на нашем датасете.
    • Внезапные изменения техзадания (окажется, что поступают кроссворды по 1337 цветов или размером 1000x1000) тоже требуют оптимизации. Для большого количество цветов можно использовать быстрый std::bitset, для размера — те же структуры данных, и так далее.

    В общем, вот такие прикольные оптимизации. "Пихание" бедного алгоритма в зависимости от условий это весело и познавательно. Можно узнать про разные крутые вещи, как встроенное декартово дерево по неявному ключу в C++ (это rope, но велосипедная писанина практически всегда работает быстрее), особые встроенные сорта деревьев и спрятанный hashtable, работающий в 3-6 раза быстрее по вставке/удалению и в 4-10 раз по записи/чтению, чем unordered_map. Не говоря уже про различные нестандартные мазафаки со стороны, например из boost.


    ROFL Omitting Format Language



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


    ROFL — рекурсивный акроним, прямо как "GNU's Not Unix". Собственно, смысл формата в том, что картинка кодируется в виде японского кроссворда, а редактор, чтобы прочитать ее, должен решить этот кроссворд. Отсюда и слово Omitting в названии — формат как бы скрывает истинное положение дел в картинке (что, кстати, может быть полезным в криптографии: можно передавать японские кроссворды с зашифрованными в нем паролями — все хакеры повесятся).


    Лучше, если формат был бы похож на Matroska — в начале файла 4 байта [52][4F][46][4C], затем в трех байтах размер картинки и количество цветов, потом описание цветов, потом цвета, каждый по 3 байта, и потом описание всех групп — длина, количество клеток и цвет.


    Формат свободный, лицензия MIT, финансовые отчисления добровольные — мне на кофе, как автору.


    Исходники


    Исходники программы лежат на GitHub. У программы есть множество возможных аргументов, генерация кроссворда из картинки, генерация картинок из кроссворда (кстати, почти все картинки в статье сгенерированы через этот код). Дополнительными библиотеками были Magick++ и args.


    В репозитории картинки для кроссворда я добавил несколько примеров, взятых из интернета (они не являются частью проекта). Спарсенные тысячи с nonograms.ru никуда не выкладывал и не буду, чтобы их не потырили, потому что они могут быть защищены правами и тырятеля могут ожидать интересные последствия.


    Особую благодарность хочу выразить автору nonograms.ru Чугунному К.А. KyberPrizrak, за создание прекрасного, лучшего и одного из крупнейших в интернете сайта по нонограммам, и за разрешение использовать материалы для статьи! Все примеры в этой статье я взял с nonograms.ru, вот список источников.


    Исходники.


    Поделиться публикацией
    AdBlock похитил этот баннер, но баннеры не зубы — отрастут

    Подробнее
    Реклама

    Комментарии 20

      0
      Круто! Когда-то давно и я пробовал писать такую программу, но скорость вашей поражает!
      Я делал перебором по строкам и затем по столбцам с одной оптимизацией: в одной строке/столбце можно стекать все группы влево, затем вправо и быстро таким способом обнаружить закрашенные клетки, они могут позднее помочь ускорить решение по перпендикулярному направлению
        +1

        Круто. Кстати, самые тупые алгоритмы перебором очень хорошо распараллеливаются.

          0
          Оранжевым помечены «нерешаемые» клетки.
          Я бы поспорил с Википедией, потому что метод частичного перебора тоже можно считать одним из способов решения.
          Например, для правой части 1 легко находится, поскольку куда её не поставь — везде ломается конфигурация следующей единицы после тройки. И это делается просто «в уме», беглым просмотром. Как это алгоритмизировать — без понятия, но я раздумывал над алгоритмом решения чёрно-белого кроссворда.

          И вопрос к автору: есть быстрая реализация для чёрно-белого варианта? Цветной-то заведомо легче.
            0

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


            Можно рассмотреть такой вариант — каждой строке и столбцу привязать два дерева Фенвика с хорошей константой и памятью за линию (лучше чем у ДО), один для черных клеток, другой для белых, и тогда CanPlace работает за log(n), но и изменение массива в одном месте работает за log(n). По моим ощущениям это может побить сложные кроссворды но почти наверняка поднимет average time за счет побочных операций. А вообще можно в программу добавить счетчик, который будет считать реальное количество тактов в CanPlace — может оказаться, что реально на входных данных оно вызывается очень мало раз, и разные структуры только портят все.


            Дерево Фенвика не подходит к цветным кроссвордам, но если бы подходило, положение наверняка было бы плохо, так как состояний у клетки больше, т.е. изменение клетки массива не значит, что мы нашли ответ для этой клетки, а просто обновили список доступных вариантов.

              0

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

              0

              Вот кстати частичные суммы (сделал для этого отдельную ветку) — https://github.com/Izaron/Nonograms/tree/blacks

              0
              Мне одному в голову пришел «дурацкий» способ упаковки изображений.
                +3

                Видимо статью до конца не осилили? Там автор сообщает, что придумал принципиально новый формат картинок, основанный на японских кроссвордах и эффекте Даннинга-Крюгера.

                0
                На www.spoj.com/problems/JCROSS можно ещё попробовать.
                Интересно, сколько всего тестов получится пройти.
                  0

                  В этих тестах может быть не единственное решение, или недетерминированный результат. Алгоритм из статьи решает 173/300 паззлов за 0.11 секунды, т.е. именно столько тестов решаются логикой.
                  Если форсированно ставить в "непонятную" клетку белый цвет и идти дальше, полностью решаются 203/300, это 0.36 секунды. Исходник.


                  Такой тип задач — дают NP-полную и решай себе с разными фокусами и фортелями. Ограничение по времени 19 секунд. Можно пойти дальше так — собрать те 97 нерешаемых тестов, запустить на них бэктрекинг/полный перебор, начиная с самого маленького. Попутно надо следить за временем с high_resolution_clock, чтобы не пробить тайм-лимит, иначе 0 баллов. Еще может помочь отсечение по времени для каждого теста, чтобы не застрять на мерзком тесте.

                    0
                    Примерно, как я и думал. Только там, видимо, всего 252 теста, решённых у первых 12-ти с одним и тем же результатом в очках, а 252-173=79 как раз «Неправильные нонограммы». Я ставил форсированно «непонятную» клетку и в белый цвет и в чёрный. Получилось 217 тестов за 0.34 сек. Но как там рекорд 252 за 0.6 сек получился — это магия какая-то.
                      0
                      Переписал свой питоновский солвер на Rust, он решает все 252 задачи за 1.76 секунд (в таблице правда написано 2.10, но там первый успешный результат указан).
                        0
                        Да, видел. И с памятью 19М тоже лучше стало. По идее, могло бы и в таблице обновиться.
                        А надо было схитрить — один или несколько тестов придержать, отладить как надо, и окончательный вариант решить полностью.
                      +1
                      Попробовал на досуге ваш код
                      с такими данными
                      1
                      50 100

                      1 1 1 1 2 3 3 4 1 7 2 7 5 3 9 3 2 2 1 0
                      2 1 2 3 1 5 1 7 2 4 5 14 3 2 1 1 3 0
                      2 1 1 3 4 7 1 7 13 1 4 2 1 1 5 3 0
                      2 2 1 1 4 7 1 18 6 6 2 3 6 3 0
                      1 1 5 7 1 8 5 5 7 6 1 1 1 5 3 0
                      1 4 7 2 5 5 4 5 7 5 1 3 4 4 0
                      1 1 5 7 1 18 2 10 9 1 1 1 4 0
                      1 4 7 2 14 6 4 18 1 1 0
                      1 1 5 7 14 32 5 0
                      4 2 1 4 7 11 8 12 4 2 1 1 5 0
                      3 1 5 7 3 8 24 3 6 1 5 0
                      3 5 5 6 3 12 5 3 5 13 1 1 2 4 0
                      2 8 1 5 6 2 6 12 15 1 1 0
                      1 10 1 5 6 7 2 3 15 2 1 1 1 1 1 0
                      4 1 2 5 6 1 1 6 7 6 2 4 6 5 2 1 2 0
                      4 3 1 6 5 1 1 5 3 1 2 3 2 1 9 3 5 2 2 1 0
                      4 2 1 6 5 1 4 2 1 2 2 1 2 12 3 6 1 2 1 0
                      4 1 6 6 1 1 1 2 1 1 1 1 2 2 1 7 10 3 2 1 0
                      5 6 6 2 1 3 1 2 1 8 9 3 8 2 1 1 1 1 0
                      1 16 11 3 1 2 1 1 2 1 1 2 6 10 2 1 3 0
                      16 5 2 1 2 1 1 1 1 3 8 5 2 2 1 2 0
                      2 14 4 5 1 2 2 3 1 2 3 6 2 3 2 1 1 3 0
                      1 3 11 4 5 2 1 2 1 1 2 1 2 7 7 2 3 1 0
                      2 4 10 4 5 1 1 4 1 6 1 2 7 3 5 2 4 1 0
                      1 6 9 8 2 3 2 2 1 1 2 1 1 1 9 2 7 1 1 0
                      1 2 3 1 8 4 1 1 4 1 1 1 1 1 1 11 7 2 3 0
                      2 1 1 2 7 7 3 1 1 1 1 1 2 3 3 5 4 3 1 0
                      4 2 1 6 5 1 1 1 1 1 1 3 1 1 4 4 1 4 1 0
                      5 1 1 2 6 1 5 2 1 1 2 3 3 10 2 4 0
                      6 1 1 1 8 8 1 1 2 1 1 8 2 9 2 3 2 0
                      7 2 1 2 1 2 1 1 1 1 1 4 1 2 2 14 3 3 2 1 1 0
                      8 2 2 1 1 1 1 5 1 1 1 1 2 2 1 7 6 3 1 2 1 0
                      3 4 1 1 3 1 3 1 1 2 2 1 1 1 3 2 3 5 5 1 1 1 0
                      5 3 4 2 3 5 1 1 1 5 2 1 3 4 5 4 1 2 0
                      5 3 4 3 7 1 1 3 1 1 1 1 3 1 6 8 3 2 2 1 2 0
                      3 1 3 4 6 3 1 1 1 1 6 1 1 7 2 6 2 1 1 0
                      3 6 1 8 3 1 1 1 1 1 5 2 4 1 4 4 1 2 1 0
                      6 3 15 2 2 2 1 5 1 2 4 1 6 3 1 2 1 1 2 0
                      9 17 10 2 1 2 4 9 2 1 1 1 1 0
                      8 2 4 4 6 2 2 4 7 1 1 3 2 3 1 0
                      5 2 1 4 4 7 3 1 3 4 3 1 7 1 0
                      5 2 5 4 4 4 5 6 3 5 1 3 2 1 1 2 1 1 0
                      9 3 4 4 2 7 6 2 2 1 10 7 7 1 1 1 1 0
                      2 6 3 2 3 3 7 4 4 2 12 18 1 0
                      8 7 3 2 2 6 2 8 14 2 4 4 1 0
                      7 5 7 2 6 2 3 1 17 2 4 2 2 1 6 0
                      6 5 5 8 7 8 19 5 3 1 7 0
                      3 1 7 5 2 4 3 2 8 4 1 2 14 1 2 1 1 1 8 0
                      3 1 7 5 4 14 6 3 2 2 2 11 2 1 5 9 0
                      4 7 5 16 1 1 3 2 8 2 7 2 4 1 3 2 2 0

                      1 1 1 4 1 25 0
                      4 4 1 1 24 0
                      1 3 4 1 16 6 0
                      1 1 6 1 1 5 2 10 1 0
                      1 1 7 1 2 4 2 13 0
                      2 1 9 1 4 5 5 0
                      1 3 4 3 1 1 5 1 2 5 3 0
                      3 3 2 3 3 1 9 4 4 0
                      1 1 1 6 7 1 7 2 4 0
                      1 2 1 3 1 3 2 1 2 3 3 4 0
                      1 2 3 1 3 3 3 4 1 4 0
                      4 4 3 6 2 1 4 0
                      1 4 1 1 3 8 3 0
                      2 5 1 3 3 1 3 0
                      3 2 5 1 1 3 3 4 0
                      2 3 10 1 1 6 6 0
                      1 1 1 13 3 1 4 3 6 0
                      1 1 16 6 4 5 0
                      2 18 1 10 4 0
                      1 23 1 1 6 1 3 0
                      23 1 5 3 0
                      11 8 1 15 0
                      9 1 5 3 4 5 0
                      7 11 5 1 1 3 3 3 2 0
                      5 15 4 1 2 2 2 3 2 0
                      2 17 4 2 3 3 1 0
                      20 3 1 2 3 2 5 0
                      1 22 2 2 3 2 5 0
                      1 11 3 4 2 1 3 3 2 4 0
                      9 1 1 4 4 1 1 4 1 2 0
                      9 1 2 7 1 1 1 1 4 2 0
                      8 2 8 2 1 1 5 2 0
                      7 2 3 6 2 1 3 1 1 4 3 0
                      5 1 2 2 4 2 1 1 5 3 0
                      4 2 4 3 1 3 2 4 4 0
                      3 10 2 1 1 1 2 6 2 0
                      2 1 4 4 3 1 2 1 1 1 1 4 2 0
                      1 14 1 2 1 1 1 4 0
                      1 12 8 1 2 1 6 0
                      11 4 1 2 1 3 1 3 1 0
                      1 2 8 4 1 1 2 3 1 2 3 0
                      1 1 9 1 1 1 1 2 1 1 2 3 2 0
                      12 1 1 1 1 2 1 1 1 8 0
                      10 1 2 1 1 2 4 1 2 8 0
                      1 8 1 2 1 1 1 2 2 1 2 1 1 4 0
                      12 1 2 1 1 2 1 1 1 1 7 0
                      5 6 2 1 1 4 1 1 2 1 7 0
                      4 4 1 1 1 5 1 1 3 1 5 0
                      2 6 1 2 1 2 1 1 4 4 1 0
                      1 7 1 1 1 1 1 2 2 1 1 1 0
                      7 1 2 1 1 2 1 5 1 1 2 0
                      7 3 3 1 4 1 1 1 2 3 0
                      5 5 1 3 1 4 3 1 2 4 0
                      4 7 1 2 1 1 3 1 1 1 3 1 0
                      12 7 1 3 1 3 4 1 0
                      1 4 5 4 2 1 1 1 1 1 4 1 0
                      11 1 1 1 2 1 2 1 1 4 2 0
                      5 7 1 1 1 1 2 1 2 4 2 0
                      6 3 1 2 1 2 1 1 2 2 6 1 0
                      4 4 4 2 1 8 1 1 5 1 0
                      2 11 5 1 1 4 8 0
                      13 2 2 1 2 7 0
                      2 8 5 2 1 3 5 0
                      2 2 1 3 3 3 3 1 3 5 2 0
                      2 8 4 4 3 9 0
                      2 10 1 3 7 3 3 7 0
                      14 5 6 3 8 0
                      2 10 5 2 2 14 5 0
                      12 6 4 2 5 8 1 4 0
                      6 3 1 1 8 3 4 3 2 1 4 0
                      9 9 3 5 3 1 2 4 0
                      1 3 4 9 6 4 1 3 4 3 0
                      1 2 3 6 1 6 3 1 1 4 3 0
                      2 6 6 2 10 5 1 2 4 0
                      1 4 3 5 5 4 4 2 8 0
                      2 4 9 7 1 1 3 3 1 8 0
                      1 2 4 8 8 3 1 3 7 0
                      1 1 10 6 4 6 4 1 5 0
                      1 1 4 11 1 4 5 3 1 0
                      2 1 14 3 10 2 2 1 1 1 0
                      14 5 3 5 2 2 1 2 0
                      1 1 2 3 3 2 1 6 4 5 2 0
                      1 2 13 2 6 11 2 1 0
                      1 1 2 18 5 4 8 2 0
                      1 1 1 3 3 2 9 2 3 2 6 1 0
                      1 1 6 1 12 1 4 2 0
                      1 1 5 3 4 2 3 3 0
                      1 1 6 1 5 2 3 1 0
                      1 8 5 1 2 0
                      2 14 1 2 0
                      3 1 5 3 2 0
                      5 2 1 1 1 5 1 1 1 2 0
                      5 1 1 2 2 12 1 1 3 0
                      5 2 3 1 1 2 1 2 1 1 4 0
                      4 1 1 1 2 1 1 1 4 0
                      1 3 2 1 1 1 1 3 5 0
                      2 4 1 1 1 2 5 0
                      5 4 1 1 1 1 2 4 0
                      11 1 1 1 2 1 1 1 5 5 0
                      1 5 4 1 2 1 1 6 5 0


                      Картинка нарисовалась, но если я правильно понял, неопределённые клетки, как и пустые, помечаются точкой. Тогда ошибки там получаются в крайних рядах.
                    +1

                    Странно, что ещё никто не написал приложение, которое разгадывает картинку, наведением камеры телефона на изображение в газете.

                    0
                    Флаг -O2 при компиляции

                    Флаг -O3 еще лучше.


                    Но вообще-то решение таких кроссвордов — это задача для SAT солвера.

                      0
                      Флаг -O3 еще лучше.

                      Не всегда. Почти все оптимизации компилятора эвристические, то есть некий дополнительный флаг (развернуть рекурсии, скрыть фреймы) в -O(x+1) по сравнению с -O(x) будет просто попыткой сделать программу более "оптимальной", но может иметь и обратный эффект. -O3 предполагается лучше -O2, как -O2 лучше -O1, это не значит что это будет правда работать быстрее на рандомной программе.


                      Довольно много вещей с -O2 работают быстрее чем с -O3, в т.ч. этот код. По-хорошему надо бы еще отдельные флаги оптимизации тестировать, чтобы затюнить быстродействие.


                      Но вообще-то решение таких кроссвордов — это задача для SAT солвера.

                      SAT солвер скушает быстро большой черно-белый кроссворд? Если не рассматривать цветные, для которых страшно генерировать конъюнктивную форму.
                      Пусть у нас есть K групп суммарно и кроссворд размером NxN, мы вычислили порядка KxN булевых переменных (для каждой группы все возможные позиции в его ряду) и состряпали формулу, она будет быстро считаться?
                      Может есть менее затратный способ составить формулу? Я к тому что если взять K=4 в среднем группы на ряд и N=50, выйдет 500'000 булевых переменных, которые будут считаться до Второго Пришествия.

                        0
                        Откуда появилась цифра 500000?
                        Кстати, это не такая большая цифра для современных SAT солверов.
                      +1
                      Делал недавно похожее на Python, в основном использовал паззлы с webpbn.com. Интересно было бы сравнить скорость с вашей версией на одинаковых задачах.

                      Только полноправные пользователи могут оставлять комментарии. Войдите, пожалуйста.

                      Самое читаемое