To efficiently explore the densest structures, we develop an unbiased random search approach based on both a piling up method to generate initial structures in an unbiased way and an iterative balance method to optimize the volume of unit cell while keeping the overlap of hard spheres minimized.