[step:Apply the chain rule along every regular forward orbit]Fix $x \in X$. First observe that the derivative determined by a local representative is independent of the representative. If $F,G$ are two $C^1$ representatives near $y\in I\setminus E$, then $H:=F-G$ is $C^1$ and vanishes on a subset of $I$ accumulating at $y$, because $I$ is a nondegenerate interval. Hence the difference quotients for $H$ along that subset are zero, and the ordinary derivative $H'(y)$ must be zero. Thus $F'(y)=G'(y)$, including at endpoints of $I$.
Next we record the local composition fact needed for the ordinary chain rule. Fix $n \in \mathbb{N}$. For each $k=0,\dots,n-1$, the point $f^k(x)$ lies in $I\setminus E$, so choose a local $C^1$ representative $F_k:U_k\to\mathbb{R}$ of $f$ near $f^k(x)$. By continuity of these finitely many maps, we may shrink open intervals $V_k\subset U_k$ around $f^k(x)$, for $k=0,\dots,n-1$, and choose an open interval $V_n$ around $f^n(x)$ so that
\begin{align*}
F_k(V_k)\subset V_{k+1}
\end{align*}
for every $k=0,\dots,n-1$. Then
\begin{align*}
G_n:=F_{n-1}\circ\cdots\circ F_0
\end{align*}
is a $C^1$ map on $V_0$. After possibly shrinking $V_0$ further along this finite construction, every $z\in V_0\cap I$ has its first $n$ iterates in the corresponding $V_k$, and the equality $G_n(z)=f^n(z)$ follows by induction from $F_k=f$ on $V_k\cap I$. Thus $G_n$ is a local $C^1$ representative of the finite iterate at $x$, and the derivative $(f^n)'(x)$ exists.
We claim that, for every $n \in \mathbb{N}$, this derivative satisfies
\begin{align*}
(f^n)'(x) = \prod_{j=0}^{n-1} f'(f^j(x)).
\end{align*}
For $n=1$, this is the identity $(f^1)'(x)=f'(x)$, and $x \notin E$ because $x \in X_{\mathrm{reg}}$. The repaired hypothesis gives a local $C^1$ representative of $f$ near $x$, and the uniqueness argument above makes the ordinary derivative well-defined.
Suppose the formula holds for some $n \in \mathbb{N}$. The local composition fact applied to $n$ gives a local $C^1$ representative of $f^n$ near $x$, and the local representative hypothesis gives a local $C^1$ representative of $f$ near $f^n(x)$. Thus the ordinary chain rule applies to the local representatives of the composition $f^{n+1}=f \circ f^n$ at $x$, and it gives
\begin{align*}
(f^{n+1})'(x) = f'(f^n(x))(f^n)'(x).
\end{align*}
Substituting the induction hypothesis yields
\begin{align*}
(f^{n+1})'(x) = f'(f^n(x))\prod_{j=0}^{n-1} f'(f^j(x)) = \prod_{j=0}^{n} f'(f^j(x)).
\end{align*}
By induction, the product formula holds for every $n \in \mathbb{N}$.[/step]