So we just have to start in the bottom layer of an all blue cube, go through all colorings of the current layer until we get to all white, then move one layer up and repeat until we reach the top layer. When we move to the next layer up, we have to color at least one vertex white or else the coloring would not change, so there are only 2^n - 1 relevant colorings in each layer. In the end we also have to add one for the all blue cube.
1, 1-1, 1-2-1, 1-3-3-1 are the layers in dimensions zero to three. The number of relevant colorings per layer are 1, 1-1, 1-3-1, 1-7-7-1, their sum plus one are 2, 3, 6, 17. So this worked for dimensions zero to two but gives the wrong answer for three which should be 20, not 17 according to the article. Where did I make the mistake? Was the visualization with the cube oversimplifying the problem?
For dimensions four I would guess the layer structure 1-4-6-4-1 which yields 96 which also does not match 168.