
We present a measurement of baryon acoustic oscillations (BAO) from Lyman-$\alpha$ (Ly$\alpha$) absorption and quasars at an effective redshift $z = 2.33$ using the complete extended Baryon Oscillation Spectroscopic Survey (eBOSS). The sixteenth and final eBOSS data release (SDSS DR16) contains all data from eBOSS and its predecessor, the Baryon Oscillation Spectroscopic Survey (BOSS), providing 210{,}005 quasars with redshift $z_q > 2.10$ used to measure Ly$\alpha$ absorption. We measure the BAO scale both from the auto-correlation of Ly$\alpha$ absorption and from its cross-correlation with 341{,}468 quasars at redshift $z_q > 1.77$. In addition to the statistical gain from the increased quasar sample and deeper observations, the main improvements over previous analyses arise from more accurate modeling of physical and instrumental correlations and from the use of new sets of mock data. Combining the BAO measurements from the auto- and cross-correlations yields constraints on the distance ratios $D_H(z = 2.33)/r_d = 8.99 \pm 0.19$ and $D_M(z = 2.33)/r_d = 37.5 \pm 1.1$, where the quoted uncertainties are statistical. These results are within $1.5\sigma$ of the prediction of the flat $\Lambda$CDM cosmology inferred from \textit{Planck} (2016). The analysis code \texttt{picca}, the catalog of flux transmission field measurements, and the $\Delta\chi^2$ surfaces are publicly available.