We generalize the K matrix formulation to non-trivial non-Abelian families of 2+1D topological orders. Given a topological order C, any topological order in the same non-Abelian family as C can be efficiently described by a = (aI ) where aI are Abelian anyons in C, together with a symmetric invertible matrix K, KIJ = kIJ − ta I ,a J where kIJ are integers, kII are even and ta I ,a J are the mutual statistics between aI , aJ . In particular, when C is a root whose rank is the smallest in the family, K becomes an integer matrix. Our results make it possible to generate the data of large numbers of topological orders instantly.