Original: S UBSAMPLING IN L ARGE G RAPHS U SING R ICCI C URVATURE
Revised: S UBSAMPLING IN L ARGE G RAPHS U SING R ICCI C URVATURE

Original: Anonymous authors

Original: Paper under double-blind review

Original: A BSTRACT
Revised: A BSTRACT

Original: In the past decades, many large graphs with millions of nodes have been collected/constructed.
Revised: In recent decades, an increasing number of large graphs with millions of nodes have been collected and constructed.

Original: The high computational cost andsignificant visualization difficulty hinder the analysis of large graphs.
Revised: Despite their utility, analyzing such graphs is hindered by high computational costs and visualization difficulties.

Original: To overcome the difficulties, researchers have developedmany graph subsampling approaches  to provide a rough sketch  that preserves  global properties.
Revised: To tackle these issues, researchers have developed graph subsampling approaches that extract representative nodes to provide a rough sketch of the graph that preserves its global properties.

Original: By selecting representative nodes, these graph subsampling methods can help researchers estimate the graph statistics, e.g., the number of communities, of the large graph from the subsample.
Revised: By selecting representative nodes, these graph subsampling methods can help researchers estimate the graph statistics, e.g., the number of communities, of the large graph from the subsample.

Original: However, the available subsampling methods, e.g., degree node sampler and random walk sampler, tend to leave out minority communities because nodes with high degrees are more likely to be sampled.
Revised: However, the available subsampling methods, such as degree node sampler and random walk sampler, tend to overlook minority communities, as they prioritize nodes with high degrees.

Original: To overcome the shortcomings of the existing methods, we are motivated to apply the community information hidden in the graph to the subsampling method.
Revised: To address this limitation, we propose leveraging community information hidden within the graph for subsampling.

Original: Though the community structure is unavailable, community structure informationcan be obtained by applying geometric methods to a graph. An analog of Ricci curvature in the manifold is defined for the graph, i.e., Ollivier Ricci curvature.
Revised: Although the community structure is typically unknown, geometric methods can reveal community structure information by defining an analog of Ricci curvature on the graph, known as the Ollivier Ricci curvature.

Original: Based on the  asymptotic results about the within-communityedge and between-communities edge’s OR curvature, we propose a subsampling algorithm based on our theoretical results, the Ollivier-Ricci curvature Gradientbased subsampling (ORG-sub) algorithm.
Revised: We introduce a new subsampling algorithm called the Ollivier Ricci curvature Gradient-based subsampling (ORG-sub) algorithm based on our asymptotic results regarding the within-community and between-communities edges’ OR curvature.

Original: Theproposed ORG-sub algorithm has two main contributions: First, ORG-sub provides a rigorous theoretical guarantee that the probability ofORG-sub taking all communities into the final subgraph converges to one.
Revised: The ORG-sub algorithm makes two significant contributions: Firstly, it provides a rigorous theoretical guarantee that the probability of taking all communities into the final subgraph converges to one.

Original: Second, extensive experiments on synthetic and benchmark datasets demonstrate the advantages of our algorithm.
Revised: Secondly, extensive experiments on synthetic and benchmark datasets demonstrate the advantages of our algorithm.

Original:  I NTRODUCTION
Revised:  I NTRODUCTION

Original: As we enter the big data era, our capacity to access large graphs (a.k.a. networks) has provided unprecedented opportunities and challenges.
Revised: As we enter the big data era, our capacity to access large graphs (a.k.a. networks) has provided unprecedented opportunities and challenges.

Original: For example, (Wang et al., 2011) presented a Twitter social network, which has more than 190 million nodes (users) who generate more than 65 million edges (tweets) every day (Wang et al., 2011).
Revised: For example, Wang et al. (2011) presented a Twitter social network, which has more than 190 million nodes (users) who generate more than 65 million edges (tweets) every day (Wang et al., 2011).

Original: Such huge networks enable researchers to tackle more complex problems.
Revised: Such huge networks enable researchers to tackle more complex problems.

Original: However, they pose great challenges to storing, visualizing, and analyzing since their sheer volumes render many computational methods infeasible.
Revised: However, they pose great challenges to storing, visualizing, and analyzing since their sheer volumes render many computational methods infeasible.

Original: Graph subsampling.
Revised: Graph subsampling.

Original: Graph subsampling is a commonly used technique to address this scalability issue because of its simplicity and efficiency.
Revised: Graph subsampling is a commonly used technique to address this scalability issue because of its simplicity and efficiency.

Original: Graph subsampling aims to take a subgraph that preserves critical features of the full graph.
Revised: Graph subsampling aims to take a subgraph that preserves critical features of the full graph.

Original: Various graph subsampling methods that preserve different graph properties have been proposed, including node sampling (Mall et al., 2013; Zeng et al., 2019), edge sampling (Krishnamurthy et al., 2005), and exploration sampling (Goodman, 1961;
Revised: Various graph subsampling methods that preserve different graph properties have been proposed, including node sampling (Mall et al., 2013; Zeng et al., 2019), edge sampling (Krishnamurthy et al., 2005), and exploration sampling (Goodman, 1961;

Original: Leskovec et al., 2005; Hübler et al., 2008; Maiya & Berger-Wolf, 2010).
Revised: Leskovec et al., 2005; Hübler et al., 2008; Maiya & Berger-Wolf, 2010).

Original: Researchers evaluate the graph subsampling approaches by measuring the similarity between the features of the original graph and those of the subgraph.
Revised: Researchers evaluate the graph subsampling approaches by measuring the similarity between the features of the original graph and those of the subgraph.

Original: The features include, e.g., degree distribution (Adamic et al., 2001), minimum cut (Hu & Lau, 2013), and the number of triangles (Seshadhri et al., 2014).
Revised: The features include, e.g., degree distribution (Adamic et al., 2001), minimum cut (Hu & Lau, 2013), and the number of triangles (Seshadhri et al., 2014).

Original: An important graph feature that has received less attention is the number of communities (denoted by M ), which plays a crucial role in identifying the community structures.
Revised: An important graph feature that has received less attention is the number of communities (denoted by M ), which plays a crucial role in identifying the community structures.

Original: Network data often have natural communities, and the identification of these communities helps answer vital questions in a variety of fields (Rohe et al., 2011).
Revised: Network data often have natural communities, and the identification of these communities helps answer vital questions in a variety of fields (Rohe et al., 2011).

Original: For example, communities in social networks may represent groups of people who share a similar interest, and communities in protein-protein interaction networks could be regulatory modules of interacting proteins (Rohe et al., 2011).
Revised: For example, communities in social networks may represent groups of people who share a similar interest, and communities in protein-protein interaction networks could be regulatory modules of interacting proteins (Rohe et al., 2011).

Original: In this paper, we focus on the setting where M is a fixed model parameter, and there is a ground truth about M .
Revised: In this paper, we focus on the setting where M is a fixed model parameter, and there is a ground truth about M .

Original: This setting is widely used in many models, e.g., stochastic block model (SBM) (Holland et al., 1983) and its variants such as degree-corrected SBM (DCBM) (Karrer & Newman, 2011).
Revised: This setting is widely used in many models, e.g., stochastic block model (SBM) (Holland et al., 1983) and its variants such as degree-corrected SBM (DCBM) (Karrer & Newman, 2011).

Original: The SBM family are arguably the most widely-used generative model for community detection from a theoretical perspective.
Revised: The SBM family is arguably the most widely-used generative model for community detection from a theoretical perspective.

Original: In fact, Vaca-Ramírez & Peixoto (2022) performed a systematic analysis of the quality of fit of the SBM for 275 real networks, and observed that “SBM is capable of providing an accurate description for the majority of networks considered”.
Revised: In fact, Vaca-Ramírez & Peixoto (2022) performed a systematic analysis of the quality of fit of the SBM for 275 real networks and observed that “SBM is capable of providing an accurate description for the majority of networks considered”.

Original: Indeed, there are other settings where the number of communities is not fixed.
Revised: Indeed, there are other settings where the number of communities is not fixed.

Original: For example, Olhede & Wolfe (2014) used SBM to approximate a nonparametric graphon model, under which case the number of communities is a hyperparameter and is not fixed.
Revised: For example, Olhede & Wolfe (2014) used SBM to approximate a nonparametric graphon model, under which case the number of communities is a hyperparameter and is not fixed.

Original: The latter hyperparameter case is beyond the scope of this paper.
Revised: The latter hyperparameter case is beyond the scope of this paper.

Original: Under the setting where M is a model parameter, many community detection methods have been proposed, such as modularity maximization (Newman, 2006;
Revised: Under the setting where M is a model parameter, many community detection methods have been proposed, such as modularity maximization (Newman, 2006;

Original: Good et al., 2010), spectral clustering (Von Luxburg, 2007; Rohe et al., 2011;
Revised: Good et al., 2010), spectral clustering (Von Luxburg, 2007; Rohe et al., 2011;

Original: Liu et al., 2018), and pseudo-likelihood based methods (Amini et al., 2013; Wang et al., 2021).
Revised: Liu et al., 2018), and pseudo-likelihood based methods (Amini et al., 2013; Wang et al., 2021).

Original: The theoretical properties of most community detection methods, such as consistency and asymptotic distributions, are built based on the assumption that
Revised: The theoretical properties of most community detection methods, such as consistency and asymptotic distributions, are built based on the assumption that

Original: M is known (Ma et al., 2021).
Revised: M is known (Ma et al., 2021).

Original: In addition, M is usually required as an input for those community detection algorithms.
Revised: In addition, M is usually required as an input for those community detection algorithms.

Original: However, in practice, we do not have the information of M , which significantly diminishes the usefulness of the aforementioned methods.
Revised: However, in practice, we do not have the information of M , which significantly diminishes the usefulness of the aforementioned methods.

Original: Existing methods for estimating M is usually very expensive.
Revised: Existing methods for estimating M are usually very expensive.

Original: For example, the cross-validation method proposed by Li et al.
Revised: For example, the cross-validation method proposed by Li et al.

Original: (2020) requires a computational cost that is cubic in the number of nodes n .
Revised: (2020) requires a computational cost that is cubic in the number of nodes n .

Original: When there are thousands or millions of nodes, the computational cost is unaffordable.
Revised: When there are thousands or millions of nodes, the computational cost is unaffordable.

Original: Thus, it is highly desirable for subsampling methods to yield subgraphs with ˜ n << n nodes preserving the number of communities M , such that we can use it to get an accurate estimation while reducing the computational cost .
Revised: Thus, it is highly desirable for subsampling methods to yield subgraphs with ˜ n << n nodes preserving the number of communities M , such that we can use it to get an accurate estimation while reducing the computational cost.

Original: Despite many successful applications, existing subsampling methods tend to leave out minority communities (i.e., a community with a smaller number of nodes) because nodes with high degrees are more likely to be sampled into subgraphs.
Revised: Despite many successful applications, existing subsampling methods tend to leave out minority communities (i.e., a community with a smaller number of nodes) because nodes with high degrees are more likely to be sampled into subgraphs.

Original: Consequently, these subsampling methods usually underestimate M , especially for graphs with imbalanced community structures.
Revised: Consequently, these subsampling methods usually underestimate M , especially for graphs with imbalanced community structures.

Original: To overcome the shortcomings of existing methods, we develop a graph subsampling method that yields subgraphs that can be used to accurately estimate M .
Revised: To overcome the shortcomings of existing methods, we develop a graph subsampling method that yields subgraphs that can be used to accurately estimate M .

Original: Achieving this goal is challenging since the community structure is hidden and unavailable.
Revised: Achieving this goal is challenging since the community structure is hidden and unavailable.

Original: Fortunately, recent studies indicate that the community structure is a geometric phenomenon by considering a graph as a Riemannian geometric object (Ni et al., 2015).
Revised: Fortunately, recent studies indicate that the community structure is a geometric phenomenon by considering a graph as a Riemannian geometric object (Ni et al., 2015).

Original: Some insights into the community information can be obtained by applying some geometric methods to a graph (Ni et al., 2015; 2019;
Revised: Some insights into the community information can be obtained by applying some geometric methods to a graph (Ni et al., 2015; 2019;

Original: Sia et al., 2019).
Revised: Sia et al., 2019).

Original: Ollivier Ricci Curvature of Graph.
Revised: Ollivier Ricci Curvature of Graph.

Original: In particular, a graph can be regarded as a discrete version of a Riemannian manifold (Ni et al., 2019).
Revised: In particular, a graph can be regarded as a discrete version of a Riemannian manifold (Ni et al., 2019).

Original: A node of the graph is analogous to a point on a
Revised: A node of the graph is analogous to a point on a

Original: Riemannian manifold, and a pair of connected nodes in a graph is analogous to two points connected by a geodesic on a manifold.
Revised: Riemannian manifold, and a pair of connected nodes in a graph is analogous to two points connected by a geodesic on a manifold.

Original: The partition of a graph into communities is analogous to the geometric decomposition of a Riemannian manifold (Ni et al., 2019).
Revised: The partition of a graph into communities is analogous to the geometric decomposition of a Riemannian manifold (Ni et al., 2019).

Original: Ricci curvature is a key tool for the geometric decomposition of a Riemannian manifold.
Revised: Ricci curvature is a key tool for the geometric decomposition of a Riemannian manifold.

Original: It measures how the Riemannian manifold deviates from the flat manifold.
Revised: It measures how the Riemannian manifold deviates from the flat manifold.

Original: Recently, (Lin et al., 2011) defines an analog of Ricci curvature for the graph, i.e., the Ollivier Ricci (abbreviated as OR) curvature.
Revised: Recently, Lin et al. (2011) defines an analog of Ricci curvature for the graph, i.e., the Ollivier Ricci (abbreviated as OR) curvature.

Original:  In this paper, we theoretically show that the OR curvatures of edges within a densely connected community are asymptotically larger than those of edges between two sparsely connected communities (Ni et al., 2015; Gosztolai & Arnaudon, 2021), as the size of the graph increases.
Revised: Previous empirical results have shown the OR curvature is related to the connectivity of the graphs (Ni et al., 2015; Gosztolai & Arnaudon, 2021). However, the relationship between OR curvature and connectivity is insufficiently explored in theoretical evidence. In this paper, we theoretically show that the OR curvatures of edges within a densely connected community are asymptotically larger than those of edges between two sparsely connected communities, as the size of the graph increases.

Original: Ollivier Ricci Curvature Gradient Based Graph Subsampling.
Revised: Ollivier Ricci Curvature Gradient Based Graph Subsampling.

Original: Based on our theoretical result, we propose an OR curvature gradient-based graph subsampling algorithm (abbreviated as ORG-sub).
Revised: Based on our theoretical result, we propose an OR curvature gradient-based graph subsampling algorithm (abbreviated as ORG-sub).

Original: Specifically, ORG-sub randomly chooses one edge as the starting point of the subgraph and calculates the OR curvature of the selected edge.
Revised: Specifically, ORG-sub randomly chooses one edge as the starting point of the subgraph and calculates the OR curvature of the selected edge.

Original: ORG-sub then gradually expands the subgraph by taking the next edge whose OR curvature shows the largest difference from the OR curvature of the previously taken edge.
Revised: ORG-sub then gradually expands the subgraph by taking the next edge whose OR curvature shows the largest difference from the OR curvature of the previously taken edge.

Original: Here, we define the difference as the OR curvature gradient (ORG).
Revised: Here, we define the difference as the OR curvature gradient (ORG).

Original: We use the ORG to guide the expansion of ORG-sub, i.e., direct ORG-sub which edge to take next.
Revised: We use the ORG to guide the expansion of ORG-sub, i.e., direct ORG-sub which edge to take next.

Original: All edges that
Revised: All edges that

Original: ORG-sub took expand into the final subgraph.
Revised: ORG-sub took expand into the final subgraph.

Original: Our proposed ORG-sub enjoys two main advantages.
Revised: Our proposed ORG-sub enjoys two main advantages.

Original: First, ORG-sub has a rigorous theoretical guarantee.
Revised: First, ORG-sub has a rigorous theoretical guarantee.

Original: In particular, under the SBM scenario, we prove that the probability of ORG-sub taking all communities into the final subgraph converges to one, indicating that even nodes in the minority community are subsampled.
Revised: In particular, under the SBM scenario, we prove that the probability of ORG-sub taking all communities into the final subgraph converges to one faster than the random walk algorithm, indicating that even nodes in the minority community are subsampled.

Original: More importantly, we theoretically show that the estimation of M by the subsampled graph converges to M of the full graph.
Revised: More importantly, we theoretically show that the estimation of M by the subsampled graph converges to M of the full graph.

