---
title: Groupwise Harmonic Centrality
url: https://www.emergentmind.com/topics/groupwise-harmonic-centrality
type: topic
---

# Groupwise Harmonic Centrality

Groupwise Harmonic Centrality generalizes classic node-centric centrality by quantifying the collective reach of a group of nodes in a network, where the importance of nodes outside the group is inversely proportional to their shortest-path distances from the group. This measure is foundational in group centrality optimization, influence maximization, and diffusion analysis in static and probabilistic networks. Rigorous theoretical and practical advances have established its computational complexity, algorithmic approximability, and generalization to stochastic influence models [2010.15435][1810.09981]. 

## 1. Formal Definition and Variants

Given a (possibly directed) graph $G = (V, E, \ell)$ with positive edge-lengths $\ell: E \to \mathbb{R}_{>0}$, let $S\subseteq V$ be a nonempty set of seed nodes. For each $v \in V \setminus S$, define the group-to-node distance as
$$d(S, v) = \min_{u\in S} d(u, v),$$
where $d(u, v)$ denotes the shortest-path length from $u$ to $v$ (set to $\infty$ if $v$ is unreachable from $S$). The groupwise harmonic centrality is
$$
\mathrm{GH}(S) = \sum_{v\in V\setminus S} \frac{1}{d(S,v)}.
$$
The reciprocal $1/d(S,v)$ is interpreted as $0$ when $d(S,v)=\infty$.

Properties include submodularity (diminishing marginal gains), normalization ($\mathrm{GH}(\emptyset)=\mathrm{GH}(V)=0$), and monotonicity in the deterministic (static) graph case [1810.09981]. The $\delta$-sphere extension counts the number of nodes within radius $\delta$, interpolating between immediately adjacent sets and full reachability.

Influence-based group harmonic centrality is defined in stochastic or cascade models (such as independent cascade or linear threshold). For an influence instance $I=(V,E,P_I)$ and $S\subseteq V$, let $(S_0=S, S_1, ..., S_{n-1})$ denote a random cascading activation sequence generated according to $P_I(S)$. For each $u\notin S$, let $d_u(S_0, ..., S_{n-1})$ be the first time $u$ is activated, or $\infty$. Define the influence-based centrality as
$$
\Psi_{\mathrm{harm}}(S; I) = \mathbb{E}_{(S_1,...,S_{n-1})\sim P_I}\left[ \sum_{u\notin S} \frac{1}{d_u(S_0,...,S_{n-1})} \right].
$$

## 2. Computational Complexity and Inapproximability

The maximization problem — finding $S\subseteq V$ with $|S|=k$ maximizing $\mathrm{GH}(S)$ — is NP-hard even on unweighted, undirected or directed graphs. This is shown via reductions from Maximum Coverage (directed) and Dominating Set (undirected) [2010.15435]. 

Approximation hardness results are as follows:

| Graph Type    | Inapproximability Bound                       |
|---------------|----------------------------------------------|
| Directed      | $1-1/e \approx 0.632$                        |
| Undirected    | $1-1/(4e) \approx 0.407$                     |

No polynomial-time algorithm achieves a better approximation unless $P=NP$ (directed: reduction via set-cover bipartite digraphs giving exact correspondence; undirected: via gap-amplified dominating set reductions).

The function is submodular but not always monotone—adding nodes may reduce total centrality due to overlap in coverage, affecting some algorithmic guarantees [2010.15435][1810.09981].

## 3. Approximation Algorithms

### Greedy Algorithm

Despite non-monotonicity, a standard greedy approach is used for maximizing $\mathrm{GH}(\cdot)$. 

**Algorithmic Steps:**
1. Initialize $S\leftarrow \emptyset$.
2. While $|S|<k$, compute for each $u\in V\setminus S$ the marginal gain $\Delta(u) := \mathrm{GH}(S\cup\{u\}) - \mathrm{GH}(S)$. 
3. Select $v=\arg\max_{u}\Delta(u)$ and set $S\gets S\cup \{v\}$.
4. Return $S$.

Let $\lambda = \frac{\ell_{\min}}{\ell_{\max}}$ where $\ell_{\min}$ and $\ell_{\max}$ are respective minimum and maximum edge lengths. The guarantees are:

| Graph Type    | Guarantee (relative to OPT)                    |
|---------------|------------------------------------------------|
| Directed      | $\lambda(1-2/e) \approx 0.264\cdot\lambda$     |
| Undirected    | $(\lambda/2)(1-1/e)\approx 0.316\cdot\lambda$  |

Scaling from unweighted to weighted graphs reduces the guarantee by $\lambda$. The analysis exploits surrogates and telescoping inequalities in the undirected case to obtain the relevant bounds [2010.15435].

### Local-Search Algorithm (Undirected)

For undirected graphs, local search (single-swap) is effective:
- Start from an initial $k$-set $S$ (e.g., greedy output).
- While there exist $u\in S$ and $v\in V\setminus S$ such that $\mathrm{GH}(S-\{u\}+\{v\}) > (1+\varepsilon)\mathrm{GH}(S)$, perform the swap.

