A calculation circuit calculates, among a plurality of neurons whose states are represented by variables each having m values (m is a positive integer of 3 or greater), two energy changes caused by a state change of a second neuron by 2n (n is an integer of 0 or greater) in positive and negative directions, based on a state change direction of a first neuron whose state has been updated and a weighting coefficient indicating magnitude of an interaction between the first and second neurons. A state transition determination circuit determines, based on magnitude relationships among a thermal excitation energy and the two energy changes, whether to allow updates of the state of the second neuron that cause the two energy change, outputs the determination results, and limits updates by which the state of the second neuron exceeds an upper limit value or falls below a lower limit value.