next up previous
Next: About this document Up: My Home Page

Among n people, a celebrity is defined as someone who is known by everyone but does not know anyone. We seek to identify the celebrity (if one is present) by asking questions of the form ``Hey, x, do you know person y?''. Show how to find the celebrity using O(n) questions. Note that there are tex2html_wrap_inline255 possible questions to ask, so we cannot ask them all.

What if we ask 1 if she knows 2, and 2 if she knows 1? If both know each other neither can be a celebrity. If neither know each other, neither can be a celebrity. If one of them knows the other, the former cannot be a celebrity.

Thus in two questions we can eliminate at least one person from celebrity status. Thus in 2(n-1) questions, we have only one possible celebrity. It is now possible to check whether the survivor is really a celebrity using n-1 additional queries, by checking whether everyone else knows them.

An Eulerian cycle in a graph visits each edge exactly once. A graph contains an Eulerian cycle iff it is connected and the degree of each vertex is even. Give an O(|E|) algorithm to find an Eulerian cycle if one exists. Observe that an cycle of edges defines a graph where each vertex is of degree 2. Thus deleting a cycle from an Eulerian graph leaves each vertex with even degree, although the graph may not be connected.

We can use depth-first search to decompose the edges of a graph into cycles. If the graph was connected, these cycles must link together. Splicing them together gives an Eulerian cycle. For example, the cycle (1,2,3,1) and (4,5,6,1,4) can be spliced together as (4,5,6,1,2,3,1,4).

Although Eulerian cycle has an efficient algorithm, the Hamiltonian cycle problem (visit each vertex exactly once) is NP-complete.

lec21-7.au The Theory of NP-Completeness

Several times this semester we have encountered problems for which we couldn't find efficient algorithms, such as the traveling salesman problem. We also couldn't prove an exponential time lower bound for the problem.  

By the early 1970s, literally hundreds of problems were stuck in this limbo. The theory of NP-Compleness, developed by Stephen Cook and Richard Karp, provided the tools to show that all of these problems were really the same problem. lec21-8.au Polynomial vs. Exponential Time   

tabular14

lec21-9.au The Main Idea

Suppose I gave you the following algorithm to solve the bandersnatch problem:   

Bandersnatch(G)
\> Convert G to an instance of the Bo-billy problem Y.
\> Call the subroutine Bo-billy on Y to solve this instance.
\> Return the answer of Bo-billy(Y) as the answer to G.

Such a translation from instances of one type of problem to instances of another type such that answers are preserved is called a reduction.  

Now suppose my reduction translates G to Y in O(P(n)):

  1. If my Bo-billy subroutine ran in O(P'(n)) I can solve the Bandersnatch problem in O(P(n)+P'(n))
  2. If I know that tex2html_wrap_inline333 is a lower-bound to compute Bandersnatch, then tex2html_wrap_inline335 must be a lower-bound to compute Bo-billy.

The second argument is the idea we use to prove problems hard! lec21-10.au

Convex Hull and Sorting

A nice example of a reduction goes from sorting numbers to the convex hull problem:   

convex-hull2.0in

We must translate each number to a point. We can map x to tex2html_wrap_inline339 .

parabola2.5in

Why? That means each integer is mapped to a point on the parabola tex2html_wrap_inline341 . lec21-11.au

Since this parabola is convex, every point is on the convex hull. Further since neighboring points on the convex hull have neighboring x values, the convex hull returns the points sorted by x-coordinate, ie. the original numbers.

Sort(S)
\> For each tex2html_wrap_inline349 , create point tex2html_wrap_inline351 .
\> Call subroutine convex-hull on this point set.
\> From the leftmost point in the hull,
\>\> read off the points from left to right.

Creating and reading off the points takes O(n) time.

What does this mean? Recall the sorting lower bound of tex2html_wrap_inline355 . If we could do convex hull in better than tex2html_wrap_inline357 , we could sort faster than tex2html_wrap_inline359 - which violates our lower bound.

