hdu 5033 Building

最新推荐文章于 2018-07-23 21:38:01 发布

会写代码的咸鱼

最新推荐文章于 2018-07-23 21:38:01 发布

阅读量343

点赞数

分类专栏：学生时代刷题归档

本文链接：https://blog.csdn.net/kep159/article/details/39484695

版权

学生时代刷题归档专栏收录该内容

146 篇文章 0 订阅

订阅专栏

网络赛的时候没有想到，智商被压制了。。。。下来之后突然有了思路。。。怒写，过之

#include<iostream>
#include<sstream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#include<map>
#include<time.h>
#include<set>
#include<string>
#include<vector>
#include<algorithm>
using namespace std;
#define inf 0x7fffffff
#define lc l,m,index<<1
#define rc m+1,r,index<<1|1
#define max_n 300005
#define mod 10000007
#define FOR(i,s,t) for(int i=(s);i<=(t);++i)
#define LL __int64
struct P
{
	double x,y;
	bool flag;
	int index;
	P(){}
	P(double x,double y):x(x),y(y){
	}
	P operator +(P p){
		return P(x+p.x,y+p.y);
	}
	P operator -(P p){
		return P(x-p.x,y-p.y);
	}
	double det(P p){
		return x*p.y-y*p.x;
	}
};
P num[max_n];
P L[max_n],R[max_n];
P qs[max_n];
const double pi=acos(-1.0);
bool cmp(P a,P b)
{
    return a.x<b.x;
}
double ANS[max_n];
int main()
{
   // printf("%lf\n",((atan2(3,-2)-atan2(1,1))/pi*180.0));
    int step=1;
    int T,n,m;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        for(int i=0;i<n;i++)
        {
            scanf("%lf%lf",&num[i].x,&num[i].y);
            num[i].flag=0;
        }
        scanf("%d",&m);
        for(int i=n;i<n+m;i++)
        {
            scanf("%lf",&num[i].x);
            num[i].y=0;
            num[i].flag=1;
            num[i].index=i-n;
        }
        sort(num,num+n+m,cmp);
        int k=0;
        int ans=0;
        for(int i=0;i<n+m;i++)
        {
            while(k>1 && (qs[k-1]-num[i]).det(qs[k-2]-qs[k-1])<=0) k--;

            if(num[i].flag==1)
            L[ans++]=qs[k-1];
            else
            qs[k++]=num[i];
        }
       // printf("%lf %lf\n",L[0].x,L[0].y);
        k=0;
        ans--;
        for(int i=n+m-1;i>=0;i--)
        {
            while(k>1 && (qs[k-1]-num[i]).det(qs[k-2]-qs[k-1])>=0) k--;

            if(num[i].flag==1)
            R[ans--]=qs[k-1];
            else
            qs[k++]=num[i];
        }
        k=0;
        for(int i=0;i<n+m;i++)
        {
            if(num[i].flag==1)
            {
                ANS[num[i].index]=(atan2(L[k].y,L[k].x-num[i].x)-atan2(R[k].y,R[k].x-num[i].x))/pi*180.0;
                k++;
            }
        }
        printf("Case #%d:\n",step++);
      // for(int i=0;i<n+m;i++)printf("%lf %lf\n",num[i].x,num[i].y);
    //   for(int i=0;i<m;i++) printf("L = %.3lf %.3lf  R = %.3lf %.3lf\n",L[i].x,L[i].y,R[i].x,R[i].y);
        for(int i=0;i<m;i++)
        printf("%.10lf\n",ANS[i]);

    }
    return 0;
}