The halo occupation distribution (HOD) describes the relation between galaxies and dark matter at the level of individual dark matter halos. The properties of galaxies residing at the centers of halos differ from those of satellite galaxies because of differences in their formation histories. Using a smoothed particle hydrodynamics (SPH) simulation and a semi-analytic (SA) galaxy formation model, we examine the separate contributions of central and satellite galaxies to the HOD, more specifically to the probability P(N|M) that a halo of virial mass M contains N galaxies of a particular class. In agreement with earlier results for dark matter subhalos, we find that the mean occupation function <N> for galaxies above a baryonic mass threshold can be approximated by a step function for central galaxies plus a power law for satellites, and that the distribution of satellite numbers is close to Poisson at fixed halo mass. For galaxy samples defined by different baryonic mass thresholds, there is a nearly linear relation between the minimum halo mass Mmin required to host a central galaxy and the mass M1 at which an average halo hosts one satellite, with M1 ~ 14 Mmin (SPH) or M1 ~ 18 Mmin (SA). The mean occupation number of young galaxies exhibits a local minimum at M ~ 10 Mmin where halos are too massive to host a young central galaxy but not massive enough to host satellites. We show that the conditional galaxy mass function at fixed halo mass cannot be described by a Schechter function because central galaxies produce a bump at high masses. We suggest parameterizations for the HOD and the conditional luminosity function that can be used to model observed galaxy clustering. Many of our predictions are in good agreement with recent results inferred from clustering in the Sloan Digital Sky Survey.