Low-lying equilibrium geometric structures of GaP (n) (n = 2-12) clusters obtained by an all-electron linear combination of atomic orbital approach, within spin-polarized density functional theory, are reported. The binding energy, dissociation energy, and stability of these clusters are studied within the local spin density approximation (LSDA) and the three-parameter hybrid generalized gradient approximation (GGA) due to Becke-Lee-Yang-Parr (B3LYP). Ionization potentials, electron affinities, hardness, and static dipole polarizabilities are calculated for the ground-state structures within the GGA. It is observed that the gallium atoms of the symmetric ground-state structures prefer to occupy the peripheral positions. It is found that the relative ordering of the isomers is influenced by the nonlocal exchange-correlation effects for small clusters. Generalized gradient approximation extends bond lengths and widens the energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), as compared to the LSDA gap. The odd-even oscillations in the dissociation energy, the second differences in energy, the HOMO-LUMO gaps, the ionization potential, the electron affinity, and the hardness are more pronounced within the GGA. The stability analysis based on the energies clearly shows the GaP5 and GaP7 clusters to be endowed with special stabilities.