Measurements of anisotropic flow coefficients (v n) and their cross-correlations using two-and multiparticle cumulant methods are reported in collisions of pp at ffiffi ffi s p ¼ 13 TeV, p-Pb at a center-ofmass energy per nucleon pair ffiffiffiffiffiffiffi ffi s NN p ¼ 5.02 TeV, Xe-Xe at ffiffiffiffiffiffiffi ffi s NN p ¼ 5.44 TeV, and Pb-Pb at ffiffiffiffiffiffiffi ffi s NN p ¼ 5.02 TeV recorded with the ALICE detector. The multiplicity dependence of v n is studied in a very wide range from 20 to 3000 particles produced in the midrapidity region jηj < 0.8 for the transverse momentum range 0.2 < p T < 3.0 GeV=c. An ordering of the coefficients v 2 > v 3 > v 4 is found in pp and p-Pb collisions, similar to that seen in large collision systems, while a weak v 2 multiplicity dependence is observed relative to nucleus-nucleus collisions in the same multiplicity range. Using a novel subevent method, v 2 measured with four-particle cumulants is found to be compatible with that from six-particle cumulants in pp and p-Pb collisions. The magnitude of the correlation between v 2 n and v 2 m , evaluated with the symmetric cumulants SCðm; nÞ is observed to be positive at all multiplicities for v 2 and v 4 , while for v 2 and v 3 it is negative and changes sign for multiplicities below 100, which may indicate a different v n fluctuation pattern in this multiplicity range. The observed long-range multiparticle azimuthal correlations in high multiplicity pp and p-Pb collisions can neither be described by PYTHIA 8 nor by impact-parameter-Glasma, MUSIC, and ultrarelativistic quantum molecular dynamics model calculations, and hence, provide new insights into the understanding of collective effects in small collision systems.