numpy.unique が使えれば、C言語の速度で計算できる