Scalable Second-order Riemannian Optimization for $K$-means Clustering