Ketepatan numerik dalam metode jumlah kotak?


13

Saya telah membaca sedikit tentang metode sum-of-square (SOS) dari survei Barak & Steurer dan catatan kuliah Barak . Dalam kedua kasus mereka menyapu masalah akurasi numerik di bawah karpet.

Dari pemahaman saya (yang diakui terbatas) tentang metode ini, hal-hal berikut ini harus benar:

Dengan adanya sistem persamaan polinomial E atas variabel bernilai riil xRn , di mana semua parameter adalah O(1) ( n , |E| , dan derajat dari setiap kendala), derajat- " 2n " ( =O(1) ) Metode SOS menemukan penugasan yang memuaskan dari variabel atau membuktikan tidak ada dalam waktu O(1) .

Pertanyaan pertama saya adalah apakah klaim di atas benar (adakah argumen naif yang tidak menggunakan SOS untuk menyelesaikan ini?). Pertanyaan kedua adalah di mana akurasi numerik cocok. Jika saya ingin mendapatkan tugas yang memenuhi semua kendala dalam akurasi aditif ε , bagaimana runtime bergantung pada 1/ε ? Secara khusus, apakah jumlahnya banyak?

Motivasi untuk ini adalah, katakanlah, menerapkan pendekatan divide-and-menaklukkan pada sistem besar sampai kasus dasar adalah sistem ukuran O(1) .

EDIT: Dari Barak-Steurer, tampak bahwa "derajat jumlah algoritma kuadrat" pada hal.9 (dan paragraf yang mengarah ke sana) semua mendefinisikan masalah untuk solusi lebih dari R , dan pada kenyataannya definisi pseudo -Distribusi di bagian 2.2 adalah lebih R . Sekarang saya melihat dari Lemma 2.2, bagaimanapun, bahwa seseorang tidak dijamin solusi / sanggahan pada derajat 2 n tanpa variabel biner.lRR2n

Jadi saya bisa memperbaiki pertanyaan saya sedikit. Jika variabel Anda bukan biner, yang dikhawatirkan adalah urutan output tidak terbatas (bahkan mungkin tidak meningkat secara monoton?). Jadi pertanyaannya adalah: apakah φ ( l ) masih meningkat? Dan jika demikian, seberapa jauh Anda harus pergi untuk mendapatkan akurasi aditif ε ?φ(l)φ(l)ε

Meskipun kemungkinan ini tidak mengubah apa pun, kebetulan saya tahu sistem saya memuaskan (tidak ada penyangkalan tingkat apa pun), jadi saya benar-benar hanya peduli tentang seberapa besar kebutuhan . Akhirnya, saya tertarik pada solusi teoritis, bukan pemecah numerik.l


Jawaban:


1

Inilah komentar Boaz Barak tentang masalah ini:

Kami melakukan akurasi numerik di bawah permadani - literatur SOS yang lebih "tradisional" dari Parrilo, Lasserre dll. Membahas masalah-masalah ini (misalnya, lihat survei Monique Laurent dan referensi di dalamnya). Diketahui bahwa hirarki adalah monoton (tidak sulit untuk melihat bahwa suatu derajat distribusi psuedo khususnya derajat l - 1 satu), dan bahwa ia akan konvergen dalam derajat yang terbatas untuk setiap set persamaan persamaan (ini adalah Positivstellensatz). Tingkat pastinya bisa bervariasi. Secara umum, jika semua koefisien dari polinomial dibatasi dan Anda mencoba untuk membedakan antara case yang ada solusinya dan case yang dalam setiap penugasan salah satu persamaan dimatikan oleh ϵ , maka seseorang dapat melarutkannya ke sebuahll1ϵ -net untuk δ terkait dengan jumlah variabel, derajat persamaan, dan ϵ , dan kemudian (dengan asumsi net cukup "bagus" dan "mirip kubus"), derajat yang diperlukan harus secara kasar mencatat ukuran net.δδϵ


Diposting sebagai jawaban untuk menghindari bot komunitas menabrak pertanyaan lagi di masa depan.
Kaveh

1

Saya pikir jawaban saya mungkin tidak cukup, tetapi tetap untuk kelengkapan (meskipun lihat komentar Boaz di bawah untuk kemungkinan jawaban yang lebih baik)

Ketika kita membatasi diri pada variabel boolean, klaim dapat dilihat ketika untuk semua i [ n ] dengan pengamatan bahwa derajat 2(xi21)Ei[n] distribusi semu adalah distribusi aktual, yaitu, misalkan Anda memiliki distribusi semu μ ( x ) lebih dari solusi x persamaan polinomial Anda E memuaskan:2nμ(x)xE

dan  x { - 1 , 1 } n μ ( x )x{1,1}nμ(x) untuk semua polinomial p dengan derajat paling banyak nx{1,1}nμ(x)p2(x)0pn

