【洛谷 P2785】 物理1（phsic1）- 磁通量【计算几何】【叉积】

解题思路

关于计算多边形面积，可以用叉积，因为叉积的几何意义就是一个四边形的面积。

设矢量$P=(x1,y1)，Q=(x2,y2)$，叉积为：$P\times Q=x1\ast y2-x2\ast y1$，其结果是一个标量，叉积的绝对值是一个平行四边形面积$

还有一个非常重要性质是：叉积可以通过它的符号（正负）判断两矢量相互之间的顺逆时针关系

即：若 $P\times Q>0$ , 则P在Q的顺时针方向。
　　若 $P\times Q<0$ , 则P在Q的逆时针方向。
　　若 $P\times Q=0$ , 则P与Q共线，但可能同向也可能反向。

最后，了解了这些再来计算多边形面积就很简单了：
拿下面这个图形举例子：

$S正方形=Soab+Sobc+Socd+Soda$ ，4个三角形中后3个是顺时针，后3个的绝对值的面积S$obcda$，但是Soab为逆时针，符号与后三个相反，4个相加刚好把多余的抵消。

我们得出规律：对多边形$P1P2P3...Pn$(顺或逆时针拐都可以)，设平面上有任意的一点P，则有：
$S(P1,P2,P3，...,Pn)=0.5\ast abs(S(P,P1,P2)+S(P,P2,P3)+、、、+S(P,Pn,P1))$

$P$是可以取任意的一点，用$(0，0)$时就有
$S=0.5\ast abs(x1\ast y2-y1\ast x2+x2\ast y3-y2\ast x3+...+xn\ast y1-yn\ast x1)$

最后这道题注意下坐标右负的就不好用（0,0）取第一个坐标当原点就好了

代码

#include<iostream>
#include<cstdio>
#include<iomanip>
#include<cmath>
#include<cstring>
#define ll long long
using namespace std;

int n;
double B,ans,x[100010],y[100010];

double work(int i,int j) {
	return (x[i]-x[1])*(y[j]-y[1])-(x[j]-x[1])*(y[i]-y[1]);
}

int main() {
	scanf("%d%lf",&n,&B);
	for(int i=1; i<=n; i++) {
		scanf("%lf%lf",&x[i],&y[i]);
	}
	for(int i=2; i<=n; i++) {
		double w=work(i,i-1);
		ans=ans+w;
	}
	ans=ans+work(n,1);
	printf("%.4lf",abs(ans)/2*B);
}