Sonlu cisimler üzerinde doğrusal dinamik sistemlerde erişilebilirliğin karmaşıklığı


10

Let A sonlu alan üzerinde bir matris F2={0,1} ve x , y alan vektörleri olduğu F2n . olup olmadığına karar vermenin hesaplama karmaşıklığıyla ilgileniyorum, tNöyle ki Atx=y , yani sonlu alanlar üzerinden doğrusal dinamik sistemler için ulaşılabilirlik probleminde.

Sorun olduğu açıktır NP (tahmin 0t<2n ve işlem At tekrar kare alma ile polinom zamanda). Ben ve arkadaşları da kanıtlamak mümkün NP var olup olmadığını kurulması ile ilgili bir sorun -completeness tN , öyle ki Atxy , componentwise eşitsizlik.

Bu problem oldukça doğal gözüküyor, ancak muhtemelen kesin terminolojinin farkında olmadığım için literatürdeki hesaplama karmaşıklığına atıfta bulunamadım. Eğer eşitlik ile sorun olup olmadığını biliyor musunuz NP -tamamlamak yoksa aslında ise P ?


3
Kişi tersinir olduğu duruma düşebilir . Görüntüleri gözlemleyin A 1 , A 2 , ... bölme odasının bir artmayacak zinciridir, ve sonuçta da sabit bir boşluk olur B (birinci içinde aslında , n ) adımları. O zaman A , W üzerinde ters çevrilebilir bir doğrusal dönüşümdür . Biri t = 1 , 2 , , n olduğunda özel durumları kolayca kontrol edebilir , daha sonra sorunu W ve x ile sınırlı A ile çözmeye devam eder.AA1,A2,WnAWt=1,2,,nAWxile ikame . Anx
Andrew Morgan

Yanıtlar:


4

Açıklık getirmek için, sorunuzu p = q = 2 özel durumu yerine (taban alanı F q ile ) karakteristiği üzerinde genelleştireceğim . P ve q'yu sabit bir sabit olarak alacağım ; Bu parametrelere tam bağımlılığın ne olduğunu anlamak için okuyucuya bırakacağım, çünkü yapılabilecek bazı ödünleşmeler var. Buradaki sonuç, probleminizin kabaca p karakteristik sonlu alanları için ayrık günlük problemine eşdeğer olmasıdır .p>0Fqp=q=2pqp

Daha spesifik olması açısından, izin sıradan farklı günlük sorunu uzantıları üzerinde , uzatma alanı verilebilir F arasında K q ve bir , b F , arayan bir tamsayı T , böylece bir = b t , ya da rapor yok varolduğunu . Let güçlü uzantıları üzerinde ayrı ayrı günlük sorun F q , verilecek F , bir , b tamsayılardır bulmak, daha önce olduğu gibi , z , m , böylece birFqFFqa,bFta=btFqF,a,bz,m tamsayı t iffveyayokolduğunu. Sonra aşağıdaki indirimler var:a=btttt=z(modm)t

  • uzantılarının ayrık oturum sorununuza belirleyici bir eşleme azalması var .Fq

  • uzantıları üzerindeki güçlü ayrık günlük problemini hesaplayan bir kehanete erişim verildiğinde probleminizi çözen etkili, deterministik bir algoritma vardır .Fq

Buna göre, birisinin sertliğinin bir kanıtı ya da yakın gelecekte probleminizin içinde olduğuna dair bir kanıt yayınlamasının olası olmadığını .PNPP

Not: uzantıları üzerinde güçlü bir ayrık günlük sorun (hâlâ sıradan ayrı günlük sorun daha görünüşte güçlü) Aşağıdaki görünüşte zayıf forma Turing azaltılabilir: bir uzantı alanı göz önüne alındığında, bölgesinin ve bulmak az, negatif olmayan bir tam sayı , böylece . Bu, sırasının bir artı en küçük negatif olmayan olduğu gerçeğinden kaynaklanır, böylece .F F q a,b F ta= b t bt b - 1 = b tFqFFqa,bFta=btbtb1=bt