Original: Second, our extensive empirical experiments based on the simulated and real-world datasets show that the estimator based on ORG-sub subgraphs accurately estimates the number of communities while greatly reducing the computation cost.
Revised: Second, our extensive empirical experiments based on the simulated and real-world datasets show that the estimator based on ORG-sub subgraphs accurately estimates M while greatly reducing the computation cost.

Original:  P RELIMINARIES 2. N OTATIONS AND D EFINITIONS
Revised:  P RELIMINARIES 2. N OTATIONS AND D EFINITIONS

Original: A graph is made up of nodes that are connected by edges.
Revised: A graph is made up of nodes that are connected by edges.

Original: Considering that directed graphs don’t approximate Riemannian manifolds in discrete form since geodesics are always bidirectional, we only focus on undirected graphs in this paper.
Revised: Considering that directed graphs don’t approximate Riemannian manifolds in discrete form since geodesics are always bidirectional, we only focus on undirected graphs in this paper.

Original: Denote a graph as G = ⟨ V, E ⟩ , where V is the node set and E is the edge set.
Revised: Denote a graph as G = ⟨ V, E ⟩ , where V is the node set and E is the edge set.

Original: Operator | · | calculates the cardinality of the set.
Revised: Operator | · | calculates the cardinality of the set.

Original: We then denote the cardinality of the node set in a graph by | V | ( | V | = n ), the number of edges in a graph by | E | , and the neighborhood set of a node v by δ ( v ) =
Revised: We then denote the cardinality of the node set in a graph by | V | ( | V | = n ), the number of edges in a graph by | E | , and the neighborhood set of a node v by δ ( v ) =

Original: { w | ( v, w ) ∈ E } .
Revised: { w | ( v, w ) ∈ E } .

Original: The degree of a node v is defined as the cardinality of its neighborhood: d v = | δ ( v ) | .
Revised: The degree of a node v is defined as the cardinality of its neighborhood: d v = | δ ( v ) | .

Original: For a subset of V , S ⊆ V , the subgraph induced by the subset S is denoted by G [ S ] =
Revised: For a subset of V , S ⊆ V , the subgraph induced by the subset S is denoted by G [ S ] =

Original: ( S, E S ) , where the edge set E S
Revised: ( S, E S ) , where the edge set E S

Original: = { ( v S , w S )
Revised: = { ( v S , w S )

Original: ∈ E | v S ∈ S, w S ∈ S } .
Revised: ∈ E | v S ∈ S, w S ∈ S } .

Original: The neighborhood of the node set S is defined by N ( S ) =
Revised: The neighborhood of the node set S is defined by N ( S ) =

Original: (cid:83) v ∈ S
Revised: (cid:83) v ∈ S

Original: δ ( v ) , and the neighboring edge set of edge e = ( u, v ) is ∆( e ) =
Revised: δ ( v ) , and the neighboring edge set of edge e = ( u, v ) is ∆( e ) =

Original: { ( x, y ) | x ∈ { u, v } , y ∈
Revised: { ( x, y ) | x ∈ { u, v } , y ∈

Original: δ ( x ) \ { u, v }} .
Revised: δ ( x ) \ { u, v }} .

Original: Let ˆ M ( G [ S ]) denote the estimator of M , obtained by using the subgraph G [ S ] for estimation.
Revised: Let (cid:99) M ( G [ S ]) denote the estimator of M , obtained by using the subgraph G [ S ] for estimation.

Original: 2. O LLIVIER R ICCI C URVATURE OF G RAPHS
Revised: 2. O LLIVIER R ICCI C URVATURE OF G RAPHS

Original: A Riemannian manifold ( M , g ) is defined as a smooth manifold M equipped with a metric g :
Revised: A Riemannian manifold ( M , g ) is defined as a smooth manifold M equipped with a metric g :

Original: T M × T M →
Revised: T M × T M →

Original: R , where T M is the tangent space of M (Lee, 2018).
Revised: R , where T M is the tangent space of M (Lee, 2018).

Original: Ricci curvature of Riemannian manifold measures the deviation of the manifold from a flat manifold (Lee, 2006;
Revised: Ricci curvature of Riemannian manifold measures the deviation of the manifold from a flat manifold (Lee, 2006;

Original: Ni et al., 2015; Samal et al., 2018).
Revised: Ni et al., 2015; Samal et al., 2018).

Original: The “flat” means the distance between two nodes is the same as the distance between the two local spaces centered around the nodes.
Revised: The “flat” means the distance between two nodes is the same as the distance between the two local spaces centered around the nodes.

Original: In the discrete version, the space defined on the node is discretely uniformly distributed in the node’s neighborhood, then the distance between two spaces can be explained as the average distance between the neighborhoods of the two nodes.
Revised: In the discrete version, the space defined on the node is discretely uniformly distributed in the node’s neighborhood, then the distance between two spaces can be explained as the average distance between the neighborhoods of the two nodes.

Original: Through such an analogy, the counterpart of
Revised: Through such an analogy, the counterpart of

Original: Ricci curvature, Ollivier Ricci curvature (Lin et al., 2011), is defined in graphs.
Revised: Ricci curvature, Ollivier Ricci curvature (Lin et al., 2011), is defined in graphs.

Original: Previous empirical results show that the OR curvature of graphs is related to the Ricci curvature of the underlying manifold: the graph with large (or small) OR curvature correspond to the underlying manifold with large (small) Ricci curvature (Ni et al., 2015; Samal et al., 2018; Ni et al., 2019;
Revised: Previous empirical results show that the OR curvature of graphs is related to the Ricci curvature of the underlying manifold: the graph with large (or small) OR curvature corresponds to the underlying manifold with large (or small)

Original:  Gosztolai & Arnaudon, 2021).
Revised: Ricci curvature (Ni et al., 2015; Samal et al., 2018; Ni et al., 2019; Gosztolai & Arnaudon, 2021).

Original: Theoretical results have shown the OR curvature of random geometric graphs converges to the Ricci curvature of the underlying Riemannian manifold (Boguna et al., 2021; van der Hoorn et al., 2021).
Revised: Theoretical results have shown the OR curvature of random geometric graphs converges to the Ricci curvature of the underlying Riemannian manifold (Boguna et al., 2021; van der Hoorn et al., 2021).

Original: Here we use a toy example to illustrate the relationship between the OR curvature of graphs and the corresponding underlying manifold.
Revised: Here we use a toy example to illustrate the relationship between the OR curvature of graphs and the corresponding underlying manifold.

Original: In Figure 1, the hyperbolic space (with negative curvature) is especially suited for tree-structured graphs (Nickel & Kiela, 2017), the Euclidean space (with zero curvature) is especially suited for grid-structured graphs, and the spherical space (with positive curvature) is especially suited for complete graphs (Ni et al., 2015).
Revised: In Figure 1, the hyperbolic space (with negative curvature) is especially suited for tree-structured graphs (Nickel & Kiela, 2017), the Euclidean space (with zero curvature) is especially suited for grid-structured graphs, and the spherical space (with positive curvature) is especially suited for complete graphs (Ni et al., 2015).

Original: We now introduce the definition of OR curvature of the graph.
Revised: We now introduce the definition of OR curvature of the graph.

Original: For α ∈ [0 , 1] and any node u with degree d u , we first define the probability distribution of node u by m αu : 
Revised: For α ∈ [0 , 1] and any node u with degree d u , we first define the probability distribution of node u by m αu : [Equation] otherwise

Original: [Equation]
Revised: [Equation]

Original: Here α is to keep the probability mass of α at node u itself and distribute the rest uniformly over the neighborhood.
Revised: Here α is to keep the probability mass of α at node u itself and distribute the rest uniformly over the neighborhood.

Original: Following Ni et al.
Revised: Following Ni et al.

Original: ; Ye et al.
Revised: ; Ye et al.

Original: (2019), we set α = 0 .
Revised: (2019), we set α = 0 .

Original: This means that each node keeps 50% of the probability mass to itself.
Revised: This means that each node keeps 50% of the probability mass to itself.

Original: Let d ( u, v ) be the geodesic distance between nodes u and v , which is the shortest path between two nodes of a graph.
Revised: Let d ( u, v ) be the geodesic distance between nodes u and v , which is the shortest path between two nodes of a graph.

Original: We shall present the Wasserstein distance of two nodes in the graph, which is defined as the transportation distance between two nodes’ probability distributions m αu and m αv .
Revised: We shall present the Wasserstein distance of two nodes in the graph, which is defined as the transportation distance between two nodes’ probability distributions m αu and m αv .

Original: Definition 2.1 (Wasserstein distance) .
Revised: Definition 2.1 (Wasserstein distance) .

Original: Let X be the metric space with two probability distributions µand µ 2 with mass 1, respectively.
Revised: Let X be the metric space with two probability distributions µand µ 2 with mass 1, respectively.

Original: A transportation plan from µ 1 to µ 2 is a mapping ξ :
Revised: A transportation plan from µ 1 to µ 2 is a mapping ξ :

Original: X × X → [0 , 1] satisfying (cid:80) y ξ ( x, y )
Revised: X × X → [0 , 1] satisfying (cid:80) y ξ ( x, y )

Original: = µ 1 ( x ) and (cid:80) x ξ ( x, y )
Revised: = µ 1 ( x ) and (cid:80) x ξ ( x, y )

Original: = µ 2 ( y ) .
Revised: = µ 2 ( y ) .

Original: Substitute µ 1 and µ 2 by probability distributions of nodes m αu and m αv , the Wasserstein distance between two nodes is,
Revised: Substitute µ 1 and µ 2 by probability distributions of nodes m αu and m αv , the Wasserstein distance between two nodes is,

Original: [Equation]
Revised: [Equation]

Original: The Wasserstein distance W ( m αu , m αv ) can be computed by linear programming Lin et al.
Revised: The Wasserstein distance W ( m αu , m αv ) can be computed by linear programming (Lin et al., 2011).

Original: Then the OR curvature (Ollivier, 2007; Lin et al., 2011) over the edge ( u, v ) is defined as follows.
Revised: Then the OR curvature (Ollivier, 2007; Lin et al., 2011) over the edge ( u, v ) is defined as follows.

Original: Definition 2.2 (Ollivier Ricci Curvature) .
Revised: Definition 2.2 (Ollivier Ricci Curvature) .

Original: Given the Wasserstein distance W ( m αu , m αv ) and the geodesic distance d ( u, v ) , we can define the Ollivier Ricci curvature as:
Revised: Given the Wasserstein distance W ( m α u , m α v ) and the geodesic distance d ( u, v ) , we can define the Ollivier Ricci curvature as:

Original: [Equation]
Revised: [Equation]

Original:  O LLIVIER R ICCI C URVATURE G RADIENT AND C OMMUNITY S TRUCTURE
Revised:  O LLIVIER R ICCI C URVATURE G RADIENT AND C OMMUNITY S TRUCTURE

Original: The OR curvature of the graph reveals the properties of the underlying Riemannian manifold and provides insights into the community structure of graphs.
Revised: The OR curvature of the graph reveals the properties of the underlying Riemannian manifold and provides insights into the community structure of graphs.

Original: Previous work has shown that the large (or small) curvature corresponds to the edge which is more (or less) connected than a grid (Ni et al., 2015; Samal et al., 2018; Ni et al., 2019; Gosztolai & Arnaudon, 2021).[Equation]
Revised: Previous work has shown that the large (or small) curvature corresponds to the edge which is more (or less) connected than a grid (Ni et al., 2015; Samal et al., 2018; Ni et al., 2019; Gosztolai & Arnaudon, 2021).

Original: As shown in Figure 2, the left is a graph with five communities, and the right one is a graph with three communities generated by the stochastic block model.
Revised: As shown in Figure 2, the left is a graph with five communities, and the right one is a graph with three communities generated by the stochastic block model.

Original: The warmer color indicates a smaller curvature, and the colder color indicates a larger curvature.
Revised: The warmer color indicates a smaller curvature, and the colder color indicates a larger curvature.

Original: Previous empirical results and the above toy example have related the OR curvature of graphs to the community structure of the graph.
Revised: Previous empirical results and the above toy example have related the OR curvature of graphs to the community structure of the graph.

Original: However, rigorous theoretical proofs that the OR curvature of the within-community edge is larger than the betweencommunities edge in a graph are not developed yet.
Revised: However, rigorous theoretical proofs that the OR curvature of the within-community edge is larger than the betweencommunities edge in a graph are not developed yet.

Original: In the Lemma 3.1, we theoretically prove the lower bound of  within-community edges’ OR curvatures are larger than the upper bound of the between-communities edges’ OR curvatures in stochastic block models (SBM).
Revised: In the Theorem 3.3, we theoretically prove the lower bound of the within-community edges’ OR curvatures are larger than the upper bound of the between-communities edges’ OR curvatures in stochastic block models (SBM).

Original: 3. A SYMPTOTIC R ESULTS OF OR C URVATURE G RADIENT IN SBM
Revised: 3. A SYMPTOTIC R ESULTS OF OR C URVATURE G RADIENT IN SBM

Original: Lemma 3.1 (Ollivier-Ricci curvature of random graphs generated by stochastic block models) .
Revised: Lemma 3.1 (Ollivier-Ricci curvature of random graphs generated by stochastic block models) .

Original: Consider a graph with M communities generated from a stochastic block model (SBM) 
Revised: Consider a graph with M communities generated from a stochastic block model (SBM) B (cid:16) [Equation]

Original: B
Revised: p

Original: [Equation]
Revised: [Equation]

Original: in
Revised: in

Original: [Equation]
Revised: [Equation]

Original:  i =1 , p out , { n i } Mi =1 (cid:17) , where p ( i ) in represents the probability of an edge connecting nodes in the same community B i , p out represents the probability of an edge connecting nodes from different communities, and { n i } Mi =1 are the sizes of each community in the SBM.
Revised: M i =1 , p out , { n i } Mi =1 (cid:17) , where p ( i ) in represents the probability of an edge connecting nodes in the same community B i , p out represents the probability of an edge connecting nodes from different communities, and { n i } Mi =1 are the sizes of each community in the SBM.

Original: We assume that, for any i = 1 , ..., M , p ( i ) in > p out , and n i n converges to r i (1 > r i > 0) , as n goes to infinity.
Revised: We assume that, for any i = 1 , ..., M , p ( i ) in > p out , and n i n converges to r i (1 > r i > 0) , as n goes to infinity.

Original: Let n ∨ = max { n i } Mi =1 , and p ∨ = max (cid:110) p ( i ) in (cid:111)
Revised: Let n ∨ = max { n i } Mi =1 , and p ∨ = max (cid:110) p ( i ) in (cid:111)

Original: M i =1 .
Revised: M i =1 .

Original: We denote the node x from the community B i by x B i , the edge ( x B i , y B i ) where x is from B i and y is from B j , and the edge’s OR curvature by κ ( x B i , y B i ) .
Revised: We denote the node x from the community B i by x B i , the edge ( x B i , y B j ) where x is from B i and y is from B j , and the edge’s OR curvature by κ ( x B i , y B j ) .

Original: The following statements hold for the OR curvature.
Revised: The following statements hold for the OR curvature.

Revised: If

Revised: [Equation]

Original: If p out ≤ (cid:113) ( n i − 1)ln n n ln n i and p ( i ) in n i + p out ( n − n i ) > 83 ln n i , for the lower bound of within-community OR curvature, we almost surely have
Revised: and [Equation] [Equation] , for the lower bound of within-community OR curvature, we almost surely have

Original: [Equation]
Revised: [Equation]

Original: If p ( i ) in 2 ni + p out
Revised: If p ( i ) [Equation] 2 n

Original: [Equation]
Revised: [Equation]

Original: n
Revised: + p

Original: [Equation]
Revised: [Equation]

Original:  ln n i , for the upper bound of within-community OR curvature, we almost surely have
Revised: 4 ln n i , for the upper bound of within-community OR curvature, we almost surely have

Original: [Equation]
Revised: [Equation]

Original: If p out >
Revised: If p [Equation] >

Original: [Equation]
Revised: [Equation]

Original: , for the upper bound of betweencommunities OR curvature, we almost surely have
Revised: , for the upper bound of betweencommunities OR curvature, we almost surely have

Original: [Equation]
Revised: [Equation]

Original: Lemma 3.2.
Revised: Lemma 3.2.

Original: Combining the conclusion of statement 1 and statement 2 in Lemma 3.1, almost surely, for any ( x ′ B i , y ′ B j ) ∈ ∆(( x B i , y B i )) , we have 
Revised: Combining the conclusion of statement 1 and statement 2 in Lemma 3.1, almost surely, for any ( x ′ B i , y ′ B i ) ∈ ∆(( x B i , y B i )) , we have [Equation]

