用鞋带公式计算std::vector存储的多边形面积,需确保顶点按顺/逆时针有序、不重复,使用double或long long避免溢出和精度错误;自交多边形需先分解为简单多边形再分别计算。

用 std::vector 存多边形顶点时,怎么算面积?
直接用「鞋带公式」(Shoelace formula),它不依赖凸凹性,只要顶点按顺时针或逆时针顺序排列就行。关键不是写多复杂,而是顺序不能错、坐标类型别溢出。
实操建议:
- 顶点存成
std::vector<:pair double>></:pair>或自定义Point结构,避免用float算面积——小数精度不够容易符号翻转 - 确保首尾不重复:输入是
[(0,0), (2,0), (2,2), (0,2)]就够了,别加个重复的(0,0)进去,否则公式里会多算一轮 - 公式本质是求和
sum += x[i] * y[i+1] - x[i+1] * y[i],最后取绝对值除以 2;下标i+1要模n,C++ 里用(i + 1) % n
double polygonArea(const std::vector<:pair double>>& pts) {
int n = pts.size();
if (n <h3>遇到自相交多边形,<code>polygonArea</code> 还准吗?</h3>
<p>不准。鞋带公式算的是「有向面积」代数和,自交会导致部分区域正负抵消。比如一个八字形,结果可能接近 0,而不是两个环的面积之和。</p>
<p>常见错误现象:</p>
<ul>
<li>明明画出来是个大矩形加个小三角,面积却比矩形还小</li>
<li>同一组点,只是输入顺序调换(比如把中间某个点提前),结果差几倍</li>
</ul>
<p>这时候不能硬套公式。得先做「多边形分解」:用算法(如单调链剖分或 ear clipping)拆成若干简单多边形,再逐个调用上面的函数。OpenCV 的 <code>cv::contourArea</code> 默认也假设输入是简单多边形;若传入自交轮廓,文档明确说结果未定义。</p><div class="aritcle_card flexRow artxards">
<div class="artcardd flexRow">
<a class="aritcle_card_img" rel="nofollow" href="/xiazai/skill4025" title="C++ 算法竞赛自动化测试数据生成与校验框架"><img
src="https://img.php.cn/upload/skill/000/000/081/178988956499722.jpg" alt="C++ 算法竞赛自动化测试数据生成与校验框架" onerror="this.onerror='';this.src='/static/lhimages/moren/morentu.png'" ></a>
<div class="aritcle_card_info flexColumn">
<a rel="nofollow" href="/xiazai/skill4025" title="C++ 算法竞赛自动化测试数据生成与校验框架" class="overflowclass">C++ 算法竞赛自动化测试数据生成与校验框架</a>
<p class="overflowclass">根据原题生成新题面、验证器及完整测试数据,自动套用 testlib 模板,用于用户要求生成测试数据时。</p>
</div>
<a rel="nofollow" href="/xiazai/skill4025" title="C++ 算法竞赛自动化测试数据生成与校验框架" class="aritcle_card_btn flexRow flexcenter"><b></b><span>下载</span>
</a>
</div>
</div>
<h3>用 <code>int</code> 坐标时,为什么面积算出来是 0?</h3>
<p>因为整数乘法中间结果溢出,或者除以 2 时截断。比如顶点都是 <code>int</code>,但 <code>x[i] * y[i+1]</code> 可能远超 <code>INT_MAX</code>,一溢出整个和就崩了。</p>
<p>解决办法很实在:</p>
<ul>
<li>强制转 <code>long long</code> 再算:写成 <code>(long long)pts[i].x * pts[j].y</code>
</li>
<li>面积结果用 <code>double</code> 或 <code>long long</code> 接,别用 <code>int</code>
</li>
<li>如果确定所有坐标 ≤ 1e4,且顶点数 ≤ 1e3,那 <code>long long</code> 足够;否则老实用 <code>double</code>,它有 52 位有效位,比 <code>int</code> 安全得多</li>
</ul>
<h3>OpenCV 的 <code>cv::contourArea</code> 和手写函数结果不一样</h3>
<p>大概率是坐标顺序或数据类型问题。OpenCV 默认把输入 contour 当作逆时针为正方向,但它的实现底层也是鞋带公式——只是做了符号归一化:返回值恒为正,且内部用了 <code>double</code> 精度累积。</p>
<p>对比时要注意:</p>
<ul>
<li>OpenCV 输入是 <code>std::vector<:point></:point></code>,其中 <code>cv::Point</code> 是 <code>int</code> 类型,但函数内部会转 <code>double</code> 计算</li>
<li>你的手写函数如果用 <code>int</code> 累加,没转 <code>long long</code>,哪怕只差一个点,中间溢出就会导致最终差几百甚至符号相反</li>
<li>OpenCV 对空 contour 或少于 3 点返回 0;你自己的函数也要加同样校验,别让 <code>n=2</code> 时还硬算</li>
</ul>
<p>真要对齐结果,最稳的方式是:统一用 <code>std::vector<:point2d></:point2d></code>(即 double 坐标)喂给 OpenCV,再拿它的输出和你手写的 double 版本比。</p></:pair>C++免费学习笔记(深入):立即使用
在学习笔记中,你将探索 C++ 的入门与实战技巧!