İlk azaltma: İddia, eşlemesinin uzantıları üzerindeki olağan ayrık günlük sorununun bu soruna azalmasıdır. Bu in çarpma aslında aşağıdaki görüntüleyip zaman doğrusal dönüşümdür bir şekilde fazla boyutlu vektör alanı . Bu nedenle bir şekilde bir soru fazla olur fazla , olan boyutlu vektörler ve bir birF 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×nmatris, her yerde . Vektör kolayca hesaplanabilir , den ve sadece temsilidir verimli yazılı olabilir, . Bu, ile bile genel ayrık günlük probleminin zor bir durumu gibi görünmektedir (ancak elbette büyüyor ). Özellikle, insanlar hala bunu ne kadar hesaplayabileceklerini görmek için yarışıyorlar.FqaaBbe1Fqnp=q=2n


İkinci azaltma: İddia, probleminizin uzantıları üzerindeki güçlü ayrık günlük sorununa . Bu azaltmanın birkaç parçası vardır, bu yüzden uzunluğu affedin. Girdi, -boyutlu vektörler ve matrisi , hepsi ; Amaç bulmak, böylece . nx,yn×nA F q ty= A t xFqnx,yn×nAFqty=Atx

Temel fikir yazmak için olarak meydana kanonik Ürdün olan test etmekte azaltabilir (JCF) bazı basit cebir güçlü ayrı günlük sorununa.y = A t xAy=Atx

Matrislerin benzerliği altında kanonik bir form kullanmanın bir nedeni, , . Bu nedenle biz dönüştürmek için şimdi, rasgele daha çok güzel bir formatında . JCF, algoritmanın geri kalanını mümkün kılan özellikle basit bir formdur. Bundan böyle, zaten JCF'de olduğunu varsayın , ancak ve uzantı alanında girişlere sahip olmasına izin .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

Not: JCF ile çalışmaktan kaynaklanan bazı incelikler vardır. Özellikle, bir kerelik adımda (ne kadar büyük olursa olsun) herhangi bir uzantısında saha işlemleri yapabileceğimizi ve verimli bir şekilde hesaplayabileceğimizi varsayacağım. a priori , bu gerçekçi değildir, çünkü JCF ile çalışmak üstel derecede bir uzama alanında (karakteristik polinomun bölünme alanı) çalışmayı gerektirebilir. Bununla birlikte, biraz dikkatli ve sınırlı bir alan üzerinde çalıştığımız gerçeğini kullanarak, bu sorunları atlatabiliriz. Özellikle, her Ürdün bloğu ile derece alanıyla en fazla üzerindenF n F q xy FFFFqFnFq böylece Ürdün bloğundaki tüm girdiler ve , karşılık gelen öğeleri içinde yaşar . alanı bloktan bloğa farklılık gösterebilir, ancak bu `` karışık gösterim' 'kullanılması, JCF'nin verimli bir şekilde tanımlanmasına izin verir, ayrıca verimli bir şekilde bulunabilir. Bu bölümün geri kalanında açıklanan algoritmanın bir kerede yalnızca bir blokla çalışması gerekir, böylece ilişkili alan içindeki alan işlemlerini yaptığı sürece algoritma etkili olacaktır. [son açıklama]xyFFF

JCF kullanımı, her bir Jordan bloğuna karşılık gelen aşağıdaki formdaki denklemleri verir:

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

Algoritma her bloğu ayrı ayrı ele alacaktır. Genel durumda, her blok için, güçlü ayrık günlük kehanetimiz için bir sorguya sahip olacağız, oracle bize bir modülerlik durumu söyleyecektir, . Biz de bir dizi alırsınız , böylece zorunluluk beklemeye . Tüm blokları işledikten sonra , tüm bu koşulların birleşimlerini karşılayan bir seçeneğinin olup olmadığını kontrol etmeliyiz . Bu, tüm setlerinde ortak bir eleman olduğundan emin olarak yapılabilir, böylece veS { 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)hepsi aynı anda tatmin olur, burada bloklar arasında değişir.j

Prosedür boyunca ortaya çıkan bazı özel durumlar da vardır. Bu durumlarda, biz formun koşullarını alırsınız bazı değeri veya form içinde bazı özel tamsayı için , belirli bloklardan ya biz bile hiçbir fark edebilirsiniz var olabilir . Bunlar sorun olmadan genel durum için mantığa dahil edilebilir.t = s s tt>t=sst

Şimdi her bir Ürdün bloğunu ele almak için alt prosedürü açıklayacağız. Böyle bir bloğu düzeltin.

