Kompleksitas jangkauan dalam sistem dinamik linier pada bidang terbatas


10

Misalkan A menjadi matriks di atas bidang hingga F2={0,1} dan x , y menjadi vektor ruang F2n . Saya tertarik pada kompleksitas komputasi untuk memutuskan apakah ada tN sedemikian rupa sehingga Atx=y , yaitu, dalam masalah keterjangkauan untuk sistem dinamika linier di atas bidang terbatas.

Masalahnya jelas dalam NP (tebak 0t<2n dan hitung At dalam waktu polinomial dengan kuadrat ulang). Saya dan rekan-rekan saya juga mampu membuktikan NP -completeness dari masalah terkait mendirikan apakah ada tN sedemikian sehingga Atxy , di mana ketimpangan componentwise.

Masalah ini tampaknya cukup alami, tetapi saya belum dapat menemukan referensi untuk kompleksitas komputasinya dalam literatur, mungkin karena saya tidak mengetahui terminologi yang tepat. Apakah Anda tahu jika masalah dengan kesetaraan adalah NP -complete atau jika itu sebenarnya dalam P ?


3
Seseorang dapat mengurangi kasus di mana dapat dibalik. Amati bahwa gambar A 1 , A 2 , ... adalah rantai subruang yang tidak bertambah, dan karenanya akhirnya menjadi ruang konstan W (pada kenyataannya dalam langkah n pertama ). Maka A adalah transformasi linear dibalik pada W . Satu dapat dengan mudah memeriksa kasus-kasus khusus ketika t = 1 , 2 , , n , setelah itu hanya untuk menyelesaikan masalah dengan A terbatas pada W dan xAA1,A2,WnAWt=1,2,,nAWxdigantikan oleh . Anx
Andrew Morgan

Jawaban:


4

Untuk kejelasan, saya akan menggeneralisasi pertanyaan Anda menjadi lebih dari karakteristik (dengan bidang dasar F q ) alih-alih kasus spesifik p = q = 2 . Saya akan mengambil p dan q sebagai konstanta tetap; Saya akan menyerahkan kepada pembaca untuk mencari tahu apa sebenarnya ketergantungan pada parameter ini, karena ada beberapa pengorbanan yang dapat dibuat. Hasil akhirnya di sini adalah bahwa masalah Anda kira-kira setara dengan masalah log diskrit untuk bidang hingga karakteristik p .p>0Fqp=q=2pqp

Untuk lebih spesifik, biarkan masalah log diskrit biasa lebih ekstensi dari akan, mengingat medan perpanjangan F dari F q , dan sebuah , b F , menemukan setiap bilangan bulat t sehingga a = b t , atau laporan bahwa tidak ada . Biarkan kuat masalah log diskrit lebih ekstensi dari F q akan, mengingat F , sebuah , b seperti sebelumnya, menemukan bilangan bulat z , m sehingga suatuFqFFqa,bFta=btFqF,a,bz,m untuk integer t iff, atau laporkan bahwa tidakada. Kemudian ada pengurangan berikut:a=btttt=z(modm)t

  • Ada pengurangan pemetaan deterministik dari ekstensi log over diskrit dari untuk masalah Anda.Fq

  • Ada algoritma yang efisien dan deterministik yang memecahkan masalah Anda ketika diberi akses ke oracle yang menghitung masalah log diskrit yang kuat pada ekstensi .Fq

Karenanya, saya menganggap tidak mungkin seseorang akan memposting bukti -kemampuan atau bukti bahwa masalah Anda ada di dalam waktu dekat.PNPP

Catatan: Masalah log diskrit yang kuat atas ekstensi dapat dikurangi menjadi bentuk yang lebih lemah berikut ini (meskipun tampaknya masih lebih kuat daripada masalah log diskret biasa): Diberikan bidang ekstensi dari , dan , temukan bilangan bulat terkecil, non-negatif sehingga . Ini mengikuti dari fakta bahwa urutan adalah satu ditambah non-negatif terkecil sehingga .F F q a,b F ta= b t bt b - 1 = b tFqFFqa,bFta=btbtb1=bt