Original: [Equation]
Revised: [Equation]

Original: Statements 1 and 2 in Lemma 3.1 together show that under certain conditions, the probability that the
Revised: Statements 1 and 2 in Lemma 3.1 together show that under certain conditions, the probability that the

Original: OR curvature of within-community edges falls in ( ( n i − 2)
Revised: OR curvature of within-community edges falls in ( ( n i − 2)

Original: p
Revised: p

Original: ( i ) in 

Revised: ( i ) in

Original: +( n − n i ) p out 2 ( n [Equation] p
Revised: +( n − n i ) p out 2 ( n i − 1) p

Original: [Equation] ) in +( n − n i ) p out , p ( i ) in ) goes to 1 as the graph size n goes to infinity.
Revised: ( i ) in +( n − n i ) p out , p ( i ) in ) goes to 1 as the graph size n goes to infinity.

Original: Statement 3 in Lemma 3.1 shows that under certain conditions, the probability that the OR curvature of between-communities edges being no larger than n − n [Equation] p [Equation]
Revised: Statement 3 in Lemma 3.1 shows that under certain conditions, the probability that the OR curvature of between-communities edges being no larger than n − n i − n j n −p out + n i − n j − 2 n −p ∨ goes to 1 as n goes to infinity.

Original:  p 
Revised: Thus, we can conclude Lemma 3.2 that the probability of the maximum difference of two within-community edges’ OR curvatures being p out ( p ( i ) in − p out ) ni

Original: [Equation]
Revised: [Equation]

Original: goes to 1 as n goes to infinity.
Revised: goes to 1 as n goes to infinity.

Original: Thus, we can conclude Lemma 3.2 that the probability of the maximum difference of two within-community edges’ OR curvatures being p out ( p ( i ) in − p out ) ni −n − ni p ( i ) in + p out goes to 1 as n goes to infinity.

Original: Since the largest OR curvature difference between within-community edges can be bounded, we can prove the difference is small enough to distinguish from the difference between κ ( x B i , y B i ) and κ ( x ∗ B i , y ∗ B j ) .
Revised: Since the largest OR curvature difference of two within-community edges can be bounded, we can prove the difference is small enough to distinguish from the difference of κ ( x B i , y B i ) and κ ( x ∗ B i , y ∗ B j ) .

Original: Theorem 3.3.
Revised: Theorem 3.3.

Original: Combining the conclusion of statement 1 and statement 3 in Lemma 3.1, if p ∨ < p ( i ) in n n i + n j , p ( i ) in < p out ( n − n i )( n + n i + n j ) n 2 − ( n i + n j )( n − 2 n i ) , and conditions given in the lemma 3.1 are satisfied, almost surely, for any ( x ∗ B i , y ∗ B j ) ∈ ∆(( x B i , y B i )) , we have
Revised: Combining the conclusion of statement 1 and statement 3 in Lemma 3.1, if p ∨ < p ( i ) in n n i + n j , p ( i ) in < p out ( n − n i )( n + n i + n j ) n 2 − ( n i + n j )( n − 2 n i ) , and conditions given in the lemma 3.1 are satisfied, almost surely, for any ( x ∗ B i , y ∗ B j ) ∈ ∆(( x B i , y B i )) , we have

Original: [Equation]
Revised: [Equation]

Original: The Theorem 3.3 shows the OR curvature of the within-community edge κ ( x B i , y B i ) is larger than the OR curvature of the between-communities edge κ ( x B i , y B j ) by the maximum of | κ ( x B i , y B i ) − κ ( x ′ B
Revised: The Theorem 3.3 shows the OR curvature of the within-community edge κ ( x B i , y B i ) is larger than the OR curvature of the between-communities edge κ ( x B i , y B j ) by the maximum of | κ ( x B i , y B i ) − κ ( x ′ B

Original: i , y ′ B i ) | .
Revised: i , y ′ B i ) | .

Original: The theorem holds under mild conditions.
Revised: The theorem holds under mild conditions.

Original: When p out is small enough, and p ∨ is not too large compared to p out , the conditions can be satisfied.
Revised: When p out is small enough, and p ∨ is not too large compared to p out , the conditions can be satisfied.

Original: When the graph size is large, e.g., thousands of or millions of nodes, the constraint of p out and p ( i ) in is loosened, and the inequality in Theorem 3.3 can be easily satisfied.
Revised: When the graph size is large, e.g., thousands of or millions of nodes, the constraint of p out and p ( i ) in is loosened, and the inequality in Theorem 3.3 can be easily satisfied.

Original: 3. OR C URVATURE G RADIENT AND C OMMUNITIES
Revised: 3. OR C URVATURE G RADIENT AND C OMMUNITIES

Original: We define the OR curvature gradient (ORG) as | κ ( x ( j ) , y ( j ) )
Revised: We define the OR curvature gradient (ORG) as | κ ( x ( j ) , y ( j ) )

Original: − κ ( x ( i ) , y ( i ) )
Revised: − κ ( x ( i ) , y ( i ) )

Original: | , where ( x ( j ) , y ( j ) ) ∈ ∆(( x ( i ) , y ( i ) )) .
Revised: | , where ( x ( j ) , y ( j ) ) ∈ ∆(( x ( i ) , y ( i ) )) .

Original: It is the difference between the OR curvature of two adjacent edges.
Revised: It is the difference between the OR curvature of two adjacent edges.

Original: According to the theoretical results, the within-community ORG | κ ( x B i , y B i ) − κ ( x ′ B i , y ′ B i ) | is significantly different from the between-communities ORG | κ ( x B i , y B i ) − κ ( x ∗ B
Revised: According to the theoretical results, the within-community ORG | κ ( x B i , y B i ) − κ ( x ′ B i , y ′ B i ) | is significantly smaller than the between-communities ORG | κ ( x B i , y B i ) − κ ( x ∗ B

Original: i , y ∗ B j ) | , this motivates us to propose a subsampling algorithm based on ORG, which extracts the community information with theoretical guarantee.
Revised: i , y ∗ B j ) | , this motivates us to propose a subsampling algorithm based on ORG, which extracts the community information with theoretical guarantee.

Original:  OR C URVATURE G RADIENT B ASED G RAPH S UBSAMPLING
Revised:  OR C URVATURE G RADIENT B ASED G RAPH S UBSAMPLING

Original: The existing graph subsampling methods usually use degree information and random walk through the graph.
Revised: The existing graph subsampling methods usually use degree information and random walk through the graph.

Original: However, the subsampled graph tends to leave out minority communities (i.e., a community with a smaller number ) since they tend to sample nodes with high degrees.
Revised: However, the subsampled graph tends to leave out minority communities (i.e., a community with a smaller number ) since they tend to sample nodes with high degrees.

Original: Due to the lack of taking advantage of the community structure information during the expansion, the available subsampling algorithms usually underestimate M .
Revised: Due to the lack of taking advantage of the community structure information during the expansion, the available subsampling algorithms usually underestimate M .

Original: From the previous work, the community structure relates to the geometric phenomenon of the underlying Riemannian Manifold.
Revised: From the previous work, the community structure relates to the geometric phenomenon of the underlying Riemannian Manifold.

Original: Thus, we develop a fast and efficient ORG-based subsampling algorithm that is able to traverse the communities of the graph instead of trapping in a single community.
Revised: Thus, we develop a fast and efficient ORG-based subsampling algorithm that is able to traverse the communities of the graph instead of trapping in a single community.

Original: We then prove the probability of the proposed ORG-sub taking all communities of the full graph into the subsample converging to one.
Revised: We then prove the probability of the proposed ORG-sub taking all communities of the full graph into the subsample converges to one as the subsample size increases.

Original: Consequently, the subsampled graph attained by the ORG-sub algorithm can get a larger proportion of the minority communities than other subsampling algorithms and helps estimate M more accurately, as shown in Figure 3.
Revised: Consequently, the subsampled graph attained by the ORG-sub algorithm can get a larger proportion of the minority communities than other subsampling algorithms and helps estimate M more accurately, as shown in Figure 3.

Original: In addition, we theoretically show the estimation of M by the subsampling algorithm converges to true M .
Revised: In addition, we theoretically show the estimation of M by the subsampling algorithm converges to true M .

Original: 4. S UBSAMPLING A LGORITHM
Revised: 4. S UBSAMPLING A LGORITHM

Original: The ORG carries the community structure information since it differentiates within-community edges and between-communities edges with theoretical guarantees in Theorem 3.3.
Revised: The ORG carries the community structure information since it differentiates within-community edges and between-communities edges with theoretical guarantees in Theorem 3.3.

Original: We can use ORG as a guide for expanding the subgraph in the direction of more communities.
Revised: We can use ORG as a guide for expanding the subgraph in the direction of more communities.

Original: First, we randomly sample an edge ( x
Revised: First, we randomly sample an edge ( x

Original: , y ) ∈ E as the starting edge.
Revised: , y ) ∈ E as the starting edge.

Original: The probability of obtaining an edge ( x, y ) atthe start is
Revised: The probability of obtaining an edge ( x, y ) at start is

Original: P (( x (1) , y (1) ))
Revised: P (( x (1) , y (1) ))

Original: =
Revised: =

Original: | E | .
Revised: | E | .

Original: To take advantage of the theoretical properties of the ORG, the ORG-based subsampler expands to the next edge whose OR curvature shows the greatest difference from the OR curvature of the previously taken edge.
Revised: To take advantage of the theoretical properties of the ORG, the ORG-based subsampler expands to the next edge whose OR curvature shows the greatest difference from the OR curvature of the previously taken edge.

Original: The edge subsampled in the i + 1 -th step ( x ( i +1) , y ( i +1) ) given the edge subsampled in the i -th step ( x ( i ) , y ( i ) ) can be expressed by:
Revised: The edge subsampled in the i + 1 -th step ( x ( i +1) , y ( i +1) ) given the edge subsampled in the i -th step ( x ( i ) , y ( i ) ) can be expressed by:

Original: [Equation]
Revised: [Equation]

Original: Since the difference between the within-community ORG is smaller than the between-community
Revised: Since the difference between the within-community ORG is smaller than the between-community

Original: ORG, the proposed subsampler will expand to another community instead of trapping in one community.
Revised: ORG, the proposed subsampler will expand to another community instead of trapping in one community.

Original: After the subsampler stops expanding or the subsampling budget is used up, we get subsampled nodes.
Revised: After the subsampler stops expanding or the subsampling budget is used up, we get subsampled nodes.

Original: Then we can get the subgraph induced by the subsampled nodes.
Revised: Then we can get the subgraph induced by the subsampled nodes.

Original: Details of the algorithm are in Algorithm 1.
Revised: Details of the algorithm are in Algorithm 1.

Original: Algorithm
Revised: Algorithm

Original: Input: Graph G

Original:  ; Number of nodes to be subsampled ˜ n ; OR curvature of the graph G calculated by parameter α ; the subsampled node set S = ∅
Revised: Input: Graph G ; Number of nodes to be subsampled ˜ n ; OR curvature of the graph G calculated by parameter α ; the subsampled node set S = ∅

Original: Initialization: Randomly
Revised: Initialization Randomly

Original: choose a node v 0 as the start of sampling.
Revised: choose a node v 0 as the start of sampling

Original: Here, t = 0 .
Revised: Here, t = 0 .

Original: Cold Start
Revised: Cold Start

Original: : Among all neighbors of v
Revised: : Among all neighbors of v

Original: , randomly 
Revised: , randomly add

Original: add v1 to S 
Revised: v to S .

Original: Here, t = 1 .
Revised: Here, t = 1 .

Original: While: | S | < ˜ n
Revised: While: | S | < ˜ n

Original: • Step 1:
Revised: • Step 1:

Original: Given the edge e t = ( v t − 1 , v t ) selected in the t -th step, get the edge curvatures of the e t ’s neighboring edges ∆( e t ) .
Revised: Given the edge e t = ( v t − 1 , v t ) selected in the t -th step, get the edge curvatures of the e t ’s neighboring edges ∆( e t ) .

Original: • Step 2: Select the edge e t +1 in the edge set ∆( e t ) that has the greatest OGR.
Revised: • Step 2: Select the edge e t +1 in the edge set ∆( e t ) that has the greatest OGR.

Original: • Step 3:
Revised: • Step 3:

Original: Add the end point of the edge e t +1 to the subsampled node set S .
Revised: Add nodes connecting the edge e t +1 to the subsampled node set S .

Original: • t = t + 1 .
Revised: • t = t + 1 .

Original: Output: Subsampled node [Figure]
Revised: Output: Subsampled node set S and the induced subgraph G [ S ]

Revised: 4. ORG- SUB FOR E STIMATING M

Original: Applying the proposed ORG-sub algorithm to the full graph, as shown in Figure 3 (A), we can get a subsample as shown in Figure 3 (B).
Revised: Applying the proposed ORG-sub algorithm to the full graph, as shown in Figure 3 (A), we can get a subsample as shown in Figure 3 (B).

Original: Compared with the subsample obtained by the degree-based subsampler in Figure 3 (C), the proposed algorithm can subsample more nodes in minority communities.
Revised: Compared with the subsample obtained by the degreebased subsampler in Figure 3 (C), the proposed algorithm can subsample more nodes in minority communities.

Original: The toy example illustrates that the proposed ORG-sub algorithm can preserve the community structure information better.
Revised: The toy example illustrates that the proposed ORG-sub algorithm preserves the community structure information better.

Original: Thus, the subsampled graph can help estimate M of the full graph.
Revised: Thus, the subsampled graph can help estimate M of the full graph.

Original: We take multiple subgraphs by applying ORG-sub r times.
Revised: We take multiple subgraphs by applying ORG-sub r times.

Original: For each subsampled graph G [ S  ] given the sample size ˜ n , we get an estimation of the number of communities, ˆ M ( G [ S  ]) , then the mean of ˆ M ( G [ S  ]) is calculated as our final estimation of M .
Revised: For each subsampled graph G [ S i ] given the sample size ˜ n at i -th time, we get an estimation of M , (cid:99) M ( G [ S i ]) , then the mean of (cid:99) M ( G [ S i ]) is calculated as our final estimation of M .

Original: As for the choice of the algorithm estimating ˆ M ( G [ S ]) , we use the state-of-the-art network cross-validation method (Li et al., 2020).Compared with other methods (Wang & Bickel, 2017; Chen & Lei, 2018), network cross-validation is more robust to complicate settings, i.e., the number of communities is large or the community structure is not conspicuous.
Revised: As for the choice of the estimation algorithm, we use the state-of-the-art network cross-validation method (Li et al., 2020). [Figure] Compared with other methods (Wang & Bickel, 2017; Chen & Lei, 2018), network crossvalidation is more robust to complicated settings, i.e., the number of communities is large or the community structure is not conspicuous.

Original: Given the graph G generated by Theorem 3.3, we apply the ORG-sub algorithm 1 to the graph G and take graph subsamples with size ˜ n denoted by G [ S ] .
Revised: Given the graph G generated by Theorem 3.3, we apply the ORG-sub algorithm 1 to the graph G and take graph subsamples with size ˜ n denoted by G [ S ] .

Original: Let v denote a node in G [ S ] , Let P ( v ∈ B i ) denote the probability of the subgraph G [ S ] carrying nodes in the community B i .
Revised: Let v denote a node in G [ S ] , and P ( v ∈ B i ) denote the probability of the subgraph G [ S ] carrying nodes in the community B i .

Original: Theorem 4.1 shows that the probability of the ORG-sub algorithm expanding to all communities converging to one as the expanding steps increase.
Revised: Theorem 4.1 shows that the probability of the ORG-sub algorithm expanding to all communities converges to one as the expanding steps increase.

Original: Theorem 4.1 (Probability of the ORG-sub’s Expansion to Communities) .
Revised: Theorem 4.1 (Probability of the ORG-sub’s Expansion to Communities) .

Original: Under the conditions in
Revised: Under the conditions in

Original: Theorem 3.3, as ˜ nn → 1 , we have is
Revised: Theorem 3.3, as ˜ nn → 1 , we have P ( v ∈ B i | v ∈ G

Original: [Equation]
Revised: [ S ])

Original:  In addition, based on Theorem 4.1, we prove the theoretical advantage of our method over a popular graph sampling method, i.e., random walk sampling (Lovász, 1996).
Revised: → 1 . In addition, based on Theorem 4.1, we prove the theoretical advantage of our method over a popular graph sampling method, i.e., random walk sampling (Lovász, 1996).

Original: Let G [ S ] denote a sampled subgraph of ˜ n nodes, obtained by the ORG-sub method, and G RW [ S ] denote a sampled subgraph obtained by the random walk sampling method.
Revised: Let G [ S ] denote a sampled subgraph of ˜ n nodes, obtained by the ORG-sub method, and G RW [ S ] denote a sampled subgraph obtained by the random walk sampling method.

