Can we find the distribution and/or expected value of $$S=\sum_{n \geq 1} \frac{1}{(\text{lcm} (n,r_n))^2}$$
where $r_n$ is a uniformly distributed random integer, $r_n \in [1,n]$ and $\text{lcm}$ is the least common multiple?
Or maybe an estimate is possible?
Some boundaries are easy to see. For example:
$$1
For two trials with $10000$ runs each and summing up to $n=20000$ I obtained the following distributions:
It seems as if there are two peaks.
For 20000 trials we obtain the following distribution:
For 100000 trials:


