Description
This set of exam problems is designed to test your capabilities in learning and understanding new algorithms, implementing them and generalizing them. Make
sure that you give yourself a lot of time in digesting the videos that I am going
to upload. You may use Python, C, C++, or Java to do the implementation.
All the work must be your own; collaboration is prohibited. Please do not search
on the Internet and try to find any help — you are going to waste your time —
I design the exam so that it is googling-proof.
0. As we have learned in cpts350 this semester, NP is the class of problems that
can be solved by nondeterministic algorithms in polynomial time. In particular,
NP-complete problems are the hardest ones in NP. Currently, it is open whether
we have efficient solutions to those problems. That is, many researchers are still
trying to find (deterministic) polynomial time algorithms to solve those NPcomplete problems, or at least to find practically efficient algorithms for them.
Such efforts would lead to successfully cracking RSA, which is known in NP
(but we do not know whether cracking RSA, i.e., factorizing a large number,
is NP-complete). There is a well-known NP-complete problem, called linear
Diophantine (written LD), which seeks nonnegative integer solutions to a linear
constraint system Q over multiple variables x1, · · · , xk, for some k; e.g.,
2×1 + x2 + 15 = 0
3×1 − 4×2 > 18
x1 + 3×2 < 27
x1 > 15
(1)
The example LD instance shown in (1) has two nonnegative integer variables x1
and x2 and four constraints. Notice that the range of the variables is unbounded;
hence, you can not assume that they are in 32-bits unsigned ! Please stare at
the example for 10 minutes and see how you would design an algorithm to find
whether it has solutions and if it has, how you would come up with a solution.
Suddenly, you find that all the knowledge you learned in the course won’t give
you any ideas that are even remotely close to designing such an algorithm. You
need a ground breaking idea. (Do not even try to solve LD using equationsolving skills (called Newton elimination) you learned in high school; variables
in those high school equations are of real values and do not apply here.)
1
1. (10pts) Formally, an LD-instance Q is given as, for some k and m,
C11x1 + · · · + C1kxk + C1 #1 0
.
.
.
Cm1x1 + · · · + Cmkxk + Cm #m 0
(2)
where all the xj ’s are nonnegative integer variables (called unknowns), and the
following are all the parameters:
• all the coefficients Cij ’s, which are integers (positive, negative, zero), and
• all the numbers Ci
’s, which are nonnegative integers, and, finally,
• all the #i
’s, each of which is in {>, <, =}.
What are the the values of the m, the k and the parameters for the example
LD instance shown in (1) ?
2. (10pts) Any algorithm that solves the LD problem is to take an instance Q
in (2) as the algorithm’s input, and to return yes (along with a solution) if the
instance has a nonnegative integer solution in the unknowns x1, · · · , xk, and to
return no if otherwise. What is the size of the input? (Hint: the size of number
4.6 billions is roughly 32. why?)
3. (10pts) Tell me why the following is NOT an algorithm to solve the LD
problem:
input an instance Q in (2)
for each nonnegative integer tuple (v1, · · · , vk)
check whether (x1 = v1, · · · , xk = vk) satisfies
all the constraints in Q.
If yes, return yes along with solution (x1 = v1, · · · , xk = vk).
The remaining of the exam asks you to implement a brilliant idea in solving the
LD problem: representing each linear constraint (such as 2x1 − 3x2 + 4 = 0)
in an instance Q using a labeled graph (i.e., a finite automaton). Hence, since
there are m constraints in the instance Q, you need then implement a graph
composition algorithm to convert the m graphs into one final graph. Then, you
perform basic graph search on the final graph to finally solve the LD problem.
4. (60pts) Watch the take-home-final-exam-help videos, where I presented the
aforementioned brilliant ideas (which were invented decades ago by some of our
brightest ancestors in algorithms) and you take notes while watching (otherwise,
you can not do problems 6 and 7). Make sure that you fully understand the
ideas. In this project, you are going to implement those ideas to solve an instance
Q in three variables and with 2 equations. Being Diophantine, the variables are
2
of nonnegative integers. As I said, you may use Python, C, C++, or Java to do
the implementation. I do not ask you to design the algorithm, instead, I present
the algorithm’s design, and you need only implement the algorithm. You need
turn-in working code and prepare for a demo.
Preparation 1.
Herein, an equation is in the form of
C1x1 + C2x2 + C3x3 + C = 0 (3)
where constants C1, C2, C3 are integers (positive,negative, zero), and constant
C is nonnegative. Hence, for instance, 3x1 + 0x2 − 4x3 − 17 = 0 is not an
equation in the above form. However, equivalently, −3x1 − 0x2 + 4x3 + 17 = 0
is in the above form.
For the equation in (3), we define
Cmax = max
d,d1,d2,d3∈{0,1}
|C1d1 + C2d2 + C3d3 + d|. (4)
(Exercise: For the aforementioned equation −3x1 − 0x2 + 4x3 + 17 = 0, what is
the value of Cmax?)
Implement a function that returns Cmax from the description of an equation in
(3), and use the above Exercise to test your function.
Preparation 2.
For a constant C ≥ 0, we use binary representation for C. For instance, if
C = 34, then in binary, C = 100010. In this case, we use
b6b5b4b3b2b1
for it, with b6 = 1, b5 = b4 = b3 = 0, b2 = 1, b1 = 0. Herein, we define KC = 6
(the number of bits needed to represent C). In particular when C = 0, we let
KC = 0.
Hence, for 1 ≤ i ≤ KC , the bi
is defined as above. However, for i = KC + 1,
the bi
is now defined as 0. (Exercise: Let C = 18. What is the value of KC ?
what is the value of b6? what is the value of b4?)
Implement a function that returns the value KC from a given constant C ≥ 0.
Implement a function that returns the value bi from a given constant C ≥ 0
and i, noticing that the i shall be in the range of 1 ≤ i ≤ KC + 1.
Algorithm 1. Equation to labeled graph (automaton)
We now construct a finite automaton M from the description of the equation
in (3).
The input alphabet of M contains exactly eight input symbols, where each
symbol is in the form of (a1, a2, a3) with a1, a2, a3 ∈ {0, 1}. (Hence, an input
symbol is a triple of three Boolean values.)
3
A state in M is a pair of values: [carry, i], where
−Cmax ≤ carry ≤ Cmax,
recalling that Cmax is defined in (4), and 1 ≤ i ≤ KC + 1.
The initial state in M is [carry = 0, i = 1]. The accepting state is [carry =
0, i = KC + 1].
For all states [carry, i] and [carry0
, i0
] and all input symbols (a1, a2, a3), the
following is true:
[carry, i]
(a1,a2,a3)
−→ [carry0
, i0
] is a transition in M (i.e., M moves from state
[carry, i] to state [carry0
, i0
] on reading input symbol (a1, a2, a3)) if and only if
the following is true:
Let R = C1a1 + C2a2 + C3a3 + bi + carry. Then, R is divisible by
2, and carry0 =
R
2
. Furthermore, if 1 ≤ i ≤ KC , then i
0 = i + 1,
else i
0 = i. (You shall use the functions implemented above to
implement this step; e.g., to compute the bi
.)
Implement a function that returns the finite automaton M from the description of an equation in (3). Notice that you shall use a graph as the data structure
of the automaton M. (The automata we constructed are all deterministic by
definition.)
Algorithm 2. Cartesian product of two labeled graphs (i.e., two automata)
Now we are given a system of 2 equations over three variables x1, x2, x3,
where we use E1(x1, x2, x3) and E2(x1, x2, x3) to indicate the two equations
and recall that the equations are in the form of (3). Suppose that, using the
Algorithm 1 presented earlier, we obtain finite automata M1 and M2 from the
two equations respectively. Then, we need construct a finite automaton M that
is the Cartesian product of M1 and M2. You need watch the aforementioned
take-home-final-exam-help videos on how to implement the Cartesian product
algorithm.
Implement a function that returns the M from the M1 and the M2 (which
are computed using the function implemented in Algorithm 1 from descriptions
of the two equations). Again, the resulting automaton M uses a graph as its
data structure.
Final step.
If there is no path from the initial to the accepting in M, then the equation
system does not have solutions, else we can find a solution by using DFS on the
graph of the M (from the initial to the accepting while collecting the sequence
of input symbols on the path). Notice that you shall reverse the sequence and
then convert it into digit before you output the solution. I will sketch a little
more detail of this step. Suppose that the following sequence of input symbols
is collected on a path from the initial to the accepting:
(1, 0, 1),(0, 1, 1),(1, 0, 0),(1, 1, 0)
4
then, what is the solution to the equation system? it is 1011, 0101, 1100 (how
did I do it? I got 1011 by picking the first bit from each input symbol in the
sequence!). After reversing them, I have 1101, 1010, 0011. Then I convert them
into digits: 13, 10, 3. Hence, the solution is actually x1 = 13, x2 = 10, x3 = 3.
Please also implement this final step and use the following two to test your
code:
T1. The following equation system does not have nonnegative integer solutions:
3x1 − 2x2 + x3 + 5 = 0
6x1 − 4x2 + 2x3 + 9 = 0
T2. The following equation system does have nonnegative integer solutions (and
manually verify that the solution found after you run your code indeed satisfies
the following constraints):
3x1 − 2x2 − x3 + 3 = 0
6x1 − 4x2 + x3 + 3 = 0
You shall turn in code and screenshots of running results on the two test
cases. You may also be requested for a demo by running your code on a new
test case provided by your TAs.
5. (10pts) Analyze the worst case time complexity for the algorithm in 4. Recall
that your complexity function is on the size of input.
6. (10pts) Sketch how to generalize the algorithm that you implemented (for
two constraints and three variables) to a general LD instance Q in (2).
7. (10pts) Argue why the generalized algorithm runs in exponential time (on
the size of input).
8. (40pts) Describe an application that your generalized algorithm in 6 above
can be used.
Grading criteria: Probs 1, 2, 3: correctness. Prob 4: if your code is perfectly
working, you get full score. Otherwise, partial credit may be given. Prob 5: the
lower (but still correct) the better. Prob 6: correctness and readability. I dont
need a detaile write-up. Prob 7: write up your reasoning. Prob 8: at least two
pages and will be graded on relevance and quality of writing.
I will not accept any hand-written solutions! Please type.
5