Original: Let v RW denote a node in G RW [ S ] .
Revised: Let v RW denote a node in G RW [ S ] .

Original: Corollary 4.1.
Revised: Corollary 4.1.1 shows that P ( v ∈ B i | v ∈ G

Original: shows that P ( v ∈ B i ) is greater than P ( v RW ∈ B i ) .
Revised: [ S ]) is greater than P ( v RW ∈ B i | v ∈ G [ S ]) .

Original: This means the ORG-sub is more likely to obtain a subgraph containing nodes from different communities.
Revised: This means the ORG-sub is more likely to obtain a subgraph containing nodes from different communities.

Original: Corollary 4.1.1. (Theoretical advantage over random walk).
Revised: Corollary 4.1.1 (Theoretical advantage over random walk) .

Original: Under the assumptions in Theorem 4.1,  we have 
Revised: Under the assumptions in Theorem 4.1, for any community B i in graph G , we have P ( v RW ∈ B i | v ∈ G [ S ]) < P ( v ∈ B i | v ∈ G

Original: [Equation]
Revised: [ S ]) .

Original: Since Theorem 4.1 theoretically guaranteesa certain proportion of each community in the full graph is subsampled, the convergence of the network cross-validation method in (Li et al., 2020) holds for the subsampled graphs.
Revised: Theorem 4.1 theoretically guarantees each community in the full graph is subsampled with high probability, thus, the convergence of the estimation of M by the ORG-sub algorithm can be theoretically guaranteed referring to Theorem 3 in (Li et al., 2020).

Original: We can then prove the convergence of the  ORG-sub-based estimation ˆ M ( G [ S i ]) , referring to Theorem 3 in (Li et al., 2020).
Revised: Since cross-validation methods tend not to give (cid:99) M > M in SBMs (Li et al., 2020; Chen & Lei, 2018), we can provide a theoretical guarantee against under-estimation of the proposed ORG-sub-based estimator (cid:99) M ( G [ S i ]) .

Original: Corollary 4.1.2.
Revised: Corollary 4.1.2.

Original: If the subsample with size ˜ n satisfies ˜ nρ log ˜ n → ∞ , where ρ is the edge density,and as the subsampling times r → ∞ , the probability the subsampled graph underestimates the true M is:
Revised: If the subsample with size ˜ n satisfies ˜ nρ log ˜ n → ∞ , where ρ is the edge density, as the subsampling times r → ∞ , the probability the subsampled graph underestimates the true M is:

Original: [Equation]
Revised: [Equation]

Original: 4. C OMPUTATIONAL A NALYSIS OF ORG- SUB
Revised: 4. C OMPUTATIONAL A NALYSIS OF ORG- SUB

Revised: Subsampling Proportions E s t i m a t i on o f t he N u m be r o f C o

Original: [Figure]
Revised: mm un i t i e s SBM DCBM 0.

Revised: Subsampling Proportions E s t i m a t i on o f t he N u m be r o f C o

Revised: mm un

Revised: i t i e s

Original: Computational complexity of Algorithm 1.
Revised: Computational complexity of Algorithm 1.

Original: We analyze the time complexity of Algorithm 1 step by step.
Revised: We analyze the time complexity of Algorithm 1 step by step.

Original: The time complexity of querying the neighboring edges and their corresponding curvatures is of order
Revised: The time complexity of querying the neighboring edges and their corresponding curvatures is of order

Original: O (  d ) .
Revised: O ( n ˜ d ) by Nys-sink algorithm (Altschuler et al., 2019), given ˜ d is the average degree.

Original: The time complexity of the sorting step to attain the most different curvature from the neighbors is O 
Revised: The time complexity of the sorting step to attaining the most different curvature from the neighbors is O ( ˜ d log ˜ d ) .

Original: ( d log d ) , where d is the degeneracy of graph G (The degeneracy of the graph is smaller than or equal to the maximum degree of the graph).

Original: Thus, the time complexity of taking one step is about O ( d log d ) .
Revised: Thus, the time complexity of taking one step is of order O ( n ˜ d ) .

Original: To sample ˜ n nodes, we need to run about ˜ n steps and get the induced subgraph from the sampled node set S .
Revised: To sample ˜ n nodes, we need to run about ˜ n steps and get the induced subgraph from the sampled node set S .

Original: Getting induced subgraph has the time complexity of order O (˜ nd ) . Thus the total complexity of the sampling procedure is O(˜ n (  d log d )) .
Revised: Since getting an induced subgraph of size ˜ n takes time complexity of order O (˜ n ˜ d ) , the total complexity of the sampling procedure is O ( n ˜ n ˜ d ) .

Original: Choice of the subsample Size.
Revised: Choice of the subsample Size.

Original: When implementing our algorithm, there are some hyper-parameters to be specified. One important hyperparameter is the size of the subsamples which determines the performance and efficiency of the algorithm.
Revised: When implementing our algorithm, we need to specify the size of the subsamples, which determines the performance and efficiency of the algorithm.

Original: Empirically, we plot the estimation of the number of communities as the function of the subsampling proportion. As shown in Figure 4,  we  select the elbow where the estimation of the number of communities starts to converge  to save computational resources and stabilize the estimation results.
Revised: Empirically, we plot the estimation of M as the function of the subsampling proportion in Figure 4, and we can select the elbow where the estimation of M starts to converge in applications to save computational resources and stabilize the estimation results.

Original:  E XPERIMENT
Revised:  E XPERIMENT

Original: We evaluate the performance of our algorithm on both synthetic datasets and real-world datasets.
Revised: We evaluate the performance of our algorithm on both synthetic datasets and real-world datasets.

Original: We use the metric, MAE = 1 r (cid:80) ri =1 (cid:12)(cid:12)(cid:12) M − ˆ M( G [ S i ])
Revised: We use the metric, MAE = 1 r (cid:80) ri =1 (cid:12)(cid:12)(cid:12) M − (cid:99) M

Original:  (cid:12)(cid:12)(cid:12) to evaluate the accuracy of subsampling based estimation.
Revised: ( G [ S i ]) (cid:12)(cid:12)(cid:12) to evaluate the accuracy of subsampling based estimation.

Original: We compare the proposed method with (1) Degree-based Node sampler (DBN) (Adamic et al., 2001), (2) Community Structure Expansion Sampler (CSE) (Maiya & Berger-Wolf, 2010) which is a community structure preserved sampler, and six benchmark exploration-based graph subsampling methods, including (3) Metropolis Hasting Random Walk Sampler (MHRW)
Revised: We compare the proposed method with (1) Degree-based Node sampler (DBN) (Adamic et al., 2001), (2) Community Structure Expansion Sampler (CSE) (Maiya & Berger-Wolf, 2010) which is a community structure preserved sampler, and six benchmark exploration-based graph subsampling methods, including (3) Metropolis Hasting Random Walk Sampler (MHRW)

Original: (Hübler et al., 2008), (4) Forest Fire Sampler (FFS) (Leskovec et al., 2005), (5) Snowball Sampler (Goodman, 1961), (6) Random Walk Sampler (RW)
Revised: (Hübler et al., 2008), (4) Forest Fire Sampler (FFS) (Leskovec et al., 2005), (5) Snowball Sampler (Goodman, 1961), (6) Random Walk Sampler (RW)

Original: (Gjoka et al., 2010), (7) Multi-dimensional Random Walk
Revised: (Gjoka et al., 2010), (7) Multi-dimensional Random Walk

Original: Sampler (MDRW) (Ribeiro & Towsley, 2010).
Revised: Sampler (MDRW) (Ribeiro & Towsley, 2010).

Original: As for other methods, we set the hyper-parameters as the default in package Little Ball of Fur (Rozemberczki et al., 2020b).
Revised: As for other methods, we set the hyper-parameters as the default in package Little Ball of Fur (Rozemberczki et al., 2020b).

Original: We replicate the experiments 30 times under each setting and compare the performance in terms of MAE for all methods.
Revised: We replicate the experiments 30 times under each setting and compare the performance in terms of MAE for all methods.

Original: All the experiments are conducted on a machine with a 40-core NVIDIA Tesla V100 GPU (3.00 GHz).
Revised: All the experiments are conducted on a machine with a 40-core NVIDIA Tesla V100 GPU (3.00 GHz).

Original: 5. D ATASETS
Revised: 5. D ATASETS

Original: Synthetic Dataset We generate synthetic datasets by the stochastic block models (SBM) and degreecorrected block models (DCBM), which can assign the community distribution to each node.
Revised: Synthetic Dataset We generate synthetic datasets by the stochastic block models (SBM) and degreecorrected block models (DCBM), which can assign the community distribution to each node.

Original: We set the community proportion as (3 / 4 , 1 / 10 , 1 / 12 , 1 / 15) with 900 nodes in total.
Revised: We set the community proportion as (3 / 4 , 1 / 10 , 1 / 12 , 1 / 15) with 900 nodes in total.

Original: The out-in-ratio (the ratio of between-communities edges over within-community edges) controls the ratio of the number of edges between the communities and within a community.
Revised: The out-in-ratio (the ratio of between-communities edges over within-community edges) controls the ratio of the number of edges between the communities and within a community.

Original: A higher ratio represents a noisier graph.
Revised: A higher ratio represents a noisier graph.

Original: The degree corrected model corrects the node degree by a power-law distribution.
Revised: The degree-corrected model corrects the node degree by a power-law distribution.

