Inspired by the multicanonical approach to simulations of first-order
phase transitions we propose for q-state Potts models a combination of
cluster updates with reweighting of the bond configurations in the Fo
rtuin-Kastelein-Swendsen-Wang representation of this model. Numerical
tests for the two-dimensional models with q = 7, 10 and 20 show that t
he autocorrelation times of this algorithm grow with the system size V
as tau proportional to V-alpha, where the exponent takes the optimal
random walk value of alpha approximate to 1.