Preciso entender como esse código funciona:

#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";
}

Mas parece que não consigo entender o comportamento dessa função:

#define foo1( a ) (a * a)

e, portanto, não entendo a saída dessas duas chamadas de função:

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

Esta é a saída do programa:

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

Não vejo por que foo1 não está se comportando como foo2. Alguém me explicaria como essa macro #define funciona?

NOTA: Não devo alterar o código, estou apenas tentando entender a saída.