I was hoping to get around that due to the use of the `initial_clusterers` parameter in the calling function. If it is None, it uses this function, which gives the "default" initial clusters as used in the reference.

The variation in k is also needed -- k varies between 10 and 30 (unless that doesn't make sense), forcing different clusters in most cases.