Original: Given the probability of an edge within a community ( p in = 0 . 8 ), we vary the probability of an edge between communities ( p out ) ( { 0 . 06 , 0 .
Revised: Given the probability of an edge within a community ( p in = 0 . 8 ), we vary the probability of an edge between communities ( p out ) ( { 0 . 06 , 0 .

Original: and subsampling proportions ( { 0 . 1 , 0 .
Revised: and subsampling proportions ( { 0 . 1 , 0 .

Original: 14 , 0 . 16 } ) from low to high for both block models to observe how our method performs as the setting becomes more challenging compared with others.
Revised: 14 , 0 . 16 } ) from low to high for both block models to observe how our method performs as the setting becomes more challenging compared with others.

Original: More details about the generation of the DCBM datasets are presented in section B.1.1 in Appendix.
Revised: More details about the generation of the DCBM datasets are presented in section B.1.1 in Appendix.

Original: Real-world Dataset We use five widely-used read-world graph datasets with labeled community structures to validate the performance of our method, including Polbooks, Facebook, Cora, Polblogs, and PubMed (Rossi & Ahmed, 2015;
Revised: Real-world Dataset We use five widely-used read-world graph datasets with labeled community structures to validate the performance of our method, including Polbooks, Facebook, Cora, Polblogs, and PubMed (Rossi & Ahmed, 2015;

Original: Rozemberczki et al., 2020a; Sen et al., 2008).
Revised: Rozemberczki et al., 2020a; Sen et al., 2008).

Original: All these graphs are considered unweighted and undirected.
Revised: All these graphs are considered unweighted and undirected.

Original: Besides, all the self-loop edges and isolated nodes are removed.
Revised: Besides, all the self-loop edges and isolated nodes are removed.

Original: Table A.4 in the Appendix summarizes the network statistics of these datasets.
Revised: Table A.4 in the Appendix summarizes the network statistics of these datasets.

Original: As we can see, the number of nodes ranges from 105 to 19,717, and the network density range from 0.0001 to 0.04.
Revised: As we can see, the number of nodes ranges from 105 to 19,717, and the network density ranges from 0.0001 to 0.04.

Original: Thus the networks we considered span a wide range.
Revised: Thus the networks we considered span a wide range.

Original: 5. R ESULTS OF S YNTHETIC D ATASET
Revised: 5. R ESULTS OF S YNTHETIC D ATASET

Original: [Figure]
Revised: [Figure]

Original: The subsampling times r for estimating the number of communities is set as 3, which is enough to get a stable estimation result Li et al.
Revised: The subsampling time r for estimating M is set as 3, which is enough to get a stable estimation result (Li et al., 2020).

Original: The results of SBM and DCBM datasets are reported in Figure 5.
Revised: The results of SBM and DCBM datasets are reported in Figure 5.

Original: The error bar is the standard deviation of 30 replications.
Revised: The error bar is the standard deviation ofreplications.

Original: We set the subsampling proportion (prop) as 0.1 when varying the probability of an edge between communities ( p out ) and set p out = 0 . 06 when varying the prop.
Revised: We set the subsampling proportion (prop) as 0.1 when varying the probability of an edge between communities (pout) and set pout=0.06 when varying the prop.

Original: Our method can get a more accurate estimation of M than other methods.
Revised: Our method can get a more accurate estimation of M than other methods.

Original: More details about the table of results, including each combination of pout and prop values, are presented in TableA.2 and A.1  in Appendix.
Revised: More details about the table of results, including each combination of pout and prop values, are presented in Table A.1 and A.2 in Appendix.

Original: We also compare the computation time for estimating the number of communities by using the whole graph and the subgraph.
Revised: We also compare the computation time for estimating M by using the whole graph and the subgraph.

Original: It turns out  the time for estimation of the subgraph together with the time for subsampling is still much shorter than the time for estimation on the whole graph.
Revised: It turns out that the time for estimating the subgraph together with the time for subsampling is still much shorter than the time for estimation on the whole graph.

Original: It is consistent with our conclusion about the complexity of the algorithm.
Revised: It is consistent with our conclusion about the complexity of the algorithm.

Original: Details about computational time are presented in Section B.1.2 in Appendix.
Revised: Details about computational time are presented in Section B.1.2 in Appendix.

Original: 5. R ESULTS OF R EAL - WORLD D ATASET
Revised: 5. R ESULTS OF R EAL - WORLD D ATASET

Original: We compare the performance of the proposed ORG-sub algorithm with subsamples generated by different algorithms by MAE defined above.
Revised: We compare the performance of the proposed ORG-sub algorithm with subsamples generated by different algorithms by MAE defined above.

Original: Table 1 records the average error of the estimation of the number of communities over 30 replications.
Revised: Table 1 records the average error of the estimation of M over 30 replications.

Original: Still, our algorithm outperforms other algorithms in most cases.
Revised: Still, our algorithm outperforms other algorithms in most cases.

Original: Let n 1 , . .
Revised: Let n 1 , . .

Original: . , n M denote the number of nodes in M communities.
Revised: . , n M denote the number of nodes in M communities.

Original: Here, we calculate the normalized Shannon Entropy of n 1 , . ..
Revised: Here, we calculate the normalized Shannon Entropy of n 1 , . .

Original:  , n M as a metric to measure the imbalance level, i.e., we calculate IM = − 1 log M (cid:80) Mi
Revised: . , n M as a metric to measure the imbalance level, i.e., we calculate IM = − 1 log M (cid:80) Mi

Original: =1 n i n log n i n and a higher IM value means that the communities are more imbalanced.
Revised: =1 n i n log n i n and a higher IM value means that the communities are more imbalanced.

Original: Table A.4 in the appendix summarizes the IM.
Revised: Table A.4 in the appendix summarizes the IM.

Original: As we can see,  PubMed and Polbooks areboth more imbalanced than Facebook.
Revised: As we can see, both PubMed and Polbooks are more imbalanced than Facebook.

Original: From Table 1, we observe that the ORG-sub performances of PubMed and Polbooks are better than that of Facebook.
Revised: From Table 1, we observe that the ORG-sub performances of PubMed and Polbooks are better than that of Facebook.

Original: This shows that ORG-sub performs better when the community of the graph is more imbalanced.
Revised: This shows that ORG-sub performs better when the community of the graph is more imbalanced.

Original: The true M of each dataset is recorded in the first column of Table 1.
Revised: The true M of each dataset is recorded in the first column of Table 1.

Original: We tested the performance of different algorithms with a broader range of sampling proportions ranging from 0.5% to 30%.
Revised: We tested the performance of different algorithms with a broader range of sampling proportions ranging from 0.5% to 30%.

Original: In particular, for small-size networks Polbooks and Facebook, we consider sampling proportions 10%, 20% and 30%.
Revised: In particular, for small-sized networks Polbooks and Facebook, we consider sampling proportions 10%, 20%, and 30%.

Original: For medium-size networks
Revised: For medium-size networks

Original: Polblogs and Cora, we consider 5%, 10%, and 20%.
Revised: Polblogs and Cora, we consider 5%, 10%, and 20%.

Original: For large-size network PubMed, we consider 0.5%, 2%, and 5%.
Revised: For large-size network PubMed, we consider 0.5%, 2%, and 5%.

Original: Similar to the synthetic dataset, we also compare the computation time for estimating the number of communities on the subgraph and the whole graph.

Original: The computation time on the subgraph (together with the subsampling time) is still much shorter than on the whole graph.
Revised: In addition, the computation time on the subgraph (together with the subsampling time) is still much shorter than on the whole graph.

Original: Details about computational time are presented in Section B.2.2 in Appendix.
Revised: Details about computational time are presented in Section B.2.2 in Appendix.

Original: [Table]
Revised: [Table]

Original: True:10% 0.

Original: [Table]

Original: Facebook

Original: True:

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: Cora

Original: True

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: Polblogs

Original: True:

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: PubMed

Original: True:0.

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original: [Table]

Original:  C ONCLUSION
Revised:  C ONCLUSION

Original: In this paper, we propose
Revised: In this paper, we propose

Original: a novel Ollivier-Ricci
Revised: a novel Ollivier-Ricci

Original: Curvature Gradient-based graph subsampling (ORGsub) method that samples a subgraph by maximizing OR curvature gradient.
Revised: Curvature Gradient-based graph subsampling (ORGsub) method that samples a subgraph by maximizing OR curvature gradient.

Original: The contribution of the ORG-sub method to the graph subsampling research line is three-fold.
Revised: The contribution of the ORG-sub method to the graph subsampling research line is three-fold.

Original: First, to the best of our knowledge, we are the first to utilize the graph’s internal topological information to subsample a large graph that preserves the number of communities.
Revised: First, to the best of our knowledge, we are the first to utilize the graph’s internal topological information to subsample a large graph that preserves the number of communities.

Original: Second, to the best of our knowledge, we are the first to bridge the gap in the consistency theory regarding  the number of community estimations.
Revised: Second, to the best of our knowledge, we are the first to bridge the gap in the consistency theory regarding subsampling algorithms for the number of community estimations in SBMs.

Original: In particular, we theoretically show that ORG-sub effectively traverse different communities and avoid trapping in one community.
Revised: In particular, we theoretically show that ORG-sub effectively traverses different communities and avoids trapping in one community.

Original: In addition, we theoretically show the advantage of ORG-sub over a popular sampling method, random walk.
Revised: In addition, we theoretically show the advantage of ORG-sub over a popular sampling method, random walk.

Original: Third, we empirically show that our method has superior performance over existing subsampling algorithms in terms of estimating the number of communities.
Revised: Third, we empirically show that our method has superior performance over existing subsampling algorithms in terms of estimating the number of communities.

Original: An interesting future direction we plan to investigate is to extend the theory of our method to SBM variants, including degree-correlated SBM, overlapping SBM, and multi-layer SBM.
Revised: An interesting future direction we plan to investigate is to extend the theory of our method to SBM variants, including degree-correlated SBM, overlapping SBM, and multi-layer SBM.

Original: Another future direction is to investigate other community-related statistics that our method has preserved.
Revised: Another future direction is to investigate other community-related statistics that our method has preserved.

Original: In fact, we have empirically observed some promising results in preserving the clustering coefficient (CC).
Revised: In fact, we have empirically observed some promising results in preserving the clustering coefficient (CC).

Original: Empirical results can be found in Section E in Appendix.
Revised: Empirical results can be found in Section E in Appendix.

Original: In the future, we will investigate the consistency theory of CC.
Revised: In the future, we will investigate the consistency theory of CC.

Revised: A CKNOWLEDGMENTS

Revised: Thanks to the partial support by NSF awards DMS-1903226, DMS-1925066, DMS- and

Revised: NIH grant R01GM1222080.

Revised: We would like to thank the referees for their valuable comments and suggestions that helped us to improve the quality of this paper.

Revised: We appreciate the time and effort the referees put into the review process and we are grateful for their contributions

Original: Appendix for “Subsampling for Large Graphs Using Ricci Curvature”
Revised: Appendix for “Subsampling for Large Graphs Using Ricci Curvature”

Original: The appendix shows the details of the proof, details of experiments, including the parameter settings of generating the synthetic datasets, a description of real-world datasets, and more experiment results on synthetic and real-world datasets.
Revised: The appendix shows the details of the proof, details of experiments, including the parameter settings of generating the synthetic datasets, a description of real-world datasets, and more experiment results on synthetic and real-world datasets.

Original: A P ROOF D ETAILS
Revised: A P ROOF D ETAILS

Original: Lemma A.1 (Ollivier-Ricci curvature of random graphs generated by stochastic block models) .
Revised: Lemma A.1 (Ollivier-Ricci curvature of random graphs generated by stochastic block models) .

Original: Consider a graph with M communities generated from a stochastic block model (SBM)
Revised: Consider a graph with M communities generated from a stochastic block model (SBM)

Original: B (cid:16) (cid:110) p ( i ) in (cid:111)
Revised: [Equation]

Original: M i =1 , p out , { n i } Mi =1 (cid:17) , where p ( i ) in represents the probability of an edge connecting nodes in the same community B i , p out represents the probability of an edge connecting nodes from different communities, and { n i } Mi =1 are the sizes of each community in the SBM.
Revised: i =1 , p out , { n i } Mi =1 (cid:17) , where p ( i ) in represents the probability of an edge connecting nodes in the same community B i , p out represents the probability of an edge connecting nodes from different communities, and { n i } Mi =1 are the sizes of each community in the SBM.

Original: We assume that, for any i = 1 , ..., M , p ( i ) in > p out , and n i n converges to r i (1 > r i > 0) , as n goes to infinity.
Revised: We assume that, for any i = 1 , ..., M , p ( i ) in > p out , and n i n converges to r i (1 > r i > 0) , as n goes to infinity.

Original: Let n ∨ = max { n i } Mi =1 , and p ∨ = max (cid:110) p ( i ) in (cid:111)
Revised: Let n ∨ = max { n i } M i =1 , and p ∨ = max (cid:110) p ( i ) in (cid:111)

Original: M i =1 .
Revised: M i =1 .

Original: We denote the node x from the community B i by x B i , the edge ( x B i , y B i ) where x is from B i and y is from B j , and the edge’s OR curvature by κ ( x B i , y B i ) .
Revised: We denote the node x from the community B i by x B i , the edge ( x B i , y B i ) where x is from B i and y is from B j , and the edge’s OR curvature by κ ( x B i , y B i ) .

Original: The following statements hold for the OR curvature.
Revised: The following statements hold for the OR curvature.

Original: If p out ≤ (cid:113) ( n i − 1)ln n n ln n i and p ( i ) in n i + p out ( n − n i ) > 83 ln n i , for the lower bound of within-community OR curvature, we almost surely have
Revised: If p out ≤ (cid:113) ( n i − 1)ln n n ln n i and p ( i ) in n i + p out ( n − n i ) > 83 ln n i , for the lower bound of within-community OR curvature, we almost surely have

Original: [Equation]
Revised: [Equation]

Original: If p ( i ) in 2 n i + p out 2 ( n − n i ) > 4 ln n i , for the upper bound of within-community OR curvature, we almost surely have
Revised: If p ( i ) in 2 n i + p out 2 ( n − n i ) > 4 ln n i , for the upper bound of within-community OR curvature, we almost surely have

Original: [Equation]
Revised: [Equation]

Original: If p out >
Revised: If p out >

Original: [Equation]
Revised: [Equation]

Original: 3 ln n ∨ )
Revised: 3 ln n ∨ )

Original: [Equation]
Revised: [Equation]

Original: p
Revised: p

Revised: [Equation]

Revised: [Equation]

Revised: p out

Original: [Equation]
Revised: [Equation]

Original: p out ( n i + n j ) + n − n i − n j
Revised: n i + n j ) + n − n i − n j

Original: [Equation]
Revised: [Equation]

Original: , for the upper bound of betweencommunities OR curvature, we almost surely have
Revised: , for the upper bound of betweencommunities OR curvature, we almost surely have

Original: [Equation]
Revised: [Equation]

Original: Lemma A.2.

Original: Combining the conclusion of statement 1 and statement 2 in Lemma A.1, almost surely, for any ( x ′ B i , y ′ B j ) ∈ ∆(( x B i , y B i )) , we have

Original: [Equation]

Original: [Equation]

Original: We will use Chernoff’s inequality to prove the lemmas
Revised: We will use Chernoff’s inequality to prove the lemmas

Original: [Equation]
Revised: [Equation]

Original: Lemma A.3.
Revised: Lemma A.2.

Original: (Chernoff’s inequality)
Revised: (Chernoff’s inequality)

Original: Let X 1 , ..., X n be the independent random variables with
Revised: Let X 1 , ..., X n be the independent random variables with

Original: [Equation]
Revised: [Equation]

Revised: ( X i

Revised: [Equation]

Revised: [Equation]

Original: [Equation]
Revised: [Equation]

Original: We consider the sum X = (cid:80) ni
Revised: We consider the sum X = (cid:80) ni

Original: =1 X i , with expectation E ( X ) =
Revised: =1 X i , with expectation E ( X ) =

Original: (cid:80) ni
Revised: (cid:80) ni

Original: =1 p i .
Revised: =1 p i .

Original: Then we have 
Revised: Then we have ( Lowertail [Equation] [Equation] ( Uppertail

Revised: [Equation]

Original: [Equation]
Revised: [Equation]

Original: WLOG, the proof of the theorem considers the situation when M = 2 .
Revised: WLOG, the proof of the theorem considers the situation when M = 2 .

Original: The conclusion can be easily generated to any arbitrary M by setting one block as block B 1 , and the other blocks as another block
Revised: The conclusion can be easily generated to any arbitrary M by setting one block as block B 1 , and the other blocks as another block

Original: B 2 .
Revised: B 2 .

Original: Before we prove the theorem, we need afew lemmas.
Revised: Before we prove the theorem, we need a

Original: Lemma A.4.
Revised: Lemma A.3.

Original: If p ( i ) in n i + p out ( n − n i ) > 83 ln n i , then the degree of all nodes in block B i fall in the range (cid:16) p ( i ) in ( n i − 1) + p out ( n − n i )) − (cid:113) 4( p ( i ) in n i + p out ( n − n i ))
Revised: If p ( i ) in n i + p out ( n − n i ) > 83 ln n i , then the degree of all nodes in block B i fall in the range (cid:16) p ( i ) in ( n i − 1) + p out ( n − n i )) − (cid:113) 4( p ( i ) in n i + p out ( n − n i ))

Original: ln n 1 , p ( i ) in ( n i − 1) + p out  n i ) + (cid:113) 6( p ( i ) in n i + p out ( n − n i ))
Revised: ln n i , p ( i ) in ( n i − 1) + p out ( n − n i ) + (cid:113) 6( p ( i ) in n i + p out ( n − n i ))

Original: ln n i (cid:17) with probability at least 1 − 2 /n
Revised: ln n i (cid:17) with probability at least 1 − 2 /n

Original: i .
Revised: i .

Original: Proof.
Revised: Proof.

Original: Without loss of generality, we just prove the first statement and the second statement can be proved in the same way.
Revised: Without loss of generality, we just prove the first statement and the second statement can be proved in the same way.

Original: For each vertex v in block B i , it is easy to show that the expected degree of vertex v is
Revised: For each vertex v in block B i , it is easy to show that the expected degree of vertex v is

Original: E ( d B i v ) =
Revised: E ( d B i v ) =

Original: p ( i ) in ( n i − 1) + p out ( n − n i ) .
Revised: p ( i ) in ( n i − 1) + p out ( n − n i ) .

Original: Applying Chernoff’s inequality with the lower tail λ
Revised: Applying Chernoff’s inequality with the lower tail λ

Original: = (cid:113) 4( p ( i ) in n i + p out ( n − n i ))
Revised: = (cid:113) 4( p ( i ) in n i + p out ( n − n i ))

Original: ln n i , we have
Revised: ln n i , we have

Revised: [Equation]

Original: [Equation]
Revised: [Equation]

Original: Applying Chernoff’s inequality with the upper tail λ
Revised: Applying Chernoff’s inequality with the upper tail λ

Original: = (cid:113) 6( p ( i ) in n i + p out ( n − n i ))
Revised: = (cid:113) 6( p ( i ) in n i + p out ( n − n i ))

Original: ln n i , we have
Revised: ln n i , we have

Original: P

Original: [Equation]
Revised: [Equation]

Original: In the last step, we used the assumption p ( i ) in( n − n i) + p out  n i  > 83 ln( n − n i) .
Revised: In the last step, we used the assumption p ( i ) in n i + p out ( n − n i ) > 83 ln n i .

Original: The probability that there is a vertex v in block B i so that d B i v
Revised: The probability that there is a vertex v in block B i so that d B i v

Original: [Equation]
Revised: [Equation]

Original: is at most n 
Revised: is at most n i (cid:16) 1 n 2 i + 1 n 2 i (cid:17)

Original: [Equation]
Revised: = 2 n i .

Original: The co-degree d xy of a pair of vertices ( x, y ) is the cardinality of the common neighborhood of x and y .
Revised: The co-degree d xy of a pair of vertices ( x, y ) is the cardinality of the common neighborhood of x and y .

Original: Roughly speaking, when p ( i ) in and p out is large, the co-degree follows the binomial distribution; when p ( i ) in and p out is small, the co-degree follows the Poisson distribution.
Revised: Roughly speaking, when p ( i ) in and p out is large, the co-degree follows the binomial distribution; when p ( i ) in and p out is small, the co-degree follows the Poisson distribution.

Original: We can expect that all co-degrees within a block and between blocks are concentrated around a small interval.
Revised: We can expect that all co-degrees within a block and between blocks are concentrated around a small interval.

Original: Lemma A.5.
Revised: Lemma A.4.

Original: If p ( i ) in 2 n i + p out 2 ( n − n i ) > 4 ln n i , then with probability at least 1 − 1 /n i , all co-degrees of a pair of vertices in block B i fall in the range (cid:16) p ( i ) in 2 ( n i − 2)+ p out
Revised: If p ( i ) in 2 n i + p out 2 ( n − n i ) > 4 ln n i , then with probability at least 1 − 1 /n i , all co-degrees of a pair of vertices in block B i fall in the range (cid:16) p ( i ) in 2 ( n i − 2)

Original: [Equation]
Revised: + p out 2 ( n − n i ) − (cid:113) 6( p ( i ) in 2 n i + p out 2 ( n − n i ))

Original:  n [Equation] n i [Equation]
Revised: ln n i , p ( i ) in 2 ( n i − 2)

Original:  (cid:113) 
Revised: + p out 2 ( n − n i ) + (cid:113) 9( p ( i ) in 2 n i + p out 2 ( n − n i ))

Original: [Equation]
Revised: ln n i (cid:17) .

Original: If p ( i ) in p out n i + p out p ( i ) in n j + p out 2 ( n − n i − n j ) >
Revised: If p ( i ) in p out n i + p out p ( i ) in n j + p out 2 ( n − n i − n j ) >

Original: 3 ln n ∨ , then with probability at least 1 − 1 /n ∨ , all co-degrees of a pair of vertices in blocks B i and B j fall in the range (cid:16) p ( i ) in p out ( n i − 1)+ p out 2 ( n − n j − 
Revised: 3 ln n ∨ , then with probability at least 1 − 1 /n ∨ , all co-degrees of a pair of vertices in blocks B i and B j fall in the range (cid:16) p ( i ) in p out ( n i − 1)+ p out 2 ( n − n j − n i )

Original: [Equation]
Revised: + p ( i ) in p out ( n j − 1) − (cid:113) 6( p ( i ) in p out n i + p out p ( i ) in n j + p out 2 ( n − n i − n j )) ln n ∨ , p ( i ) in p out ( n i − 1)+ p out 2 ( n − n j − n i )+ p ( i ) in p out ( n j − 1)+ (cid:113) 9( p ( i ) in p out n i + p out p ( i ) in n j + p out 2 ( n − n i − n j )) ln n ∨ (cid:17) .

Original: Proof.
Revised: Proof.

Original: For a pair of vertices x ∈ B i and y ∈ B j , the co-degree d x Bi y Bj is the sum of n − 2 independent random variables X 1 , ..., X n − 2 with expectation E ( d x Bi y Bi )
Revised: For a pair of vertices x ∈ B i and y ∈ B i , the co-degree d x Bi y Bi is the sum of n − 2 independent random variables X 1 , ..., X n − 2 with expectation E ( d x Bi y Bi )

Original: = p ( i ) in 2 ( n i − 2)
Revised: = p ( i ) in 2 ( n i − 2)

Original: + p out 2 ( n − n j ) .
Revised: + p out 2 ( n − n i ) .

Original: Applying Chernoff’s inequality with the lower tail λ 
Revised: Applying Chernoff’s inequality with the lower tail λ =

Original: = (cid:113) 6( p ( i ) in 2 n i + p out 2 ( n − n j ))
Revised: [Equation] [Equation] p [Equation]

Original: ln n i , we have
Revised: i [Equation] in

Revised: [Equation]

Revised: n i + p

Original: [Equation]
Revised: [Equation]

Original:  If p ( i ) in 2 n i + p out 2 ( n − n i ) > 4 ln n i , we apply Chernoff’s inequality with the upper tail λ
Revised: ) ln n i , we have [Equation] If p ( i ) in 2 n i + p out [Equation]

Original: = (cid:113) 9( p ( i ) in 2 n i + p out 2 ( n − n 
Revised: n − n [Equation] > 4 ln n i [Equation]

Revised: we apply Chernoff’s inequality with the upper tail λ =

Original: [Equation]
Revised: [Equation]

Original: Now the number of pairs is at most n 2 i / 2 .
Revised: Now the number of pairs is at most n 2 i / 2 .

Original: the sum of the probabilities of small events is at most n 
Revised: the sum of the probabilities of small events is at most n 2 i 2 (cid:16) 1 n [Equation] +

Original: [Equation]
Revised: [Equation]

Original: i
Revised: n

Original: [Equation]
Revised: [Equation]

Revised: (cid:17) = 1 n i .

Original: For a pair of vertices x ∈ B i and y ∈ B j , the co-degree d x Bi y
Revised: For a pair of vertices x ∈ B i and y ∈ B j , the co-degree d x Bi y

Original: Bj is the sum of n − 2 independent random variables X 1 , ..., X n − 2 with expectation E ( d x Bi y
Revised: Bj is the sum of n − 2 independent random variables X 1 , ..., X n − 2 with expectation E ( d x Bi y

Original: Bj )
Revised: Bj )

Original: = p ( i ) in p out ( n i − 1) + p out 2 ( n − n j − n i )
Revised: = p ( i ) in p out ( n i − 1) + p out 2 ( n − n j − n i )

Original: + p ( i ) in p out ( n j − 1) .
Revised: + p ( i ) in p out ( n j − 1) .

Original: Applying Chernoff’s inequality with the lower tail λ
Revised: Applying Chernoff’s inequality with the lower tail λ

Original: = (cid:113) 6( p ( i ) in p out n i + p out p ( i ) in n j + p out 2 ( n − n i − n j ))
Revised: = (cid:113) 6( p ( i ) in p out n i + p out p ( i ) in n j + p out 2 ( n − n i − n j ))

Original: ln n ∨ , we have
Revised: ln n ∨ , we have

Original: [Equation]
Revised: [Equation]

Revised: [Equation]

Revised: [Equation]

Original: [Equation]
Revised: [Equation]

Original: If p ( i ) in p out n i + p out p ( i ) in n j + p out 2 ( n − n i − n j ) >
Revised: If p ( i ) in p out n i + p out p ( i ) in n j + p out 2 ( n − n i − n j ) >

Original: 3 ln n ∨ , we apply Chernoff’s inequality with the upper tail λ
Revised: 4 ln n ∨ , we apply Chernoff’s inequality with the upper tail λ

Original: = (cid:113)9( p ( i ) in p out n i + p out p [Equation] in n
Revised: = (cid:113)

Original: [Equation]
Revised: [Equation]

Original: p
Revised: p

Original: [Equation]
Revised: [Equation]

Original: ln n ∨

Original: [Equation]

Original: [Equation]

Original: [Equation]

Original: 33) <

Original: n
Revised: n

Original: [Equation]
Revised: [Equation]

Original: The number of pairs is at most n ∨ 2 .
Revised: The number of pairs is at most n ∨ 2 .

Original: the sum of the probabilities of small events is at most ( n ∨ 2 ) (cid:16) 1 n ∨ 3 (cid:17)
Revised: the sum of the probabilities of small events is at most [Equation]

Original: Referring to the lemma mentioned in Lin et al.
Revised: Referring to the lemma mentioned in Lin et al.

Original: (2011), we can construct the lower and upper bound of the Ricci curvature of graphs generated by the stochastic block model.
Revised: (2011), we can construct the lower and upper bound of the Ricci curvature of graphs generated by the stochastic block model.

Original: Lemma A.6.
Revised: Lemma A.5.

Original: Suppose that ϕ : Γ( x ) \ N ( y ) → Γ( y ) \ N ( x ) is an injective mapping.
Revised: Suppose that ϕ : Γ( x ) \ N ( y ) → Γ( y ) \ N ( x ) is an injective mapping.

Original: Then we have
Revised: Then we have

Revised: [Equation]

Original: [Equation]
Revised: [Equation]

Original: Lemma A.7.
Revised: Lemma A.6.

Original: For the lower bound, we will construct a matching M from Γ( x ) \ N ( y ) to Γ( y ) \ N ( x ) as follows.
Revised: For the lower bound, we will construct a matching M from Γ( x ) \ N ( y ) to Γ( y ) \ N ( x ) as follows.

Original: Let U 

Revised: Let U

Original: = Γ( x ) \ N ( y ) and V 0 =
Revised: = Γ( x ) \ N ( y ) and V 0 =

Original: Γ( y ) \ N ( x ) .
Revised: Γ( y ) \ N ( x ) .

Original: Pick up a vertex u 1 ∈ U 0 .
Revised: Pick up a vertex u 1 ∈ U 0 .

Original: Reveal the neighborhood of u 1 in V 0 .
Revised: Reveal the neighborhood of u 1 in V 0 .

Original: Pick a vertex in the neighborhood, and denote it by v 1 .
Revised: Pick a vertex in the neighborhood, and denote it by v 1 .

Original: Let U 

Revised: Let U

Original: = U 0 \{ u 1 } and V 1 = V 0 \{ v 1 } and continue this process.
Revised: = U 0 \{ u 1 } and V 1 = V 0 \{ v 1 } and continue this process.

Original: The process ends when Γ( u i +1 )
Revised: The process ends when Γ( u i +1 )

Original: ∩ V i = ∅ .
Revised: ∩ V i = ∅ .

Original: The probability that the maximum matching between U 0 and V 0 is at most k is less thank (cid:88)
Revised: The probability that the maximum matching between U 0 and V 0 is at most k is less than

Original: [Equation]
Revised: [Equation]

Original: Choose k
Revised: Choose k

Original: = ⌊| V 0 | − (3 ln n ) /p ⌋ .
Revised: = ⌊| V 0 | − (3 ln n ) /p ⌋ .

Original: With probability at least 1 − 1 /n 2 , there is a matching M of size k between Γ( x ) \ N ( y ) and Γ( y ) \ N ( x ) .
Revised: With probability at least 1 − 1 /n 2 , there is a matching M of size k between Γ( x ) \ N ( y ) and Γ( y ) \ N ( x ) .

Original: Now we extend the matching M to an injective mapping ϕ : Γ( x ) \ N ( y )
Revised: Now we extend the matching M to an injective mapping ϕ : Γ( x ) \ N ( y )

Original: → Γ( y ) \ N ( x ) arbitrarily.
Revised: → Γ( y ) \ N ( x ) arbitrarily.

Original: Applying Lemmas A.7, with probability at least 1 − 4 /n , we have
Revised: Applying Lemmas A.6, with probability at least 1 − 4 /n , we have

Original: [Equation]
Revised: [Equation]

Original: Proof.
Revised: Proof.

Original: Given the range of the degrees and co-degrees, lower bound and upper bound of Ricci curvatures with respect to d x and d xy , we can prove the theorem ?? as follows.
Revised: Given the range of the degrees and co-degrees, lower bound and upper bound of Ricci curvatures with respect to d x and d xy , we can prove the remark 3 of Lemma A.1 as follows.

Original: For vertex x ∈ B i , y ∈ B i , given p ( i ) in n i + p out ( n − n i ) > 83 ln n i , we have
Revised: For vertex x ∈ B i , y ∈ B i , given p ( i ) in n i + p out ( n − n i ) > 83 ln n i , we have

Original: [Equation]
Revised: [Equation]

Original: For vertex x ∈ B i , y ∈ B i , given p ( i ) in 2 n i + p out 2 ( n − n i ) > 4 ln n i , we have
Revised: For vertex x ∈ B i , y ∈ B i , given p ( i ) in 2 n i + p out 2 ( n − n i ) > 4 ln n i , we have

Original: [Equation]
Revised: [Equation]

Original: If p[Equation] 2  p ( i ) in ≤
Revised: ( ln n p 2 out n ) . For vertex x ∈ B i , y ∈ B j , let n ∨ = max { n i } Mi =1 , and p ∨ = max (cid:110) p ( i ) in (cid:111)

Original: [Equation]
Revised: M i =1 .

Original:  ln n  ) 2 n ln n i  , we haveκ ( x B i , y B i [Equation] = ( n
Revised: If p out > (cid:115) (3 ln n ∨ ) (cid:46) (cid:20) p ( i ) in p out ( n i + n j ) + n − n i − n j (cid:21) , we have

Original: [Equation]
Revised: [Equation]

Original: O ( ln n p 2 out n ) .
Revised: Lemma A.7.

Original: For vertex x ∈ B i , y ∈ B j , let n ∨ = max { n i} Mi =1 , and p ∨ = max (cid:110) p ( i ) in (cid:111) M i =1 . If p out > (cid:115) (3 ln n
Revised: Combining the conclusion of statement 1 and statement 2 in Lemma A.1, almost surely, for any ( x ′ B i , y ′ B j ) ∈ ∆(( x B i , y B i )) , we have

Original: [Equation]
Revised: [Equation]

Original: ) (cid:46) (cid:20) p i) in p out ( n i + n j ) + n − n i − n j (cid:21) , we have
Revised: i

Original: [Equation]
Revised: [Equation]

Revised: B

Original: [Equation]
Revised: [Equation]

Original: Proof.
Revised: Proof.

Original: We can prove the Lemma A.2 by the lower bound and upper bound of the within-community edges’ OR curvature from Lemma A.1.
Revised: We can prove the Lemma A.7 by the lower bound and upper bound of the within-community edges’ OR curvature from Lemma A.1.

Original: Theorem A.8.
Revised: Theorem A.8.

Original: Combining the conclusion of statement 1 and statement 3 in Lemma A.1, if p ∨ < p ( i ) in n n i + n j , p ( i ) in < p out ( n − n i )( n + n i + n j ) n 2 − ( n i + n j )( n − 2 n i ) and conditions given in Lemma A.1, almost surely, for any ( x ∗ B i , y ∗ B j ) ∈ ∆(( x B i , y B i )) , we have
Revised: Combining the conclusion of statement 1 and statement 3 in Lemma A.1, if p ∨ < p [Equation] in n n [Equation]

Revised: , p

Original: [Equation]
Revised: [Equation]

Revised: i ) in <

Original: [Equation]
Revised: [Equation]

Revised: i ) and conditions given in Lemma A.1, almost surely, for any ( x ∗ B i ,

Revised: y ∗ B j ) ∈ ∆(( x B i , y B i )) , we have

Revised: [Equation]

Original: Proof.
Revised: Proof.

Original: We first get the lower bound of the OR curvature difference of within community edge κ ( x B i , y B i ) and between community edge κ ( x ∗ B i , y ∗ B j ) , which can be obtained by the lower bound of κ ( x B i , y B i ) and upper bound of κ ( x ∗ B
Revised: We first get the lower bound of the OR curvature difference of within community edge κ ( x B i , y B i ) and between community edge κ ( x ∗ B i , y ∗ B j ) , which can be obtained by the lower bound of κ ( x B i , y B i ) and upper bound of κ ( x ∗ B

Original: i , y ∗ B j ) from the Lemma A.1.
Revised: i , y ∗ B j ) from the Lemma A.1.

Original: Given the upper bound of the OR curvature difference of edges in the same community [Equation] we can prove under certain conditions [Equation]
Revised: Given the upper bound of the OR curvature difference of edges in the same community, we can prove under certain conditions, the upper bound is smaller than the lower bound we obtained in the last step.

Original: the upper bound is smaller than the lower bound we obtained in the last step.

Original: Denote r i = n i /n and r j = n j /n .
Revised: Denote r i = n i /n and r j = n j /n .

Original: [Equation]
Revised: [Equation]

Original: Theorem A.9 (Probability of the ORG-sub’s Expansion to Communities) .
Revised: Theorem A.9 (Probability of the ORG-sub’s Expansion to Communities) .

Original: Given the graph G generated by Theorem A.8, we apply the ORG-sub Algorithm 1 to the graph G and take graph subsamples with size ˜ n denoted by G [ S ] .
Revised: Given the graph G generated by Theorem A.8, we apply the ORG-sub Algorithm 1 to the graph G and take graph subsamples with size ˜ n denoted by G [ S ] .

Original: Under the conditions in Theorem A.8, and if ˜ n/n → 1 , we can prove that, the probability of the final subgraph G [ S ] taken by ORG-sub containing nodes in the community B i is:
Revised: Under the conditions in Theorem A. [Equation] and if ˜ n/n → 1 , we can prove that, the probability of the final subgraph G [ S ] taken by ORG-sub containing nodes in the community B i is:

Original: [Equation]
Revised: [Equation]

Original: Proof.
Revised: Proof.

Original: Without loss of generality, we show the theorem holds for a two-block stochastic block model (SBM).
Revised: Without loss of generality, we show the theorem holds for a two-block stochastic block model (SBM).

Original: We assume the graph is generated from a two-block stochastic block model (SBM), where the probability of an edge within a block is p in , and the probability of an edge between the blocks is p out .
Revised: We assume the graph is generated from a two-block stochastic block model (SBM), where the probability of an edge within a block is p in , and the probability of an edge between the blocks is p out .

Original: We denote the subsampled graph from the original graph G by G [ S ] , where S is a subset of the full vertices set V of graph G , and G [ S ] is the subgraph induced by the subset S .
Revised: We denote the subsampled graph from the original graph G by G [ S ] , where S is a subset of the full vertices set V of graph G , and G [ S ] is the subgraph induced by the subset S .

Original: Denote the probability matrix of generating the adjacency matrix A by P :
Revised: Denote the probability matrix of generating the adjacency matrix A by P :

Original: [Equation]

Original: [Equation]
Revised: [Equation]

Original: where the size of the first community is n 1 , the size of the second community is n 2 , and the number of nodes of the graph is n = n 1 + n 2 .
Revised: where the size of the first community is n 1 , the size of the second community is n 2 , and the number of nodes of the graph is n = n 1 + n 2 .

Original: The node set belonging to the community 1 is B 1 = { v 1 , ..., v n 1 } , and the node set belonging to the community 2 is B 2 = { v n 1 +1 ,...,n } .
Revised: The node set belonging to the community 1 is B 1 = { v 1 , ..., v n 1 } , and the node set belonging to the community 2 is B 2 = { v n 1 +1 ,...,n } .

Original: Denote each entry of the probability matrix P by p ij , then the observed A is generated from:
Revised: Denote each entry of the probability matrix P by p ij , then the observed A is generated from:

Original: [Equation]
Revised: [Equation]

Original: [Equation]

Original: Define the block of the matrix A as A o 1 o 2 = { e ij =
Revised: Define the block of the matrix A as A o 1 o 2 = { e ij =

Original: ( v i , v j ) | v i ∈ B o 1 , v j ∈ B o 2 } , where o 1 , o 2 ∈ { 1 , 2 } .
Revised: ( v i , v j ) | v i ∈ B o 1 , v j ∈ B o 2 } , where o 1 , o 2 ∈ { 1 , 2 } .

Original: We denote the edge subsampled at the l -th step by e ( l ) ij , and the neighboring edge set of edge e ij by ∆( e ij ) .
Revised: We denote the edge subsampled at the l -th step by e ( l ) ij , and the neighboring edge set of edge e ij by ∆( e ij ) .

Original: For step 1, given the  edgesubsampled in the initialization step e (0) ij , we have:
Revised: For each step, the probability of subsampling an edge in a block is a constant.

Original: [Equation]
Revised: After subsampling one edge in each step, the nodes connected by the subsampled edges are subsampled into node set S .

Original: Iteratively, we get the probability of subsampling edge e ij at the  l -th step: [Equation] ∈ ∆( e ( l ) ij ) ∩
Revised: We denote the probability of adding node ˜ v at the step l belonging to community B o , where o ∈ { 1 , 2 } , by P (˜ v ( l ) ∈ B o ) .

Original: ( A 12 ∪ A 21 ) P ( e ( l ) ij ) P ( e ( l − 1) ij ) , [Equation] By the above equations, for each step
Revised: After ˜ n steps, the probability of the obtained subsample G [ S ] containing node ˜ v from the node set B 1 is:

Original: [Equation]
Revised: [Equation]

Original: the probability of subsampling an edge in a block is a constant. After subsampling one edge in each step
Revised: Analogously

Original: [Equation]
Revised: [Equation]

Original: thenodes connected by the subsampled edges are subsampled into node set S . We denote the probability of adding node ˜ vat the step l belonging tocommunity B o , where o ∈ { 1 , 2 } , by P (˜ v ( l ) ∈ B o ) . After ˜ n steps, the probability of the obtained subsample G [ S ] containing node ˜ v from the node set B 1 is:
Revised: the probability of subsampling node ˜ v belonging to the node set B 2 after n steps is:

Original: [Equation]
Revised: [Equation]

Original: ∈ A 22 ). Analogously, the probability of subsampling node ˜ v belonging to the node set B 2 after n steps is:
Revised: Since (1 − p out ) n 1 (1 − (1 − p in ) n 2 −

Original: [Equation]
Revised: [Equation]

Original: Since (1 − p out ) n 1(1 − (1 − p in ) n 2 − 2 ) ( (1 − p out ) n 2 (1 − (1 − p in ) n 1 − 2 ) ) can be small when n( n 2 ) and ˜ n is large, the proposed ORG-sub subsampler can traverse to new communities with large probability.
Revised: 1 − p out ) n 2 (1 − (1 − p in ) n 1 − 2 ) ) can be small when n( n 2 ) and ˜ n is large, the proposed ORG-sub subsampler can traverse to new communities with high probability.

