在python2.7的scikit-learn之外，有没有K-means++的实现？

import numpy as np

def kmeans_plusplus(X, K):
    centroids = [X[np.random.randint(X.shape[0])]]
    while len(centroids) < K:
        distances = np.array([np.min(np.sum((x - centroids)**2, axis=1)) for x in X])
        probabilities = distances / np.sum(distances)
        new_centroid = X[np.random.choice(X.shape[0], p=probabilities)]
        centroids.append(new_centroid)
    return np.array(centroids)

# 示例用法
X = np.array([[1, 2], [1, 4], [1, 0], [4, 2], [4, 4], [4, 0]])
K = 2
centroids = kmeans_plusplus(X, K)
print(centroids)