Znalezienie najrzadszego rozwiązania układu równań liniowych

13

Jak trudno jest znaleźć najrzadsze rozwiązanie układu równań liniowych?

Bardziej formalnie, rozważ następujący problem decyzyjny:

Instancja: układ równań liniowych o współczynnikach całkowitych i liczbie c .

Pytanie: Czy istnieje rozwiązanie dla systemu z co najmniej c zmiennymi przypisanymi do zera?

Próbuję również ustalić, na czym polega zależność od c . To znaczy, może problemem jest FPT z parametrem c .

Wszelkie pomysły lub referencje są bardzo mile widziane.

Michael Wehar
źródło

Odpowiedzi:

12

Rozważ problem MAX-LIN(R) polegający na maksymalizacji liczby spełnionych równań liniowych w pewnym pierścieniu R , który często jest NP-twardy, na przykład w przypadku R=Z

Weźmy przykład tego problemu, gdzie A jest macierzą n × m . Niech k = m + 1 . Skonstruuj nowy układ liniowy ˜ A ˜ x = ˜ b , gdzie ˜ A jest macierzą k n × ( k n + m ) , ˜ x jest teraz wektorem wymiarowym ( k n + m ) , a ˜ bAx=bAn×mk=m+1A~x~=b~A~kn×(kn+m)x~(kn+m)b~jest wektorem wymiarowym :kn

gdzieInjestmacierzą tożsamościn×n.

A~=[AInInInInInInIn],b~=[b00]
Inn×n

Należy pamiętać, że system ten jest zawsze spełniony przez wektor . W rzeczywistości pierwsze m wpisów ˜ x może być dowolne i istnieje jakiś wektor rozwiązania z tym prefiksem.x~=(0bbb)Tmx~

Teraz twierdzą, że część równania A x = b są spe IFF istnieje rzadki roztwór ~ A ~ x = ~ b , która ma co najmniej hemibursztynianu n k zera. Jest tak, ponieważ każdy spełniony wiersz A x = b daje k potencjalnych zer, gdy x jest rozszerzone do ˜ xδAx=bA~x~=b~δnkAx=bkxx~

Zatem jeśli znajdziemy rzadkość najrzadszego rozwiązania , zmaksymalizujemy również δ , dzieląc rzadkość przez k .A~x~=b~δk

Dlatego uważam, że twój problem jest trudny NP.

Joe Bebel
źródło
1
Chłodny! Dziękuję za podzielenie się. Jak myślisz, jaka jest zależność od c? Czy uważasz, że możemy rozwiązać to w mniej niż gdzie jest rozmiarem wejściowym? npoly(n)(nc)n
Michael Wehar,
1
Jasne: jeśli założymy, że podano ci, które elementów są równe zero, możesz po prostu usunąć te elementy z aby uzyskać niższy wymiar a także usunąć odpowiednie kolumny z aby uzyskać . Następnie użyj eliminacji gaussowskiej, aby zdecydować, czy możliwy jest zredukowany układ ; jeśli tak, to znalazłeś rzadkie rozwiązanie. Następnie wypróbuj wszystkie możliwe . x x x A A A x = b ( ncxxxAAAx=b Ax(nc)Ax
Joe Bebel,
1
@MichaelWehar Nie wiem, czy ten problem dotyczy FPT, czy nie
Joe Bebel,
6

Problemem jest NP-zupełny, przez redukcję z następującym problemem: podany w macierzy wpisami całkowite i całkowitej wektora z pozycji, nie istnieje do 0-1 wektor z ?A b n x A x = bm×nAbnxAx=b

Dla każdej współrzędnej wektora , xxix

  • wprowadzić nowych równań z , i 100(n+m)xi+yi,k=0k=1,,100(n+m)
  • wprowadzić nowych równań z . 100(n+m)xi+zi,k=1k=1,,100(n+m)

Ponadto użyj starego układu równań .Ax=b

Istnieje rozwiązanie 0-1 oryginalnego układu , tylko wtedy, gdy nowy system ma rozwiązanie, w którym co najmniej zmiennych wynosi zero.100 ( n + m ) nAx=b100(n+m)n

Gamow
źródło
4

Ten problem jest trudny w różnych ustawieniach. Jak stwierdzono w innych odpowiedziach na to pytanie, problem jest NP-zupełny w stosunku do liczb całkowitych.

W przetwarzaniu sygnału macierz i wektory mają racjonalne wpisy, a problem ten jest czasem nazywany problemem rzadkiej rekonstrukcji . W tym ustawieniu problem jest NP-zupełny (patrz Twierdzenie 1).

W teorii kodowania wpisy pochodzą z pola skończonego, a problem ten jest czasem nazywany problemem dekodowania o najwyższym prawdopodobieństwie . W tym ustawieniu problem jest NP-zupełny, a nie w czasie podwykonawczym , przy założeniu wykładniczej hipotezy o czasie. Ponadto, zgodnie z poprzednią wersją artykułu na arXiv (patrz Lemat C.2 w wersji 1 artykułu), problem polega na ukończeniu W [1].

argentpepper
źródło
Twoje odniesienie do kompletności W [1] nie wydaje się mieć „Lemmy C.2”.
@RickyDemer W wersji 1 artykułu, do którego prowadzi link, znajduje się Lemma C.2. Jednak wersja 2 wydaje się mieć inny tytuł i została bardzo niedawno zmieniona.
Michael Wehar,
Ten lemat wykorzystuje inną parametryzację niż PO.
Och, nie zdawałem sobie sprawy, że istnieje zaktualizowana wersja, przyjrzę się jej i odpowiednio zaktualizuję swoją odpowiedź.
argentpepper
Jak wspomniałem w poprzednim komentarzu, „lemat używa innej parametryzacji niż OP”, więc nawet jeśli założymy, że wynik jest prawdziwy (pomimo usunięcia z wersji 2), pytanie OP dotyczące sparametryzowanej złożoności nadal będzie otwarte.