I assume that $1<\alpha<2$. If not, as remarked by the questioner, $f(x)$ is Lipschitz on $(0,1)$ and the problem is simpler.
You can use the following approach. Take a small $\delta>0$. You wish to bound $V=\sum_i |f(x_i)-f(y_i)|$ whenever $S=\{(x_1,y_1),\dots,(x_n,y_n)\}$ is a finite set of pairwise disjoint intervals contained in $(0,1)$ with $\sum_i |x_i-y_i|<\delta$.
Let $N$ be the smallest integer such that $N>(2\pi \delta)^{-1}$, and split the interval into two pieces, $(0,1/(2\pi N))$ and $(1/(2\pi N),1)$. By splitting an interval in $S$ if necessary, which does not decrease $V$, you can assume that each member of $S$ is in either $(0,1/(2\pi N))$ or $(1/(2\pi N),1)$.
Let $V_1$ be the portion of $V$ coming from the intervals in $(0, 1/(2\pi N))$. Here, the function $f(x)$ has alternating local maxima and minima. The maxima are close to the values $x=2/((4n+1)\pi)$; let $M_n$ be the value of $x$ close to $x=2/((4n+1)\pi)$ where $f(x)$ has a local maximum. Similarly, the minima are close to the values $x=2/((4n+3)\pi)$; let $m_n$ be the value of $x$ close to $x=2/((4n+3)\pi)$ where $f(x)$ has a local minimum.
Argue that
$$V_1\le |f(M_N)|+|f(M_N)-f(m_N)|+|f(m_N)-f(M_{N+1})|+\cdots\ \ \ (1)$$
and that
$$
f(M_n)=f(m_n)=O(n^{-\alpha}).\qquad (2)
$$
Then, combining $(1)$ and $(2)$, conclude that $V_1=O(N^{1-\alpha})$.
Let $V_2$ be the portion of $V$ coming from the intervals in $(1/(2\pi N),1)$. In this portion of $(0,1)$, $|f'(x)|$ is bounded above, so the function $f(x)$ is Lipschitz. Argue that the Lipschitz constant is $O(\delta^{\alpha-2})$. Therefore, $V_2$ is $O(\delta^{\alpha-1})$.
Since $V=V_1+V_2$, adding the above estimates together should prove that $V\to 0$ as $\delta\to 0$.