Original: When p out is small, that is, the edges density between communities is relatively small, the proposed ORG-sub needs more steps to expand to another community.
Revised: When p out is small, that is, the edges density between communities is relatively small, the proposed ORG-sub needs more steps to expand to another community.

Original: As the traveling steps of the proposed subsampler ORG-sub increase, we can prove that the probability of the proposed
Revised: As the traveling steps of the proposed subsampler ORG-sub increase, we can prove that the probability of the proposed

Original: ORG-sub algorithm taking all communities into the final subgraph converges to one under the two-stochastic-block model.
Revised: ORG-sub algorithm taking all communities into the final subgraph converges to one under the two-stochastic-block model.

Original: In general cases, the graph is generated by SBM with multiple communities, we can always think of the community that the subsampler has never been to as one block, and other communities as the other block.
Revised: In general cases, the graph is generated by SBM with multiple communities, we can always think of the community that the subsampler has never been to as one block, and other communities as the other block.

Original: Given the conclusion under the two-stochastic-block model, the ORG-sub can gradually expand the subsampled graph to the block we have never been to.
Revised: Given the conclusion under the two-stochastic-block model, the ORG-sub can gradually expand the subsampled graph to the block we have never been to.

Original: Thus, we can prove the ORG-sub can expand to any community under a multiple-stochastic-block model.
Revised: Thus, we can prove the ORG-sub can expand to any community under a multiple-stochastic-block model.

