let $f=\sum_{n=0}^\infty a_nx^n$ and $g =\sum_{n=0}^\infty b_nx^n$. Suppose the gcd of $a_n$ and $b_n$ are both 1. Then show the gcd of the coefficients of $fg$ is also one.
I know that since R is in particular a domain, then $R[[x]]$ is also domain. However, I'm not really sure how to proceed further.
I like the proof of @JohnBrevik fine, but here’s my version, which you may find more visual:
For a prime $p$ of $R$, the ring $R/(p)$ is an integral domain. If we have a series $f\in R[[x]]$ let $\tilde f$ be the corresponding series in $R/(p)$. In case $p$ does not divide all coefficients of $f$, the series $\tilde f$ will have an initial degree $n$ so that $\tilde f(x)=x^n\varphi(x)$ with $\varphi(0)\ne0$, i.e. $\varphi$ has nonzero constant term.
Now under your hypotheses on $f$ and $g$, if $p$ is any prime of $R$, then we have $\tilde f(x)=x^n\varphi(x)$ and $\tilde g=x^m\psi(x)$ where as above, $\varphi(x)$ and $\psi(x)$ have nonzero constant terms. Since $R/(p)$ is an integral domain, $\varphi\psi$ will also have nonzero constant term. And $\widetilde{fg}(x)=\tilde f(x)\tilde g(x)=x^{m+n}\varphi(x)\psi(x)$, whose $x^{m+n}$-coefficient is nonzero. Thus $fg$ has a coefficient indivisible by $p$. Do this for all $p$, and get your result.