  The Quantum Inverse Scattering Method is a scheme for solving integrable
models in $1+1$ dimensions, building on an $R$-matrix that satisfies the
Yang--Baxter equation and in terms of which one constructs a commuting family
of transfer matrices. In the standard formulation, this $R$-matrix acts on a
tensor product of vector spaces. Here, we relax this tensorial property and
develop a framework for describing and analysing integrable models based on
planar algebras, allowing non-separable \textit{$R$-operators} satisfying
\textit{generalised} Yang--Baxter equations. We also re-evaluate the notion of
integrals of motion and characterise when an (algebraic) \textit{transfer
operator} is polynomial in a single integral of motion. We refer to such models
as {\em polynomially integrable}. In an eight-vertex model, we demonstrate that
the corresponding transfer operator is polynomial in the natural hamiltonian.
In the Temperley--Lieb loop model with loop fugacity $\beta\in\mathbb{C}$, we
likewise find that, for all but finitely many $\beta$-values, the transfer
operator is polynomial in the usual hamiltonian element of the Temperley--Lieb
algebra $\mathrm{TL}_n(\beta)$, at least for $n\leq17$. Moreover, we find that
this model admits a second canonical hamiltonian, and that this hamiltonian
also acts as a polynomial integrability generator for small $n$ and all but
finitely many $\beta$-values.
