We present K-band imaging of two \sim 30 ^ { \prime } \times 30 ^ { \prime } fields covered by the MUSYC Wide NIR Survey . The 1030 and 1255 fields were imaged with ISPI on the 4m Blanco telescope at CTIO to a 5 \sigma point-source limiting depth of K \sim 20 ( Vega ) . Combining this data with the MUSYC Optical UBVRIz imaging , we created multi-band K-selected source catalogs for both fields . These catalogs , together with the MUSYC K-band catalog of the ECDF-S field , were used to select K < 20 BzK galaxies over an area of 0.71 deg ^ { 2 } . This is the largest area ever surveyed for BzK galaxies . We present number counts , redshift distributions and stellar masses for our sample of 3261 BzK galaxies ( 2502 star-forming ( sBzK ) and 759 passively evolving ( pBzK ) ) , as well as reddening and star formation rate estimates for the star-forming BzK systems . We also present 2-point angular correlation functions and spatial correlation lengths for both sBzK and pBzK galaxies and show that previous estimates of the correlation function of these galaxies were affected by cosmic variance due to the small areas surveyed . We have measured correlation lengths r _ { 0 } of 8.89 \pm 2.03 Mpc and 10.82 \pm 1.72 Mpc for sBzK and pBzK galaxies respectively . This is the first reported measurement of the spatial correlation function of passive BzK galaxies . In the \Lambda CDM scenario of galaxy formation , these correlation lengths at z \sim 2 translate into minimum masses of \sim 4 \times 10 ^ { 12 } M _ { \odot } and \sim 9 \times 10 ^ { 12 } M _ { \odot } for the dark matter ( DM ) halos hosting sBzK and pBzK galaxies respectively . The clustering properties of the galaxies in our sample are consistent with them being the descendants of bright LBG at z \sim 3 , and the progenitors of present-day > 1 L ^ { * } galaxies .