Without using the prime number theorem, we obtain the asymptotics of the $r$th largest prime divisor of a harmonically distributed random positive integer $N$; harmonic asymptotics are obtained from asymptotics of the zeta distribution via Tauberian methods. (Knuth and Trabb-Pardo need a strong form of the prime number theorem to obtain the distributions when $N$ is uniformly distributed.) A trick brings in Poisson variates, and then we can use the methods developed for the fractional length of the $r$th longest cycle in a random permutation.