Численная точность в методе суммы квадратов?

13

Я немного читал о методе суммы квадратов (SOS) из опроса Barak & Steurer и лекционных заметок Barak . В обоих случаях они охватывают вопросы числовой точности под ковриком.

Из моего (по общему признанию ограниченного) понимания метода следует следующее:

Для любой системы полиномиальных равенств над вещественными переменными x R n , где все параметры O ( 1 ) ( n , | E | и степень каждого ограничения), степень - « 2 n » ( = O ( 1 ) ) Метод SOS находит удовлетворительное назначение переменных или доказывает, что в O ( 1 ) времени не существует. ExRnO(1)n|E|2n=O(1)O(1)

Мой первый вопрос: верно ли приведенное выше утверждение (есть ли наивный аргумент, который не использует SOS для решения этой проблемы?). Второй вопрос, где числовая точность подходит. Если я хочу получить назначение, которое удовлетворяет всем ограничениям с точностью до аддитивной , как время выполнения зависит от 1 / ε ? В частности, это полином?ε1/ε

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

EDIT: От Barak-Steurer, оказывается , что «степень сумм квадратов алгоритм» на (и п.9 пунктов , ведущих к нему) все определения задач для решения более R , а на самом деле определение псевдо -распределение в разделе 2.2 находится над R . Однако теперь из леммы 2.2 я вижу, что не гарантируется решение / опровержение на степени 2 n без бинарных переменных.lRR2n

Поэтому я могу немного уточнить свой вопрос. Если ваши переменные не являются двоичными, беспокойство заключается в том, что последовательность выходов не является конечной (может быть, даже не монотонно возрастающей?). Таким образом, вопрос: φ ( l ) все еще увеличивается? И если да, то как далеко вы должны пройти, чтобы получить аддитивную точность ε ?φ(l)φ(l)ε

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

Джереми Кун
источник

Ответы:

1

Вот комментарий Боаза Барака по этому вопросу:

Мы подметаем числовую точность под ковром - более «традиционная» литература SOS о Паррило, Лассере и т. Д. Имеет дело с этими проблемами (например, см. Обзоры Моник Лоран и ссылки в них). Известно, что эта иерархия является монотонной (нетрудно видеть, что псевдораспределение степени является, в частности, степенью l - 1 ), и что она будет сходиться в конечной степени для любого фиксированного набора уравнений (это Positivstellensatz). Точная степень может варьироваться. Как правило, если все коэффициенты полиномов ограничены, и вы пытаетесь различить случай, когда существует решение, и случай, когда в любом назначении одно из уравнений выключено на ϵ , то можно дискретизировать это доll1ϵ для δ, связанной с числом переменных, степенью уравнений и ϵδδϵ , а затем (при условии, что сеть достаточно «хороша» и «подобна кубу») требуемая степень должна быть приблизительно равна log-размеру сети.

каве
источник
Размещено в качестве ответа, чтобы не допустить повторения вопроса ботом сообщества в будущем.
Каве
1

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

Когда мы ограничиваемся булевыми переменными, утверждение можно увидеть, когда для всех i [ n ] с наблюдением, что псевдораспределения степени 2 n являются действительными, то есть предположим, что вы имеете псевдораспределение µ ( x ) по решениям x ваших полиномиальных равенств E, удовлетворяющих:(xi21)Ei[n]2nμ(x)xE

и Σ х { - 1 , 1 } п μ ( х ) р 2 ( х ) 0 для всех полиномов р со степенью в большей пx{1,1}nμ(x)x{1,1}nμ(x)p2(x)0pn

Но полиномы степени включают индикаторный полином (например, x 1 = 1 , x 2 = - 1 , x 3 = 1 имеет 2 - 3 ( 1 + x 1 ) ( 1 - x 2 ) ( 1 + x 3 ), который все ноль в другом месте и 1 в этом назначении). Таким образом, µ ( x ) 0 для всех x {nx1=1,x2=1,x3=123(1+x1)(1x2)(1+x3)μ(x)0 , поэтому мы приходимвыводу μ является фактическим распределением по решениям Е . Степень л псевдо-распределения могут быть найдены с помощью полуопределенного программированиячтобы найти соответствующую степень л оператор псевдо-ожидание в п О ( ) время, поэтому мы можем найти фактическое распределение М во время п О ( п ) с помощью этого псевдо- ожидание (теперь фактическое ожидание), чтобы найти все моменты μ .x{1,1}nμEnO()μnO(n)μ

Так что, если , то можно найти распределение решений E за O ( 1 ) по времени. Конечно, поиск методом перебора гарантирует то же самое.|E|=O(1)EO(1)

Однако, если решения не обязательно являются булевыми, то псевдо-ожиданий степени недостаточно, чтобы найти распределение по решениям. Как видно выше, доказательство того, что псевдораспределения степени 2 n являются действительными, зависит от того факта, что полиномов степени n достаточно, чтобы «выбрать» отдельные присвоения, что в общем случае неверно. Еще один способ просмотра состоит в том, что булево-переменные полиномы считаются2n2nn , поэтому степень каждого монома не превосходит n .mod(xi2)n

Например, можно было бы рассмотреть возможность замены каждой бинарной переменной с 4-кратной переменной, скажем, в том числе . Тогда вам потребуется псевдо-ожидание степени 4 n , чтобы гарантировать восстановление распределения по решениям.(xi21)(xi24)E4n

Теперь, для теоретических гарантий, похоже, что аппроксимация корня системы полиномов также известна как 17-я проблема Смейла, и, по-видимому, существует рандомизированный (Лас-Вегас) алгоритм полиномиального времени, который решает эту проблему - см. Http://arxiv.org /pdf/1211.1528v1.pdf . Обратите внимание, что это, кажется, в модели Блюм-Шуб-Смейл, поэтому реальные операции являются примитивными. Я не уверен, дает ли это гарантию, что вам нужно.

Джо Бебель
источник
xiRO(2n)=O(1)
Ой, моя ошибка! Да, это применимо к более общим настройкам, хотя часто мы просто предполагаем, что находимся в гиперкубе. Я обновил свой ответ, хотя мой ответ будет менее ясным, чем я надеялся.
Джо Бебель
10
Мы подметаем числовую точность под ковром - более «традиционная» литература SOS о Паррило, Лассере и т. Д. Имеет дело с этими проблемами (например, см. Обзоры Моник Лоран и ссылки в них). Известно, что эта иерархия является монотонной (нетрудно увидеть, что степень псевдо-распределение - это, в частности, степень -11) и что он будет сходиться в конечной степени для любого фиксированного набора уравнений (это Positivstellensatz).
Вооз Барак
9
Точная степень может варьироваться. Как правило, если все коэффициенты полиномов ограничены, и вы пытаетесь различить случай, когда существует решение, и случай, когда в любом назначении одно из уравнений выключеноεтогда можно было бы дискретизировать это до δсеть для δ связанные с числом переменных, степенью уравнений и εи затем (при условии, что сеть достаточно «хороша» и «подобна кубу»), требуемая степень должна быть примерно равна логу размера сети.
Вооз Барак
4
@BoazBarak может быть, это может быть ответом?
Суреш Венкат