BZOJ1035: [ZJOI2008]Risk（最小左转法+点定位）

传送门

题意：
给一个平面图和图上的一些点，问这些点与哪些点在相邻的区域。

题解：

首先用最小左转法确定每一条线左边的区域。
接下来用扫描线做最小点定位。
最后统计在相邻区域的点对即可。

有一个比较难处理的地方是区域相互包含，一开始不知道怎么做就去膜了Claris题解，发现对于每个不同的联通块，确定最左上方的点，做一遍扫描线，之后内层向外层连一条边就好了。

还有凸包的比较函数比较难调（我调了两天。。），注意边界。

#include<bits/stdc++.h>
using namespace std;
inline int read(){
    char ch=getchar();int i=0,f=1;
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch)){i=(i<<1)+(i<<3)+ch-'0';ch=getchar();}
    return i*f;
}
inline void W(int x){
    static int buf[50];
    if(!x){putchar('0');return;}
    if(x<0){putchar('-');x=-x;}
    while(x)buf[++buf[0]]=x%10,x/=10;
    while(buf[0])putchar(buf[buf[0]--]+'0');
}

const int Maxn=1e4+50;
const double PI=acos(-1.0);
const double eps=1e-9;
inline int sgn(double x){return (int)(x>eps)-(x<-eps);}
struct P{
    int x,y;
    P(){}
    P(int x,int y):x(x),y(y){}
    friend inline P operator -(const P &a,const P &b){return P(a.x-b.x,a.y-b.y);}
    friend inline int operator *(const P &a,const P &b){return a.x*b.y-a.y*b.x;}
    inline double slope(){return atan2(y,x);}
}p[Maxn];
struct E{
    int a,b;
    E(){}
    E(int a,int b):a(a),b(b){}
}e[Maxn];

int n,m,a[Maxn],tot,g[Maxn],v[Maxn],nt[Maxn],ec=1,bl[Maxn],bl_p[Maxn],Pcnt;
inline void add(int x,int y){nt[++ec]=g[x];g[x]=ec;v[ec]=y;}
map< int,int >S_p[Maxn<<1];
inline int getid(int x,int y){
    x+=10000;
    if(S_p[x].count(y))return S_p[x][y];
    else return (S_p[x].insert(make_pair(y,++tot)),p[tot]=P(x-10000,y),tot);
}
namespace Area{

int vis[Maxn];
struct cmp{
    inline bool operator ()(const pair<double,int> &a,const pair<double,int> &b){
        int t=sgn(a.first-b.first);
        return (t<0)||(t==0&&a.second<b.second);
    }
};
set< pair<double,int>,cmp >S_l[Maxn];
int st[Maxn],top;
inline int area(){
    int res=0;
    for(int i=1;i<=top;i++)res+=p[v[st[i]^1]]*p[v[st[i]]];
    return res;
}
inline int getsuf(int now){
    double ang=(p[v[now]]-p[v[now^1]]).slope();
    int t=v[now];
    set< pair<double,int> >::iterator suf=S_l[t].lower_bound(make_pair((ang<=0)?ang+PI:ang-PI,0));
    if(suf==S_l[t].begin())return (*--S_l[t].end()).second;
    return (*--suf).second;
}
inline void dfs(int x){
    vis[x]=1;st[++top]=x;
    int suf=getsuf(x);
    if(vis[suf]){
        int ar=area();
        int nowP=(ar>0?++Pcnt:0);
        for(int i=1;i<=top;i++)bl[st[i]]=nowP;          
        for(;top>=1;--top){
            int a=v[st[top]^1],b=v[st[top]];
            S_l[a].erase(S_l[a].lower_bound(make_pair((p[b]-p[a]).slope(),st[top])));
        }
    }else dfs(suf);
}
inline void getarea(){
    for(int i=1;i<=tot;i++){
        for(int j=g[i];j;j=nt[j])
        S_l[i].insert(make_pair((p[v[j]]-p[i]).slope(),j));
    }
    for(int i=2;i<=ec;i++)
        if(!vis[i])dfs(i);
}

}

