A discrete random variable has uniform distribution with parameter $\theta$ and we have a sample $X$ of size $n$.
$$P_{\theta}(\xi=k)=\theta^{-1},\theta\in\{1,2\dots\},k\in\{1,2,\dots\theta\}$$
I've got required expectation as a function from generalized harmonic numbers.
\begin{gather*} P(X_{(n)}=k)=\\ =\sum\limits_{i=1}^nP(\text{exactly $i$ observations equals $k$ and other are less})=\\ =\sum\limits_{i=1}^n\binom{i}{n}\left(\frac{1}{\theta}\right)^i\left(\frac{k-1}{\theta}\right)^{n-i}=\left(\frac{k}{\theta}\right)^n.\\ \mathbb{E}[X_{(n)}]= \sum\limits_{k=1}^\theta k\left(\frac{k}{\theta}\right)^n= \frac{1}{\theta^n}\sum\limits_{k=1}^\theta k^{n+1}= \frac{1}{\theta^n}H_{\theta,-(n+1)}. \end{gather*}
Is there a way to represent expectation of maximum from sample $X$ with a simpler expression?