<p>When the size of a magneto-electro-elastic cylindrical shell (MEECS) is reduced to micro-/nano-scale, the size-dependent flexomagnetic effect (FME) and flexoelectric effect (FEE) significantly influence their multi-physical coupling behaviors. To investigate these effects on the post-buckling behaviors of an MEECS, a nonlinear post-buckling model is developed based on the higher-order shear deformation theory (HSDT) and magneto-electro-elastic (MEE) constitutive relations with the FME and FEE. The equilibrium path and the corresponding shell deformations are obtained with a set of newly developed generalized displacement functions within the framework of the Galerkin approach. These displacement functions are established based on the trigonometric series expansions, which accurately satisfy the clamped boundary conditions (BCs). The effects of geometry, flexomagnetic/flexoelectric coefficients, and external electromagnetic fields on the post-buckling behaviors of an MEECS with the FME and FEE are analyzed. Numerical results indicate that the FME decreases the upper critical load of an MEECS, whereas the FEE exhibits an opposite effect by increasing it.</p>