(2006), there are two terms, and the second is easily computable given the feature expectations obtained by forward/backward and the entropy for the sequence.
unlike the previous method, here the first term can be efficiently calculated as well.