Статистики на отрезках. Корневая эвристика — различия между версиями

Материал из Викиконспекты
Перейти к: навигация, поиск
м (Обработка запроса)
м (rollbackEdits.php mass rollback)
 
(не показано 69 промежуточных версий 9 участников)
Строка 1: Строка 1:
{{Определение
 
|definition=
 
'''Корневая эвристика (Sqrt-декомпозиция)''' — это метод, или структура данных, которая позволяет выполнять некоторые ассоциативные операции над отрезками (суммирование элементов подмассива, нахождение минимума/максимума и т.д.) за <tex> O(\sqrt n)</tex>. }}
 
  
== Предпосчет ==
+
'''Корневая эвристика (Sqrt-декомпозиция)''' {{---}} это подход к реализации ассоциативных операций (например, суммирование элементов, нахождение минимума/максимума и т.д.) над идущими подряд элементами некоторого множества размера <tex>n</tex> за <tex> O(\sqrt n)</tex>.
[[Файл:sqrt.png|right|540px]]
 
Пусть нам дан массив <tex>A</tex> размерности <tex>n</tex>. Cделаем следующий предпосчет:
 
* разделим массив <tex>A</tex> на блоки длины <tex>len = \lfloor \sqrt{n} \rfloor</tex> ;
 
* в каждом блоке заранее предпосчитаем необходимую нам операцию;
 
* результаты предпосчёта запишем в массив <tex>B</tex> размерности <tex>cnt</tex>, где <tex>cnt = \left\lceil \frac{n}{len} \right\rceil</tex> — количество блоков.
 
  
 +
== Построение ==
 +
Пусть дан массив <tex>A</tex> размерности <tex>n</tex>. Cделаем следующие действия:
 +
* разделим массив <tex>A</tex> на блоки длины <tex>len = \lfloor \sqrt{n} \rfloor</tex> ,
 +
* в каждом блоке заранее посчитаем необходимую операцию,
 +
* результаты подсчета запишем в массив <tex>B</tex> размерности <tex>cnt</tex>, где <tex>cnt = \left\lceil \dfrac{n}{len} \right\rceil</tex> {{---}} количество блоков.
  
Пример реализации предпосчета для запроса "подсчет суммы":
+
[[Файл:sqrt.png|358px]]
<pre>
+
 
precalc()
+
Пример реализации построения массива <tex>B</tex> для операции <tex> \circ </tex>:
    for i = 0 to n - 1
+
<code>
        B[i / len] += A[i]
+
'''void''' build():
</pre>
+
    '''for''' i = 0 ... cnt
 +
        B[i] = neutral                <font color=green>// neutral {{---}} нейтральный элемент для операции <tex> \circ </tex> </font>
 +
    '''for''' i = 0 ... n - 1
 +
        B[i / len] = B[i / len] <tex> \circ </tex> A[i]
 +
</code>
  
  
Пердпосчет, очевидно, происходит за <tex>O(n)</tex> времени.
+
Построение, очевидно, происходит за <tex>O(n)</tex> времени.
  
 
== Обработка запроса ==
 
== Обработка запроса ==
[[Файл:sqrt(sum).png|right|520px]]
+
Пусть получен запрос на выполнение операции на отрезке <tex>[l, r]</tex>. Отрезок может охватить некоторые блоки массива <tex>B</tex> полностью, а так же не более двух блоков (начальный и конечный) {{---}} не полностью.
Пусть мы получили запрос на выполнение операции на отрезке <tex>[l, r]</tex>. Отрезок может охватить некоторые блоки массива <tex>B</tex> полностью, а так же не более двух блоков (начальный и конечный) - не полностью.
 
  
Таким образом, для того чтобы найти результат операции на отрезке <tex>[l, r]</tex> нам необходимо вручную выполнить ее на "хвостах", а потом выполнить ее для полученного результата и полных блоков, предпосчет которых мы сделали заранее.
+
Таким образом, для того чтобы найти результат операции на отрезке <tex>[l, r]</tex> необходимо вручную выполнить ее на "хвостах", а потом выполнить ее для полученного результата и полных блоков, значения которых мы посчитали заранее.
  
 +
[[Файл:sqrt(sum).png|358px]]
  
Пример реализации обработки запроса "подсчет суммы на отрезке <tex>[l, r]</tex> " :
+
Пример реализации обработки запроса:
  
<pre>
+
<tex> \circ </tex> {{---}} операция, для которой было сделано построение.
request(l, r)
 
    left = l / len
 
    right = r / len
 
    end = (left + 1) * len - 1
 
    sum = 0
 
  
    if left == right
+
<code>
        for i = l to r
+
'''T''' query('''int''' l, '''int''' r):
    sum += A[i]
+
    left = l / len
    else
+
    right = r / len
        for i = l to end
+
    end = (left + 1) * len - 1
            sum += A[i]