Bloktaki sadece son koordinat üzerine odaklanarak başlayın. koşulu , olmasını gerektirir . Başka bir deyişle, bazı alan uzantılarındaki ayrık günlük sorununun bir . Daha sonra, bunu çözmek için bir çözüm kullanıyoruz, bu da çözüm getirmiyor ya da üzerinde modülerlik koşulu veriyor . "Çözüm yok" döndürülürse, bunu belirterek geri döneriz. Aksi takdirde, ile eşdeğer bir .y k = λ t x k F q t t = zy=Atxyk=λtxkFqty k = λ t x kt=z(modm)yk=λtxk

Diğer koordinatları ele almak için aşağıdaki formülle başlarız (bakınız, örneğin, burada ): 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 } Öncelikle, şu davaya bakalım:xk=0 . anlamına gelen modülerlik koşuluna zaten sahip olduğumuz için, da varsayabiliriz . Ama sonra sadece ve ilk girişlerine ve Ürdün bloğunun sol üst matrisine odaklanmayı azaltabiliriz. Bundan sonra, olduğunu .yk=λtxkyk=0k1xy(k1)×(k1)xk0

İkincisi, olan davayı ele alacağız . Bu durumda, Ürdün bloğunun yetkileri özel bir forma sahiptir ve başka bir koşul olmaksızın bazı için veya başka bir zorlar . Davaları şaşırtmayacağım, ancak her birinin verimli bir şekilde kontrol edilebileceğini söylemek için yeterli. (Alternatif olarak, tersinir olduğu duruma düşebiliriz ; soru hakkındaki yorumuma bakın.)t = z z k t > k Aλ=0t=zzkt>kA

Son olarak, genel davaya ulaşıyoruz. anlamına gelen modülerlik koşuluna zaten sahip olduğumuzdan, bu varsayabilir ve için bir stand-in olarak kullanabiliriz . Daha genel olarak, ı temsil etmek için kullanabiliriz . Bu nedenle, aşağıdaki sistemin bazı seçenekleri için uygun olup olmadığını kontrol etmeliyiz : 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} Bunu gözlemleyin denklemin mevcut olup olmadığı sadece ; çünkü bağımlılık sadece polinomdur, t(modp)tt bir tamsayı olmalı ve yukarıdaki denklemler karakteristiği alanının üzerindedir . Bu nedenle her bir değerini ayrı ayrı deneyebiliriz . Geri döneceğimiz , sadece sistemin tatmin olduğu seçimleridir .pt{0,1,,p1}St

Bu nedenle, şimdi, bazı özel durumlar dışında, blok başına alt modülerlik koşulu ve bir set böylece biri bazı . Bu koşullar, bu özel Jordan bloğundaki ile eşdeğerdir . Bu yüzden bunları alt prosedürden iade ediyoruz. Özel durumlar ya bulunamayacağı sonucuna varır (bu durumda alt prosedür derhal bunun bir göstergesini döndürür) ya da modülerlik koşulu ve tamsayı için gibi bazı özel koşullarımız var veya bazı tamsayılar içint=a(modm)St=s(modp)sSy=Atxtt=a(modm)t=sst> . Her durumda, söz konusu koşulların tümü bu Jordan bloğundaki eşdeğerdir . Dolayısıyla, yukarıda belirtildiği gibi, alt prosedür sadece bu koşulları döndürür.y=Atx

Bu, blok başına alt prosedürün ve bir bütün olarak algoritmanın özelliklerini sonuçlandırır. Doğruluk ve verimlilik önceki tartışmadan sonra gelir.


İkinci azaltımda JCF kullanan incelikler : İkinci azaltmada belirtildiği gibi, JCF ile çalışmaktan kaynaklanan bazı incelikler vardır. Bu sorunları hafifletmek için birkaç gözlem var:

  • Sonlu cisimlerin uzantıları normaldir . Bu, üzerinden indirgenemez bir polinomsa , kökünü içeren uzantısının tüm köklerini içerdiği anlamına gelir . Diğer bir deyişle, bir indirgenemez polinom bölme alan derecesi derecesi, sadece var fazla .PFqFqPPPddFq

  • Ürdün kanonik formunun, birincil rasyonel kanonik form (PRCF) adı verilen ve alan uzantılarının yazılmasını gerektirmeyen bir genellemesi vardır . Özellikle, eğer girdilerle matrisidir , o zaman geç olabilir bir matris için girdilerle , ayrıca PRCF'de. Buna ek olarak, girişlerin iddia ise alanı canlı uzanan tüm öz içeren , ardındanAFqA=P1QPP,QFqQAFFqAQaslında JCF'de olacak. Böylece JCF'sinin hesaplanmasını PRCF'nin hesaplanması için özel bir durum olarak görebiliriz.A

  • PRCF şeklini kullanarak, bir JCF işlem faktör olabilir olarakA

    1. arasında PRCF işlem fazlaAFq

    2. her bir blok bir PRCF işlem arasında PRCF içinde (Ara makalesinden gösterimde borçlanma) bir uzantı alanı üzerinde, , tüm öz içeren seçilirCAFFC

    Bu çarpanlara temel avantajı blok karakteristik polinomlar olmasıdır bütün olacak indirgenemez ve dolayısıyla ilk gözlem, biz seçebilir derecesi büyüklüğü için , en fazla olduğu ( ) . Dezavantajı, şimdi JCF'nin her bir bloğunu temsil etmek için farklı genişletme alanları kullanmamız gerektiğinden temsil, atipik ve karmaşıktır.CFCnFq

