Natknąłem się na następujący interesujący problem: niech będą wielomianami na polu liczb rzeczywistych i załóżmy, że wszystkie ich współczynniki są liczbami całkowitymi (tzn. Istnieje dokładna reprezentacja tych wielomianów). W razie potrzeby możemy założyć, że stopień obu wielomianów jest równy. Oznaczmy przez (resp. X_q ) największa wartość bezwzględna niektórych (rzeczywistej lub zespolonej) pierwiastka wielomianu p (resp. Q ). Czy można rozstrzygać o właściwości x_p = x_q ?x p x qq x p = x q
Jeśli nie, to czy ta właściwość obowiązuje dla niektórych ograniczonych rodzin wielomianów? W kontekście, z którego wynika ten problem, wielomiany są charakterystycznymi wielomianami macierzy, a ich pierwiastki są wartościami własnymi.
Zdaję sobie sprawę z niektórych algorytmów numerycznych do obliczania pierwiastków wielomianów / wartości własnych, jednak wydaje się, że tutaj one nie są przydatne, ponieważ dane wyjściowe tych algorytmów są jedynie przybliżone. Wydaje mi się, że algebra komputerowa może się tu przydać, niestety nie mam prawie żadnej wiedzy w tej dziedzinie.
Nie szukam szczegółowego rozwiązania tego problemu, jednak każda intuicja i pomysły na znalezienie rozwiązania byłyby pomocne.
Z góry dziękuję.
Odpowiedzi:
Nie mam również wiedzy w tej dziedzinie, ale myślę, że mogę udzielić niekonstruktywnej odpowiedzi.
Teoria pierwszego rzędu rzeczywistych pól zamkniętych jest rozstrzygalna. Twój problem można określić jako układ równań algebraicznych i nierówności względem rzeczywistych liczb algebraicznych. Rozważ zmienne . Chcesz wiedzieć, czy następujący system jest zadowalający: x 1 , … , x deg P , y 1 , … , y deg P , x ′ 1 , … , x ′ deg P , y ′ 1 , … , y ′ deg P \ begin { wyrównaj *} P (x_j + i \, y_j) & = 0 i \ text {for \ (1 \ le j \ le \ deg P \)} \\ Q (x'_k + i \, y'_k) & = 0 i \ text {dla \ (1 \ le k \ le \ deg Q \)} \\2 ( degP.+ degQ ) x1, … , XdegP., y1, … , YdegP., x′1, … , X′degP., y′1,…,y′degP.
P(x_j+i\,y_j) &= 0 & \text{for \(1 \le j \le \deg P\)} \\
Q(x'_k+i\,y'_k) &= 0 & \text{for \(1 \le k \le \deg Q\)} \\
x_j^2 + y_j^k &\le x_1^2 + x_2^2 & \text{for \(2 \le j \le \deg P\)} \\
x'_j^2 + y'_j^k &\le x'_1^2 + x'_2^2 & \text{for \(2 \le k \le \deg Q\)} \\
x_1^2 + y_1^2 = x'_1^2 + y'_1^2 \\
\end{align*}
Pierwsze dwie rodziny równań wyrażają, że i są pierwiastkami wielomianów, kolejne dwie rodziny nierówności wyrażają, że i mają największą wartość bezwzględną, a ostatnie równanie porównuje te największe wartości bezwzględne.x ′ k + ixjot+ iyjot x 1 + ix′k+ iy′k x ′ 1 + ix1+ iy1 x′1+ iy′1
Można zadecydować, czy ten system jest zadowalający: problem można rozstrzygnąć. Jednak to stwierdzenie nie jest prawdopodobnie najskuteczniejszym sposobem na obejście tego.
Bardziej użyteczna odpowiedź prawdopodobnie dotyczy teorii baz Gröbnera . Jeśli próbujesz rozwiązać ten problem sam, myślę, że przeczytanie kilku pierwszych rozdziałów jakiejkolwiek książki algebry obliczeniowej da ci niezbędne tło. Jeśli chcesz tylko rozwiązać swój podstawowy problem, prawdopodobnie istnieje gotowy algorytm do zaimplementowania.
źródło
Mogę się mylić w tej kwestii: nie mam zbyt dużej wiedzy w tej dziedzinie (gdzie są eksperci !?), ale uważam, że mam dość szybki algorytm do tego, o co pytasz.
Przypuszczam, dla uproszczenia, że wszystkie korzenie są prawdziwe. Znajdź przedział związany z pierwiastkiem o najwyższej wartości bezwzględnej (tj. Przedział taki, że i dla wszystkich innych pierwiastków ). Taki odstęp można znaleźć dzięki połączeniu dychotomii i twierdzenia Sturma . Teraz obliczenie wielomianowe GCD z i . Sprawdź, czy ma pierwiastek w (ponownie z twierdzeniem Sturma).I x P ∈ I x ′ P ∉ I P R P Q R IP. ja xP.∈ Ja x′P.∉ Ja P. R P. Q R ja
Jeśli nie myli, jest taki jak korzeń tylko wtedy, gdy i mają wspólny pierwiastek , która z kolei jest możliwe tylko wtedy, gdy jest pierwiastkiem . Zarówno zastosowanie twierdzenia Sturma, jak i GCD jest raczej szybkie (w rzeczywistości nie większe niż kwadratowe wielkości wielomianów).P Q I x P QR P. Q ja xP. Q
To tylko szkic, ale przekształcenie go w prawdziwy algorytm nie wymaga wiele , w rzeczywistości podejrzewam, że użycie Maple lub Mathematica uczyniłoby to trywialnym.
źródło