มีอัลกอริธึมเวลาพหุนามเพื่อพิจารณาว่าช่วงของชุดเมทริกซ์มีเมทริกซ์การเปลี่ยนแปลงหรือไม่?


30

ฉันต้องการค้นหาอัลกอริทึมเวลาพหุนามที่กำหนดว่าช่วงของชุดเมทริกซ์ที่กำหนดมีเมทริกซ์การเปลี่ยนแปลง

หากมีใครรู้ว่าปัญหานี้มีระดับความซับซ้อนที่แตกต่างกันนั่นจะเป็นประโยชน์


แก้ไข: ฉันติดแท็กคำถามนี้ด้วยการเขียนโปรแกรมเชิงเส้นเพราะฉันมีความสงสัยอย่างมากว่าหากวิธีการแก้ปัญหาดังกล่าวมีอยู่ก็จะเป็นขั้นตอนวิธีการเขียนโปรแกรมเชิงเส้น เหตุผลที่ฉันเชื่อว่าเป็นเพราะจุดที่สูงที่สุดของBirkhoff polytopeเป็นเมทริกซ์การเปลี่ยนรูปที่แม่นยำ หากคุณสามารถหาฟังก์ชั่นวัตถุประสงค์ที่มีการขยายหรือย่อให้เล็กสุดที่จุดยอดของ Birkhoff polytope คุณสามารถ จำกัด การทำงานของคุณกับการแยกของ polytope และพื้นที่ย่อยเวกเตอร์ของคุณจากนั้นขยายให้ใหญ่ที่สุดในเวลาพหุนาม ถ้าค่านี้เป็นเมทริกซ์การเปลี่ยนแปลงคุณจะรู้ว่าเซตนั้นมีการเปลี่ยนแปลง นี่คือความคิดของฉันในเรื่องนี้


แก้ไข 2: หลังจากคิดเพิ่มอีกดูเหมือนว่าการเปลี่ยนรูปเป็นองค์ประกอบของ Birkhoff Polytope กับ Euclidean norm , เราพิจารณา Birkhoff polytope เป็นตัวนูนของเมทริกซ์การเปลี่ยนแปลง บางทีนั่นอาจจะสำคัญ n×nnn×n


แก้ไข 3: ฉันเพิ่มแท็กการเขียนโปรแกรม semidefinite เพราะหลังจากความคิดเห็นก่อนหน้านี้ฉันเริ่มคิดว่าโซลูชันการเขียนโปรแกรม semidefinite อาจเป็นไปได้


2
เมทริกซ์อินพุตมีชนิดของรายการใด?

รายการอาจอยู่ในสาขาใดมีอิสระในการตั้งค่าเมทริกซ์; อย่างไรก็ตามคุณต้องการฟิลด์ที่มีขนาดใหญ่พอ (ดังนั้นฟิลด์ของคุณลักษณะ 2 จะไม่ดีสำหรับตัวอย่าง)
Nick

สามารถอธิบายว่าชุดของเมทริกซ์คืออะไร?
Mohammad Al-Turkistany

โมฮัมหมัด: ฉันคิดว่ามันเป็นการผสมผสานเชิงเส้นของเซตเมทริกซ์
Vivek Bagaria

4
@ DavidRicherby ฉันคิดว่าความสับสนของโมฮัมหมัดมาจากความจริงที่ว่าเรามักคิดว่าเมทริกซ์เป็นตัวแทนของแผนที่เชิงเส้นและบางครั้งช่วงของแผนที่เชิงเส้นบางครั้งก็ใช้เป็นอีกคำหนึ่งสำหรับช่วงของมัน แต่นั่นไม่เข้าท่าเลยดังนั้นฉันคิดว่าเราจะคิดว่าเมทริกซ์เป็นองค์ประกอบของปริภูมิเวกเตอร์
Sasho Nikolov

คำตอบ:


5

ทฤษฎีบท. ปัญหาในการโพสต์คือ NP-hard โดยการลดลงจากเซตย่อยรวม

