Supersymmetric elementary string states in the compactified heterotic string theory are described by small black holes that have zero area event horizon. In this paper we compute the supersymmetric index of such elementary string states using gravitational path integral. The dominant contribution to the path integral comes from an Euclidean rotating black hole solution of the supergravity theory with a finite area event horizon, but the logarithm of the index, computed from the saddle point, vanishes. Nevertheless we show that the solution is singular on certain subspaces of the horizon where higher derivative corrections can be important, and once the higher derivative corrections are taken into account the solution could yield a finite result for the logarithm of the index whose form agrees with the microscopic results up to an overall numerical constant. While the numerical constant is not determined in our analysis, we show that it is independent of the details of the compactification and even the number of non-compact dimensions, in agreement with the microscopic results.