$4\mid 100\:$ & $\:4\mid \color{#90f}{3^{333}}\!+1\,$ (by $\bmod 4\!:\ 3^{333}\!+1\equiv (-1)^{333}\!+1\equiv -1+1\equiv 0),\,$ so
using $\ ab\bmod ac = a\,(b\bmod c) = $ mod Distributive Law to factor out $\,a = 4$
$3^{\large 333}\!+\!1\bmod 100 = 4\left[\dfrac{\color{#90f}{3^{\large 333}}\!+\!1}4\bmod 25\right] = 4\left[\dfrac{\color{#90f}{23}\!+\!1}4\bmod 25\right] = 4[6],\ $ by
$\!\bmod 25\!:\,\ \color{#c00}{3^{\large 20}\equiv 1}\ \Rightarrow\ \color{#90f}{3^{\large 333}}\equiv \underbrace{3(\color{#0a0}{3^{\large 3}})^{\large 4}(\color{#c00}{3^{\large 20}})^{\large 16}}_{\textstyle 3\,(\color{#0a0}2)^{\large 4}\,\color{#c00}1^{\large 16}}\equiv \color{#90f}{23},\,$ by $\rm\color{#c00}{Euler\ \phi(5^2)=4\cdot 5}$
Remark $ $ The computation of $\,\color{#90f}{3^{333}}\bmod 25\,$ in the prior line is a special case of the method of modular order reduction - the key to simplifying the computation of large modular powers.
As explained in the first link, the mod Distributive Law is an operational form of CRT. We can use the standard CRT formula along with the above computation mod $25$ to obtain the same result. In fact already $23\equiv -1\pmod{\!4}$ so we're done - we don't even need to apply the CRT formula.
Another quick way: $\, 3^{\large 333} = 3(-1+10)^{\large 166} =\,\ldots\,$ so expanding by Binomial Theorem shows only the first $2$ terms survive $\!\bmod 100 = 10^2$.