Thus convex hull must take tex2html_wrap_inline361 as well!!!

Observe that any tex2html_wrap_inline363 convex hull algorithm also gives us a complicated but correct tex2html_wrap_inline365 sorting algorithm as well.

lec22-2.au What is a problem?

A problem is a general question, with parameters for the input and conditions on what is a satisfactory answer or solution.   

An instance is a problem with the input parameters specified.

Example: The Traveling Salesman

Problem: Given a weighted graph G, what tour tex2html_wrap_inline369 minimizes tex2html_wrap_inline371 .

Instance: tex2html_wrap_inline373 , tex2html_wrap_inline375 , tex2html_wrap_inline377 , tex2html_wrap_inline379 , tex2html_wrap_inline381 , tex2html_wrap_inline383

f673in

Solution: tex2html_wrap_inline385 cost= 27

A problem with answers restricted to yes and no is called a decision problem. Most interesting optimization problems can be phrased as decision problems which capture the essence of the computation.  

lec22-3.au Example: The Traveling Salesman Decision Problem.  

Given a weighted graph G and integer k, does there exist a traveling salesman tour with cost tex2html_wrap_inline397 k?

Using binary search and the decision version of the problem we can find the optimal TSP solution.

For convenience, from now on we will talk only about decision problems.

Note that there are many possible ways to encode the input graph: adjacency matrices, edge lists, etc. All reasonable encodings will be within polynomial size of each other.

The fact that we can ignore minor differences in encoding is important. We are concerned with the difference between algorithms which are polynomial and exponential in the size of the input. lec22-4.au Satisfiability

Consider the following logic problem:   

Instance: A set V of variables and a set of clauses C over V.

Question: Does there exist a satisfying truth assignment for C?

Example 1: tex2html_wrap_inline409 and tex2html_wrap_inline411

A clause is satisfied when at least one literal in it is TRUE. C is satisfied when tex2html_wrap_inline415 TRUE.

Example 2: tex2html_wrap_inline417 ,

Although you try, and you try, and you try and you try, you can get no satisfaction.  

There is no satisfying assigment since tex2html_wrap_inline419 must be FALSE (third clause), so tex2html_wrap_inline421 must be FALSE (second clause), but then the first clause is unsatisfiable!

For various reasons, it is known that satisfiability is a hard problem. Every top-notch algorithm expert in the world (and countless other, lesser lights) have tried to come up with a fast algorithm to test whether a given set of clauses is satisfiable, but all have failed. lec22-5.au Further, many strange and impossible-to-believe things have been shown to be true if someone in fact did find a fast satisfiability algorithm.

Clearly, Satisfiability is in NP, since we can guess an assignment of TRUE, FALSE to the literals and check it in polynomial time. lec22-10.au P versus NP

The precise distinction between whether a problem is in P or NP is somewhat technical, requiring formal language theory and Turing machines to state correctly.  

However, intuitively a problem is in P, (ie. polynomial) if it can be solved in time polynomial in the size of the input.

A problem is in NP if, given the answer, it is possible to verify that the answer is correct within time polynomial in the size of the input.  

Example P - Is there a path from s to t in G of length less than k.

Example NP - Is there a TSP tour in G of length less than k. Given the tour, it is easy to add up the costs and convince me it is correct.

Example not NP - How many TSP tours are there in G of length less than k. Since there can be an exponential number of them, we cannot count them all in polynomial time.

Don't let this issue confuse you - the important idea here is of reductions as a way of proving hardness. lec22-7.au 3-Satisfiability

Instance: A collection of clause C where each clause contains exactly 3 literals, boolean variable v.  

Question: Is there a truth assignment to v so that each clause is satisfied?

Note that this is a more restricted problem than SAT. If 3-SAT is NP-complete, it implies SAT is NP-complete but not visa-versa, perhaps long clauses are what makes SAT difficult?!

After all, 1-Sat is trivial!

Theorem: 3-SAT is NP-Complete