Reduksi pertama: Klaimnya adalah bahwa masalah log diskrit biasa atas ekstensi dari pemetaan-mengurangi masalah ini Ini mengikuti fakta bahwa perkalian dalam adalah transformasi linear ketika kita melihat sebagai ruang vektor dimensi di atas . Karenanya pertanyaan dari bentuk atas menjadi atas , di mana adalah -dimensi vektor, dan adalahF q n F q n n F q a= b t F q n a = B t e F q a ,e nBn×n F q a aBbe 1 F q n p=q=2nFqFqnFqnnFqa=btFqna=BteFqa,enBn×nmatrix, all . Vektor dapat dengan mudah dihitung dari , dari , dan hanyalah representasi dari , yang dapat ditulis secara efisien . Ini tampaknya masih menjadi kasus sulit dari masalah log diskrit umum, bahkan dengan (tetapi tumbuh , tentu saja). Secara khusus, orang masih berlomba untuk melihat seberapa jauh mereka dapat menghitungnya.FqaaBbe1Fqnp=q=2n


Reduksi kedua: Klaimnya adalah bahwa masalah Anda berkurang menjadi masalah log diskrit yang kuat atas ekstensi . Pengurangan ini memiliki beberapa bagian untuk itu, jadi maafkan panjangnya. Biarkan input menjadi berdimensi vektor dan matriks , seluruh ; tujuannya adalah menemukan sehingga . nx,yn×nA F q ty= A t xFqnx,yn×nAFqty=Atx

Ide dasarnya adalah untuk menulis di Jordan canonical form (JCF), dari mana kita dapat mengurangi pengujian untuk masalah log diskrit yang kuat dengan beberapa aljabar langsung.y = A t xAy=Atx

Salah satu alasan untuk menggunakan bentuk kanonik di bawah kesamaan matriks adalah bahwa jika , maka . Karenanya kita dapat mentransformasikan ke , di mana sekarang berada dalam format yang jauh lebih bagus daripada arbitrer . JCF adalah bentuk yang sangat sederhana, yang memungkinkan sisa dari algoritma. Jadi mulai sekarang, asumsikan bahwa sudah ada di JCF, tetapi juga memungkinkan dan dapat memiliki entri dalam bidang ekstensi .A t = P - 1 J t P y = A t x ( P y ) = J t ( P x ) J A A x , y , A F qA=P1JPAt=P1JtPy=Atx(Py)=Jt(Px)JAAx,y,AFq

Catatan: Ada beberapa seluk-beluk yang timbul dari bekerja dengan JCF. Khususnya, saya akan berasumsi bahwa kita dapat melakukan operasi lapangan dalam setiap ekstensi (tidak peduli seberapa besar) dalam satu langkah waktu, dan bahwa kita dapat menghitung JCF secara efisien. a priori , ini tidak realistis, karena bekerja dengan JCF mungkin memerlukan bekerja di bidang ekstensi (bidang pemisahan polinom karakteristik) dari derajat eksponensial. Namun, dengan hati-hati, dan menggunakan fakta bahwa kami sedang mengerjakan bidang yang terbatas, kami dapat menghindari masalah ini. Secara khusus, kami akan mengaitkan dengan setiap Jordan memblokir bidang paling banyak lebih dariF n F q xy FFFFqFnFq sehingga semua entri di blok Jordan dan elemen-elemen terkait , semuanya tinggal di dalam . Kolom mungkin berbeda dari blok ke blok, tetapi menggunakan `` representasi campuran' 'ini memungkinkan untuk deskripsi JCF yang efisien, yang juga dapat ditemukan secara efisien. Algoritme yang dijelaskan dalam sisa bagian ini hanya perlu bekerja dengan satu blok pada satu waktu, jadi selama ia melakukan operasi lapangannya dalam bidang terkait , algoritma akan efisien. [komentar akhir]xyFFF

