Cow Uncle 学习了叉积的一点运用，叉积真的不错

Cow Uncle

Time Limit: 4000/2000MS (Java/Others) Memory Limit: 128000/64000KB (Java/Others)

SubmitStatus

Problem Description

South China Algorithm University (SCAU) 是一个生态环境优秀的校园。走在校道上，你可以看见牛羊猫狗鸡鸭鹅，等等。

牛在校道上走当然不是没人管的。放牛大叔通常会带N头大牛小牛去到一片宽阔的草地上吃草，草地上有M块大石头，放牛大叔会在石头上坐着，看着这群大牛小牛。当然，并非所有石头的位置都是那么好，放牛大叔要看着所有的牛，所以挑选的位置P的视野必须小于180度（也就是在他面前，最左边的牛A和最右边的牛B与P形成的角度∠APB不能大于等于180度）。

现在给出N头大牛小牛的位置，然后给出M个石头的位置，请你分别求出每个石头是否可以看着全部牛。如果可以，求出相应位置的视野度数。（可以假设大牛小牛不会在石头上吃草，而且大牛小牛都很聪明，不会在吃同一位置的草）

Input

有多组数据输入。

每组数据第一个数是牛的个数N，之后N行每行有两个数(Xi, Yi)，表示牛的位置。

紧接着的是石头的个数M，之后M行有每行有两个数(Xj, Yj)，表示石头的位置。

数据范围：

3≤M,N≤1000

|Xi|,|Yi|≤1000000

输入数据保证不会全部牛都在同一条直线上。

Output

对于每一块石头输出一行。如果这个位置的视野可以看到全部的牛，那么输出这个位置的视野的度数是多少（保留两位小数）。否则，输出“Bad Position”。

Sample Input

4
0 0
1 0
1 1
0 1
3
0.5 0.5
1 0.5
2 0

Sample Output

Bad Position
Bad Position
45.00

 #include<iostream>
 #include<algorithm>
 #include<math.h>
 #include<stdio.h>
 using namespace std;
 #define PI 3.1415926535898
 struct point
 {
     double x,y;
 };
 point p[],res[];
 int n,top;
 double Dist(const point &arg1, const point &arg2)
 {
     return sqrt( 1.0*(arg1.x - arg2.x)*(arg1.x - arg2.x) + (arg1.y - arg2.y)*(arg1.y - arg2.y) );
 }
 bool multi(point p0,point p1,point p2)
 {
     return (p1.x-p0.x)*(p2.y-p0.y)>=(p2.x-p0.x)*(p1.y-p0.y);
 }
 bool cmp(const point &a,const point &b)
 {
     point temp=p[];
     double xmt=(a.x-temp.x)*(b.y-temp.y)-(b.x-temp.x)*(a.y-temp.y);
     if(xmt)                             //向量不共线就按逆时针旋转
         return xmt>;
     return Dist(a,temp)>Dist(b,temp);//向量共线取最长的。
 }
 void graham()//p[0]是左下角的元素
 {
     res[]=p[];
     sort(p+,p+n,cmp);//按照极角排序
     res[]=p[];
     res[]=p[];
     top=;
     for(int i=; i<n; i++)
     {
         while(multi(p[i],res[top],res[top-]))
             top--;
         res[++top]=p[i];
     }
 }
 bool check(point temp)
 {
     int i,m=multi(temp,res[],res[]);
     for(i=;i<top;i++)
     {
         if(multi(temp,res[i],res[i+])!=m)return ;
     }
     if(multi(temp,res[i],res[])!=m)return ;
     return ;
 }
 void work(point temp)
 {
     int bi,en;
     bi=en=;
     for(int i=;i<=top;i++)
     {
         if(multi(temp,res[i],res[bi]))bi=i;
         if(multi(temp,res[en],res[i]))en=i;
     }
     double ans=atan2(res[en].y-temp.y,res[en].x-temp.x)-atan2(res[bi].y-temp.y,res[bi].x-temp.x);
     if(ans<0.0)ans+=*PI;
     ans=ans*/PI;
     printf("%.2lf\n",ans);
 }
 int main()
 {
     int i,mina,m;
     point temp;
     while(~scanf("%d",&n))
     {
         scanf("%lf%lf",&p[].x,&p[].y);
         mina=;
         for(i=; i<n; i++)
         {
             scanf("%lf%lf",&p[i].x,&p[i].y);
             if(p[i].y<p[mina].y||(p[i].y==p[mina].y&&p[i].x<p[mina].x))
                 mina=i;
         }
         swap(p[].x,p[mina].x),swap(p[].y,p[mina].y);
         graham();
         scanf("%d",&m);
         while(m--)
         {
             scanf("%lf%lf",&temp.x,&temp.y);
             if(check(temp))
             {
                 printf("Bad Position\n");
             }
             else
             {
                 work(temp);
             }
         }
     }
 }