Tetapi derajat polinomial mencakup polinomial indikator (misalnya, x 1 = 1 , x 2 = - 1 , x 3 = 1 memiliki 2 - 3 ( 1 + x 1 ) ( 1 - x 2 ) ( 1 - x 2 ) ( 1 + xnx1=1,x2=1,x3=1 yang semua nol di tempat lain dan 1 pada tugas itu). Jadi μ ( x ) 0 untuk semua x {23(1+x1)(1x2)(1+x3)μ(x)0 , sehingga kami menyimpulkan μ adalah distribusi yang sebenarnya atas solusi dari E . Derajat distribusi semu dapat ditemukan dengan menggunakan pemrograman semidefinit untuk menemukan derajat terkait operator ekspektasi semu dalam waktu n O ( ) , sehingga kita dapat menemukan distribusi aktual μ dalam waktu n O ( n ) dengan menggunakan pseudo- ekspektasi (sekarang ekspektasi aktual) untuk menemukan semua momen μ .x{1,1}nμEnO()μnO(n)μ

Jadi, jika , maka Anda dapat menemukan distribusi solusi ke E dalam O ( 1 ) waktu. Tentu saja, pencarian brute force menjamin hal yang sama.|E|=O(1)EO(1)

Namun, jika solusi belum tentu boolean, maka derajat- harapan semu tidak cukup untuk menemukan distribusi lebih dari solusi. Seperti dapat dilihat di atas, bukti bahwa derajat 2 n pseudo-distribusi adalah distribusi aktual tergantung pada kenyataan bahwa derajat dan polinomial cukup untuk 'memilih' tugas individu, yang tidak benar secara umum. Cara lain untuk melihatnya adalah bahwa polinomial variabel boolean dipertimbangkan2n2nn, sehingga derajat setiap monomial paling banyak adalahn.mod(xi2)n

Sebagai contoh, seseorang dapat mempertimbangkan mengganti setiap variabel biner dengan variabel 4-ary, katakan dengan memasukkan . Maka Anda harus memilikipseudo-harapangelar 4 n untuk menjamin pemulihan distribusi atas solusi.(xi21)(xi24)E4n

Sekarang, untuk jaminan teoretis, sepertinya mendekati akar sistem polinomal juga dikenal sebagai masalah ke-17 Smale, dan tampaknya ada algoritma waktu polinomial acak (Las Vegas) yang memecahkan masalah ini - lihat http://arxiv.org /pdf/1211.1528v1.pdf . Perhatikan bahwa ini tampaknya berada dalam model Blum-Shub-Smale, jadi operasi nyata adalah primitif. Saya tidak yakin apakah ini memberikan jaminan yang Anda butuhkan.


Saya pikir saya mungkin tidak membuat ini jelas: variabel saya berada di R , karena kalau tidak saya bisa melakukan pencarian O ( 2 n ) = O ( 1 ) yang sepele . Saya telah memperbarui pertanyaan untuk mencerminkan ini. SDP / SOS juga berlaku untuk masalah optimisasi input nyata, bukan? xiRO(2n)=O(1)
Jeremy Kun

Ups, kesalahan saya! Ya, itu memang berlaku untuk pengaturan yang lebih umum, walaupun seringkali kami menganggap kami berada di hypercube. Saya telah memperbarui jawaban saya, meskipun jawaban saya akan kurang jelas dari yang saya harapkan.
Joe Bebel

10
Kami melakukan akurasi numerik di bawah permadani - literatur SOS yang lebih "tradisional" dari Parrilo, Lasserre dll. Membahas masalah-masalah ini (misalnya, lihat survei Monique Laurent dan referensi di dalamnya). Diketahui bahwa hierarki adalah monoton (tidak sulit untuk melihat bahwa suatu derajat distribusi psuedo khususnya derajat - 1 satu), dan bahwa hierarki akan konvergen dalam derajat terbatas untuk setiap set persamaan persamaan (ini adalah Positivstellensatz). 1
Boaz Barak

9
..Tingkat tepatnya dapat bervariasi. Secara umum, jika semua koefisien dari polinomial dibatasi dan Anda mencoba untuk membedakan antara case yang ada solusinya dan case yang dalam setiap penugasan salah satu persamaan dimatikan oleh , maka seseorang dapat melarutkannya ke δ - net untuk δ terkait dengan jumlah variabel, derajat persamaan, dan ϵ , dan kemudian (dengan asumsi net cukup "bagus" dan "mirip kubus"), derajat yang diperlukan harus secara kasar mencatat ukuran net. ϵδδϵ
Boaz Barak

4
@BoazBarak mungkin ini bisa menjadi jawaban?
Suresh Venkat
Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.