Solution (source code)

= Solution

With $r_B^2=\hbar/(qB)$, direct substitution gives $[a,a^\dagger]=[b,b^\dagger]=1$.