pta 7-10 几几几何

几几几何

题目描述：
按照逆时针绕向给出一个凸多边形的 n 个顶点 P0,P1,··· ,Pn−1，再给出凸多边形内部（含边界） 一点 Q。现在要将这个凸多边形在地上无滑动地滚动一周，初始时 P0P1 边与地面接触，假设当前是 PiP(i+1) mod n 边与地面接触，那么滚动一下之后则是 P(i+1) mod nP(i+2) mod n 边与地面接触。不难发现， 从初始状态滚动 n 下之后 P0P1 边再一次与地面接触，这时认为凸多边形已经滚动了一周。Keven现在想知道凸多边形滚动一周之后点 Q 经过的轨迹长度，于是请教了你，你能帮助他解决这个问题吗？

输入格式:
第一行包含一个整数 T（T < 10），表示测试数据的组数。 接下来依次描述 T 组测试数据。对于每组测试数据： 第一行包含一个整数 n，表示凸多边形的顶点数。 接下来 n 行，每行包含两个整数 xi,yi，按照逆时针的顺序给出凸多边形 n 个顶点 P1,P2,··· ,Pn 的坐标。 最后一行包含两个整数 xQ,yQ，表示点 Q 的坐标。 保证点 Q 在凸多边形内部（含边界）。

输出格式:
对于每组测试数据，输出一行信息 “Case #x: y”（不含引号），其中 x 表示这是第 x 组测试数据， y 表示凸多边形滚动一周之后点 Q 经过的轨迹长度，四舍五入精确到小数点后 3 位，数据保证答案的小 数点后第 4 位不是 4 或 5。

输入样例:
在这里给出一组输入。例如：

2
4
0 0
0 1
1 1
1 0
0 0
4
0 0
0 2
2 2
2 0
1 1
输出样例:
在这里给出相应的输出。例如：

Case #1: 5.363
Case #2: 8.886

首先，我们将每个顶点和Q点相连，我们通过观察发现：Q的轨迹可以被分为n段，每一段对应Q点绕点i（ i 从1到n）旋转的弧长，这一段距离用弧长公式l=r×a(弧长=半径×角度)可求，半径r用两点间距离公式易求，角度我们用acos（）函数求:
对于i点，i-1点，i+1点，我们可以求出i到i+1和i到i-1的两个向量n1,n2，通过公式cos(theta)=(n1×n2)/(|n1|×|n2|),再用acos（）函数求出theta即可，然后累加得答案。

附上AC代码：
#include<cstdio>
#include<cstring>
#include<queue>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxn=100000;
int T,n,a[maxn],b[maxn],X,Y;
double get_d(int x,int y)
{
	return sqrt((x-X)*(x-X)+(y-Y)*(y-Y));
}
double get(int x1,int y1,int x2,int y2,int tmp)
{
	double t1=(x1*x2+y1*y2);
	double t2=sqrt(x1*x1+y1*y1)*sqrt(x2*x2+y2*y2);
	double d=get_d(a[tmp],b[tmp]);
	double theta= acos(-t1/t2);
	return d*theta;
}
int main()
{
	scanf("%d",&T);
	for(int kk=1;kk<=T;kk++)
	{
		double ans=0;
		scanf("%d",&n);
		for(int i=1;i<=n;i++) scanf("%d%d",&a[i],&b[i]);
		scanf("%d%d",&X,&Y);
		ans+=get(a[n]-a[1],b[n]-b[1],a[2]-a[1],b[2]-b[1],1);
		ans+=get(a[n-1]-a[n],b[n-1]-b[n],a[1]-a[n],b[1]-b[n],n);
		for(int i=2;i<n;i++) 
			ans+=get(a[i+1]-a[i],b[i+1]-b[i],a[i-1]-a[i],b[i-1]-b[i],i);
		printf("Case #%d: %.3lf\n",kk,ans);
	}
	return 0;
 }