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)}