Original: Corollary A.9.1. (Theoretical advantage over random walk).
Revised: Corollary A.9.1 (Theoretical advantage over random walk) .

Original: Under the assumptions in Theorem 4.1,  we have[Equation]
Revised: Under the assumptions in Theorem 4.1, for any community B i in graph G , we have

Original: [Equation]
Revised: [Equation]

Original: Proof.
Revised: Proof.

Original: For the proposed method, Theorem A.9 has proved that after walking ˜ n steps, the probability of the obtained subsample G [ S ] by the proposed ORG-sub containing vertexes v from community
Revised: For the proposed method, Theorem A.9 has proved that after walking ˜ n steps, the probability of the obtained subsample G [ S ] by the proposed ORG-sub containing vertexes v from community

Original: B i is:
Revised: B i is:

Original: [Equation]
Revised: [Equation]

Original: After ˜ n steps, the probability of the obtained subsample G RW [ S ] by the random walk subsampler containing vertexes in community B i is
Revised: After ˜ n steps [Equation] the probability of the obtained subsample G RW [ S ] by the random walk subsampler containing vertexes in community B i is

Original: [Equation]
Revised: [Equation]

Original: Corollary A.9.2.
Revised: Corollary A.9.2.

Original: If the subsample with size ˜ n satisfies ˜ nρ log ˜ n → ∞ , where ρ controls the expected node degree of subsampled graph, and as the subsampling times r → ∞ , the probability the subsampled graph underestimates the true M is:
Revised: If the subsample with size ˜ n satisfies ˜ nρ log ˜ n → ∞ , where ρ is the edge density, as the subsampling times r → ∞ , the probability the subsampled graph underestimates the true M is:

Original: [Equation]
Revised: [Equation]

Original: Proof.
Revised: Proof.

Original: From Theorem A.9, we know that the true number of communities M ( G [ S ]) of subgraph G [ S ] is the same as the true number of communities M of the original graph G .
Revised: From Theorem A.9, we know that, if we take enough steps for subsampling, the true number of communities M ( G [ S ]) of subgraph G [ S ] is the same as the true number of communities M of the original graph G .

Original: Theorem 3 of Li et al.
Revised: Theorem 3 of Li et al.

Original: (2020) proves the one-sided consistency of estimating the number of communities of the stochastic block model under the following assumptions.
Revised: (2020) proves the one-sided consistency of estimating the number of communities of the stochastic block model under the following assumptions.

Original: First, the expected node degree λ n = nρ has the order  n log n , where n is the number of nodes of the graph G .
Revised: First, the expected node degree λ n = nρ has the order λ n log n , where n is the number of nodes of the graph G .

Original: Second, there exists a constant γ such that min k n k ≥ γN , where n k is the size of the k -th community.
Revised: Second, there exists a constant γ such that min k n k ≥ γN , where n k is the size of the k -th community.

Original: Since our subsampled graph satisfies ˜ nρ log ˜ n  , we can easily verify that the subsampled subgraph satisfies the assumptions in Theorem 3 of Li et al.
Revised: Since our subsampled graph satisfies ˜ nρ log ˜ n → ∞ , we can easily verify that the subsampled subgraph satisfies the assumptions in Theorem 3 of Li et al.

Original: The estimate of the number of communities ˆ M  i  corresponding to S i , the i -th subsample, satisfies:
Revised: The estimate of the number of communities (cid:99) M ( G [ S i ]) corresponding to S i , the i -th subsample, satisfies:

Original: [Equation]
Revised: [Equation]

Original: As the replications of the subsamples increase, it’s trivial to show:
Revised: As the replications of the subsamples increase, it’s trivial to show:

Original: [Equation]
Revised: [Equation]

Original: B D ETAILS OF E XPERIMENTS
Revised: B D ETAILS OF E XPERIMENTS

Original: We evaluate the performance of our algorithm on both synthetic datasets and real-world datasets using the metrics presented in the main context.
Revised: We evaluate the performance of our algorithm on both synthetic datasets and real-world datasets using the metrics presented in the main context.

Original: We also compare the proposed method with seven benchmark exploration-based graph subsampling methods, including Metropolis Hasting Random
Revised: We also compare the proposed method with seven benchmark exploration-based graph subsampling methods, including Metropolis Hasting Random

Original: Walk Sampler (MHRW), Forest Fire Sampler (FFS), Snowball Sampler, Community Structure
Revised: Walk Sampler (MHRW), Forest Fire Sampler (FFS), Snowball Sampler, Community Structure

Original: Expansion Sampler (CSE), degree-based node sampler (DBN), Random Walk Sampler (RW) and
Revised: Expansion Sampler (CSE), degree-based node sampler (DBN), Random Walk Sampler (RW) and

Original: Multi-dimensional Random Walk Sampler (MDRW).
Revised: Multi-dimensional Random Walk Sampler (MDRW).

Original: We set the hyper-parameters of the seven methods as the default in the package Little Ball of Fur .
Revised: We set the hyper-parameters of the seven methods as the default in the package Little Ball of Fur .

Original: The rejection constraint of the MHRW method is set as 1; the burning probability of FFS is 0.4; the bound on the degree of Snowball Sampler is set as 50.
Revised: The rejection constraint of the MHRW method is set as 1; the burning probability of FFS is 0.4; the bound on the degree of Snowball Sampler is set as 50.

Original: The hyper-parameter used to calculate the Ricci curvature of graphs is set as α = 0 .
Revised: The hyper-parameter used to calculate the Ricci curvature of graphs is set as α = 0 .

Original: The times of subsampling r in Algorithm 2 for estimation of the number of communities is set as 3, which is enough to get a stable estimation result.
Revised: The times of subsampling r in Algorithm 2 for estimation of the number of communities is set as 3, which is enough to get a stable estimation result.

Original: All the experiments are conducted on a workstation with a 40-core NVIDIA Tesla V100 GPU (3.00 GHz).
Revised: All the experiments are conducted on a workstation with a 40-core NVIDIA Tesla V100 GPU (3.00 GHz).

Original: B. E XPERIMENTS ON S YNTHETIC D ATASETS B.1. P ARAMETERS S ETTINGS
Revised: B. E XPERIMENTS ON S YNTHETIC D ATASETS B.1. P ARAMETERS S ETTINGS

Original: We use stochastic block models (SBM), and degree corrected block models (DCBM), which can assign the community distribution to each node.
Revised: We use stochastic block models (SBM), and degree corrected block models (DCBM), which can assign the community distribution to each node.

Original: To create graphs with unbalanced communities, we set the community proportion as (3 / 4 , 1 / 10 , 1 / 12 , 1 / 15) with 900 nodes.
Revised: To create graphs with unbalanced communities, we set the community proportion as (3 / 4 , 1 / 10 , 1 / 12 , 1 / 15) with 900 nodes.

Original: The out-in-ratio controls the ratio of an edge’s probability of between-communities and within-community.
Revised: The out-in-ratio controls the ratio of an edge’s probability of between-communities and within-community.

Original: A higher out-inratio represents a noisier graph, i.e., harder to distinguish the community from the graph.
Revised: A higher out-inratio represents a noisier graph, i.e., harder to distinguish the community from the graph.

