We investigate the effect of mixing of radioactive nickel (\(^{56}\)Ni) on the early-time color evolution of Type Ib and Ic supernovae (SNe Ib/Ic) using multi-group radiation hydrodynamics simulations. We consider both helium-rich and helium-poor progenitors. Mixing of \(^{56}\)Ni is parameterized using a Gaussian distribution function. We find that the early-time color evolution with a weak \(^{56}\)Ni mixing is characterized by three different phases: initial rapid reddening, blueward evolution due to the delayed effect of \(^{56}\)Ni heating, and redward evolution thereafter until the transition to the nebular phase. With a strong \(^{56}\)Ni mixing, the second phase disappears. We compare our models with the early-time color evolution of several SNe Ib/Ic (SN1999ex, SN 2008D, SN 2009jf, iPTF13bvn, SN 1994I, SN 2007gr, SN 2013ge, and 2017ein) and find signatures of relatively weak and strong \(^{56}\)Ni mixing for SNe Ib and SNe Ic, respectively. This suggests that SNe Ib progenitors are distinct from SN Ic progenitors in terms of helium content and that \(^{56}\)Ni mixing is generally stronger in the carbon-oxygen core and weaker in the helium-rich envelope. We conclude that the early-time color evolution is a powerful probe of \(^{56}\)Ni mixing in SNe Ib/Ic.