HDU--1080--计算几何--精度问题

这次xhd面临的问题是这样的：在一个平面内有两个点，求两个点分别和原点的连线的夹角的大小。

注：夹角的范围[0，180]，两个点不会在圆心出现。

Input

输入数据的第一行是一个数据T，表示有T组数据。
每组数据有四个实数x1,y1,x2,y2分别表示两个点的坐标,这些实数的范围是[-10000,10000]。

Output

对于每组输入数据，输出夹角的大小精确到小数点后两位。

Sample Input
2
1 1 2 2
1 1 1 0
Sample Output
0.00
45.00

思路：方法不一，我采用余弦定理；

注意：1.double y1；这里标识符编译通不过，改为了yy1;

2.实例二中算出的答案应该为acos(-1)/4,所以采用余弦定理计算出的答案应该ans/acso(-1.0)*180;

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<stack>
using namespace std; 
double x1,yy1;
double x2,y2;
double pf(double x){
	return x*x;
}
int main(){
	int t;
	scanf("%d",&t);
	while(t--){
		scanf("%lf%lf%lf%lf",&x1,&yy1,&x2,&y2);
		double a=sqrt(pf(x1)+pf(yy1)),b=sqrt(pf(x2)+pf(y2)),c=sqrt(pf(x1-x2)+pf(yy1-y2));
		double ans=acos((pf(a)+pf(b)-pf(c))/(2.0*a*b))/acos(-1.0)*180;
		printf("%.2lf\n",ans);
	} 
	return 0;
}