Necesito entender cómo funciona este código:

#define foo1( a ) (a * a)              // How does this work?
inline int foo2( int a ) { return (a * a); }
int goo1( int x ) { return foo1(foo2(x)); }
int goo2( int& x ) { return x = foo2(foo1(x)); }
int goo3( int& x, int y ) { return foo2(foo1(y + 1)); }
void goo4( int& x, int y ) { x = foo1(foo2(y + 1)); }

int main(){
int i = 2, j = 1, a = 2+3;

cout << "foo1 = " << foo1( 1+2 ) << "\n";  // How does this work?
cout << "foo2 = " << foo2( 2 + 1 ) << "\n";
cout << "goo1 = " << goo1( i ) << "\n";
cout << "goo2 = " << goo2( j ) << "\n";
cout << "goo3 = " << goo3( i, j ) << "\n";  // How does this work?
goo4( i, j );
cout << " i = " << i << " j = " << j << "\n";
}

Pero no parece que pueda entender el comportamiento de esta función:

#define foo1( a ) (a * a)

y por lo tanto no entiendo el resultado de estas dos llamadas a funciones:

foo1( 1+2 )
goo3( i, j )

Este es el resultado del programa:

foo1 = 5
foo2 = 9
goo1 = 16
goo2 = 1
goo3 = 9
 i = 16 j = 1

No puedo ver por qué foo1 no se comporta como foo2. ¿Alguien me explicaría cómo funciona esta macro #define?

NOTA: No debo cambiar el código, solo estoy tratando de entender el resultado.