A fully polynomial-time approximation scheme (FPTAS) is an algorithm for finding approximate solutions to function problems, especially optimization problems. An FPTAS takes as input an instance of the problem and a parameter ε > 0. It returns as output a value which is at least 1 − ε {\displaystyle 1-\varepsilon } times the correct value, and at most 1 + ε {\displaystyle 1+\varepsilon } times the correct value. In the context of optimization problems, the correct value is understood to be the value of the optimal solution, and it is often implied that an FPTAS should produce a valid solution (and not just the value of the solution). Returning a value and finding a solution with that value are equivalent assuming that the problem possesses self reducibility. Importantly, the run-time of an FPTAS is polynomial in the problem size and in 1/ε. This is in contrast to a general polynomial-time approximation scheme (PTAS). The run-time of a general PTAS is polynomial in the problem size for each specific ε, but might be exponential in 1/ε. The term FPTAS may also be used to refer to the class of problems that have an FPTAS. FPTAS is a subset of PTAS, and unless P = NP, it is a strict subset.
Relation to other complexity classes All problems in FPTAS are fixed-parameter tractable with respect to the standard parameterization. Any strongly NP-hard optimization problem with a polynomially bounded objective function cannot have an FPTAS unless P=NP. However, the converse fails: e.g. if P does not equal NP, knapsack with two constraints is not strongly NP-hard, but has no FPTAS even when the optimal objective is polynomially bounded.
Converting a dynamic program to an FPTAS Woeginger presented a general scheme for converting a certain class of dynamic programs to an FPTAS.
Input The scheme handles optimization problems in which the input is defined as follows:
The input is made of n vectors, x1,...,xn. Each input vector is made of some a {\displaystyle a} non-negative integers, where a {\displaystyle a} may depend on the input. All components of the input vectors are encoded in binary. So the size of the problem is O(n+log(X)), where X is the sum of all components in all vectors.
Extremely-simple dynamic program It is assumed that the problem has a dynamic-programming (DP) algorithm using states. Each state is a vector made of some b {\displaystyle b} non-negative integers, where b {\displaystyle b} is independent of the input. The DP works in n steps. At each step i, it processes the input xi, and constructs a set of states Si. Each state encodes a partial solution to the problem, using inputs x1,...,xi. The components of the DP are:
A set S0 of initial states. A set F of transition functions. Each function f in F maps a pair (state,input) to a new state. An objective function g, mapping a state to its value. The algorithm of the DP is:
Let S0 := the set of initial states. For k = 1 to n do: Let Sk := {f(s,xk) | f in F, s in Sk−1} Output min/max {g(s) | s in Sn}. The run-time of the DP is linear in the number of possible states. In general, this number can be exponential in the size of the input problem: it can be in O(n Vb), where V is the largest integer than can appear in a state. If V is in O(X), then the run-time is in O(n Xb), which is only pseudo-polynomial time, since it is exponential in the problem size which is in O(log X). The way to make it polynomial is to trim the state-space: instead of keeping all possible states in each step, keep only a subset of the states; remove states that are "sufficiently close" to other states. Under certain conditions, this trimming can be done in a way that does not change the objective value by too much. To formalize this, we assume that the problem at hand has a non-negative integer vector d = (d1,...,db), called the degree vector of the problem. For every real number r>1, we say that two state-vectors s1,s2 are (d,r)-close if, for each coordinate j in 1,...,b: r − d j ⋅ s 1 , j ≤ s 2 , j ≤ r d j ⋅ s 1 , j {\displaystyle r^{-d_{j}}\cdot s_{1,j}\leq s_{2,j}\leq r^{d_{j}}\cdot s_{1,j}} (in particular, if dj=0 for some j, then s 1 , j = s 2 , j {\displaystyle s_{1,j}=s_{2,j}} ). A problem is called extremely-benevolent if it satisfies the following three conditions:
… excerpt ends here. Continue reading the full article.