This yields a $(1+\varepsilon)/6$-approximation, supported by submodularity and matroid exchange properties [2010.15435].

### Influence-Based Extension

For influence models, the reverse-reachable (RR) set framework efficiently approximates influence-based group centrality. For each sampled RR set $R_v$ from node $v$, the marginal contribution is $1/d_{R_v}(S, v)$, and estimates are accumulated over many such sets for unbiased estimators. Relative-error guarantees can be made with appropriate number of samples, adapting the IMM framework [1810.09981]. 

## 4. Shapley-Value Decomposition and Uniqueness

Groupwise harmonic centrality can be interpreted as a cooperative-game value function on $2^{V}\to\mathbb{R}$. The Shapley value embedding decomposes the group score into individual node scores:
$$
\Phi_{\mathrm{harm}}(v; I) = \sum_{T\subseteq V\setminus\{v\}} \frac{|T|!(n-|T|-1)!}{n!}\Bigl[\Psi_{\mathrm{harm}}(T\cup\{v\};I)-\Psi_{\mathrm{harm}}(T;I)\Bigr].
$$

Equivalently, averaging over random orderings,
$$
\Phi_{\mathrm{harm}}(v; I) = \mathbb{E}_{\pi}\left[\Psi_{\mathrm{harm}}(S_{\pi,v} \cup \{v\}; I) - \Psi_{\mathrm{harm}}(S_{\pi,v}; I)\right],
$$
where $S_{\pi,v}$ is the set of nodes preceding $v$ in permutation $\pi$.

A characterization theorem establishes that influence-based group harmonic centrality is the unique Bayesian (linear in profile), permutation-invariant extension conforming to its deterministic (static-graph) definition. Layered-graph BFS instances form a basis for this uniqueness result [1810.09981].

## 5. Engineering and Scalability Considerations

Efficient maximization for large-scale graphs is achieved via several engineering techniques [2010.15435]: 

- **Lazy Marginal Evaluation**: Maintain a heap of overestimates $\hat{\mathrm{GH}}(S, u)\geq \Delta(u)$ to accelerate marginal computation; recompute exact values only when necessary (following the “lazy greedy” paradigm).
- **Pruned SSSP**: For candidate $u$, a BFS or Dijkstra search from $u$ is terminated at $v$ when $d(u,v)\geq d(S,v)$. Partial sums and upper bounds enable early pruning and avoid traversing the full graph.
- **Distance Labels**: Precompute and update $d_S[v]$ via multi-source BFS/Dijkstra per iteration. For local search, maintain the nearest and second-nearest group members for each $v$ to speed up swap marginal calculations.
- **Parallelization**: Both greedy and local-search procedures benefit from multi-threaded computation. Marginals for different candidates, and swap tests in local search, can be evaluated in parallel. Reported scaling to 16–32 cores yields $6$--$10\times$ speedups.

In practice, these optimizations reduce per-iteration work from $O(n+m)$ to empirically $O(n^{0.6}$–$n^{0.8})$ [2010.15435].

## 6. Experimental Results and Quality Benchmarks

Comprehensive experiments were conducted on [2010.15435]:
- **Small (n$\leq 2\,000$)**: KONECT, DIMACS, and set-cover–derived graphs for comparison to ILP optima.
- **Medium/Large (n=$10^4$–$2\cdot 10^5$, m$\leq 5\cdot 10^6$)**: Complex networks.
- **High-diameter weighted (n$\leq 10^4$)**: Road networks.

Findings indicate:
- On small graphs, greedy and local-search are within $0.5\%$ and $0.3\%$ of optimum, respectively; random baselines achieve $70-80\%$ of OPT.
- On large unweighted graphs, greedy yields a $40–60\%$ improvement over random, local search adds $1–2\%$ further quality at $5–20\times$ higher runtime.
- For high-diameter weighted graphs, greedy achieves $2.4–2.6\times$ the random baseline, with local-search again gaining $0.5–1\%$ but at $50$–$400\times$ the runtime.
- Both methods demonstrate excellent parallel scaling.

## 7. Practical Guidance and Theoretical Significance

- **Greedy algorithms** are recommended for large-scale problems ($n > 10^5$), balancing runtime with near-optimal solution quality ($>99\%$ of optimum), with provable guarantees parameterized by $\lambda$.
- **Local search** is favored when highest solution quality is essential and computational resources are less constrained; most beneficial for moderate graphs ($n\leq 10^4$) or small groups ($k \leq 20$).
- **Pruned search and lazy marginal updates** dramatically improve scalability.
- **Parallel implementation** is highly effective and straightforward for both greedy and local search, justifying use on modern multicore architecture.

Groupwise harmonic centrality encapsulates a rigorous group-based notion of accessibility and influence, with robust theoretical foundations, unique influence-based characterization, and scalable maximization algorithms applicable in real-world network analytics [2010.15435][1810.09981].

Source: https://www.emergentmind.com/topics/groupwise-harmonic-centrality