There are likely more-elementary approaches, but here is an argument using some basic function-field arithmetic. Let's work a bit more generally than in the original question, considering a base field $F$, not of characteristic $3$. Let $R$ be the ring $F[x,y] / (x^3 + y^3 - 1)$, and let $L$ be the fraction field of $R$. Let's assume furthermore that $F$ contains an element $\omega$ that is a primitive cube root of $1$. We will show that $R$ is not a UFD by showing that the ideal $(x, y - 1)$ is not principal.
The ring $R$ has three evident units, $\eta_1 = x + y$, $\eta_2 = \omega x + y$, and $\eta_3 = \omega^2 x + y$, satisfying $\eta_1 \eta_2 \eta_3 = 1$. This rank-$2$ group of units complicates the arithmetic of $R$, as discussion in comments elsewhere noted. The field $L$ is something like a totally real cubic number field. (If $F$ did not contain a primitive cube root of $1$, then we would only have a rank-$1$ group of units, and the arithmetic would be simpler. The field $L$ would be something like a cubic number field with one real and one complex infinite place.)
The function field $L$ over $F$ has three infinite places (where the valuation of $x$ is negative), all of which have degree $1$. We can describe these places explicitly using the coordinates $u = 1/x$ and $v = y/x$. The infinite places are $P = (u, v + 1)$, $Q = (u, v + \omega)$, and $R = (u, v + \omega^2)$. Using these places, we find that the divisor of $\eta_1$ is $2P - Q - R$; that the divisor of $\eta_2$ is $-P + 2Q - R$; and that the divisor of $\eta_3$ is $-P-Q+2R$. As a consequence of these divisor calculations, we have $2P - Q \sim R$ and $3P \sim 3Q$, where $\sim$ denotes linear equivalence of divisors in the function field $L$ over $F$.
Let $T$ be a degree-$1$ place of $L$, and suppose that we have $T \sim aP + bQ + cR$ for some integers $a, b, c$, where $\sim$ denotes linear equivalence of divisors in $L$ as above. A bit of algebra with the divisors of the $\eta_i$ shows that $T$ is equivalent to one of $P, Q, R$. Indeed, using $-P + 2Q \sim R$, we can eliminate $R$ to get $T \sim a'P + b'Q$. Using $3P \sim 3Q$, we then have $T \sim a'' P + b''Q$, where $b''$ is $0, 1,$ or $-1$. Since the total degree of $a''P + b''Q$ is $1$, in the first case we have $T \sim P$. In the same way, the second case gives $T \sim Q$, and, using $2P - Q \sim R$, the third case gives $T \sim R$.
Let $\mathfrak{p}$ be a prime ideal of $R$ whose residue field is $F$, such as $(x, y - 1)$. We show that $\mathfrak{p}$ is not principal. Suppose for contradiction that $\mathfrak{p}$ were principal, generated by an element $f$ of $R$. The divisor of $f$ would then be $T - (aP + bQ + cR)$, where $T$ is the degree-$1$ place of $L$ corresponding to $\mathfrak{p}$, and where $a,b,c$ are suitable integers. By the argument above, we could then conclude that one of $T\sim P$, $T \sim Q$, or $T\sim R$ holds. Since $T$ is a finite place, it is distinct from the three infinite places, and so these linear equivalences are impossible, as the genus of $L$ is $1$.
(There are various ways to see that these linear equivalences are impossible. Suppose, for example, that $g$ were an element of $L$ with divisor $T - P$. Then the $1$-form $g \frac{dx}{y^2}$ would have a single singularity, a simple pole at $P$ with non-zero residue. But this is impossible, since the sum of residues of a meromorphic differential for $L$ must be $0$.)