Haskell, 278 ตัวอักษร
(∈)=elem
r v[][]=[(>>=(++" ").show.fromEnum.(∈v))]
r v[]c@(a:b:_)=r(a:v)c[]++r(-a:v)c[]++[const"UNSOLVABLE"]
r v(a:b:c)d|a∈v||b∈v=r v c d|(-a)∈v=i b|(-b)∈v=i a|1<3=r v c(a:b:d)where i w|(-w)∈v=[]|1<3=r(w:v)(c++d)[]
t(n:_:c)=(r[][]c!!0)[1..n]++"\n"
main=interact$t.map read.words
ไม่ดุร้าย ทำงานในเวลาพหุนาม แก้ไขปัญหาที่ยาก (60 ตัวแปร 99 ข้อ) ได้อย่างรวดเร็ว:
> time (runhaskell 1933-2Sat.hs < 1933-hard2sat.txt)
1 1 1 0 0 0 0 0 0 1 1 0 0 1 0 1 1 1 0 1 1 0 0 1 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 1 0 0 0 0 1 0 1 1 1 1 0
real 0m0.593s
user 0m0.502s
sys 0m0.074s
และจริงๆแล้วเวลาส่วนใหญ่นั้นใช้ในการรวบรวมรหัส!
แฟ้มแหล่งที่มาเต็มกับกรณีทดสอบและการทดสอบอย่างรวดเร็วตรวจสอบที่มีอยู่
Ungolf'd:
-- | A variable or its negation
-- Note that applying unary negation (-) to a term inverts it.
type Term = Int
-- | A set of terms taken to be true.
-- Should only contain a variable or its negation, never both.
type TruthAssignment = [Term]
-- | Special value indicating that no consistent truth assignment is possible.
unsolvable :: TruthAssignment
unsolvable = [0]
-- | Clauses are a list of terms, taken in pairs.
-- Each pair is a disjunction (or), the list as a whole the conjuction (and)
-- of the pairs.
type Clauses = [Term]
-- | Test to see if a term is in an assignment
(∈) :: Term -> TruthAssignment -> Bool
a∈v = a `elem` v;
-- | Satisfy a set of clauses, from a starting assignment.
-- Returns a non-exhaustive list of possible assignments, followed by
-- unsolvable. If unsolvable is first, there is no possible assignment.
satisfy :: TruthAssignment -> Clauses -> [TruthAssignment]
satisfy v c@(a:b:_) = reduce (a:v) c ++ reduce (-a:v) c ++ [unsolvable]
-- pick a term from the first clause, either it or its negation must be true;
-- if neither produces a viable result, then the clauses are unsolvable
satisfy v [] = [v]
-- if there are no clauses, then the starting assignment is a solution!
-- | Reduce a set of clauses, given a starting assignment, then solve that
reduce :: TruthAssignment -> Clauses -> [TruthAssignment]
reduce v c = reduce' v c []
where
reduce' v (a:b:c) d
| a∈v || b∈v = reduce' v c d
-- if the clause is already satisfied, then just drop it
| (-a)∈v = imply b
| (-b)∈v = imply a
-- if either term is not true, the other term must be true
| otherwise = reduce' v c (a:b:d)
-- this clause is still undetermined, save it for later
where
imply w
| (-w)∈v = [] -- if w is also false, there is no possible solution
| otherwise = reduce (w:v) (c++d)
-- otherwise, set w true, and reduce again
reduce' v [] d = satisfy v d
-- once all caluses have been reduced, satisfy the remaining
-- | Format a solution. Terms not assigned are choosen to be false
format :: Int -> TruthAssignment -> String
format n v
| v == unsolvable = "UNSOLVABLE"
| otherwise = unwords . map (bit.(∈v)) $ [1..n]
where
bit False = "0"
bit True = "1"
main = interact $ run . map read . words
where
run (n:_:c) = (format n $ head $ satisfy [] c) ++ "\n"
-- first number of input is number of variables
-- second number of input is number of claues, ignored
-- remaining numbers are the clauses, taken two at a time
ในเวอร์ชัน golf'd satisfyและformatถูกนำไปกลิ้งreduceแม้ว่าเพื่อหลีกเลี่ยงการส่งผ่านnให้reduceส่งคืนฟังก์ชันจากรายการของตัวแปร ( [1..n]) ไปยังผลลัพธ์ของสตริง
- แก้ไข: (330 -> 323) สร้าง
sตัวดำเนินการจัดการกับ newline ได้ดีขึ้น
- แก้ไข: (323 -> 313) องค์ประกอบแรกจากรายการผลลัพธ์ที่ขี้เกียจมีขนาดเล็กกว่าตัวดำเนินการลัดวงจรแบบกำหนดเอง เปลี่ยนชื่อฟังก์ชั่นแก้ปัญหาหลักเพราะฉันชอบใช้
∮เป็นผู้ประกอบการ!
- แก้ไข: (313 -> 296) เก็บคำสั่งเป็นรายการเดียวไม่ใช่รายการ ประมวลผลสององค์ประกอบพร้อมกัน
- แก้ไข: (296 -> 291) ผสานทั้งสองฟังก์ชั่นวนซ้ำซึ่งกันและกัน มันราคาถูกกว่าแบบอินไลน์
★ดังนั้นตอนนี้เปลี่ยนชื่อทดสอบ∈
- แก้ไข: (291 -> 278) การจัดรูปแบบเอาต์พุตแบบอินไลน์เข้าสู่การสร้างผลลัพธ์