I wanted to show that the running maximum say $\max_{t\in [0,1]}W_{t}$ has continuous distribution without taking help from the fact that it is absolutely continuous and has the distribution of $|W_{1}|$. To clarify, I want to just use the definition of Brownian motion and to prove this. (I am always talking about the standard brownian motion here).
Let us restrict to the set (of probability $1$ ) where $W_{t}$ is continuous.
Then I want to show that $P(\{\max_{t\in[0,1]}W_{t}=x\})=0$ for each $x\in\Bbb{R}$ .
So let us consider an enumeration of rationals $\{r_{n}\}_{n\in\Bbb{N}}$ in $[0,1]$ and define $Q_{n}=\{r_{1},...,r_{n}\}$ .
We have that $\max_{t\in[0,1]}W_{t}=\sup_{r\in\Bbb{Q}\cap[0,1]}W_{r}$ .
Also $\sup_{r\in\Bbb{Q}\cap[0,1]}W_{r}=x\implies W_{r}\leq x\,,\forall r\in\Bbb{Q}\cap[0,1]$ and for each $m\in\Bbb{N}$ , there exists $r\in\Bbb{Q}\cap[0,1]$ such that $W_{r}\in (x,x-\frac{1}{m}]$ .
So I want to write $\displaystyle P(\sup_{r\in[0,1]\cap\Bbb{Q}} W_{r}=x)\leq P\big(\bigcap_{m\in\Bbb{N}}\bigcup_{n\in\Bbb{N}}\{\exists r\in Q_{n}\,: W_{r}\in (x-\frac{1}{m},x]\}\big)=\lim_{m\to\infty}P(\bigcup_{n\in\Bbb{N}}\{\exists r\in Q_{n}: W_{r}\in(x-\frac{1}{m},x])$
$$\leq \lim_{m\to\infty}P\bigg(\bigcup_{n\in\Bbb{N}}\big\{\exists r\in Q_{n}:W_{r}\in (x-\frac{1}{m2^{n}},x]\big\}\bigg)\leq \lim_{m\to\infty}\lim_{n\to\infty}\sum_{k=1}^{n}P(W_{r_{k}}\in(x-\frac{1}{m2^{n}},x])$$
$$\leq \lim_{m\to\infty}\lim_{n\to\infty}\sum_{k=1}^{n}\frac{n}{m2^{n}}=0$$
However I am having trouble making this precise and I obviously am making a mistake in the step of considering $\lim_{m\to\infty}\lim_{n\to\infty}\sum_{k=1}^{n}P(W_{r_{k}}\in(x-\frac{1}{m2^{n}},x])$ . Can anyone tell me how do I make it correct?. Also is there an easier way to do this? (I mean without using the fact that it has a density).
Another way I tried to think of is that if I can show that $\max_{t\in[0,1]} W_{t}=\max_{r\in\Bbb{Q}\cap[0,1]}W_{r}$ (note that I mean $\max$ and not $\sup$) , then I can do it easily as for any finitely many points $r_{1},...,r_{n}$ , we have $P(\max_{r\in Q_{n}}W_{r}=x)\leq \sum_{k=1}^{n} P(W_{r_{k}}=x)=0$ for all $n$ . So I can proceed with the countable union and conclude what I want. But this is not necessarily true for continuous functions. I mean for example , $\sin(2x)$ achieves it's maximum at $x=\frac{\pi}{4}\in[0,1]$ and not at a rational point.
EDIT:
I also tried the following but I am unsure of it. Consider the sets $A_{n}=\bigcap_{m\in\Bbb{N}}\bigg\{\exists r\in Q_{n}:W_{r}\in(x-\frac{1}{m},x]\bigg\}$ . Then $A_{n}$'s increase to $\bigcap_{m\in\Bbb{N}}\bigg\{\exists r\in \Bbb{Q}\cap[0,1]:W_{r}\in(x-\frac{1}{m},x]\bigg\}$ .
But $\displaystyle P(A_{n})=P\bigg(\bigcap_{m\in\Bbb{N}}\bigg\{\bigcup_{r\in Q_{n}}W_{r}\in (x-\frac{1}{m},x]\bigg\}\bigg)=\lim_{m\to\infty}P\bigg(\bigcup_{r\in Q_{n}}\bigg\{W_{r}\in(x-\frac{1}{m},x]\bigg\}\bigg)$
(We could switch the intersection over $m$ to limit because for each fixed $n$, the sets $\{\exists r\in Q_{n}:W_{r}\in(x-\frac{1}{m},x]\}$ is a decreasing sequence)
$$\leq \lim_{m\to\infty} \sum_{k=1}^{n}\frac{1}{m}\frac{1}{\sqrt{r_{k}}}=0$$ (By using a simple bound for standard gaussian that $P(X\in(a,b))\leq (b-a)$ for $X\sim N(0,1))$
Also to note that the maximum is always $\geq 0$ a.s. . So we can ignore $r=0$ and consider $r\in\Bbb{Q}\cap(0,1]$ as $W_{0}=0$ a.s.
So $P(A_{n})=0\,,\forall n\in\Bbb{N}$ and hence $P(\bigcup_{n\in\Bbb{N}}A_{n})=0$ which is what is required.
We will prove (with our hands tied behind our back, as requested) that for each $x>0$, we have $P(\{\max_{t\in[0,1]}W_{t}=x\})=0$ (The case $x=0$ follows from time-inversion, for instance.)
It suffices to show that for each $\epsilon \in (0,1)$, the event $A^{\epsilon}_x:= \{\max_{t\in[\epsilon,1]}W_{t}=x\}$ has $P(A^{\epsilon}_x)=0$, since for $x>0$, $$\{\max_{t\in[0,1]}W_{t}=x\}=\cup_{m=2}^\infty A^{1/m}_x\,.$$ Let $H^{\epsilon}_z:=\{\max_{t\in[0,1-\epsilon]}W_{t}=z\}$ and let $\mu$ be the distribution of $W_\epsilon$. Then by independence of increments of Brownian motion, we have $$P(A^{\epsilon}_x)= \int_{\mathbb R} P(H^{\epsilon}_{x-y})\, d\mu(y) \,. \tag{*} $$
Since $P(H_{x-y}^\epsilon)$ could be positive for at most countably many $y$, and $\mu$ is a continuous measure, the integral in $(*)$ must vanish.