We present a new framework for analysis and visualization of complex networks based on structural information retrieved from their distance k-graphs and B-matrices. The construction of B-matrices for graphs with more than 1 million edges requires massive Breadth-First Search (BFS) computations and is facilitated using new software prepared for distributed environments. Our framework benefits from data parallelism inherent to all-pair shortest-path problem and extends Cassovary, an open-source in-memory graph processing engine, to enable multinode computation of distance k-graphs and related graph descriptors. We also introduce a new type of B-matrix, constructed using clustering coefficient vertex invariant, which can be generated with a computational effort comparable with the one required for a previously known degree B-matrix, while delivering an additional set of information about graph structure. Our approach enables efficient generation of expressive, multidimensional descriptors useful in graph embedding and graph mining tasks. The experiments showed that the new framework is scalable and for specific all-pair shortest-path task provides better performance than existing generic graph processing frameworks. We further present how the developed tools helped in the analysis and visualization of real-world graphs from Stanford Large Network Dataset Collection. Figure 2. Results of degree B-matrix computation performance tests for the three graphs from SNAP database: Epinions (75,879 vertices, 508,837 edges), Slashdot0902 (82,168 vertices, 948,464 edges), and Web Notre Dame (325,729 vertices, 1,497,134 edges). (a) Dependency of computation time on a number of processors for Epinions and Slashdot0902 graphs, (b) Dependency of computation time on a number of processors for Web Notre Dame graph, (c) Speedup vs. number of processors for Epinions and Slashdot0902 graphs, (d) Speedup vs. number of processors for Web Notre Dame graph.Figure 3. Time of loading and indexing list of edges by Cassovary versus graph size (number of vertices + number of nodes).In this section, we examine discriminating capabilities of B-matrices in the unsupervised learning task. To this end, we constructed feature vectors representing graphs from different groups by coarse graining and row packing of the B-matrices. Next, the dimensionality reduction techniques such as principal component analysis (PCA) and t-distributed stochastic neighbor embedding [48] were used to obtain 2D embedding of graphs for visual comparison of cluster proximity. Let Y stand for V or E symbol, so that B Y;F denotes B V;F or B E;F matrix of a graph G. The set F contains degree or clustering coefficient vertex-domain functions defined for consecutive distance k-graphs. The number of categories n was chosen so that n < jV .G/j. The size of category s b is the same for all categories (bins) in a test setup. The parameter s b reflects granularity of sampling