Consider the standard one-dimensional Brownian motion $B_t$. Set $$X_t = \int_0^t\mathbf{1}_{[0, 1]}(B_s)ds.$$ Here $\mathbf{1}_{[0, 1]}(\cdot)$ is the indicator function of the interval $[0, 1]$. $X_t$ denotes the total time spent by Brownian motion in $[0, 1]$ before time $t$, which can also be written using the Brownian local time as $$X_t = \int_0^1L^a(t)da.$$
Question:
Can we compute the distribution of $X_t$? If we can't get the specific expression of its density, is there any method to get the decay estimate (lower and upper bound) of its density?