Original: We set the probability of an edge within a block as 0.8 and vary the probability of an edge between blocks (pout) for both SBM and DCBM ( { 0 . 06 , 0 .
Revised: We set the probability of an edge within a block as 0.8 and vary the probability of an edge between blocks (pout) for both SBM and DCBM ( { 0 . 06 , 0 .

Original: The higher the pout value (higher the out-in-ratio) is, the noisier the subsampled graph is.
Revised: The higher the pout value (higher the out-in-ratio) is, the noisier the subsampled graph is.

Original: We also change the subsampling proportion ( { 0 . 1 , 0 .
Revised: We also change the subsampling proportion ( { 0 . 1 , 0 .

Original: 14 , 0 . 16 } ) from low to high for both block models to observe how our method performs compared with others as the settings change.
Revised: 14 , 0 . 16 } ) from low to high for both block models to observe how our method performs compared with others as the settings change.

Original: The degree corrected model corrects the node degree using a power-law distribution.
Revised: The degree corrected model corrects the node degree using a power-law distribution.

Original: More parameters need to be set before generating the graphs following
Revised: More parameters need to be set before generating the graphs following

Original: DCBM.
Revised: DCBM.

Original: We set the average node degree as 40 , with is consistent with the setting in Li et al.
Revised: We set the average node degree as 40 , with is consistent with the setting in Li et al.

Original: (2020), in which the assumptions in our algorithm hold.
Revised: (2020), in which the assumptions in our algorithm hold.

Original: The node degree follows a power-law distribution with the lower bound as 1 and the scaling parameter as 5.
Revised: The node degree follows a power-law distribution with the lower bound as 1 and the scaling parameter as 5.

Original: B.1. A DDITIONAL R ESULTS OF S YNTHETIC D ATASETS
Revised: B.1. A DDITIONAL R ESULTS OF S YNTHETIC D ATASETS

Original: [Table]
Revised: [Table]

Original: [Table] 
Revised: [Table] [Table] [Table] [Table]

Original: We use synthetic datasets with hidden community structures to evaluate the performance of our method.
Revised: [Table] [Table] [Table] [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Revised: [Table]

Original:  We replicate the experiments 30 times under each setting and compare the performance of the error of estimation of the number of communities and the computation time.
Revised: [Table] [Table] We use synthetic datasets with hidden community structures to evaluate the performance of our method. We replicate the experiments 30 times under each setting and compare the performance of the error of estimation of the number of communities and the computation time.

Original: The average of 30 replications of each setting is recorded in Table A.1 and Table A.2.
Revised: The average of 30 replications of each setting is recorded in Table A.1 and Table A.2.

Original: In these tables, Column
Revised: In these tables, Column

Original: Prop denotes the subsampling proportion, and column Pout denotes the probability of an edge between communities.
Revised: Prop denotes the subsampling proportion, and column Pout denotes the probability of an edge between communities.

Original: The rest of the columns are the results of different subsampling methods.
Revised: The rest of the columns are the results of different subsampling methods.

Original: The performance of our algorithm is better than that of other methods.
Revised: The performance of our algorithm is better than that of other methods.

Original: As for the results for SBM in Table
Revised: As for the results for SBM in Table

Original: A.2, the error of estimating the number of communities decreases as the subsampling proportion increases, and the error increases as the observed graphs are noisier (pout increases).
Revised: A.2, the error of estimating the number of communities decreases as the subsampling proportion increases, and the error increases as the observed graphs are noisier (pout increases).

Original: As for the computational time recorded in Table A.3 for SBM, we observe that the estimation time of the full sample is two orders of magnitude larger than the time of subsampling and estimation. 
Revised: As for the computational time recorded in Table A.3 for SBM, we observe that the estimation time of the full sample is two orders of magnitude larger than the time of subsampling and estimation. As the

Original: As the subsampling proportion increases, the computational time of subsampling methods increases.
Revised: [Table] subsampling proportion increases, the computational time of subsampling methods increases.

Original: Though our method’s computational time is slightly larger than other subsampling methods, our method has better accuracy.
Revised: Though our method’s computational time is slightly larger than other subsampling methods

Original: [Table]
Revised: [Table]

Revised: our method has better accuracy.

Original: [Table]
Revised: [Table]

Original: B. E XPERIMENTS ON R EAL - WORLD D ATASETS B.2. D ESCRIPTION OF THE D ATASETS
Revised: B. E XPERIMENTS ON R EAL - WORLD D ATASETS B.2. D ESCRIPTION OF THE D ATASETS

Original: We use five real-world datasets to evaluate the performance of our method as well.
Revised: We use five real-world datasets to evaluate the performance of our method as well.

Original: These datasets are widely used as benchmark datasets in the study of community detection for graphs.
Revised: These datasets are widely used as benchmark datasets in the study of community detection for graphs.

Original: They can be easily found at KONECT, UCI network data repository and Stanford Network Analysis Project
Revised: They can be easily found at KONECT, UCI network data repository and Stanford Network Analysis Project

Original: Rossi & Ahmed (2015); Rozemberczki et al.
Revised: Rossi & Ahmed (2015); Rozemberczki et al.

Original: (2020a);
Revised: (2020a);

Original: Sen et al.
Revised: Sen et al.

Original: All datasets are considered undirected and unweighted graphs, and all self-loop edges and isolated nodes are removed.
Revised: All datasets are considered undirected and unweighted graphs, and all self-loop edges and isolated nodes are removed.

Original: More details about the datasets are as follow.
Revised: More details about the datasets are as follow

Original: Polblogs : The political blogs network dataset was collected by Adamic and Glance in 2005.
Revised: Polblogs : The political blogs network dataset was collected by Adamic and Glance in 2005.

Original: A network is constructed among all the posts published by liberal or conservative bloggers (two communities).
Revised: A network is constructed among all the posts published by liberal or conservative bloggers (two communities).

Original: Each node represents one post.
Revised: Each node represents one post.

Original: An edge connects two nodes if one of them is cited by the other.
Revised: An edge connects two nodes if one of them is cited by the other.

Original: Polbooks : This is a network of books about US politics published around the 2004 presidential election and sold by online bookseller Amazon.com .
Revised: Polbooks : This is a network of books about US politics published around the 2004 presidential election and sold by online bookseller Amazon.com .

Original: All the books are divided into four communities by NI-LPA.
Revised: All the books are divided into four communities by NI-LPA.

Original: Edges between books represent frequent co-purchasing of books by the same buyers.
Revised: Edges between books represent frequent co-purchasing of books by the same buyers.

Original: Facebook : This is an ego-network dataset of the “friend circles” of one anonymous user on Facebook.
Revised: Facebook : This is an ego-network dataset of the “friend circles” of one anonymous user on Facebook.

Original: The network forms friend circles such as family members, high school friends or other friends that are “hand-labeled” by the user.
Revised: The network forms friend circles such as family members, high school friends or other friends that are “hand-labeled” by the user.

Original: Cora : The Cora dataset describes the citation relationship among scientific publications classified into seven classes.
Revised: Cora : The Cora dataset describes the citation relationship among scientific publications classified into seven classes.

Original: After preprocessing, there remain 2,485 nodes and 5,069 links.
Revised: After preprocessing, there remain 2,485 nodes and 5,069 links.

Original: PubMed : The PubMed dataset describes the citation relationship among scientific publications classified into seven classes.
Revised: PubMed : The PubMed dataset describes the citation relationship among scientific publications classified into seven classes.

Original: After preprocessing, there remain 19,717 nodes and 44,338 links.
Revised: After preprocessing, there remain 19,717 nodes and 44,338 links.

Original: Table A.4 presents the summary statistics for those real networks we used in this paper.
Revised: Table A.4 presents the summary statistics for those real networks we used in this paper.

Original: As we can see, these datasets cover different levels of network size, density, number fo communities, clustering coefficient (CC), and imbalance.
Revised: As we can see, these datasets cover different levels of network size, density, number fo communities, clustering coefficient (CC), and imbalance.

Original: The degree distributions of all the datasets are presented in Figure
Revised: The degree distributions of all the datasets are presented in Figure

Original: A.1.
Revised: A.1.

Original: [Table]
Revised: [Table]

Original: Dataset

Original: Full

Original: Ricci

Original: Prop

Original: ORG-sub

Original: MHRW

Original: CSE FFS SnowBall DBN
Revised: [Table]

Original: RW MDRW Polbooks 1. Facebook 29.14 0.
Revised: [Table]

Original: Cora 191.06 0. Polblogs 48. PubMed NA 41.
Revised: B.2. C OMPUTATIONAL T IME

Original: [Table]

Original: The comparisons of errors in estimating the number of communities are presented in the main context.
Revised: The comparisons of errors in estimating the number of communities are presented in the main context.

Original: The computation time is also obtained in A.5.
Revised: The computation time is also obtained in A.5.

Original: Similar to the synthetic datasets, the estimation time of the full sample is two orders of magnitude larger than the time of subsampling (including the time for computing OR curvature for each edge in the graph) and estimation.
Revised: Similar to the synthetic datasets, the estimation time of the full sample is two orders of magnitude larger than the time of subsampling (including the time for computing OR curvature for each edge in the graph) and estimation.

Original: The result also shows that the computation time for estimating the number of communities using the subsampling method is much shorter than that using the full dataset.
Revised: The result also shows that the computation time for estimating the number of communities using the subsampling method is much shorter than that using the full dataset.

Original: The complexity of all methods is influenced by the size of the graph, especially by the number of nodes.
Revised: The complexity of all methods is influenced by the size of the graph, especially by the number of nodes.

Original: It is no wonder that the computation time of all methods in estimating the number of communities for dataset Cora is much longer than other datasets.
Revised: It is no wonder that the computation time of all methods in estimating the number of communities for dataset Cora is much longer than other datasets.

Original: B. A DDITIONAL EMPIRICAL RESULTS
Revised: B. A DDITIONAL EMPIRICAL RESULTS

Original: In this section, we provide additional empirical results.
Revised: In this section, we provide additional empirical results.

Original: In real-world networks, the ground truth is not available.
Revised: In real-world networks, the ground truth is not available.

Original: We used some data sets that are widely used for community detection tasks, where the communities are labeled with domain knowledge.
Revised: We used some data sets that are widely used for community detection tasks, where the communities are labeled with domain knowledge.

Original: The manually labeled community structures correspond to different underlying aspects of the nodes.
Revised: The manually labeled community structures correspond to different underlying aspects of the nodes.

Original: For example, in the Cora citation network (where a node represents a paper), the community label corresponds to the paper’s subject (e.g., neural network).
Revised: For example, in the Cora citation network (where a node represents a paper), the community label corresponds to the paper’s subject (e.g., neural network).

Original: We agree that the manually labeled number of communities might not be the ground truth, which is a limitation for real networks.
Revised: We agree that the manually labeled number of communities might not be the ground truth, which is a limitation for real networks.

Original: To overcome this limitation, we also employ another evaluation metric to assess the performance of our methods.
Revised: To overcome this limitation, we also employ another evaluation metric to assess the performance of our methods.

Original: In particular, we do not compare our estimation with the manually labeled number of communities.
Revised: In particular, we do not compare our estimation with the manually labeled number of communities.

Original: We compare our estimation with the estimation obtained from the full data.
Revised: We compare our estimation with the estimation obtained from the full data.

Original: These comparison results can be used to evaluate whether the subsampled subgraph can be a good surrogate to carry out computations of interest for the full data.
Revised: These comparison results can be used to evaluate whether the subsampled subgraph can be a good surrogate to carry out computations of interest for the full data.

Original: Table A.6 shows the estimation difference between the subsampled and full graphs for four datasets, i.e., Polbooks, Facebook, Cora, and
Revised: Table A.6 shows the estimation difference between the subsampled and full graphs for four datasets, i.e., Polbooks, Facebook, Cora, and

Original: Polblogs.
Revised: Polblogs.

Original: Here we do not show the results for PubMed, since the estimation results of PubMed full graph are not available due to prohibitive computation.
Revised: Here we do not show the results for PubMed, since the estimation results of PubMed full graph are not available due to prohibitive computation.

Original: As we can see, using this new metric, our method ORG-sub still outperforms other subsampling methods.[Table] Prop. ORG-sub
Revised: As we can see, using this new metric, our method ORG-sub still outperforms other subsampling methods.

Original: MHRW

Original: CSE

Original: FFS

Original: SnowBall

Original: DBN

Original: RW

Original: MDRW

Original: Polbooks

Original: [Table]
Revised: [Table]

Original: C D ISCUSSION ON EDGE DENSITY
Revised: C D ISCUSSION ON EDGE DENSITY

Original: Indeed, the performance of our method depends on the edge density from both theory and empirical results.
Revised: Indeed, the performance of our method depends on the edge density from both theory and empirical results.

Original: On one hand, Corollary 4.1.2 basically assumes that ρ > log (˜ n ) ˜ n , where ˜ n is the number of sampled nodes, and ρ is the edge density.
Revised: On one hand, Corollary 4.1.2 basically assumes that ρ > log (˜ n ) ˜ n , where ˜ n is the number of sampled nodes, and ρ is the edge density.

Original: When ˜ n = 100 , we require ρ > 0 . 046 , and when ˜ n = 1000 , we only require ρ > 0 .
Revised: When ˜ n = 100 , we require ρ > 0 . 046 , and when ˜ n = 1000 , we only require ρ > 0 .

Original: As the subsample size becomes larger, the constraint imposed on the edge density becomes weaker.
Revised: As the subsample size becomes larger, the constraint imposed on the edge density becomes weaker.

Original: On the other hand, we include additional simulation studies to show the empirical performance of our method with different network densities.
Revised: On the other hand, we include additional simulation studies to show the empirical performance of our method with different network densities.

Original: Particularly, we use the same
Revised: Particularly, we use the same

Original: SBM setting in Section 5.1 to generate synthetic data, except that we let p in = 0 .
Revised: SBM setting in Section 5.1 to generate synthetic data, except that we let p in = 0 .

Original: 8 λ and p out = 0 .
Revised: 8 λ and p out = 0 .

Original: 1 λ , where λ ∈ { 0 . 2 , 0 . 4 , . .
Revised: 1 λ , where λ ∈ { 0 . 2 , 0 . 4 , . .

Original: , 1 } controls the edge density level.
Revised: , 1 } controls the edge density level.

Original: Here, we fix the sampling proportion as 10%, that is, we sample 90 nodes.
Revised: Here, we fix the sampling proportion as 10%, that is, we sample 90 nodes.

Original: Table A.7 shows the edge density and average estimation error under different λ .
Revised: Table A.7 shows the edge density and average estimation error under different λ .

Original: As we can see, an increasing edge density comes with a lower error.
Revised: As we can see, an increasing edge density comes with a lower error.

Original: [Table]
Revised: [Table]

Original: D R ELATED WORK ON GRAPH SUBSAMPLING
Revised: D R ELATED WORK ON GRAPH SUBSAMPLING

Original: Existing graph subsampling techniques can be categorized into three main groups: node-based, edge-based, and exploration-based.
Revised: Existing graph subsampling techniques can be categorized into three main groups: node-based, edge-based, and exploration-based.

Original: Node-based sampling.
Revised: Node-based sampling.

Original: Random Node (RN) sampling (Leskovec & Faloutsos, 2006) and Degreebased Node sample (DBN) (Leskovec & Faloutsos, 2006) sampling are two most common node-base sampling methods.
Revised: Random Node (RN) sampling (Leskovec & Faloutsos, 2006) and Degreebased Node sample (DBN) (Leskovec & Faloutsos, 2006) sampling are two most common node-base sampling methods.

Original: RN selects a set of nodes uniformly at random from the graph, while DBN selects a node with a probability that is proportional to its degree.
Revised: RN selects a set of nodes uniformly at random from the graph, while DBN selects a node with a probability that is proportional to its degree.

Original: DBN has been shown to favor high-degree nodes.
Revised: DBN has been shown to favor high-degree nodes.

Original: Edge-based sampling.
Revised: Edge-based sampling.

Original: Random Edge (RE) sampling (Leskovec & Faloutsos, 2006) generates an induced subgraph by selecting edges uniformly at random.
Revised: Random Edge (RE) sampling (Leskovec & Faloutsos, 2006) generates an induced subgraph by selecting edges uniformly at random.

Original: Some variants of RE have been proposed, such as Random Node-Edge (RNE) sampling (Rafiei, 2005) that randomly selects a node and then randomly chooses an adjacent edge.
Revised: Some variants of RE have been proposed, such as Random Node-Edge (RNE) sampling (Rafiei, 2005) that randomly selects a node and then randomly chooses an adjacent edge.

Original: Previous studies (Leskovec & Faloutsos, 2006) have demonstrated that neither RE nor RNE preserves community structures because the resulting sampled graphs are often sparsely connected.
Revised: Previous studies (Leskovec & Faloutsos, 2006) have demonstrated that neither RE nor RNE preserves community structures because the resulting sampled graphs are often sparsely connected.

Original: Meanwhile, both RE and REN slightly favor high-degree nodes because the probability of selecting a node increases with its degree.
Revised: Meanwhile, both RE and REN slightly favor high-degree nodes because the probability of selecting a node increases with its degree.

Original: Exploration-based sampling.
Revised: Exploration-based sampling.

Original: SnowBall sampling (Goodman, 1961) and Forest Fire sampling (FFS) (Leskovec et al., 2005) are two basic exploration-based sampling methods, selecting a fixed fraction of neighbors visited at each iteration.
Revised: SnowBall sampling (Goodman, 1961) and Forest Fire sampling (FFS) (Leskovec et al., 2005) are two basic exploration-based sampling methods, selecting a fixed fraction of neighbors visited at each iteration.

Original: Random walk (RW)
Revised: Random walk (RW)

Original: (Lovász, 1996) is also a popular explorationbased graph sampling method, which selects the next node at random from the neighbors of the currently selected node.
Revised: (Lovász, 1996) is also a popular explorationbased graph sampling method, which selects the next node at random from the neighbors of the currently selected node.

Original: One of the major limitations of RW is that RW is inherently biased towards visiting high-degree nodes.
Revised: One of the major limitations of RW is that RW is inherently biased towards visiting high-degree nodes.

Original: To overcome the drawback of RW, researchers proposed Metropolis
Revised: To overcome the drawback of RW, researchers proposed Metropolis

Original: Hasting Random Walk Sampler (MHRW)
Revised: Hasting Random Walk Sampler (MHRW)

Original: (Hübler et al., 2008) and Multi-dimensional Random Walk
Revised: (Hübler et al., 2008) and Multi-dimensional Random Walk

Original: Sampler (MDRW) (Ribeiro & Towsley, 2010).
Revised: Sampler (MDRW) (Ribeiro & Towsley, 2010).

Original: Most of the aforementioned literature aims at getting a subgraph that preserves some network summary statistics, such as degree distribution.
Revised: Most of the aforementioned literature aims at getting a subgraph that preserves some network summary statistics, such as degree distribution.

Original: Only a few works aim at obtaining a subgraph that preserves the community information.
Revised: Only a few works aim at obtaining a subgraph that preserves the community information.

Original: (Maiya & Berger-Wolf, 2010) proposed a local greedy search-based community structure expansion sampling (CSE) method to optimize the preservation of community structures.
Revised: (Maiya & Berger-Wolf, 2010) proposed a local greedy search-based community structure expansion sampling (CSE) method to optimize the preservation of community structures.

Original: Despite many successful applications of existing graph sampling methods.
Revised: Despite many successful applications of existing graph sampling methods.

Original: They all have various limitations.
Revised: They all have various limitations.

Original: In particular, node-based and edge-based sampling methods sample nodes or edges independently, ignoring the neighborhoods of seed nodes.
Revised: In particular, node-based and edge-based sampling methods sample nodes or edges independently, ignoring the neighborhoods of seed nodes.

Original: They might obtain a disconnected subgraph from a connected graph (Wu et al., 2016).
Revised: They might obtain a disconnected subgraph from a connected graph (Wu et al., 2016).

Original: For RW and FF, it has been shown that they could get trapped inside the communities and leave other communities out of the sample entirely (Wu et al., 2016).
Revised: For RW and FF, it has been shown that they could get trapped inside the communities and leave other communities out of the sample entirely (Wu et al., 2016).

Original: While the MHRW algorithm ensures that the subgraph preserves degree distribution, the performance of this algorithm is dependent on its sample acceptance ratio.
Revised: While the MHRW algorithm ensures that the subgraph preserves degree distribution, the performance of this algorithm is dependent on its sample acceptance ratio.

Original: It has been shown that the acceptance ratio of MHRW is typically very low in real-world networks .
Revised: It has been shown that the acceptance ratio of MHRW is typically very low in real-world networks .

Original: Therefore, MHRW generally suffers from the sample rejection problem, which clearly degrades the performance of
Revised: Therefore, MHRW generally suffers from the sample rejection problem, which clearly degrades the performance of

Original: MHRW.
Revised: MHRW.

Original: For CSE, the selection procedure of the next node is not based on a mathematical framework and one cannot compute the probability of visiting sampled nodes in CSE (Salehi et al., 2012).
Revised: For CSE, the selection procedure of the next node is not based on a mathematical framework and one cannot compute the probability of visiting sampled nodes in CSE (Salehi et al., 2012).

Original: E E MPIRICAL INVESTIGATION OF FUTURE WORK
Revised: E E MPIRICAL INVESTIGATION OF FUTURE WORK

Original: Theoretically, the OR curvature has been proven to be related to some network summaries concerning communities, such as the eigenvalues of the graph Laplacian Bauer et al.
Revised: Theoretically, the OR curvature has been proven to be related to some network summaries concerning communities, such as the eigenvalues of the graph Laplacian Bauer et al.

Original: (2011) and the clustering coefficient (CC) Jost & Liu (2014).
Revised: (2011) and the clustering coefficient (CC) Jost & Liu (2014).

Original: Empirically, we calculated the CC value of the full graph and that of the subsampled graph for two real datasets: Facebook and Cora.
Revised: Empirically, we calculated the CC value of the full graph and that of the subsampled graph for two real datasets: Facebook and Cora.

Original: We set the subsampling proportion to 10%.
Revised: We set the subsampling proportion to 10%.

Original: We replicated the sampling procedure 100 times and thus obtained 100 subsamples.
Revised: We replicated the sampling procedure 100 times and thus obtained 100 subsamples.

Original: Using these subsamples, we obtained the 95% confidence interval (CI) by calculating [Table] [Table] [Table]
Revised: Using these subsamples, we obtained the 95% confidence interval (CI) by calculating mean ± 1.96*sd (standard deviation).

Original: [Table]
Revised: For Facebook data, the CI of the CC value is [0.420, 0.537], and the CC value obtained from the full graph is 0.476.

Original:  covers 
Revised: For Cora data, the CI is [0.234, 0.290], and the CC value of the full graph is 0.238. As we can see, the 95% CI covers the CC value of the full graph for both real data sets.

Original: These observations suggest that our subsampling algorithm could preserve CC to some extent.
Revised: These observations suggest that our subsampling algorithm could preserve CC to some extent.

Original: The consistency theory on CC using our method is under investigation and the results will be reported in future publications.
Revised: The consistency theory on CC using our method is under investigation, and the results will be reported in future publications.