计算最终结果时强制类型转化为 long long

ans += (p[j+1] - p[j]) * (q[i+1] - q[i]);

改为

ans += (long long)(p[j+1] - p[j]) * (long long)(q[i+1] - q[i]);

其中 $p$, $q$ 分别存储离散化的 $x$, $y$ 坐标