namespace Scanline{

int vis[Maxn],tp,tcnt;
struct T{
    int a,b,u,id;
    T(){}
    T(int a,int b,int u):a(a),b(b),u(u){}
    friend inline bool operator <(const T &a,const T &b){
        if(a.a==b.a&&a.u!=b.u)return a.u>b.u;
        else if(a.a==b.a)return (p[a.b]-p[a.a])*(p[b.b]-p[a.a])<0;
        else if(p[a.a].x!=p[b.a].x)return p[a.a].x<p[b.a].x;
        else return p[a.a].y>p[b.a].y;
    }
}t[Maxn];
struct cmp{
    inline bool operator ()(const T &a,const T &b){
        if(a.a==a.b)return (p[b.b]-p[a.a])*(p[b.a]-p[a.a])<0;
        if(b.a==b.b)return (p[a.b]-p[b.a])*(p[a.a]-p[b.a])>0;
        int t1=(p[b.b]-p[a.a])*(p[b.a]-p[a.a]),t2=(p[b.b]-p[a.b])*(p[b.a]-p[a.b]);
        int t3=(p[a.b]-p[b.a])*(p[a.a]-p[b.a]),t4=(p[a.b]-p[b.b])*(p[a.a]-p[b.b]);
        return (t1<0&&t2<=0)||(t2<0&&t1<=0)||(t3>0&&t4>=0)||(t4>0&&t3>=0);
    }
};
set<T,cmp>S;
void dfs(int x){
    vis[x]=1;if(p[x].y>p[tp].y)tp=x;
    for(int j=g[x];j;j=nt[j])if(!vis[v[j]])dfs(v[j]);
}
inline void connect(){
    tcnt=0;S.clear();
    for(int i=1,j=tot;i<=j;i++)
        if(g[i]&&!vis[i]){
            dfs(tp=i);
            t[++tcnt]=T(tp,tp,1);
        }
    for(int i=1;i<=tot;i++)
        for(int j=g[i];j;j=nt[j])
            if(p[v[j]].x>p[i].x){
                t[++tcnt]=T(i,v[j],0);
                t[++tcnt]=T(v[j],i,2);
            }
    sort(t+1,t+tcnt+1);
    for(int i=1;i<=tcnt;i++){
        if(t[i].u==0)S.insert(t[i]);
        else if(t[i].u==2)S.erase(S.find(T(t[i].b,t[i].a,0)));
        else{
            set<T>::iterator it=S.lower_bound(t[i]);
            if(it==S.begin())continue;
            --it;add(it->a,t[i].a);add(t[i].a,it->a);
        }
    }
}
inline void getarea(){
    tcnt=0;S.clear();
    for(int i=1;i<=n;i++){
        t[++tcnt]=T(a[i],a[i],1);
    }
    for(int i=1;i<=tot;i++){
        for(int j=g[i];j;j=nt[j]){
            if(p[v[j]].x>p[i].x){
                t[++tcnt]=T(i,v[j],0);
                t[tcnt].id=j^1;
                t[++tcnt]=T(v[j],i,2);
            }
        }
    }
    sort(t+1,t+tcnt+1);
    for(int i=1;i<=tcnt;i++){
        if(t[i].u==0)S.insert(t[i]);
        else if(t[i].u==2)S.erase(S.find(T(t[i].b,t[i].a,0)));
        else{
            set<T>::iterator it=S.lower_bound(t[i]);
            --it;bl_p[t[i].a]=bl[(it)->id];
        }
    }
}
}
bool isadj[Maxn][Maxn];
int main(){
    n=read(),m=read();
    for(int i=1;i<=n;i++){
        int x=read(),y=read();
        a[i]=getid(x,y);
    }
    for(int i=1;i<=m;i++){
        int x=read(),y=read();
        e[i].a=getid(x,y);
        x=read(),y=read();
        e[i].b=getid(x,y);
        add(e[i].a,e[i].b);
        add(e[i].b,e[i].a);
    }
    Scanline::connect();
    Area::getarea();
    Scanline::getarea();
    for(int i=2;i<=ec;i++)
        if(bl[i]&&bl[i^1]){
            isadj[bl[i]][bl[i^1]]=1;
        } 
    for(int i=1;i<=n;i++){
        static int st[Maxn],top;
        top=0;
        for(int j=1;j<=n;j++)if(i!=j&&isadj[bl_p[a[i]]][bl_p[a[j]]])st[++top]=j;
        W(top);
        for(int j=1;j<=top;j++)putchar(' '),W(st[j]);
        if(i!=n)putchar('\n');
    }
}