Algoritma Berlekamp–Rabin

Foto Algoritma Berlekamp–Rabin

Dalam teori bilangan, algoritma Berlekamp–Rabin adalah metode probabilistik pencarian akar dari polinomial pada medan $\mathbb {Z} _{p}$. Metode ini ditemukan oleh Elwyn Berlekamp pada tahun 1970 sebagai tambahan algoritma untuk faktorisasi polinomial pada medan berhingga. Algoritma kemudian dimodifikasi oleh Rabin untuk medan berhingga sebarang pada tahun 1979. Sebelum Berlekamp, metode ini juga ditemukan secara terpisah oleh peneliti lain.

Sejarah

Metode ini diusulkan oleh Elwyn Berlekamp dalam karyanya tahun 1970 tentang faktorisasi polinomial pada medan berhingga. Karya aslinya tidak memiliki bukti kebenaran formal dan kemudian disempurnakan dan dimodifikasi untuk medan berhingga sebarang oleh Michael Rabin. Pada tahun 1986, René Peralta mengusulkan algoritma serupa untuk mencari akar kuadrat $\mathbb {Z} _{p}$. Perumuman metode Peralita untuk persamaan kubik dibuat pada tahun 2000.

Pernyataan masalah

Misalkan $p$ adalah bilangan prima ganjil, dan misalkan pula polinomial ${\textstyle f(x)=a_{0}+a_{1}x+\cdots +a_{n}x^{n}}$ atas medan $\mathbb {Z} _{p}$ dari modulo sisa $p$. Tujuan pernyataan masalah ini adalah bahwa algoritma Berlekamp harus menemukan semua $\lambda$ di $\mathbb {Z} _{p}$ sehingga ${\textstyle f(\lambda )=0}$ di $\mathbb {Z} _{p}$.

Algoritma

Pengacakan

Misalkan ${\textstyle f(x)=(x-\lambda _{1})(x-\lambda _{2})\cdots (x-\lambda _{n})}$. Menemukan semua akar polinomial ini sama saja dengan mencari faktorisasinya menjadi faktor linear. Untuk menemukan faktorisasi tersebut, cukup membagi polinomial menjadi dua pembagi non-trivial dan memfaktorkannya secara rekursif. Lebih lanjut, misalkan polinomial ${\textstyle f_{z}(x)=f(x-z)=(x-\lambda _{1}-z)(x-\lambda _{2}-z)\cdots (x-\lambda _{n}-z)}$, untuk $z$ setiap anggota $\mathbb {Z} _{p}$. Jika polinomial tersebut dinyatakan sebagai hasilkali $f_{z}(x)=p_{0}(x)p_{1}(x)$, maka dalam bentuk polinomial awal, mengartikan bahwa $f(x)=p_{0}(x+z)p_{1}(x+z)$. Fungsi yang merupakan menyediakan faktorisasi dibutuhkan dari $f(x)$.

Klasifikasi elemen $\mathbb {Z} _{p}$

Menurut kriteria Euler, untuk setiap monomial $(x-\lambda )$, berlaku tepat salah satu dari sifat-sifat berikut:

1. Monomial sama dengan $x$ jika $\lambda =0$,
2. Pembagian monomial ${\textstyle g_{0}(x)=(x^{(p-1)/2}-1)}$ jika $\lambda$ adalah residu kuadrat modulo $p$,
3. Pembagian monomial ${\textstyle g_{1}(x)=(x^{(p-1)/2}+1)}$ jika $\lambda$ adalah non-residul kuadrat modulo $p$.

Jadi, jika $f_{z}(x)$ tidak habis dibagi $x$, yang dapat diperiksa secara terpisah, maka $f_{z}(x)$ sama dengan hasilkali dari faktor persekutuan terbesar $\operatorname {FPB} (f_{z}(x);g_{0}(x))$ dan $\operatorname {FPB} (f_{z}(x);g_{1}(x))$.

Metode Berlekamp

Sifat-sifat di atas mengarah pada algoritma berikut:

1. Hitung secara eksplisit koefisien $f_{z}(x)=f(x-z)$,
2. Hitung sisa ${\textstyle x,x^{2},x^{2^{2}},x^{2^{3}},x^{2^{4}},\ldots ,x^{2^{\lfloor \log _{2}p\rfloor }}}$ modulo $f_{z}(x)$ dengan mengkuadratkan polinomial dan mengambil sisa modulo $f_{z}(x)$,
3. Menggunakan eksponensiasi yang dikuadratkan dan polinomial yang dihitung pada langkah sebelumnya, hitung sisa ${\textstyle x^{(p-1)/2}}$ modulo ${\textstyle f_{z}(x)}$,
4. Jika ${\textstyle x^{(p-1)/2}\not \equiv \pm 1{\pmod {f_{z}(x)}}}$, maka $\operatorname {FPB}$ yang disebutkan diatas memberikan faktorisasi non-trivial dari $f_{z}(x)$,
5. Jika tidak, maka semua akar $f_{z}(x)$ adalah residu atau non-residu secara bersamaan dan harus memilih $z$ yang lain.

Jika $f(x)$ habis dibagi oleh setiap polinomial primitif $g(x)$ atas $\mathbb {Z} _{p}$ maka saat menghitung $\operatorname {FPB}$ dengan $g_{0}(x)$ dan $g_{1}(x)$ akan diperoleh faktorisasi non-trivial dari $f_{z}(x)/g_{z}(x)$, sehingga algoritma memungkinkan untuk menemukan semua akar polinomial sebarang atas $\mathbb {Z} _{p}$.

