Day 21 is my favourite puzzle so far. Part 2 has a nice solution by observing that the number f(i) of reachable fields after i iterations is a quadratic function of the number of iterations i for i large enough, i.e. `f(i)=a*i^2+b*i+c`. The coefficients a,b,c can be derived from three values of f obtained by simulation of the process for few iterations.

