refactor
connect components in grapgh
find digital root
solve linear equation system by forward substitution on lower triangular matrix
using boost accumulator for statistics