Penggunaan JCF memberi kita persamaan bentuk berikut, dengan setiap persamaan yang sesuai dengan blok Jordan:

[y1y2y3yk1yk]=[λ1λ1λ1λ1λ]t[x1x2x3xk1xk]

Algoritma akan menangani setiap blok secara terpisah. Dalam kasus umum, untuk setiap blok, kami akan memiliki pertanyaan untuk oracle log diskrit kami yang kuat, dari mana oracle akan memberi tahu kami kondisi modularitas, . Kami juga akan mendapatkan set sehingga harus tahan . Setelah memproses semua blok, kita perlu memeriksa bahwa ada pilihan yang memenuhi konjungsi dari semua kondisi ini. Hal ini dapat dilakukan dengan memastikan ada unsur umum di seluruh set sehingga persamaan danS { 0 , 1 , , p - 1 } s S [ t = st=z(modm)S{0,1,,p1} t s S t = ssS[t=s(modp)]tsSt = z jt=s(modp)jt=zj(modmj)semua puas secara bersamaan, di mana berkisar di atas blok.j

Ada juga beberapa kasus khusus yang muncul sepanjang prosedur. Dalam kasus ini, kita akan mendapatkan kondisi dalam bentuk untuk beberapa nilai , atau dari bentuk untuk beberapa tertentu bilangan bulat , dari blok tertentu, atau kita mungkin bahkan menemukan bahwa ada bisa eksis . Ini dapat dimasukkan ke dalam logika untuk kasus umum tanpa masalah.t = s s tt>t=sst

Kami sekarang menjelaskan prosedur untuk menangani setiap blok Jordan. Perbaiki blok seperti itu.

Mulailah dengan fokus hanya pada koordinat terakhir di blok. Kondisi mengharuskan . Dengan kata lain, ini adalah contoh dari masalah log diskrit dalam beberapa ekstensi bidang . Kami kemudian menggunakan oracle untuk menyelesaikannya, yang menghasilkan tidak ada solusi, atau memberikan kondisi modularitas pada . Jika "tidak ada solusi" dikembalikan, kami kembali menunjukkan itu. Kalau tidak, kita mendapatkan kondisi , yang setara dengan .y k = λ t x k F q t t = zy=Atxyk=λtxkFqty k = λ t x kt=z(modm)yk=λtxk

Untuk menangani koordinat lain, kita mulai dengan rumus berikut (lihat, misalnya, di sini ): xk=0yk=λtxkyk=0k-1xy(k-1)×(k-1)xk0

[λ1λ1λ1λ1λ]t=[λt(t1)λt1(t2)λt2(tk1)λtk+1λt(t1)λt1(tk2)λtk+2λt(t1)λt1λt]
\ ddots & \ ddots & \ vdots \\ & & & & & lambda ^ t & \ binom {t} {1} \ lambda ^ {t-1} \\ & & & & & & & lambda ^ t \ end {bmatrix } Pertama, mari kita selesaikan kasus inixk=0 . Karena kita sudah memiliki kondisi modularitas yang menyiratkan , kita dapat mengasumsikan bahwa juga. Tapi kemudian kita bisa mengurangi untuk fokus pada entri pertama dan , dan submatrix kiri dari blok Jordan. Jadi mulai sekarang, asumsikan .yk=λtxkyk=0k1xy(k1)×(k1)xk0

Kedua, kami akan menangani kasus di mana . Dalam hal ini, kekuatan blok Jordan memiliki bentuk khusus, dan memaksa baik untuk beberapa , atau , tanpa syarat lain. Saya tidak akan mengulangi kasus, tetapi cukup untuk mengatakan bahwa masing-masing dapat diperiksa secara efisien. (Atau, kita dapat mengurangi kasus di mana tidak dapat dibalik; lihat komentar saya pada pertanyaan.)t = z z k t > k Aλ=0t=zzkt>kA

Akhirnya, kami sampai pada kasus umum. Karena kita sudah memiliki kondisi modularitas yang menyiratkan bahwa , kita dapat mengasumsikan kondisi itu berlaku, dan menggunakan sebagai . Secara umum, kita dapat menggunakan untuk mewakili . Jadi kita perlu memeriksa apakah sistem berikut ini berlaku untuk beberapa pilihan : yk=λtxkykxk1λtykxk1λzλtzt

[y1y2y3yk1yk]=[ykxk1(t1)ykxk1λ1(t2)ykxk1λ2(tk1)ykxk1λ(k1)ykxk1(t1)ykxk1λ1(tk2)ykxk1λ(k2)ykxk1(t1)ykxk1λ1ykxk1][x1x2x3xk1xk]
y_kx_k ^ {- 1} \ end {bmatrix} \ begin {bmatrix} x_1 \\ x_2 \\ x_3 \\ \ vdots \\ x_ {k-1} \\ x_ {k} \\ \ end {bmatrix} Perhatikan itu apakah persamaan memegang hanya bergantung pada ; ini karena ketergantungan pada hanya polinomial, t(modp)tt harus bilangan bulat, dan persamaan di atas adalah bidang karakteristik . Karenanya kita bisa mencoba setiap nilai secara terpisah. Himpunan akan kami kembalikan hanyalah pilihan yang membuat sistem puas.pt{0,1,,p1}St

Jadi sekarang, kecuali untuk beberapa kasus khusus, per-blok telah menemukan kondisi modularitas , dan himpunan sehingga salah satu dari harus tahan untuk beberapa . Kondisi ini setara dengan dalam blok Jordan spesifik ini. Jadi kami mengembalikan ini dari sub-prosedur. Kasus-kasus khusus menyimpulkan bahwa tidak ada (di mana kasus prosedur segera mengembalikan indikasi itu), atau kita memiliki kondisi modularitas dan beberapa kondisi khusus seperti untuk integer , atau untuk bilangan bulatt=a(modm)St=s(modp)sSy=Atxtt=a(modm)t=sst> . Bagaimanapun, kondisi yang terlibat semuanya setara dengan dalam blok Jordan ini. Jadi, seperti yang disebutkan di atas, prosedur hanya mengembalikan kondisi ini.y=Atx

Ini menyimpulkan spesifikasi dari sub-prosedur per-blok, dan algoritma secara keseluruhan. Itu benar dan efisiensi mengikuti dari diskusi sebelumnya.


Seluk-beluk dengan menggunakan JCF dalam reduksi kedua: Seperti disebutkan dalam reduksi kedua, ada beberapa kehalusan yang muncul dari bekerja dengan JCF. Ada beberapa pengamatan untuk mengurangi masalah ini:

  • Ekstensi bidang terbatas adalah normal . Ini berarti bahwa jika adalah sebuah tereduksi polinomial lebih , maka setiap perpanjangan mengandung akar mengandung semua akar . Dengan kata lain, bidang pemisah dari polinomial dengan derajat dapat direduksi hanya memiliki derajat atas .PFqFqPPPddFq

  • Ada generalisasi bentuk kanonik Jordan, yang disebut bentuk kanonik rasional primer (PRCF), yang tidak memerlukan ekstensi lapangan untuk dituliskan. Secara khusus, jika adalah matriks dengan entri dalam , maka kita dapat menulis untuk beberapa matriks dengan entri dalam , di mana lebih lanjut ada di PRCF. Selain itu, jika kita berpura-pura bahwa entri langsung di bidang extending yang berisi semua nilai eigen , makaAFqA=P1QPP,QFqQAFFqAQsebenarnya akan berada di JCF. Dengan demikian kita dapat melihat komputasi JCF sebagai kasus khusus menghitung PRCF.A

  • Dengan menggunakan bentuk PRCF, kita dapat menghitung faktor JCF sebagaiA

    1. menghitung PRCF dari overAFq

    2. menghitung PRCF dari setiap blok (meminjam notasi dari artikel Wikipedia) di RRCF , pada bidang ekstensi , di mana dipilih untuk memuat semua nilai eigenCAFFC

    Keuntungan utama dengan faktorisasi ini adalah bahwa polinomial karakteristik dari blok semuanya akan dapat direduksi , dan karenanya, dengan pengamatan pertama kami, kami dapat memilih untuk memiliki derajat ukuran (yang paling banyak ) lebih dari . Kelemahannya adalah bahwa sekarang kita harus menggunakan bidang ekstensi yang berbeda untuk mewakili setiap blok JCF, sehingga representasi tidak khas dan rumit.CFCnFq

Dengan demikian, diberikan kemampuan untuk menghitung PRCF secara efisien, kita bisa menghitung encoding cocok dari JCF efisien, dan encoding ini adalah agar bekerja dengan setiap blok tertentu dari JCF dapat dilakukan dalam bidang perluasan tingkat paling lebih .nFn

Sedangkan untuk menghitung PRCF secara efisien, makalah " A Canonical Form Algorithm " (KR Matthews, Math. Bohemica 117 (1992) 315-324) memberikan algoritma yang efisien untuk menghitung PRCF ketika faktorisasi polinomial karakteristik dari diketahui. . Untuk karakteristik tetap (seperti yang kita miliki), memfaktorkan polinomial univariat pada bidang terbatas dapat dilakukan dalam waktu polinomial deterministik (lihat misalnya " Pada Algoritma Faktorisasi Baru untuk Polinom atas Bidang Hingga " (H. Niederreitter dan R. Gottfert, Math. Of Komputasi 64 (1995) 347-353).), Sehingga PRCF dapat dihitung secara efisien.A


Bisakah JCF dihitung secara efisien? Bagaimanapun, keberadaannya mungkin perlu memperluas lapangan.
Emil Jeřábek

@ EmilJeřábek Terima kasih- Saya kira saya telah bekerja di bawah asumsi implisit bahwa itu mudah, tetapi saya tidak benar-benar tahu secara spesifik. Tampaknya sangat terkait dengan memfaktorkan polinomial univariat atas bidang yang terbatas, yang dapat dilakukan cukup efisien untuk tujuan di atas, setidaknya menurut Wikipedia . ...
Andrew Morgan

Jadi tebakan saya adalah bahwa JCF dapat ditemukan secara efisien, tetapi saya tidak sepenuhnya yakin. Anda menyebutkan harus memperluas bidang-ini diperlukan untuk pengurangan (log diskrit atas bidang hingga ukuran konstan adalah mudah), jadi seharusnya tidak mengejutkan. Saya khawatir tentang tingkat ekstensi meskipun- sementara nilai eigen hanya memiliki tingkat , nilai dan adalah kombinasi linear dari nilai daya eigen, jadi mereka mungkin perlu hidup dalam bidang ukuran. Saya akan mencatat kemungkinan jebakan ini dalam jawaban saya, meskipun saya pikir itu masih menyumbang ide yang cukup untuk tinggal di sekitar. nxiyin!
Andrew Morgan

Baik. Elemen-elemen hidup dalam bidang pemisahan dari polinomial karakteristik dari matriks, yang mungkin merupakan polinomial sembarang dari derajat n, sehingga bidang pemisah mungkin memiliki tingkat setinggi kira-kira (jika perhitungan saya benar). Tapi mungkin ini bisa dielakkan entah bagaimana. Mari kita pisahkan faktor char (bahkan faktorisasi derajat yang berbeda sudah cukup). Bisakah kita entah bagaimana mengidentifikasi ruang eigens yang sesuai dengan akar setiap faktor? Artinya, alih-alih JCF penuh, kami akan mendapatkan matriks diagonal blok di atas bidang asli, di mana setiap blok akan memiliki nilai eigen ...exp(nlogn)
Emil Jeřábek

... dalam perpanjangan gelar paling banyak . Maka kami mungkin dapat memproses setiap blok secara terpisah. (Ini hanya gagasan yang kabur, saya belum mencoba menyelesaikannya.)n
Emil Jeilábek
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.