region method 2

compute the order of skelton by BFS