We first present a catalogue of photometric redshifts for 14.68 million galaxies derived from the 7-band photometric data of Hyper Suprime-Cam Subaru Strategic Program and the Wide-field Infrared Survey Explorer using the nearest-neighbour algorithm. The redshift uncertainty is about 0.024 for galaxies of z<0.7, and steadily increases with redshift to about 0.11 at z~2. From such a large data set, we identify 21,661 clusters of galaxies, among which 5537 clusters have redshifts z>1 and 642 clusters have z>1.5, significantly enlarging the high redshift sample of galaxy clusters. Cluster richness and mass are estimated, and these clusters have an equivalent mass of M_{500} > 0.7*10^{14} Msun. We find that the stellar mass of the brightest cluster galaxies (BCGs) in each richness bin does not significantly evolve with redshift. The fraction of star-forming BCGs increases with redshift, but does not depend on cluster mass.