Partir de la definición de Polinomios de Bernoulli en función de su función de generación
$$\frac{t e^{xt}}{e^t - 1} = \sum_{n=1}^\infty \frac{t^n}{n!} B_n(x)$$
Si se multiplican ambos lados por $e^{-2\pi i kx}$ para cualquier $k \in \mathbb{Z}$ e integramos, obtenemos
$$ \sum_{n=1}^\infty \frac{t^n}{n!} \int_0^1 B_n(x) e^{-2\pi i k x} dx = \frac{t}{e^{t}-1} \int_0^1 e^{(t-2\pi i k) x} dx = \frac{t}{t-2\pi i k} $$ Comparando los coeficientes de $t^n$ de dos lados, encontramos para $n > 0$ ,
$$\int_0^1 B_n(x) e^{-2\pi ik x} dx = \begin{cases} 0, & k = 0,\\ - \frac{n!}{(2\pi ik)^n}, & k \ne 0 \end{cases} $$ En consecuencia, para $n > 0$ el polinomio de Bernoulli $B_n(x)$ tiene lo siguiente Expansión en serie de Fourier sobre $(0,1)$ .
$$B_n(\{x\}) = -\frac{n!}{(2\pi i)^n} \sum_{|k|>0} \frac{e^{2\pi i kx}}{k^n} \tag{*1}$$
Para grandes $n$ el LHS de $(*1)$ está dominado por los dos términos con $|k| = 1$ . Esto significa que para cualquier $x \in (0,1)$ , si se fija la paridad de $n$ y enviar $n$ a $\infty$ tendremos
$$(-1)^{\lfloor n/2 \rfloor - 1} \frac{(2\pi)^n}{2(n!)} B_n(\{x\}) \quad\to\quad \begin{cases} \cos(2\pi x), &n \text{ even}\\ \sin(2\pi x), &n \text{ odd} \end{cases} \tag{*2} $$
Pasemos a la evaluación de la integral $u_n$ para $n > 0$ .
Sea $u = 1/x$ tenemos
$$\begin{align} &\int_0^1 B_n\left(\left\{ \frac1x \right\}\right) \frac{dx}{x} = \lim_{N\to\infty} \int_{1/N}^1 B_n\left(\left\{ \frac1x \right\}\right) \frac{dx}{x}\\ =& \lim_{N\to\infty} \int_1^N B_n(\{u\}) \frac{du}{u} = \lim_{N\to\infty} \sum_{j=1}^{N-1} \int_j^{j+1} B_n(\{u\}) \frac{du}{u}\\ =& \lim_{N\to\infty} \int_0^1 B_n(u) \sum_{j=1}^{N-1} \frac{1}{u+j} du = \lim_{N\to\infty} \int_0^1 B_n(u) \sum_{j=1}^{N-1} \left( \frac{1}{u+j} - \frac{1}{j} \right) du\\ =& \int_0^1 B_n(u) \sum_{j=1}^\infty \left( \frac{1}{u+j} - \frac{1}{j} \right) du \end{align} $$ Compara la serie de la última integral con la siguiente expansión de función digamma
$$\psi(1+z) = -\gamma + \sum_{k=1}^\infty \left(\frac{1}{k} - \frac{1}{k+z}\right)$$
encontramos para $n > 0$ , $$\int_0^1 B_n\left(\left\{ \frac1x \right\}\right) \frac{dx}{x} = -\int_0^1 B_n(u) \psi(1+u) du\tag{*3}$$
Combine $(*2)$ y $(*3)$ para grandes $k$
$$u_{2k+1} \sim (-1)^{k-1}\frac{2C(2k+1)!}{(2\pi)^{2k+1}} \quad\implies\quad |u_{2k+1}| \to \infty \;\;\text{ as }\;\; k \to \infty $$
donde $C$ es una constante definida por una integral $$\begin{align} C & = -\int_0^1 \sin(2\pi x) \psi(1+x) dx = \int_1^\infty \frac{\sin(2\pi x)}{x} dx = \frac{\pi}{2} - \text{Si}(2\pi)\\ & \approx 0.1526447506622681689855415293400020129561... \end{align}$$