We include it as a regularization term with a parameter λ that controls the importance of the interactions.
the M-step is as follows: Since there are no closed form solutions for those parameters, we use gradient ascent to solve them parameter-by-parameter.