Boost-based implementation of Bairstow's (YES, I'm aware that it's not 'Baristow's'  as in code, fix if you like) method of polynomial decomposition.