Compare commits
11 Commits
| Author | SHA1 | Date | |
|---|---|---|---|
| 35e00b921b | |||
| ff1db55ca0 | |||
| 9aedda643f | |||
| a276cd658f | |||
| 763e86f307 | |||
| 2f79036d1f | |||
| 1d09b82a41 | |||
| 70672ae603 | |||
| 4f7a6f2383 | |||
| 8edb8f8a45 | |||
| 99ef773ecb |
2
.gitignore
vendored
2
.gitignore
vendored
@@ -1,7 +1,7 @@
|
||||
|
||||
.vscode/
|
||||
*.sage.py
|
||||
./*.pdf
|
||||
*.pdf
|
||||
## Core latex/pdflatex auxiliary files:
|
||||
*.aux
|
||||
*.lof
|
||||
|
||||
14
algo.sty
Normal file
14
algo.sty
Normal file
@@ -0,0 +1,14 @@
|
||||
\def\begin@lg{\begin{minipage}{1in}\begin{tabbing}
|
||||
\quad\=\qquad\=\qquad\=\qquad\=\qquad\=\qquad\=\qquad\=\kill}
|
||||
\def\end@lg{\end{tabbing}\end{minipage}}
|
||||
|
||||
\newenvironment{algorithm}
|
||||
{\begin{tabular}{|l|}\hline\begin@lg}
|
||||
{\end@lg\\\hline\end{tabular}}
|
||||
|
||||
\newenvironment{algo}
|
||||
{\begin{center}\begin{algorithm}}
|
||||
{\end{algorithm}\end{center}}
|
||||
|
||||
\def\argmax{\operatornamewithlimits{arg\,max}}
|
||||
\def\argmin{\operatornamewithlimits{arg\,min}}
|
||||
2
chao.sty
2
chao.sty
@@ -44,8 +44,8 @@
|
||||
\RequirePackage{microtype}
|
||||
\RequirePackage[mathletters]{ucs} % allow Unicode in .tex file
|
||||
\RequirePackage[utf8]{inputenc}
|
||||
|
||||
\RequirePackage[charter]{mathdesign} % change fonts
|
||||
\RequirePackage{XCharter} % change fonts
|
||||
\RequirePackage{berasans, beramono}
|
||||
\RequirePackage{eucal}
|
||||
\RequirePackage[nocompress]{cite} % other convenient stuff
|
||||
|
||||
191
main.tex
191
main.tex
@@ -1,6 +1,6 @@
|
||||
\documentclass[12pt]{article}
|
||||
\usepackage{chao}
|
||||
|
||||
\usepackage{algo}
|
||||
|
||||
\geometry{a4paper,margin=2cm}
|
||||
% \usepackage[breakable, theorems, skins]{tcolorbox}
|
||||
@@ -23,7 +23,7 @@
|
||||
% \end{tcolorbox}
|
||||
% }
|
||||
|
||||
\title{connectivity interdiction}
|
||||
\title{Connectivity Interdiction Notes}
|
||||
\author{}
|
||||
\date{}
|
||||
|
||||
@@ -31,26 +31,20 @@
|
||||
\DeclareMathOperator*{\len}{len}
|
||||
|
||||
\begin{document}
|
||||
|
||||
\maketitle
|
||||
|
||||
\section{``Cut-free'' Proof}
|
||||
\begin{problem}[b-free knapsack]\label{bfreeknap}
|
||||
Consider a set of elements $E$ and two weights $w:E\to \Z_+$ and
|
||||
$c:E\to \Z_+$ and a budget $b\in \Z_+$. Given a feasible set $\mathcal
|
||||
F\subset 2^E$, find $\min_{X
|
||||
\in \mathcal F, F\subset E} w(X\setminus F)$ such that $c(F)\leq b$.
|
||||
Consider a set of elements $E$ with weights $w:E\to \Z_+$ and cost $c:E\to \Z_+$ and a budget $b\in \Z_+$. Given a feasible set $\mathcal F\subset 2^E$, find $\min_{X\in \mathcal F, F\subset E} w(X\setminus F)$ such that $c(F)\leq b$.
|
||||
\end{problem}
|
||||
Always remember that $\mathcal F$ is usually not explicitly given.
|
||||
Note that $\mathcal F$ is usually not explicitly given.
|
||||
|
||||
\begin{problem}[Normalized knapsack]\label{nknap}
|
||||
Given the same input as \autoref{bfreeknap}, find $\min \limits_{X
|
||||
\in \mathcal F, F\subset E} \frac{w(X\setminus F)}{B-c(F)}$ such that
|
||||
$c(F)\leq b$.
|
||||
Given the same input as \autoref{bfreeknap}, find $\min \limits_{X\in \mathcal F, F\subset E} \frac{w(X\setminus F)}{B-c(F)}$ such that $c(F)\leq b$.
|
||||
\end{problem}
|
||||
In \cite{vygen_fptas_2024} the normalized min-cut problem use $B=b+1$. Here we use any integer $B>b$ and see how their method works.
|
||||
|
||||
Denote by $\tau$ the optimum of \autoref{nknap}. Define a new weight
|
||||
$w_\tau:\E\to \R$,
|
||||
Let $\tau$ be the optimum of \autoref{nknap}. Define a new weight $w_\tau:\E\to \R$,
|
||||
|
||||
\[
|
||||
w_\tau(e)=\begin{cases}
|
||||
@@ -63,7 +57,7 @@ w_\tau(e)=\begin{cases}
|
||||
Let $(X^N,F^N)$ be the optimal solution to \autoref{nknap}.
|
||||
Every element in $F^N$ is heavy.
|
||||
\end{lemma}
|
||||
The proof is exactly the same as \cite[Lemma 1]{vygen_fptas_2024}.
|
||||
The proof is the same as \cite[Lemma 1]{vygen_fptas_2024}.
|
||||
|
||||
The following two lemmas show (a general version of) that the optimal cut $C^N$ to normalized min-cut is exactly the minimum cut under weights $w_\tau$.
|
||||
|
||||
@@ -85,9 +79,9 @@ The following two lemmas show (a general version of) that the optimal cut $C^N$
|
||||
Thus by \autoref{lem:lb}, $X^N$ gets the minimum.
|
||||
\end{proof}
|
||||
|
||||
Now we show the counter part of \cite[Theorem 5]{vygen_fptas_2024}, which states the optimal solution to \autoref{bfreeknap} is a $\alpha$-approximate solution to $\min_{F\in \mathcal{F}} w_\tau(F)$.
|
||||
% Now we show the counter part of \cite[Theorem 5]{vygen_fptas_2024}, which states the optimal solution to \autoref{bfreeknap} is a $\alpha$-approximate solution to $\min_{F\in \mathcal{F}} w_\tau(F)$.
|
||||
|
||||
\begin{lemma}[Lemma 4 in \cite{vygen_fptas_2024}]\label{lem:conditionalLB}
|
||||
\begin{lemma}\label{lem:conditionalLB}
|
||||
Let $(X^*,F^*)$ be the optimal solution to \autoref{bfreeknap}.
|
||||
$X^*$ is either an $\alpha$-approximate solution to $\min_{X\in\mathcal F}
|
||||
w_\tau(X)$ for some $\alpha>1$, or $w(X^*\setminus F^*)\geq
|
||||
@@ -96,7 +90,7 @@ Now we show the counter part of \cite[Theorem 5]{vygen_fptas_2024}, which states
|
||||
|
||||
% In fact, corollary 1 and theorem 5 are also the same as those in
|
||||
% \cite{vygen_fptas_2024}.
|
||||
Then following the argument of Corollary 1 in \cite{vygen_fptas_2024}, assume that $X^*$ is not an $\alpha$-approximate solution to $\min_{X\in\mathcal F}
|
||||
Then following arguments in \cite[Corollary 1]{vygen_fptas_2024}, assume that $X^*$ is not an $\alpha$-approximate solution to $\min_{X\in\mathcal F}
|
||||
w_\tau(X)$ for some $\alpha>1$. We have
|
||||
\[
|
||||
\frac{w(C^N\setminus F^N)}{w(C^*\setminus F^*)}\leq \frac{\tau(B-c(F^N))}{\tau(\alpha B-b)}\leq \frac{B}{\alpha B-b},
|
||||
@@ -104,16 +98,14 @@ Then following the argument of Corollary 1 in \cite{vygen_fptas_2024}, assume th
|
||||
where the second inequality uses \autoref{lem:conditionalLB}.
|
||||
One can see that if $\alpha>2$, $\frac{w(C^N\setminus F^N)}{w(C^*\setminus F^*)}\leq \frac{B}{\alpha B-b} <1$ which implies $(C^*,F^*)$ is not optimal. Thus for $\alpha >2$, $X^*$ must be a $2$-approximate solution to $\min_{X\in\mathcal F} w_\tau(X)$.
|
||||
|
||||
Finally we get a knapsack version of Theorem 4:
|
||||
\begin{theorem}[Theorem 4 in \cite{vygen_fptas_2024}]\label{thm:main}
|
||||
Finally we get a general version of \cite[Theorem 4]{vygen_fptas_2024}:
|
||||
\begin{theorem}\label{thm:main}
|
||||
Let $X^{\min}$ be the optimal solution to $\min_{X\in\mathcal F} w_\tau(X)$.
|
||||
The optimal set $X^*$ in \autoref{bfreeknap} is a
|
||||
2-approximation to $X^{\min}$.
|
||||
\end{theorem}
|
||||
|
||||
Thus to obtain a FPTAS for \autoref{bfreeknap}, one need to design a FPTAS for
|
||||
\autoref{nknap} and a polynomial time alg for finding all 2-approximations to
|
||||
$\min_{X\in\mathcal F} w_\tau(X)$.
|
||||
Thus to obtain a FPTAS for \autoref{bfreeknap}, one need to design a FPTAS for \autoref{nknap} and a polynomial time algorithm for finding all 2-approximations to $\min_{X\in\mathcal F} w_\tau(X)$.
|
||||
|
||||
\paragraph{FPTAS for \autoref{nknap} in \cite{vygen_fptas_2024}} (The name
|
||||
``FPTAS'' here is not precise since we do not have a approximation scheme but
|
||||
@@ -134,20 +126,21 @@ have a non-zero lowerbound on $\opt^N$.
|
||||
|
||||
\begin{conjecture}
|
||||
Let $(C,F)$ be the optimal solution to connectivity interdiction. The optimum
|
||||
cut $C$ can be computed in polynomial time. In other words, connectivity
|
||||
interdiction is almost as easy as knapsack.
|
||||
cut $C$ can be computed in polynomial time.
|
||||
\end{conjecture}
|
||||
|
||||
Note that there is a FPTAS algorithm for finding $C$ in \cite{vygen_fptas_2024}.
|
||||
|
||||
\section{Connections}
|
||||
For unit costs, connectivity interdiction with budget $b=k-1$ is the same
|
||||
For unit weight and cost, connectivity interdiction with budget $b=k-1$ is the same
|
||||
problem as finding the minimum weighted edge set whose removal breaks $k$-edge
|
||||
connectivity.
|
||||
|
||||
It turns out that \autoref{nknap} is just a necessary ingredient for MWU.
|
||||
Authors of \cite{vygen_fptas_2024} $\subset$ authors of
|
||||
\cite{chalermsook_approximating_2022}.
|
||||
\autoref{nknap} may come from an intermediate problem of MWU methods for positive covering LPs.
|
||||
% Authors of \cite{vygen_fptas_2024} $\subset$ authors of
|
||||
% \cite{chalermsook_approximating_2022}.
|
||||
|
||||
How to derive normalized min cut for connectivity interdiction?
|
||||
Can we get an FPTAS using LP methods?
|
||||
|
||||
|
||||
\begin{equation*}
|
||||
@@ -175,9 +168,9 @@ s.t.& & \sum_{e\in T} x_e&\geq 1 & &\forall T\quad \text{($x$ is a c
|
||||
|
||||
Now this LP looks similar to the normalized min-cut problem.
|
||||
|
||||
A further reformulation (the new $x$ is $x-y$) gives us the following,
|
||||
A further reformulation ($x\from x-y$) gives us the following integer program,
|
||||
|
||||
\begin{equation}\label{lp:conninterdict}
|
||||
\begin{equation}\label{ip:interdiction}
|
||||
\begin{aligned}
|
||||
\min& & \sum_{e} x_e w(e) & & \\
|
||||
s.t.& & \sum_{e\in T} x_e+y_e&\geq 1 & &\forall T\quad \text{($x+y$ is a cut)}\\
|
||||
@@ -189,23 +182,21 @@ s.t.& & \sum_{e\in T} x_e+y_e&\geq 1 & &\forall T\quad \text{($x+y$
|
||||
|
||||
Note that now this is almost a positive covering LP. Let $L(\lambda)= \min \{ w(C\setminus F)-\lambda(b-c(F)) | \forall \text{cut $C$}\;\forall F\subset C
|
||||
% \land c(F)\leq b
|
||||
\}$.
|
||||
Consider the Lagrangian dual,
|
||||
\begin{equation*}
|
||||
\}$
|
||||
and consider the Lagrangian dual,
|
||||
\[
|
||||
\max_{\lambda\geq 0} L(\lambda)= \max_{\lambda\geq 0} \min \left\{ w(C\setminus F)-\lambda(b-c(F)) | \forall \text{cut $C$}\;\forall F\subset C
|
||||
% \land c(F)\leq b
|
||||
\right\}
|
||||
\end{equation*}
|
||||
\right\}.
|
||||
\]
|
||||
|
||||
|
||||
% At this point, it becomes clear how the normalized min-cut is implicated in \cite{vygen_fptas_2024}. The optimum of normalized min-cut is exactly the value of $\lambda$ when $L(\lambda)$ is 0.
|
||||
We have shown that the budget $B$ in normalized min-cut does not really matter as long as $B>b$. Note that $L(\lambda)$ and the normalized min-cut look similar to the principal sequence of partitions of a graph and the graph strength problem.
|
||||
|
||||
\subsection{graph strength}
|
||||
% Assume that the graph $G$ is connected.
|
||||
Given a graph $G=(V,E)$ and a cost function $c:V\to \Z_+$, the strength $\sigma(G)$ is defined as $\sigma(G)=\min_{\Pi}\frac{c(\delta(\Pi))}{|\Pi|-1}$, where $\Pi$ is any partition of $V$, $|\Pi|$ is the number of parts in the partition and $\delta(\Pi)$ is the set of edges between parts. Note that an alternative formulation of strength (using graphic matroid rank) is $\sigma(G)=\min_{F\subset E} \frac{|E-F|}{r(E)-r(F)}$, which in general is the fractional optimum of matroid base packing.
|
||||
For a graph $G=(V,E)$ with edge capacity $c:V\to \Z_+$, the strength $\sigma(G)$ is defined as $\sigma(G)=\min_{\Pi}\frac{c(\delta(\Pi))}{|\Pi|-1}$, where $\Pi$ is any partition of $V$, $|\Pi|$ is the number of parts in the partition and $\delta(\Pi)$ is the set of edges between parts. Note that an alternative formulation of strength (using graphic matroid rank function) is $\sigma(G)=\min_{F\subset E} \frac{c(E-F)}{r(E)-r(F)}$, which in general is the fractional optimum of matroid base packing.
|
||||
|
||||
The principal sequence of partitions of $G$ is a piecewise linear concave curve $L(\lambda)=\min_\Pi c(\delta(\Pi))-\lambda |\Pi|$. ($L(\lambda)=\min_{F\in E}c(E\setminus F)-\lambda(r(E)-r(F)+1)$) Cunningham used principal partition to computed graph strength\cite{cunningham_optimal_1985}. There is a list of good properties mentioned in \cite[Section 6]{chekuri_lp_2020}(implicated stated in \cite{cunningham_optimal_1985}).
|
||||
The principal sequence of partitions of $G$ is a pwl concave curve $L(\lambda)= \min_\Pi c(\delta(\Pi))-\lambda |\Pi|$. (alternatively, $L(\lambda)=\min_{F\in E}c(E\setminus F)-\lambda(r(E)-r(F)+1)$) Cunningham used principal partition to computed graph strength\cite{cunningham_optimal_1985}. There is a list of good properties mentioned in \cite[Section 6]{chekuri_lp_2020}(implicated stated in \cite{cunningham_optimal_1985}).
|
||||
\begin{enumerate}
|
||||
\item We can assume $G$ is connected and deal with the smallest strength component. One can see this by fractional base packing on the direct sum of matroids. Note that on disconnected graphs we should use the edge set definition instead of partitions.
|
||||
\item $L(\lambda)$ is piecewise linear concave since it is the lower envelope of some line arrangement.
|
||||
@@ -227,6 +218,16 @@ Now we focus on $L(\lambda)=\min \{w(C\setminus F)-\lambda(b-c(F)) | \forall \te
|
||||
% Now we can see that $B$ in the denominator $B-c(F)$ should be the cost of mincut in $G$.
|
||||
Let $B$ be the minimum cost of cuts in $G$.
|
||||
We have $-\lambda(b-B)\leq w(C\setminus F)-\lambda(b-c(F))$ for any cut $C$ and $F\subsetneq C$. Thus the upperbound is $\e=\min \frac{w(C\setminus F)}{B-c(F)}$.
|
||||
Note that $B-c(F)$ is always positive since otherwise the curve is not concave.
|
||||
|
||||
Now we depict the curve. The first segment is determined by $C=F\in \argmin_{X\subset V} c(\delta(X))$. The first breakpoint is $\e=\min \frac{w(C\setminus F)}{B-c(F)}$ and for the second segment we have $(C,F)\in \argmin\frac{w(C\setminus F)}{B-c(F)}$ with minimum $c(F)$.
|
||||
Now there are two cases:
|
||||
\begin{enumerate}
|
||||
\item if $c(F)\leq b$, then the first breakpoint is the optimal $\lambda$.
|
||||
\item if $c(F)\in (b,B)$, we have to look at more segments.
|
||||
\end{enumerate}
|
||||
Unfortunately, the seconds case is possible (consider a path with parallel edges) and the number of segments can be exponential.
|
||||
|
||||
|
||||
% \note{It remains to show that the optimal solution at $\e$ guarantees $c(F)\leq b$? or maybe we don't need this for normalized mincut. I think normalized min-cut should not require $c(F)\leq b$. Further checks are needed. What we can guarantee is that $c(F)\leq B$.}
|
||||
|
||||
@@ -235,9 +236,9 @@ Consider $L(\lambda)$ for cut problem. One can see that the optimal $\lambda$ is
|
||||
|
||||
|
||||
\subsection{integrality gap}
|
||||
I guess the 2-approximate min-cut enumeration algorithm implies an integrality gap of 2 for cut interdiction problem.
|
||||
I guess the 2-approximate min-cut enumeration algorithm implies an integrality gap of 2 for cut interdiction problem. \textcolor{red}{Which is wrong}.
|
||||
|
||||
First consider the dual of linear relaxation of \autoref{lp:conninterdict}.
|
||||
First consider the dual of linear relaxation of \autoref{ip:interdiction}.
|
||||
|
||||
\begin{equation}\label{lp:dualcutint}
|
||||
\begin{aligned}
|
||||
@@ -247,14 +248,14 @@ s.t.& & \sum_{T\ni e} z_T &\leq w(e) & &\forall e\in E\\
|
||||
& & z_T,\lambda &\geq 0 & &
|
||||
\end{aligned}
|
||||
\end{equation}
|
||||
We want to prove something like tree packing theorem for \autoref{lp:dualcutint}.
|
||||
\begin{conjecture}
|
||||
The optimum of \autoref{lp:dualcutint} is $\min \set{\frac{w(C\setminus F)}{B-c(F)}| \forall \text{cut $C$}, c(F)\leq b}$, where $B$ is the cost of mincut in $G$ and $b$ is the budget.
|
||||
\end{conjecture}
|
||||
% We want to prove something like tree packing theorem for \autoref{lp:dualcutint}.
|
||||
% \begin{conjecture}
|
||||
% The optimum of \autoref{lp:dualcutint} is $\min \set{\frac{w(C\setminus F)}{B-c(F)}| \forall \text{cut $C$}, c(F)\leq b}$, where $B$ is the cost of mincut in $G$ and $b$ is the budget.
|
||||
% \end{conjecture}
|
||||
|
||||
I believe the previous conjecture is not likely to be true.
|
||||
% I believe the previous conjecture is not likely to be true.
|
||||
|
||||
\paragraph{Weight truncation} Assuming we know the optimal $\lambda$ to the LP dual, \autoref{lp:dualcutint} in fact gives the idea of weight truncation. The capacity of each edge $e$ in the ``tree packing'' is $\min\{c(e)\lambda,w(e)\}$. Therefore, the optimum of \autoref{lp:dualcutint} is $\Lambda_{w_\tau}^{fr}-b\lambda$, where $\Lambda_{w_\tau}^{fr}$ is the fractional mincut on $G$ with weights $w_\tau$.
|
||||
\paragraph{Weight truncation} Assuming that the optimal $\lambda$ to the LP dual is known, \autoref{lp:dualcutint} in fact gives the idea of weight truncation. The capacity of each edge $e$ in the ``tree packing'' is $\min\{c(e)\lambda,w(e)\}$. Therefore, the optimum of \autoref{lp:dualcutint} is $\Lambda_{w_\tau}^{fr}-b\lambda$, where $\Lambda_{w_\tau}^{fr}$ is the fractional mincut on $G$ with weights $w_\tau$.
|
||||
|
||||
\paragraph{The optimal $\lambda$} Denote by $\lambda^*$ the optimal $\lambda$ that maximizes $L(\lambda)$. From the previous argument on the first segment of $L(\lambda)$ we know that $\lambda^* \geq \min \frac{w(C\setminus F)}{B-c(F)}$. Now assume $\lambda^* > \min_{c(F)\leq b} \frac{w(C\setminus F)}{b-c(F)}$. We have $\min w(C\setminus F)-\lambda^*(b-c(F))<w(C\setminus F)-\min_{c(F)\leq b} \frac{w(C\setminus F)}{b-c(F)}(b-c(F))=0$ since the optimum must be achieved by $F$ such that $0\leq b-c(F)$(the slope). The negative optimum of $L(\lambda)$ contradicts the fact that $L(0)=0$ and $L$ is concave. Hence, the optimal solution $\lambda^*$ is in the range $[\min\frac{w(C\setminus F)}{B-c(F)},\min_{c(F)\leq b}\frac{w(C\setminus F)}{b-c(F)}]$.
|
||||
|
||||
@@ -265,7 +266,7 @@ $\lambda_i=\min \frac{w(C\setminus F)-w(C_{i-1}\setminus F_{i-1})}{c(F_{i-1})-c(
|
||||
The proof is using the argument for showing $\lambda_1=\min \frac{w(C\setminus F)}{B-c(F)}$ and induction. $\lambda_i$ looks similar to normalized mincut but is related to the slope and vertical intercept of a previous segment.
|
||||
|
||||
\begin{conjecture}\label{conj:intgap}
|
||||
\autoref{lp:conninterdict} has an integrality gap of 4.
|
||||
\autoref{ip:interdiction} has an integrality gap of 4.
|
||||
\end{conjecture}
|
||||
|
||||
However, \autoref{conj:intgap} is wrong. The integrality gap is unbounded. Consider a cycle $C_n$ of $n$ vertices with two special edges $e_1,e_2$. Let $L$ be a large number.
|
||||
@@ -286,9 +287,92 @@ For LP, we assign $x=0$ and $y_e=\frac{1}{n-2}$ for every edge except $e_1$. The
|
||||
|
||||
What is the gap if we only relax $\lambda$ in the Lagrangian dual?
|
||||
|
||||
\section{Random Stuff}
|
||||
\section{LP method}
|
||||
Consider the following IP for connectivity interdiction:
|
||||
|
||||
\subsection{remove box constraints}
|
||||
\begin{equation*}
|
||||
\begin{aligned}
|
||||
\min& & \sum_{e} x_e w(e) & & \\
|
||||
s.t.& & \sum_{e\in T} x_e+y_e&\geq 1 & &\forall T\quad \text{($x+y$ is a cut)}\\
|
||||
& & \sum_{e} y_e c(e) &\leq b & &\text{(budget for $F$)}\\
|
||||
% & & x_e&\geq y_e & &\forall e\quad(F\subset C)\\
|
||||
& & y_e,x_e&\in\{0,1\} & &\forall e
|
||||
\end{aligned}
|
||||
\end{equation*}
|
||||
|
||||
We write its Lagrangian dual.
|
||||
\[
|
||||
\max_{\lambda\geq 0} \min_{\text{cut $C$ and $F\subseteq C$}} w(C-F)-\lambda(b-c(F))
|
||||
\]
|
||||
|
||||
For fixed $\lambda$, we define $L(\lambda)=\min\limits_{\text{cut $C$ and $F\subseteq C$}} w(C)-w(F)+\lambda c(F)$. Note that the Lagrangian dual is $\max_\lambda L(\lambda) - \lambda b$.
|
||||
|
||||
If the optimal $C^*$ of $L(\lambda)$ is known, any edge in $C^*$ with $w(e)\geq c(e)\lambda$ will be added to $F^*$. So now the hard part is to find $C^*$. We reweight edges.
|
||||
|
||||
\[
|
||||
w_\lambda(e)=\begin{cases}
|
||||
w(e) & \text{if $w(e)< \lambda\cdot c(e)$ (light elem)}\\
|
||||
\lambda\cdot c(e) & \text{otherwise (heavy elem)}
|
||||
\end{cases}
|
||||
\]
|
||||
|
||||
Then clearly edges in $F^*$ are heavy and edges in $C^*-F^*$ are light.
|
||||
|
||||
Now consider the min cut under capacity $w_\lambda$. For any cut $C$ in $G$, its capacity $w_\lambda(C)$ can be write as the sum of heavy elements ($\lambda c(F)$) and light elements ($w(C-F)$) which is at least $L(\lambda)$.
|
||||
On the other hand, $C^*$ is a feasible cut and the capacity of $C^*$ is exactly $L(\lambda)$. Thus $C^*$ is a mincut in the reweighted graph and we can easily find it.
|
||||
|
||||
\subsection{Approximation}
|
||||
|
||||
How is $C^*$ related to the optimal solution to IP?
|
||||
|
||||
% The argument is almost identical to lyft optimization
|
||||
|
||||
One can see that the Lagrangian dual (denoted by LD) is at most the optimum of IP. So we have $\opt(LD) \leq \opt(IP)$.
|
||||
|
||||
We can assume WLOG that the optimal $\lambda^*$ in LD is the intersection of two lines with positive and negative slopes. Then there exists an optimal solution $(\lambda^*, C^{LD},F^{LD})$ to LD such that $c(F^{LD})\leq b$. Then we have
|
||||
\begin{equation}\label{eq:ub}
|
||||
L(\lambda^*)\geq \opt(LD)+\lambda^* b -\lambda c(F^{LD}) = w(C^{LD}-F^{LD}) \geq \opt(IP)=w(C^*-F^*),
|
||||
\end{equation}
|
||||
since $\opt(IP)$ is the smallest $b$-free min cut.
|
||||
|
||||
We have $L(\lambda^*)\leq w_{\lambda^*}(C^*)$ since $L(\lambda^*)$ is the value of the minimum cut in $(G,w_{\lambda^*})$. Now we prove $L(\lambda^*)+b\lambda \geq w_{\lambda^*}(C^*)$.
|
||||
|
||||
\begin{equation*}
|
||||
\begin{aligned}
|
||||
L(\lambda^*)+b\lambda^* &\geq \opt(IP)+\lambda^* c(F^*)\\ &=w(C^*-F^*)+\lambda^* c(F^*)\\ &\geq w_{\lambda^*}(C^*)
|
||||
\end{aligned}
|
||||
\end{equation*}
|
||||
The first line uses \autoref{eq:ub} and the fact that $c(F^*)\leq b$. The last line follows from the definition of $w_{\lambda}$.
|
||||
Finally the approximation guarantee $w_{\lambda^*}(C^*)\in [L(\lambda^*),2L(\lambda^*))$ easily follows since $L(\lambda^*)-\lambda^* b = LD > 0$.
|
||||
|
||||
\subsection{complexity}
|
||||
|
||||
\begin{algo}
|
||||
\underbar{\textsc{$b$-free MinCut}}($G,w,c,b$):\\
|
||||
compute $\lambda^*$ using parametric search\\
|
||||
reweight $G$ with $w_\lambda$\\
|
||||
for each 2-approx mincut $C$ in $(G,w_\lambda)$:\\
|
||||
\quad run FPTAS for knapsack to compute $\min \set{w(C-F)| F\subset C, c(F)\leq b}$\\
|
||||
return the optimal $(C,F)$
|
||||
\end{algo}
|
||||
|
||||
To compute $\lambda^*$ we need to use parametric search.
|
||||
\begin{lemma}[\cite{salowe_parametric}]
|
||||
Let $S(n)$ be the complexity of solving the Lagrangian dual problem for fixed $\lambda$ (where $n$ is the size of the input), then one can compute $\lambda^*$ using parametric search in $O(S(n)^2)$ time.
|
||||
\end{lemma}
|
||||
It follows directly from the preceding lemma that $\lambda^*$ can be computed in $\tilde O(m^2)$ time.
|
||||
|
||||
Reweighting the graph takes linear time.
|
||||
Finding $<2$-approx mincut takes $\tilde O(n^3)$.
|
||||
An $1+\e$ approximate solution to knapsack can be found in time $\tilde O(m+\frac{1}{\e^2})$ \cite{10.1145/3618260.3649730}.
|
||||
The total complexity is $\tilde O(mn^3+\frac{n^3}{\e^2})$.
|
||||
|
||||
\bibliographystyle{plain}
|
||||
\bibliography{ref}
|
||||
|
||||
\appendix
|
||||
|
||||
\section{remove box constraints}
|
||||
Given a positive covering LP,
|
||||
|
||||
\begin{equation*}
|
||||
@@ -317,11 +401,10 @@ suppose that $c(f)<x_f\in x^*$. Consider all constraints $\sum_{e\in T\setminus
|
||||
x_e\geq k-c(f)$ on $T\ni f$. For any such constraint, we have $\sum_{e\in T} c(e)
|
||||
x_e>k$ since we assume $x_f>c(f)$, which means we can decrease $x_f$ without violating any constraint. Thus it contradicts the assumption that $x^*$ is optimal. Then we can add redundant constraints $x_e\leq c(e) \;\forall e$ to LP2 and see that LP1 and LP2 have the same optimum.
|
||||
|
||||
This applies to \cite{chalermsook_approximating_2022} but cannot get an improvement on their algorithm.(MWU does not care the number of constraints.) So does this trick apply to connectivity interdiction?
|
||||
This applies to \cite{chalermsook_approximating_2022} but cannot get an improvement on their algorithm.(MWU does not care the number of constraints.)
|
||||
|
||||
\[
|
||||
\min_{\text{cut C}, f\in C}\frac{\sum_{e\in C\setminus\set{f}}w(e)x_e}{k-c(f)}
|
||||
\]
|
||||
\bibliographystyle{plain}
|
||||
\bibliography{ref}
|
||||
|
||||
\end{document}
|
||||
|
||||
28
ref.bib
28
ref.bib
@@ -81,5 +81,33 @@ abstract = { This paper considers the problem of designing fast, approximate, co
|
||||
pages = {1334--1353},
|
||||
}
|
||||
|
||||
@inproceedings{10.1145/3618260.3649730,
|
||||
author = {Chen, Lin and Lian, Jiayi and Mao, Yuchen and Zhang, Guochuan},
|
||||
title = {A Nearly Quadratic-Time FPTAS for Knapsack},
|
||||
year = {2024},
|
||||
isbn = {9798400703836},
|
||||
publisher = {Association for Computing Machinery},
|
||||
address = {New York, NY, USA},
|
||||
url = {https://doi.org/10.1145/3618260.3649730},
|
||||
doi = {10.1145/3618260.3649730},
|
||||
abstract = {We investigate the classic Knapsack problem and propose a fully polynomial-time approximation scheme (FPTAS) that runs in O(n + (1/)2) time. Prior to our work, the best running time is O(n + (1/)11/5) [Deng, Jin, and Mao’23]. Our algorithm is the best possible (up to a polylogarithmic factor), as Knapsack has no O((n + 1/)2−δ)-time FPTAS for any constant δ > 0, conditioned on the conjecture that (min, +)-convolution has no truly subquadratic-time algorithm.},
|
||||
booktitle = {Proceedings of the 56th Annual ACM Symposium on Theory of Computing},
|
||||
pages = {283–294},
|
||||
numpages = {12},
|
||||
keywords = {Approximation scheme, Knapsack},
|
||||
location = {Vancouver, BC, Canada},
|
||||
series = {STOC 2024}
|
||||
}
|
||||
|
||||
|
||||
@incollection{salowe_parametric,
|
||||
author = {Salowe, Jeffrey S.},
|
||||
title = {Parametric search},
|
||||
year = {1997},
|
||||
isbn = {0849385245},
|
||||
publisher = {CRC Press, Inc.},
|
||||
address = {USA},
|
||||
booktitle = {Handbook of Discrete and Computational Geometry},
|
||||
pages = {683–695},
|
||||
numpages = {13}
|
||||
}
|
||||
Reference in New Issue
Block a user