Branching fractions and CP-violating asymmetries of charmless B 0 s → P P P decays with P = K, π are studied using a simple model based on the factorization approach. The penguin-dominated modes B 0 s → K 0 K + π − and K 0 K − π + have the largest rates among the three-body B s decays and they are dominated by the K * 0 (1430) resonances and nonresonant contributions. The branching fraction of B 0 s → K 0 K + K − is predicted to be of order 1.4 × 10 −6 , which is in the middle of the measured result (0.2 − 3.4) × 10 −6 obtained by LHCb. We find sizable CP violation in K 0 π + π − , K 0 π 0 π 0 , K 0 K + K − and K S K S K S channels. Just as the B u sector, the CP asymmetries of B 0 s → K 0 π + π − and B 0 s → K 0 K + K − have similar magnitudes but are opposite in sign. Several modelindependent relations based on U -spin symmetry for CP asymmetries are derived. Although they cannot be tested by the present available data, they can be checked by the dynamical calculations. Since the U -spin symmetry which relates various matrix elements of scalar densities, for example, K 0 K|dq|0 and K 0 π|sq|0 , is badly broken, the derived U -spin relations are generally not well respected.