Akar kuadrat modular

Misalkan persamaan ${\textstyle x^{2}\equiv a{\pmod {p}}}$ mempunyai anggota $\beta$ dan $-\beta$ sebagai akarnya. Solusi persamaan ini ekuivalen dengan faktorisasi polinomial ${\textstyle f(x)=x^{2}-a=(x-\beta )(x+\beta )}$ atas $\mathbb {Z} _{p}$. Dalam masalah kasus istimewa ini, cukup hitung $\operatorname {FPB} (f_{z}(x);g_{0}(x))$ saja. Untuk polinomial tersebut, ada tepat satu dari sifat-sifat yang berlaku sebagai berikut:

1. FPB sama dengan $1$, yang berarti bahwa $z+\beta$ dan $z-\beta$ merupakan non-residu kuadratik,
2. FPB sama dengan $f_{z}(x)$, yang berarti kedua bilangan tersebut merupakan residu kuadratik,
3. FPB sama dengan $(x-t)$, yang berarti tepat salah satu dari bilangan tersebut merupakan residu kuadrat.

Pada kasus ketiga, FPB sama dengan $(x-z-\beta )$ atau $(x-z+\beta )$. Hal ini memungkinkan untuk menulis solusi sebagai ${\textstyle \beta =(t-z){\pmod {p}}}$.

Contoh

Asumsi bahwa persamaan ${\textstyle x^{2}\equiv 5{\pmod {11}}}$ dapat diselesaikan. Caranya adalah dengan memfaktorkan $f(x)=x^{2}-5=(x-\beta )(x+\beta )$. Nilai $z$ dapat dibagi menjadi beberapa kasus dengan memisalkannya sebagai berikut:

1. Misalkan $z=3$, maka $f_{z}(x)=(x-3)^{2}-5=x^{2}-6x+4$. Jadi, $\operatorname {FPB} (x^{2}-6x+4;x^{5}-1)=1$. Karena bilangan $3\pm \beta$ merupakan non-residu kuadrat, maka perlu dicari nilai $z$ lain.
2. Misalkan $z=2$, maka $f_{z}(x)=(x-2)^{2}-5=x^{2}-4x-1$. Jadi, ${\textstyle \operatorname {FPB} (x^{2}-4x-1;x^{5}-1)\equiv x-9{\pmod {11}}}$. Karena ${\textstyle x-9=x-2-\beta }$, maka $\beta \equiv 7{\pmod {11}}$ dan ${\textstyle -\beta \equiv -7\equiv 4{\pmod {11}}}$.

Hasil pemeriksaan yang dilakukan secara manual memperlihatkan bahwa ${\textstyle 7^{2}\equiv 49\equiv 5{\pmod {11}}}$ dan ${\textstyle 4^{2}\equiv 16\equiv 5{\pmod {11}}}$.

Bukti kebenaran

Algoritma menemukan faktorisasi $f_{z}(x)$ dalam semua kasus, kecuali ketika semua bilangan $z+\lambda _{1},z+\lambda _{2},\ldots ,z+\lambda _{n}$ adalah residu kuadrat atau non-residu secara bersamaan. Menurut teori siklotomi, peluang kejadiannya untuk kasus ketika $\lambda _{1},\ldots ,\lambda _{n}$ juga merupakan semua residu atau non-residu secara bersamaan (yaitu, ketika $z=0$ tidak berhasil) diestimasi sebagai $2^{-k}$, untuk $k$ adalah jumlah nilai yang berbeda dalam $\lambda _{1},\ldots ,\lambda _{n}$. Dengan cara ini, bahkan untuk kasus galat $k=1$ dan $f(x)=(x-\lambda )^{n}$, maka estimasi peluang galatnya adalah $1/2$, dan untuk kasus akar kuadrat modular, probabilitas galat setidaknya bernilai $1/4$.

Kompleksitas

Misalkan suatu polinomial merupakan polinomial berderajat $n$, maka kompleksitas algoritma dapat diturunkan sebagai berikut:

1. Karena menurut teorema binomial mengatakan bahwa ${\textstyle (x-z)^{k}=\sum \limits _{i=0}^{k}{\binom {k}{i}}(-z)^{k-i}x^{i}}$, maka dapat dilakukan transisi dari $f(x)$ ke $f(x-z)$ dalam waktu $O(n^{2})$.
2. Perkalian polinomial dan mengambil sisa dari satu polinomial yang modulo dengan yang lain dapat dilakukan di ${\textstyle O(n^{2})}$. Jadi, perhitungan ${\textstyle x^{2^{k}}{\bmod {f}}_{z}(x)}$ dilakukan di ${\textstyle O(n^{2}\log p)}$.
3. Eksponensial biner bekerja di $O(n^{2}\log p)$.
4. Mengambil $\operatorname {FPB}$ dari dua polinomial melalui algoritma Euklides yang bekerja di $O(n^{2})$.

Jadi, seluruh prosedur dapat dilakukan di $O(n^{2}\log p)$. Dengan menggunakan transformasi Fourier cepat dan algoritma setengah-FPB, maka kompleksitas algoritmanya dapat ditingkatkan menjadi $O(n\log n\log pn)$. Untuk kasus akar kuadrat modular, derajatnya adalah $n=2$, sehingga seluruh kompleksitas algoritma dalam kasus tersebut dibatasi dengan $O(\log p)$ per iterasi.

Referensi

Templat:Algoritma teori bilangan