Bueno, se podría simplificar como sigue:\begin{align}
\sum_{{\Large 1 \leq a,b \leq L} \atop {\Large \gcd(a,b)=1}}{(L+1-a)(L+1-b)}& =L^2+2 \cdot \left(\sum_{{\Large 1 \leq b<a \leq L} \atop {\Large \gcd(a,b)=1}}{(L+1-a)(L+1-b)}\right) \\
& =L^2+2 \cdot \left(\sum_{\Large 2 \leq a \leq L}{(L+1-a)\sum_{{\Large 1 \leq b<a} \atop {\Large \gcd(a,b)=1}}{(L+1-b)}}\right) \\
& =L^2+2 \cdot \left(\sum_{\Large 2 \leq a \leq L}{(L+1-a)\left((L+1)\varphi{(a)}-\frac{a\varphi{(a)}}{2}\right)}\right) \\
& =L^2+2\cdot \left(\sum_{\Large 2 \leq a \leq L}{(L+1-a)\left(L+1-\frac{a}{2}\right)\varphi{(a)}} \right)
\end {Alinee el}
Ahora basta con sumar $a$. ¿Es que lo suficientemente rápido como para usted?
Edición: Una forma simple de calcular de forma recursiva sería $\varphi{(a)}$ $2 \leq a \leq L$ generar la secuencias de Farey $F_n$ hasta $n=L$ y tomar la diferencia del número de términos de secuencias sucesivas en el camino: $|F_a|-|F_{a-1}|=\varphi{(a)}$