The fi rst untagged decay-time-integrated amplitude analysis of B 0 s ! K 0 S K decays is performed using a sample corresponding to 3: 0 fb of pp collision data recorded with the LHCb detector during 2011 and 2012. The data are described with an amplitude model that contains contributions from the intermediate resonances K 9892) 0;+, K 2 91430) 0;+ and K 0 91430) 0;+, and their charge conjugates. Measurements of the branching fractions of the decay modes B 0 s ! K 9892) K and B 0 s ! K 9892) 0 K 0 are in agreement with, and more precise than, previous results. The decays B 0 s ! K 0 91430) K and B 0 s ! K 0 91430) 0 K 0 are observed for the fi rst time, each with signi fi cance over 10 standard deviations.