library(cluster)

input_data = cbind(rnorm(200),
                   rnorm(200))
k = 6

clara_clusters1 = clara(input_data, k, medoids.x = FALSE)$clustering
clara_clusters2 = clara(input_data, k, medoids.x = FALSE)$clustering

pam_clusters1 = pam(input_data, k, cluster.only = TRUE)
pam_clusters2 = pam(input_data, k, cluster.only = TRUE)

print(clara_clusters1 - clara_clusters2)
print(pam_clusters1 - pam_clusters2)