Proof: 3-SAT is NP - given an assignment, just check that each clause is covered. To prove it is complete, a reduction from tex2html_wrap_inline467 must be provided. We will transform each clause independantly based on its length.

Suppose the clause tex2html_wrap_inline469 contains k literals.

If none of the original variables in a clause are TRUE, there is no way to satisfy all of them using the additional variable:

But if any literal is TRUE, we have n-3 free variables and n-3 remaining 3-clauses, so we can satisfy each of them. tex2html_wrap_inline523

Since any SAT solution will also satisfy the 3-SAT instance and any 3-SAT solution sets variables giving a SAT solution - the problems are equivallent. If there were n clauses and m total literals in the SAT instance, this transform takes O(m) time, so SAT and 3-SAT.

Note that a slight modification to this construction would prove 4-SAT, or 5-SAT,... also NP-complete. However, it breaks down when we try to use it for 2-SAT, since there is no way to stuff anything into the chain of clauses. It turns out that resolution gives a polynomial time algorithm for 2-SAT.

lec22-9.au Having at least 3-literals per clause is what makes the problem difficult. Now that we have shown 3-SAT is NP-complete, we may use it for further reductions. Since the set of 3-SAT instances is smaller and more regular than the SAT instances, it will be easier to use 3-SAT for future reductions. Remember the direction to reduction!

lec24-2.au A Perpetual Point of Confusion

Note carefully the direction of the reduction.  

We must transform every instance of a known NP-complete problem to an instance of the problem we are interested in. If we do the reduction the other way, all we get is a slow way to solve x, by using a subroutine which probably will take exponential time.

This always is confusing at first - it seems bass-ackwards. Make sure you understand the direction of reduction now - and think back to this when you get confused. lec24-3.au Integer Programming

Instance: A set v of integer variables, a set of inequalities over these variables, a function f(v) to maximize, and integer B.  

Question: Does there exist an assignment of integers to v such that all inequalities are true and tex2html_wrap_inline565 ?

Example:

A solution to this is tex2html_wrap_inline567 , tex2html_wrap_inline569 .

Example:

Since the maximum value of f(v) given the constraints is tex2html_wrap_inline573 , there is no solution.

Theorem: Integer Programming is NP-Hard

Proof: By reduction from Satisfiability

Any set instance has boolean variables and clauses. Our Integer programming problem will have twice as many variables as the SAT instance, one for each variable and its compliment, as well as the following inequalities:

lec24-4.au For each variable tex2html_wrap_inline575 in the set problem, we will add the following constraints:

Our maximization function and bound are relatively unimportant: tex2html_wrap_inline589 B=0.

Clearly this reduction can be done in polynomial time.

lec24-5.au We must show:

  1. Any SAT solution gives a solution to the IP problem.

    In any SAT solution, a TRUE literal corresponds to a 1 in the IP, since if the expression is SATISFIED, at least one literal per clause in TRUE, so the sum in the inequality is tex2html_wrap_inline595 1.

  2. Any IP solution gives a SAT solution.

    Given a solution to this IP instance, all variables will be 0 or 1. Set the literals correspondly to 1 variable TRUE and the 0 to FALSE. No boolean variable and its complement will both be true, so it is a legal assignment with also must satisfy the clauses.

Neat, sweet, and NP-complete!

lec24-6.au

Things to Notice

  1. The reduction preserved the structure of the problem. Note that the reduction did not solve the problem - it just put it in a different format.
  2. The possible IP instances which result are a small subset of the possible IP instances, but since some of them are hard, the problem in general must be hard.
  3. The transformation captures the essence of why IP is hard - it has nothing to do with big coefficients or big ranges on variables; for restricting to 0/1 is enough. A careful study of what properties we do need for our reduction tells us a lot about the problem.
  4. It is not obvious that IP tex2html_wrap_inline607 NP, since the numbers assigned to the variables may be too large to write in polynomial time - don't be too hasty!




next up previous
Next: About this document Up: My Home Page

Steve Skiena
Tue Sep 15 17:28:55 EDT 1998