+
    res = neutral                      <font color=green> // neutral {{---}} нейтральный элемент для операции <tex> \circ </tex> </font>
        for i = left + 1 to right - 1
+
    '''if''' left == right
            sum += B[i]
+
        '''for''' i = l ... r
        for i = right * len to r
+
            res = res <tex> \circ </tex> A[i]
            sum += A[i]
+
    '''else'''
</pre>
+
        '''for''' i = l ... end
 +
            res = res <tex> \circ </tex> A[i]
 +
        '''for''' i = left + 1 ... right - 1
 +
            res = res <tex> \circ </tex> B[i]
 +
        '''for''' i = right * len ... r
 +
            res = res <tex> \circ </tex> A[i]
 +
</code>
  
  
Размер каждого из "хвостов", очевидно, не превосходит длины блока <tex>len</tex>, а количество блоков не превосходит <tex>cnt</tex>. Поскольку и <tex>len</tex>, и <tex>cnt</tex> мы выбирали <tex>~ ~ \approx \sqrt{n}</tex>, то для выполнения операции на отрезке <tex>[l, r]</tex> нам понадобится <tex>O(\sqrt{n})</tex> времени.
+
Размер каждого из "хвостов", очевидно, не превосходит длины блока <tex>len</tex>, а количество блоков не превосходит <tex>cnt</tex>. Поскольку <tex>len</tex> было выбрано равным <tex>\lfloor \sqrt{n} \rfloor</tex> , а <tex>cnt</tex> было выбрано равным  <tex>\left\lceil \dfrac{n}{len} \right\rceil</tex> , то для выполнения операции на отрезке <tex>[l, r]</tex> понадобится <tex>O(\sqrt{n})</tex> времени.
  
 
== Запрос на изменение элемента ==
 
== Запрос на изменение элемента ==
[[Файл:sqrt(+delta).png|right|264px]]
+
Реализация данного запроса будет зависеть от того, имеет ли операция, для которой сделано построение, обратную операцию и обладает ли она свойством коммутативности.
Реализации данного запроса будет зависеть от того, имеет ли операция, для которой мы сделали предпосчет, обратную операцию и обладает ли она свойством коммутативности.
+
* если оба условия выполняются, то запрос на изменение элемента можно сделать за <tex>O(1)</tex> времени,
* если обратная операция существует, и выполняется свойство коммутативности, то нам не придется заново пересчитывать значение для блока;
+
* если хотя бы одно из условий не выполняется, то запрос на изменение элемента можно сделать за <tex>O(\sqrt{n})</tex> времени.
* если хотя бы одно из условий не выполняется, то нам придется заново пересчитать значение для блока, к которому принадлежит элемент указанный в запросе, и записать полученный результат в соответствующий элемент массива <tex>B</tex>.
 
  
 +
[[Файл:sqrt(+delta).png|358px]]
  
 
Примеры реализации:
 
Примеры реализации:
  
<tex>p</tex> - номер элемента из массива <tex>A</tex>, который необходимо заменить; <tex>newValue</tex> - новое значение для данного элемента.
+
* <tex>p</tex> {{---}} номер элемента из массива <tex>A</tex>, который необходимо заменить,
 +
* <tex>\mathtt{newValue}</tex> {{---}} новое значение для данного элемента.
  
 
Запрос на изменение элемента для операции, у которой есть обратная операция, и выполняется свойство коммутативности:
 
Запрос на изменение элемента для операции, у которой есть обратная операция, и выполняется свойство коммутативности:
  
 
<code>
 
<code>
  tmp = B[p / len] <tex> \circ </tex> inverse(A[p]) // <tex> \circ </tex> - операция, для которой был сделан предпосчет; inverse(A[p]) - обратный элемент
+
  '''function''' set('''int''' p, '''T''' newValue):
A[p] = newValue
+
    tmp = B[p / len] <tex> \circ </tex> inverse(A[p])   <font color=green>// inverse(A[p]) {{---}} обратный элемент</font>
B[p / len] = tmp <tex> \circ </tex> newValue
+
    A[p] = newValue
 +
    B[p / len] = tmp <tex> \circ </tex> newValue
 
</code>
 
</code>
  
Запрос на изменение элемента для поиска минимума (выполняется свойство коммутативности, но нет обратной операции):
+
'''Замечание:''' важность наличия свойства коммутативности подчеркивает следующий контрпример. Известно, что умножение матриц не коммутативно. Возьмем блок <tex> b_0 </tex>, как показано на иллюстрации выше, со следующими значениями:  
  
<pre>
+
<tex> b_0 = \begin{pmatrix} 27 & 32 \\ 42 & 50 \end{pmatrix} </tex> ,
index = len * p / cnt
 
A[p] = newValue
 
B[p / len] = newValue
 
for i = index to index + len
 
    B[p / len] = min(A[i], B[p / len])
 
