Let $u = 1 + \sqrt{x}$. Then $\sqrt{x} = u - 1$, so $x = (u-1)^2$.
Differentiating to find $dx$: $dx = 2(u-1) du$.
Now, change the limits of integration:
- When $x = 0$, $u = 1 + \sqrt{0} = 1$.
- When $x = 1$, $u = 1 + \sqrt{1} = 2$.
Substitute into the integral:
$$ \int_0^1 (1 + \sqrt{x})^8 dx = \int_1^2 u^8 \cdot 2(u-1) du = 2 \int_1^2 (u^9 - u^8) du $$
Now integrate:
$$ = 2 \left[ \frac{u^{10}}{10} - \frac{u^9}{9} \right]_1^2 $$
Evaluate at the limits:
$$ = 2 \left[ \left(\frac{2^{10}}{10} - \frac{2^9}{9}\right) - \left(\frac{1^{10}}{10} - \frac{1^9}{9}\right) \right] $$
$$ = 2 \left[ \left(\frac{1024}{10} - \frac{512}{9}\right) - \left(\frac{1}{10} - \frac{1}{9}\right) \right] $$
$$ = 2 \left[ \frac{1023}{10} - \frac{511}{9} \right] = 2 \left[ \frac{1023 \cdot 9 - 511 \cdot 10}{90} \right] $$
$$ = 2 \left[ \frac{9207 - 5110}{90} \right] = 2 \left[ \frac{4097}{90} \right] = \frac{4097}{45} $$