Not sure about the reference, but the "probability proof" is already nearly trivial.
Just notice that you can run your assignment step by step, at each step either placing some particular weighted edge between a uniformly chosen random pair of not yet joined vertices, or putting a uniformly chosen random remaining edge between a fixed pair of not yet joined vertices (which option to choose is entirely up to you at each step and you can even make your choice between the two depending on what has been constructed already). Now observe that after placing the shortest edge, you should place the second shortest one disjointly. Then you can fill the rest of the four vertex graph in whatever way, but after that the shortest remaining edge has to be placed disjointly with the first two. Then you fill the 6 vertex graph, after which the shortest remaining edge has to be placed disjointly with the first three, and so on. This immediately leads to the formula
$$
p(n)=\prod_{k=1}^{n-1}\frac{{2n-2k\choose 2}}{{2n-2k\choose 2}+2k(2n-2k)}\,,
$$
which after simple algebraic manipulations turns into your expression.
Probably, this is equivalent to what you know yourself, but I'm a bit perplexed about why you would want anything simpler or different.