I just wanted to write down an answer that didn't use the Mayer-Vietoris sequence explicitly, although it is certainly just this underneath. I however found it useful to try to get some deeper insight on the inner workings of the theorem. The general lesson that learnt from this, which is of course a simple corollary of the Mayer-Vietoris sequence, is that if a $k$-form $\omega$ can be written as an exact form $\text{d}\alpha$ on $U$ and $\text{d}\beta$ on $V$, and $U\cap V$ is such that every closed $(k-1)$-form is exact, then $\omega$ can in fact be written as an exact form on $U\cup V$.
The procedure in order to do this is the following. First, note that the difference between the two expression of $\omega$ on $U\cap V$ is closed, since $\text{d}(\alpha-\beta)=\omega-\omega=0$. The assumption on $U\cap V$ implies then that this difference is exact. In other words, there is a $(k-2)$-form $\gamma$ such that $\alpha=\beta+\text{d}\gamma$ on the intersection. Of course, we would like this to extend to all of $U\cup V$. In order to do so, we choose a partition of unity $\{\rho_U,\rho_V\}$ subordinate to $\{U,V\}$. We can then define the $(k-2)$-forms $\rho_U\gamma$ on $V$ and $\rho_V\gamma$ on $U$. We can then slightly modify our $(k-1)$-forms to $\alpha'=\alpha-\text{d}(\rho_V\gamma)$ and $\beta'=\beta+\text{d}(\rho_U\gamma)$. They are still representations of $\omega$ on their respective domains since they differ from the original ones by closed forms but now have a very nice property: they agree on the intersection
$$\alpha'-\beta'=\alpha-\beta-\text{d}((\rho_V+\rho_U)\gamma)=\alpha-\beta-\text{d}\gamma=\alpha-\alpha=0.$$
They can thus be extended to a single $(k-1)$-form $\mu$. Moreover $\omega=\text{d}\mu$.
The intuition here is clear. The function $\rho_V$ vanishes far away from $V$, such that $\alpha'$ coincides with $\alpha$ there. However, as we are about to leave $U$ and enter $V$ it becomes $1$, such that $\alpha'=\beta$ there. In other words, $\alpha'$ smoothly interpolates between $\alpha$ and $\beta$. Moreover, it does so in a very special way which is where the hypothesis enters. Note that we could've still done this smoothing directly on the difference of our two forms $\rho_V\text{d}\gamma$. However, unlike $\text{d}\gamma$, this is not closed.
One thus only need to consider a cover of $\mathbb{R}^2\setminus\{p,q\}$ by two half-spaces, each surrounding only one of the removed points and intersecting on a strip, which would satisfy the conditions above $H^1(U\cap V)=0$. The proof that any 2-form is exact on each of these opens can be found in Ted Shiffrin's answer in how to compute the de Rham cohomology of the punctured plane just by Calculus?. Alternatively, one can further decompose each of these opens following Bruno Joyal's excellent answer.