Böylece, PRCF'yi verimli bir şekilde hesaplama yeteneği göz önüne alındığında, JCF'nin uygun bir kodlamasını verimli bir şekilde hesaplayabiliriz ve bu kodlama, JCF'nin herhangi bir belirli bloğu ile çalışmanın, en fazla derece üzerinde bir genişleme alanında yapılabileceği .nFn

Etkin bir şekilde, kağıt "PRCF hesaplamak için de bir rasyonel kanonik formu algoritması karakteristik polinomunun çarpanlara zaman" (KR Matthews, 117 (1992 Math. Bohemica) 315-324) PRCF hesaplamak için etkili bir algoritma verir bilinir . Sabit karakteristikler için (sahip olduğumuz gibi), tek değişkenli polinomları sonlu alanlar üzerinde çarpanlara ayırma, deterministik polinom zamanında yapılabilir (bkz. Örneğin " Sonlu Alanlar Üzerindeki Polinomlar İçin Yeni Bir Çarpanlara Ayırma Algoritması " (H. Niederreitter ve R. Gottfert, Math. Hesaplama 64 (1995) 347-353).), Böylece PRCF verimli bir şekilde hesaplanabilir.A


JCF verimli bir şekilde hesaplanabilir mi? Her neyse, varlığı alanın genişletilmesini gerektirebilir.
Emil Jeřábek

@ EmilJeřábek Teşekkürler- Sanırım kolay olduğu örtülü bir varsayımla çalıştım, ama aslında detayları bilmiyorum. En azından Wikipedia'ya göre , yukarıdaki amaçlar için yeterince verimli bir şekilde yapılabilen sınırlı alanlar üzerinde tek değişkenli polinomları çarpanlara ayırmakla güçlü bir şekilde ilişkili gibi görünmektedir . ...
Andrew Morgan

Benim tahminim JCF'nin verimli bir şekilde bulunabileceğidir, ama tam olarak emin değilim. Alanı genişletmek zorunda olduğunuzdan bahsediyorsunuz - bu küçültme için gereklidir (sabit boyutlu sonlu alanlar üzerindeki ayrık kütük kolaydır), bu yüzden sürpriz olmamalıdır. derecesi hakkında endişeleniyorum - özdeğerler sadece derecesine sahipken, ve değerleri özdeğerlerin güçlerinin doğrusal kombinasyonlarıdır, bu yüzden boyutunda bir alanda yaşamak gerekebilir. Cevabımda bu olası tuzakları not edeceğim, ancak hala etrafta kalmak için yeterli fikirlere katkıda bulunduğunu düşünüyorum. nxiyin!
Andrew Morgan

Sağ. Elemanlar, n derecesinin keyfi bir polinomu olabilen matrisin karakteristik polinomunun yarma alanında yaşarlar, bu nedenle yarılma alanı kabaca kadar yüksek olabilir (eğer hesaplamam doğru). Ama belki bu bir şekilde atlatılabilir. Karakter poli'sini çarpanlarına ayıralım (farklı derece çarpanlarına ayırma bile yeterli olmalıdır). Bir şekilde her faktörün köklerine karşılık gelen öz uzamayı tanımlayabilir miyiz? Yani, tam JCF yerine, orijinal alanın üzerinde her bloğun öz değerlerine sahip olacağı bir blok çapraz matris elde ederiz ...exp(nlogn)
Emil Jeřábek

... en fazla derecesinin uzantısında . Sonra her bloğu ayrı ayrı işleyebiliriz. (Bu sadece belirsiz bir fikir, bunu n
çözmeye
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.