Consider $(1)$, $H_n$ is the nth-harmonic number
$$\sum_{n=1}^{\infty}\left({1\over (nx)^2-1}+{2\over (nx)^2-2^2}+{3\over (nx)^2-3^2}+\cdots+{x-1\over (nx)^2-(x-1)^{2}}\right)=S\tag1$$ $x\ge2$
How does one show that $$\color{blue}{S={H_{x-1}\over 2}}?$$
An attempt:[Edited]
Because $(1)$ is the difference of two squares form, so it becomes
$$\sum_{n=1}^{\infty}\left[\left({1\over nx-1}+{1\over nx-2}+{1\over nx-3}+\cdots+{x-1\over nx-(x-1)}\right)-\left({1\over nx+1}+{1\over nx+2}+{1\over nx+3}+\cdots+{x-1\over nx+(x-1)}\right)\right]=2S\tag2$$
Not sure how to proceed next
From Cotangent identity (formula 18) and Digamma reflection formula: $$ \begin{align} \pi &\,\cot(\pi z)=\frac{1}{z}+2z\,\sum_{n=1}^{\infty}\,\frac{1}{z^2-n^2}=\psi(1-z)-\psi(z) \\[3mm] S &= \sum_{n=1}^{\infty}\,\sum_{m=1}^{x-1}\,\frac{m}{(nx)^2-m^2} = \sum_{m=1}^{x-1}\,\frac{m}{x^2}\sum_{n=1}^{\infty}\,\frac{1}{n^2-(m/x)^2} \\[3mm] &= \frac{1}{2x}\,\sum_{m=1}^{x-1}\,\left[-2\left(\frac{m}{x}\right)\sum_{n=1}^{\infty}\,\frac{1}{\left(\frac{m}{x}\right)^2-n^2}\right] \\[3mm] &= \frac{1}{2x}\,\sum_{m=1}^{x-1}\,\left[\frac{x}{m}-\pi\,\cot\left(\pi\,\frac{m}{x}\right)\right] \\[3mm] &= \frac{1}{2x}\,\sum_{m=1}^{x-1}\,\left[\frac{x}{m}+\psi\left(\frac{m}{x}\right)-\psi\left(1-\frac{m}{x}\right)\right] \\[3mm] &= \frac{1}{2}\,\sum_{m=1}^{x-1}\frac{1}{m}+\frac{1}{2x}\left[\sum_{m=1}^{x-1}\psi\left(\frac{m}{x}\right)-\sum_{m=1}^{x-1}\psi\left(\frac{x-m}{x}\right)\right] \\[3mm] &= \frac{1}{2}\,\sum_{m=1}^{x-1}\frac{1}{m} = \color{red}{\frac{H_{x-1}}{2}} \end{align} $$ Where: $$ \color{blue}{\sum_{m=1}^{x-1}\psi\left(\frac{x-m}{x}\right)}=\psi\left(\frac{x-1}{x}\right)+\psi\left(\frac{x-2}{x}\right)+\cdots+\psi\left(\frac{2}{x}\right)+\psi\left(\frac{1}{x}\right)=\color{blue}{\sum_{m=1}^{x-1}\psi\left(\frac{m}{x}\right)} $$