</pre>
 
  
 +
<tex> a_0 = \begin{pmatrix} 1 & 1 \\ 1 & 2 \end{pmatrix} </tex> ,
  
Таким образом, запрос на изменение элемента происходит не более чем за длину блока <tex>len</tex>, т.е. не более чем за <tex>O(\sqrt{n})</tex> времени.
+
<tex> a_1 = \begin{pmatrix} 2 & 2 \\ 2 & 3 \end{pmatrix} </tex> ,
  
==Источники==
+
<tex> a_2 = \begin{pmatrix} 3 & 3 \\ 3 & 4 \end{pmatrix} </tex>.
* [http://www.e-maxx.ru/algo/sqrt_decomposition Maximal:: algo:: Sqrt - декомпозиция]
+
 
 +
Пусть необходимо изменить значение матрицы <tex> a_1 </tex> на следующее:
 +
 
 +
<tex> \mathtt{newValue} = a_1 = \begin{pmatrix} 4 & 4 \\ 4 & 5 \end{pmatrix} </tex>.
 +
 
 +
Тогда значения <tex> a_1^{-1} </tex>, <tex> tmp  </tex> и новое значение <tex> a_1 </tex> таковы :
 +
 +
 
 +
<tex> a_1^{-1} = \begin{pmatrix} 1,5 & -1 \\ -1 & 1 \end{pmatrix} </tex>,
 +
 
 +
<tex> tmp = b \cdot a_1^{-1} = \begin{pmatrix} 8,5 & 5 \\ 13 & 8 \end{pmatrix} </tex> ,
 +
 
 +
<tex> a_1 = \begin{pmatrix} 4 & 4 \\ 4 & 5 \end{pmatrix} </tex>.
 +
 
 +
Тогда новое значение <tex> b_0  </tex> следующее:
 +
 +
 
 +
<tex> b_0 = \begin{pmatrix} 54 & 59 \\ 84 & 92 \end{pmatrix} </tex>.
 +
 
 +
Хотя правильный результат: <tex> b_0 = \begin{pmatrix} 51 & 60 \\ 78 & 92 \end{pmatrix} </tex>.
 +
 
 +
Запрос на изменение элемента для операции, у которой хотя бы одно из условий не выполняется:
 +
 
 +
<code>
 +
'''function''' set('''int''' p, '''T''' newValue):
 +
    index = len * (p / len)
 +
    A[p] = newValue
 +
    B[p / len] = neutral            <font color = green> // neutral {{---}} нейтральный элемент для операции <tex> \circ </tex> </font>
 +
    '''for''' i = index ... index + len - 1
 +
        B[p / len] = B[p / len] <tex> \circ </tex> A[i]
 +
</code>
 +
 
 +
==См. также==
 +
* [[Дерево отрезков. Построение]]
 +
* [[Многомерное дерево отрезков]]
 +
 
 +
==Источники информации==
 +
* [http://www.e-maxx.ru/algo/sqrt_decomposition Maximal:: algo:: Sqrt-декомпозиция]
 
* [http://habrahabr.ru/post/138946/#habracut Sqrt-декомпозиция (корневая оптимизация)]
 
* [http://habrahabr.ru/post/138946/#habracut Sqrt-декомпозиция (корневая оптимизация)]
 +
  
 
[[Категория: Дискретная математика и алгоритмы]]
 
[[Категория: Дискретная математика и алгоритмы]]
 
[[Категория: Дерево отрезков]]
 
[[Категория: Дерево отрезков]]

Текущая версия на 19:30, 4 сентября 2022

Корневая эвристика (Sqrt-декомпозиция) — это подход к реализации ассоциативных операций (например, суммирование элементов, нахождение минимума/максимума и т.д.) над идущими подряд элементами некоторого множества размера [math]n[/math] за [math] O(\sqrt n)[/math].

Построение

Пусть дан массив [math]A[/math] размерности [math]n[/math]. Cделаем следующие действия:

  • разделим массив [math]A[/math] на блоки длины [math]len = \lfloor \sqrt{n} \rfloor[/math] ,
  • в каждом блоке заранее посчитаем необходимую операцию,
  • результаты подсчета запишем в массив [math]B[/math] размерности [math]cnt[/math], где [math]cnt = \left\lceil \dfrac{n}{len} \right\rceil[/math] — количество блоков.

Sqrt.png

Пример реализации построения массива [math]B[/math] для операции [math] \circ [/math]:

void build():
    for i = 0 ... cnt
        B[i] = neutral                // neutral — нейтральный элемент для операции [math] \circ [/math] 
    for i = 0 ... n - 1
        B[i / len] = B[i / len] [math] \circ [/math] A[i]


Построение, очевидно, происходит за [math]O(n)[/math] времени.

Обработка запроса

Пусть получен запрос на выполнение операции на отрезке [math][l, r][/math]. Отрезок может охватить некоторые блоки массива [math]B[/math] полностью, а так же не более двух блоков (начальный и конечный) — не полностью.

Таким образом, для того чтобы найти результат операции на отрезке [math][l, r][/math] необходимо вручную выполнить ее на "хвостах", а потом выполнить ее для полученного результата и полных блоков, значения которых мы посчитали заранее.

Sqrt(sum).png

Пример реализации обработки запроса:

[math] \circ [/math] — операция, для которой было сделано построение.

T query(int l, int r):
    left = l / len
    right = r / len
    end = (left + 1) * len - 1
    res = neutral                       // neutral — нейтральный элемент для операции [math] \circ [/math] 
    if left == right
        for i = l ... r
            res = res [math] \circ [/math] A[i]
    else
        for i = l ... end
            res = res [math] \circ [/math] A[i]
        for i = left + 1 ... right - 1
            res = res [math] \circ [/math] B[i]
        for i = right * len ... r
            res = res [math] \circ [/math] A[i]


Размер каждого из "хвостов", очевидно, не превосходит длины блока [math]len[/math], а количество блоков не превосходит [math]cnt[/math]. Поскольку [math]len[/math] было выбрано равным [math]\lfloor \sqrt{n} \rfloor[/math] , а [math]cnt[/math] было выбрано равным [math]\left\lceil \dfrac{n}{len} \right\rceil[/math] , то для выполнения операции на отрезке [math][l, r][/math] понадобится [math]O(\sqrt{n})[/math] времени.

Запрос на изменение элемента

Реализация данного запроса будет зависеть от того, имеет ли операция, для которой сделано построение, обратную операцию и обладает ли она свойством коммутативности.

  • если оба условия выполняются, то запрос на изменение элемента можно сделать за [math]O(1)[/math] времени,
  • если хотя бы одно из условий не выполняется, то запрос на изменение элемента можно сделать за [math]O(\sqrt{n})[/math] времени.

Sqrt(+delta).png

Примеры реализации:

  • [math]p[/math] — номер элемента из массива [math]A[/math], который необходимо заменить,
  • [math]\mathtt{newValue}[/math] — новое значение для данного элемента.

Запрос на изменение элемента для операции, у которой есть обратная операция, и выполняется свойство коммутативности:

function set(int p, T newValue):
    tmp = B[p / len] [math] \circ [/math] inverse(A[p])   // inverse(A[p]) — обратный элемент
    A[p] = newValue
    B[p / len] = tmp [math] \circ [/math] newValue

Замечание: важность наличия свойства коммутативности подчеркивает следующий контрпример. Известно, что умножение матриц не коммутативно. Возьмем блок [math] b_0 [/math], как показано на иллюстрации выше, со следующими значениями:

[math] b_0 = \begin{pmatrix} 27 & 32 \\ 42 & 50 \end{pmatrix} [/math] ,

[math] a_0 = \begin{pmatrix} 1 & 1 \\ 1 & 2 \end{pmatrix} [/math] ,

[math] a_1 = \begin{pmatrix} 2 & 2 \\ 2 & 3 \end{pmatrix} [/math] ,

[math] a_2 = \begin{pmatrix} 3 & 3 \\ 3 & 4 \end{pmatrix} [/math].

Пусть необходимо изменить значение матрицы [math] a_1 [/math] на следующее:

[math] \mathtt{newValue} = a_1 = \begin{pmatrix} 4 & 4 \\ 4 & 5 \end{pmatrix} [/math].

Тогда значения [math] a_1^{-1} [/math], [math] tmp [/math] и новое значение [math] a_1 [/math] таковы :


[math] a_1^{-1} = \begin{pmatrix} 1,5 & -1 \\ -1 & 1 \end{pmatrix} [/math],

[math] tmp = b \cdot a_1^{-1} = \begin{pmatrix} 8,5 & 5 \\ 13 & 8 \end{pmatrix} [/math] ,

[math] a_1 = \begin{pmatrix} 4 & 4 \\ 4 & 5 \end{pmatrix} [/math].

Тогда новое значение [math] b_0 [/math] следующее:


[math] b_0 = \begin{pmatrix} 54 & 59 \\ 84 & 92 \end{pmatrix} [/math].

Хотя правильный результат: [math] b_0 = \begin{pmatrix} 51 & 60 \\ 78 & 92 \end{pmatrix} [/math].

Запрос на изменение элемента для операции, у которой хотя бы одно из условий не выполняется:

function set(int p, T newValue):
    index = len * (p / len)
    A[p] = newValue
    B[p / len] = neutral              // neutral — нейтральный элемент для операции [math] \circ [/math] 
    for i = index ... index + len - 1
        B[p / len] = B[p / len] [math] \circ [/math] A[i]

См. также

Источники информации