แน่นอนว่าปัญหาดังกล่าวไม่น่าจะมีอัลกอริทึมโพลีเวลาตามที่ร้องขอโดย op


นี่คือสัญชาตญาณ ปัญหาในการโพสต์คือ

  • มีเมทริกซ์การเปลี่ยนแปลงในช่วงของชุดเมทริกซ์ที่กำหนดหรือไม่?

นี่คือหลักเช่นเดียวกับ

  • มีเมทริกซ์การเปลี่ยนแปลงที่ (คิดว่าเมทริกซ์เป็นเวกเตอร์) เป็นไปตามข้อ จำกัด เชิงเส้นที่กำหนดหรือไม่?

นี่คือสิ่งเดียวกัน

  • มีการจับคู่ที่สมบูรณ์แบบ (ในกราฟ bipartite ที่สมบูรณ์) ซึ่งการเกิดเวกเตอร์เป็นไปตามข้อ จำกัด เชิงเส้นที่กำหนดหรือไม่?

การลดชุดย่อย - ผลรวมของปัญหาหลังคือแบบฝึกหัดมาตรฐาน

นี่คือหลักฐานรายละเอียด


กำหนดปัญหากลางต่อไปนี้:

จับคู่ผลรวม:

การป้อนข้อมูล: สมบูรณ์ฝ่ายกราฟกับที่ไม่ใช่เชิงลบน้ำหนักขอบจำนวนเต็มและที่ไม่ใช่เชิงลบจำนวนเต็มเป้าหมายTG=(U,V,E)T

เอาท์พุท: ไม่มีการจับคู่ที่สมบูรณ์แบบของน้ำหนักว่าT ?GT


แทรก 1 เซตย่อยรวมผลรวมของโพลีไทม์จะลดลงเป็น Matching-Sum

พิสูจน์สิ่งนี้ว่าเป็นแบบฝึกหัดการบ้านมาตรฐาน หลักฐานสิ้นสุดแล้ว

บทที่ 2 การ จับคู่แบบรวมโพลีไทม์ช่วยลดปัญหาในการโพสต์

