Jaccard-constrained dense subgraph discovery

Abstract

AbstractFinding dense subgraphs is a core problem in graph mining with many applications in diverse domains. At the same time many real-world networks vary over time, that is, the dataset can be represented as a sequence of graph snapshots. Hence, it is natural to consider the question of finding dense subgraphs in a temporal network that are allowed to vary over time to a certain degree. In this paper, we search for dense subgraphs that have large pairwise Jaccard similarity coefficients. More formally, given a set of graph snapshots and input parameter $$\alpha$$ α , we find a collection of dense subgraphs, with pairwise Jaccard index at least $$\alpha$$ α , such that the sum of densities of the induced subgraphs is maximized. We prove that this problem is NP-hard and we present a greedy, iterative algorithm which runs in $${\mathcal {O}} \mathopen {} \left( nk^2 + m\right)$$ O n k 2 + m time per single iteration, where k is the length of the graph sequence and n and m denote number of vertices and total number of edges respectively. We also consider an alternative problem where subgraphs with large pairwise Jaccard indices are rewarded. We do this by incorporating the indices directly into the objective function. More formally, given a set of graph snapshots and a weight $$\lambda$$ λ , we find a collection of dense subgraphs such that the sum of densities of the induced subgraphs plus the sum of Jaccard indices, weighted by $$\lambda$$ λ , is maximized. We prove that this problem is NP-hard. To discover dense subgraphs with good objective value, we present an iterative algorithm which runs in $${\mathcal {O}} \mathopen {}\left( n^2k^2 + m \log n + k^3 n\right)$$ O n 2 k 2 + m log n + k 3 n time per single iteration, and a greedy algorithm which runs in $${\mathcal {O}} \mathopen {}\left( n^2k^2 + m \log n + k^3 n\right)$$ O n 2 k 2 + m log n + k 3 n time. We show experimentally that our algorithms are efficient, they can find ground truth in synthetic datasets and provide good results from real-world datasets. Finally, we present two case studies that show the usefulness of our problem.