A GPU-based Large-scale Monte Carlo Simulation Method for Systems with Long-range Interactions