บทพิสูจน์ของเลมม่า 2แก้ไขอินพุตจับคู่ผลรวม: กราฟสองฝ่ายที่สมบูรณ์มีค่าน้ำหนักจำนวนเต็มแบบไม่ลบค่าลบw : U × V N +และเป้าหมายT N +โดยที่U = { U 1 , ... , U n }และV = { โวลต์1 , ... , V n } สำหรับแต่ละฉันG=(U,V,E)w:U×VN+TN+U={u1,,un}V={โวลต์1,...,โวลต์n} , กำหนด M ( i j )ให้เป็นเมทริกซ์ใน R ( n + 1 ) × ( n + 1 )โดยที่ M ( i j ) i j = T , และ M ( ฉันj ) n + 1 , n + 1 = w ( uผม,J{1,2,...,n}M(ผมJ)R(n+1)×(n+1)MผมJ(ผมJ)=T , และรายการอื่น ๆ ทั้งหมดเป็นศูนย์ การลด outputs ชุดต่อไปของการฝึกอบรม: { M ( ฉันเจ) : ฉัน, เจ{ 1 , ... , n } } สิ่งนี้กำหนดการลดMn+1,n+1(ผมJ)=W(ยูผม,โวลต์J)

{M(ผมJ):ผม,J{1,...,n}}.

ข้อเรียกร้อง ช่วงของชุดเมทริกซ์นี้ประกอบด้วยเมทริกซ์เหล่านั้นเป็นไปตามข้อ จำกัด เชิงเส้นM h , n + 1 = M n + 1 , h = 0สำหรับh n ทั้งหมดและ ข้อ จำกัด เชิงเส้นMR(n+1)×(n+1)Mชั่วโมง,n+1=Mn+1,ชั่วโมง=0ชั่วโมงn

Σผม=1nΣJ=1nMผมJW(ยูผม,โวลต์J)=TMn+1,n+1.

(การพิสูจน์การอ้างสิทธิ์)โดยการตรวจสอบแต่ละเมทริกซ์ในชุดตอบสนองข้อ จำกัด เหล่านี้ดังนั้นการรวมกันเชิงเส้นของเมทริกซ์เหล่านั้นจะตรงกันข้ามถ้าM R ( n + 1 ) × ( n + 1 )เป็นไปตามข้อ จำกัด จากนั้นMเท่ากับการรวมกันเชิงเส้นM = n ฉัน= 1n j = 1 α ฉันj M ( ฉันjM(ผมJ)MR(n+1)×(n+1)Mของการฝึกอบรมที่ α ฉันJ = M ฉันJ / M ( ฉันเจ) ฉันJ = M ฉันJ / T หมายเหตุโดยเฉพาะอย่างยิ่งว่าโดยคำจำกัดความต่างๆและข้อ จำกัด เชิงเส้น M ' n + 1 , n + 1 = Σ ฉันเจ α ฉันเจ W ( U ฉัน , วีเจ ) = Σ ฉันเจเอ็มM'=Σผม=1nΣJ=1nαผมJM(ผมJ)αผมJ=MผมJ/MผมJ(ผมJ)=MผมJ/T สิ่งนี้พิสูจน์ข้อเรียกร้อง)

Mn+1,n+1'=ΣผมJαผมJW(ยูผม,โวลต์J)=ΣผมJMผมJW(ยูผม,โวลต์J)/T=(TMn+1,n+1)/T=Mn+1,n+1.

ตอนนี้เราแสดงให้เห็นว่าการลดลงนั้นถูกต้อง นั่นคือรูปแบบของกราฟให้มีน้ำหนักTการจับคู่และถ้าหากชุดของการฝึกอบรมครอบคลุมการเปลี่ยนแปลงเมทริกซ์GT

( เฉพาะในกรณีที่. ) ครั้งแรกคิดว่ากราฟที่กำหนดมีน้ำหนักTที่สมบูรณ์แบบการจับคู่M ' ให้M { 0 , 1 } ( n + 1 ) × ( n + 1 )เป็นเมทริกซ์การเปลี่ยนรูปแบบn × n ที่สอดคล้องกันโดยเพิ่มแถวและคอลัมน์เพิ่มเติมเช่นM n + 1 , n + 1 = 1และM h , n +GTM'M{0,1}(n+1)×(n+1)n×nMn+1,n+1=1สำหรับทุกชั่วโมงn จากนั้น Σ n ฉัน= 1 Σ n J = 1 M ฉันเจ W( U ฉัน ,วีเจ )เป็นน้ำหนักของ M ' , ที่อยู่,Tและ M n + 1 , n + 1 =1Mชั่วโมง,n+1=Mn+1,ชั่วโมง=0ชั่วโมงnΣผม=1nΣJ=1nMผมJW(ยูผม,โวลต์J)M'TMn+1,n+1=1ดังนั้นข้อ จำกัด เชิงเส้นในการระงับการเรียกร้องและช่วงของการตั้งค่าที่กำหนดของการฝึกอบรมที่มีการเปลี่ยนแปลงเมทริกซ์MM

( ถ้า. ) ตรงกันข้ามสมมติว่าช่วงที่มีการเปลี่ยนแปลงใด ๆ เมทริกซ์Mโดยเรียกร้องเพียงคนเดียวไม่ใช่ศูนย์รายการในแถวn + 1หรือคอลัมน์n + 1คือM n + 1 , n + 1ดังนั้น (เป็นMเป็นเมทริกซ์การเปลี่ยนแปลง) มันจะต้องเป็นที่M n + 1 , n + 1 = 1 ดังนั้นการลบแถวสุดท้ายและคอลัมน์ให้n × nเปลี่ยนแปลงเมทริกซ์ ให้M เป็นคู่ที่สมบูรณ์แบบที่สุดMn+1n+1Mn+1,n+1MMn+1,n+1=1n×nM'สอดคล้องกับที่ n × nเปลี่ยนแปลงเมทริกซ์ น้ำหนักของ M 'คือ Σ n ฉัน= 1 Σ n J = 1 M ฉันเจ W ( U ฉัน , วีเจ )ซึ่ง (โดยเรียกร้อง) เป็น T M n + 1 , n + 1 = T ดังนั้นกราฟที่ให้มีการจับคู่น้ำหนัก T , พิสูจน์เล็มม่า 2Gn×nM'Σผม=1nΣJ=1nMผมJW(ยูผม,โวลต์J)TMn+1,n+1=TT  

นี่คือหลักฐานการล่าช้าของเลมม่า 1:

หลักฐานของเล็มม่า 1.กำหนดอินสแตนซ์ย่อย - ผลรวม , การลดเอาท์พุตอินสแตนซ์จับคู่ผลรวม( G = ( U , V , E ) , T )โดยที่U = { u 1 , u 2 , , u 2 n } , V = { v 1 , v 2 ,(W,T)ยังไม่มีข้อความ+n×ยังไม่มีข้อความ+(G=(ยู,V,E),T)ยู={ยู1,ยู2,...,ยู2n}สำหรับแต่ละฉัน{ 1 , ... , n }ขอบ ( U ฉัน , วีฉัน )มีน้ำหนัก W ฉันและขอบที่เหลืออยู่ทั้งหมดมีศูนย์น้ำหนักV={โวลต์1,โวลต์2,...,โวลต์2n}ผม{1,...,n}(ยูผม,โวลต์ผม)Wผม

สำหรับการจับคู่ที่สมบูรณ์แบบกับน้ำหนักขอบรวมกับชุดS = { i : ( u i , v i ) M , i n }เป็นวิธีแก้ปัญหาสำหรับอินสแตนซ์ย่อย - Sum ที่กำหนด (เนื่องจากสิ่งเหล่านี้ไม่ใช่ ขอบน้ำหนักเป็นศูนย์ในM )TS={ผม:(ยูผม,โวลต์ผม)M,ผมn}M

ในทางกลับกันหากมีวิธีแก้ปัญหาสำหรับอินสแตนซ์ย่อย - ผลรวมพูดว่าด้วยi S w i = T , ชุดของขอบ{ ( u i , v i ) : i S }คือ จับคู่บางส่วนที่มีน้ำหนักTและจะขยายได้อย่างง่ายดายเพื่อน้ำหนักTการจับคู่ที่สมบูรณ์แบบโดยการเพิ่มตัวอย่างเช่นไปตั้ง (ศูนย์น้ำหนัก) ขอบ:S{1,...,n}ΣผมSWผม=T{(ยูผม,โวลต์ผม):ผมS}TT

{(ยูผม+n,โวลต์ผม+n):ผมS}ผม{1,...,n}S{(ยูผม,โวลต์ผม+n),(ยูผม+n,โวลต์ผม)}.

สิ่งนี้พิสูจน์เล็มม่า 1. ทฤษฎีบทต่อจากเล็มมาส 1 และ 2    


ป.ล. ตามคำตอบนี้ข้อ จำกัด ของการจับคู่ผลรวมกับอินสแตนซ์ที่มีขอบน้ำหนักแบบพหุนามมีอยู่ใน P แต่ฉันแน่ใจว่าการ จำกัด ปัญหาในโพสต์ถึงเมทริกซ์ด้วยโพลีโนเมียล ) รายการยังคง NP ยาก


2
ดูเหมือนว่าคุณจะใช้เปลือกนูนของเมทริกซ์แทนการขยาย ช่วงของเมทริกซ์ที่คุณอธิบายคือพื้นที่เต็มของเมทริกซ์ หรือว่าฉันขาดอะไรไป?
วาเนสซ่า

@ Squark คุณถูกต้อง - ฉันตีความ "span" ผิด ขอบคุณ ฉันแก้ไขหลักฐานเพื่อใช้คำจำกัดความที่ถูกต้องของ span (เหมือนกับการรวมกันเชิงเส้นของเมทริกซ์)
Neal Young

ดี! ฉันคิดว่ามันจะเป็นการดีกว่าที่จะคูณนิยามของด้วยw ( u i , v j )ดังนั้นเราจึงไม่ต้องหารด้วยบางสิ่งที่อาจเป็น 0 นอกจากนี้ดูเหมือนว่าข้อพิสูจน์สามารถลดความซับซ้อนลงได้โดยการรวมสองการลดลงโดยไม่มีปัญหาระดับกลาง M(ผมJ)W(ยูผม,โวลต์J)
วาเนสซ่า

จุดดีเกี่ยวกับการหารด้วยศูนย์ ฉันจะแก้ไขมัน ฉันจะแยกการลดสองอย่างออกจากกัน แต่สำหรับฉันแล้วมันเป็นวิธีที่ง่ายกว่ามาก
Neal Young

3

เกี่ยวกับปัญหาการคำนวณเส้นผ่านศูนย์กลางของ polytope ที่แสดงเป็นจุดตัดของ halfspaces ปัญหาคือ NP-hard โดยทั่วไปและ NP-hard เพื่อประมาณภายในปัจจัยคงที่ใด ๆ ดูที่กระดาษของ Briedenและการอ้างอิงในนั้น ฉันคิดว่าสำหรับโพลีโอพรีแบบสมมาตรส่วนกลาง, SDP จะให้ประมาณโดยที่mคือจำนวนของความไม่เท่าเทียมกันที่กำหนด polytope ฉันร่างภาพนี้ใต้เส้นO(logm)ม.

ในกรณีของคุณ Birkhoff polytope ไม่ได้เป็นสมมาตรจากส่วนกลาง แต่ทำงานร่วมกับตัวเรือนูนของพอเพียงPและ- Pเพื่อจุดประสงค์ของคุณ ฉันคิดว่า polytope "สมมาตร Birkhoff" นี้สามารถแสดงเป็นชุดของเมทริกซ์สี่เหลี่ยมจัตุรัสM ทั้งหมดที่ตอบสนอง:PPPM

ผม* * * *,J* * * *:ΣผมMผมJ* * * *=ΣJMผม* * * *J=

ผม,J:-1MผมJ1

-11

หากนี่เป็นการแสดงที่ถูกต้อง (ไม่แน่ใจ) คุณสามารถเพิ่มข้อ จำกัด ที่ จำกัด polytope นี้ไปยังพื้นที่ย่อยที่คุณกำหนด ไม่ยากที่จะปรับ SDP ด้านล่างบรรทัดให้ตรงกับการเป็นตัวแทนนี้ แต่ฉันเลือกที่จะไม่ผ่านมันเพื่อให้การจัดการสัญกรณ์

ฉันไม่แน่ใจว่าเส้นผ่านศูนย์กลางโดยประมาณสำหรับปัญหาของคุณ: อาจช่วยให้คุณตัดสินใจได้ว่าพื้นที่ย่อยที่ให้นั้นอยู่ใกล้กับเมทริกซ์การเปลี่ยนแปลงหรือห่างไกลจากทั้งหมด แต่ไม่ได้ทำการคำนวณ


P={x:-Ax}Aม.×n

α2=สูงสุดΣผม=1nโวลต์ผม22

ภายใต้:

1im:j=1nAijvj22bi2

vinPαP

αO(logm)diam(P)(vi)i=1nαxPαO(logm)nggix~i=gTvi

E x~22=α2
im:E |(Ax~)i|2bi2    E maxi=1m|(Ax~)i|biClogm.
Cm

xxPx221Clogmα12Clogmx~Px~212α

โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.