I have an incomplete answer, but maybe you can finish it.
Let v' be the reflection of v over some axis. If you try to draw a vector s such that v = v' + s, using the triangle method to add the vectors v and s, you'll see that s is orthogonal to the reflection axis, and thus colinear to a.
The drawing should show an isosceles triangle, bisected by the reflection axis, where the two equally-sized sides are |v| and |v'|, the third side being of size |s|. From this triangle you should be able to calculate the magnitude of s as equal to 2 . |v| . cos theta, where "theta" is the angle between v (or v') and the normal to the reflection axis. (Try it on the drawing!)
Now, the expression 2 . |v| . cos theta is equal in magnitude to 2 . (va) / |a|. Multiplying this quantity by an unitary vector in the direction of a, namely a / |a|, you should get your vector s... or almost. There is a +/-1 factor not accounted yet for: has s the same direction of a or the opposite?
(This sign is the missing part that I was lazy to complete - I'd guess it has to do with v being or not on the same side of the axis than a; examining the sign of (va) should give a clue.)
Doing the multiplication above (and correcting the sign), you should get s = 2 (va) . a / |a|^2 = 2 (va) . a / (aa), and when replacing this expression for s into v' = v - s, you get your reflection formula.
Edit: Oh, I was pwned again. :)