Solving large p-median problems by a multistage hybrid approach using demand points aggregation and variable neighbourhood search.