Inspired by [1]}, the constraint solely relies on the graph topology as follows:
\( C^{cont} = 1 - \frac{Tr \left(P^T (\overrightarrow{A} - diag(A)) P \right)}{Tr \left( P^T P - \frac{1}{{1}_{P_i^2 \ge \epsilon }} P^T P \right)}\) 
