In the kinetic theory of the discrete Boltzmann models (DBMs) the standard similarity shock waves solutions are, like the Broadwell shock waves, solutions of scalar Riccati equations. The corresponding microscopic densities, as well as the macroscopic conservative quantities: mass, momentum and energy, are monotonic. In thefirst part of the paper we construct new exact similarity shockwavessolutions of systems of coupled Riccati equations, which arise for the DBMs with n conservation laws andn+kindependent densities (here withn= 3,k= 2). The former standard shock waves are rational functions of two linear polynomials with an exponential variable while for the new waves the polynomials are quadratic and consequently nonmonotonic behaviours are possible. We stress that from the mathematical point of view we enter into a new domain of exact self-similar solutions associated to the systems of coupled Riccati equations. Due to the fact that, for the considered models, there exists a linear differential relation between two microscopic densities we obtain two different classes of solutions. For one class these two densities are necessarily of the former standard type while for the other all densities are of the new type of self-similar solutions. Furthermore we prove an important result: Let us consider microscopic densities associated to opposite velocities along the shock-axis. To any solution, there exists another one, called partner solution, such that those densities are exchanged with opposite similarity variable. This means that from any solution with mass, energy and pressure given, there exists a partner solution where these macroscopic quantities have opposite similarity variable.