With the advent of location-based social networks, users can tag their daily activities in different locations through check-ins. These check-in locations signify user preferences for various socio-spatial activities and can be used to improve the quality of services in some applications such as recommendation systems, advertising, and group formation. To support such applications, in this paper, we formulate a new problem of <italic>identifying top-k</italic> <italic>S</italic><italic>ocio-</italic><italic>S</italic><italic>patial co-engaged</italic> <italic>L</italic><italic>ocation</italic> <italic>S<italic/>election</italic> (<italic>SSLS</italic>) for users in a social graph, that selects the best set of <inline-formula><tex-math notation="LaTeX">$k$</tex-math><alternatives><mml:math><mml:mi>k</mml:mi></mml:math><inline-graphic xlink:href="li-ieq2-3151095.gif"/></alternatives></inline-formula> locations from a large number of location candidates relating to the user and her friends. The selected locations should be (i) <italic>spatially and socially relevant</italic> to the user and her friends, and (ii) <italic>diversified both spatially and socially</italic> to maximize the coverage of friends in the socio-spatial space. This problem has been proved as NP-hard. To address such a challenging problem, we first develop an <monospace>Exact</monospace> solution by designing some pruning strategies based on derived bounds on diversity. To make the solution scalable for large datasets, we also develop an approximate solution by deriving relaxed bounds and advanced termination rules to filter out insignificant intermediate results. To further accelerate the efficiency, we present one fast exact approach and a meta-heuristic approximate approach by avoiding the repeated computation of diversity at the running time. Finally, we have performed extensive experiments to evaluate the performance of our proposed algorithms against three adapted existing methods using four large real-world datasets.