HDU 5033

单调栈（序列）分析待补，正好区域赛前可以重温一下。

/*
ID:esxgx1
LANG:C++
PROG:B
*/
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cmath>
using namespace std;

const int maxn = 100007;


const double PI = acos(-1.0);

typedef pair<double, double> po;
#define px        first
#define py        second

const double EPS = 1e-9;
inline int sgn(double x){return x < -EPS ? -1 : x > EPS;}
template<typename T> T det(const T &x0, const T &y0, const T &x1, const T &y1) {return x0*y1 - x1*y0;}
double det(const po &p0, const po &p1) {return det(p0.px, p0.py, p1.px, p1.py);}
po operator -(const po &a, const po &b) { return make_pair(a.px - b.px, a.py - b.py);}
template<typename T> int clkws(const T &p, const T &p1, const T &p2){ return sgn(det(p1 - p, p2 - p));}

struct obss {
    po p;
    int id;
} obs[maxn];

po bl[maxn];
double ll[maxn], rr[maxn];

po st[maxn];
int sttf;

int cmp(const po &a, const po &b)
{
  return a.px < b.px;
}

int cmp2(const obss &a, const obss &b)
{
  return a.p.px < b.p.px;
}

int main(void)
{
    #ifndef ONLINE_JUDGE
    freopen("in.txt", "r", stdin);
    #endif

    int tt;
    scanf("%d", &tt);
    for(int t=1; t<=tt; ++t) {
        int nn;
        scanf("%d", &nn);
        for(int i=0; i<nn; ++i)
            scanf("%lf%lf", &bl[i].px, &bl[i].py);
        sort(bl, bl+nn, cmp);

        int qq;
        scanf("%d", &qq);
        for(int i=0; i<qq; ++i) {
            scanf("%lf", &obs[i].p.px);
            obs[i].id = i;
        }
        sort(obs, obs+qq, cmp2);

        int att;

        sttf = 0, att = 0;
        for(int i=0; i<qq; ++i) {
            for(; bl[att].px < obs[i].p.px; ++att) {
                //while(sttf>0 && st[sttf].py <= bl[att].py) --sttf;
                while(sttf>1 && clkws(st[sttf-1], st[sttf], bl[att]) >= 0) --sttf;
                st[++sttf] = bl[att];

            }
            while(sttf>1 && clkws(st[sttf-1], st[sttf], obs[i].p) >= 0) --sttf;
            ll[obs[i].id] = PI - atan2(st[sttf].py - obs[i].p.py, st[sttf].px - obs[i].p.px);
            //printf("ll[%d] : %f st(%f %f) obs(%f %f)
", i, ll[i], st[sttf].px, st[sttf].py, obs[i].px, obs[i].py);
        }

        sttf = 0, att = nn-1;

        for(int i=qq; i--; ) {
            for(; bl[att].px > obs[i].p.px; --att) {
                //while(sttf>0 && st[sttf].py <= bl[att].py) --sttf;
                while(sttf>1 && clkws(st[sttf-1], st[sttf], bl[att]) <= 0) --sttf;
                st[++sttf]= bl[att];
                //printf("bl[%d] : %f %f st[%d](%f %f) obs(%f %f)
", att, bl[att].px, bl[att].py, rr[i], sttf, st[sttf].px, st[sttf].py, obs[i].px, obs[i].py);
            }
            while(sttf>1 && clkws(st[sttf-1], st[sttf], obs[i].p) <= 0) --sttf;
            rr[obs[i].id] = atan2(st[sttf].py - obs[i].p.py, st[sttf].px - obs[i].p.px);
            //printf("rr[%d] : %f st[%d](%f %f) obs(%f %f)
", i, rr[i], sttf, st[sttf].px, st[sttf].py, obs[i].px, obs[i].py);
        }
        printf("Case #%d:
", t);
        for(int i=0; i<qq; ++i)
            printf("%.9f
", (PI - ll[i] - rr[i])/PI * 180);
    }
    return 0;
}

2014-09-28 01:27:10

Accepted

5033

718MS

7364K

3167 B

G++