实现 solution(S0, K, r, sigma, T, n_paths, n_steps, seed) -> tuple[float, float]。在风险中性 GBM 下用蒙特卡洛对终值欧式看涨期权定价，并用 几何均值亚式期权 作为控制变量来缩小标准误。返回 (price_estimate, standard_error)。