Kinetics of heterogeneous catalytic reactions are often complicated by various factors, and from the perspective of statistical physics the development of the corresponding models is frequently challenging especially in the case of practically important alloy catalysts. To extend the basics in this area, I use a generic kinetic model of N_{2} formation inthe NO reaction with such species as CO or H_{2}. The focus is onthe reaction on the surface of a bimetallic alloy with low integral fraction of one of the metals so that the alloy is formed only at the surface. The main goal is to illustrate the specifics of the reaction kinetics and to clarify the accuracy of various approximations which are often inevitable for the analysis of such systems. The key results are as follows. (i) In the baseline one-metal case with the Langmuir-type equations, the model predicts the existence of the optimal adsorbate binding energy corresponding to the maximal reaction rate. (ii) With suitable substitution of the symbols, the equations employed for a uniform surface [item (i)] can be used for the mean-field description of reaction on the surface of a random alloy. In particular, the maximum in the dependence of the reaction rate on the binding energy is converted to the maximum with respect to the alloy composition. (iii) The reaction kinetics on the random-alloy surface have been described exactly. The corresponding maximum in the reaction rate is lower and somewhat smeared compared to that predicted in the mean-field approximation for an alloy or for the optimal monometallic catalyst. (iv) The reaction kinetics have been scrutinized also in the case of two-dimensional (2D) segregation of metal atoms at the surface in the limits of slow and rapid adsorbate diffusion between the spots. The kinetics are shown to be qualitatively different in these limits and also compared to the random-alloy case. (v) The thermodynamic criteria for 2D segregation of metal atoms at the surface have been derived as well.