【BZOJ2391】Cirno的忧郁（三角剖分）（平衡树维护极角序扫描线）

传送门

---

题解：

首先题目似乎没有说不存在三点共线的情况。。。

于是，我随便找了5篇题解，全部可以被以下两组数据中的一组或两组hack。

3 1
0 2
2 2
2 0
1 2 1
2
3
1 2 3
3
3 2 1

3 1
0 2
2 2
2 0
1 1 1
2
3
1 2 3
3
3 2 1

先讲一下思路：

考虑类似面积的处理，我们从原点向多边形顶点连线，将原多边形剖成若干三角形。计算面积直接这些三角形的有向面积加起来就行了。这里计算权值显然也是三角形内部加加减减就行了。

问题就是怎么处理出这些三角形内部点的权值和。

简单来说就是：极角序扫描线，同时平衡树按照极角序维护权值和。

注意一下三点共线的细节。

我的代码计算几何部分没有用任何浮点数。

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

namespace IO{
	inline char gc(){
		static cs int Rlen=1<<22|1;static char buf[Rlen],*p1,*p2;
		return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;
	}template<typename T>T get(){
		char c;bool f=false;while(!isdigit(c=gc()))f=c=='-';
		T num=c^48;while(isdigit(c=gc()))num=((num+(num<<2))<<1)+(c^48);
		return f?-num:num;
	}inline int gi(){return get<int>();}
}using namespace IO;

using std::cerr;
using std::cout;

struct Pnt{
	int x,y;Pnt(){}Pnt(int _x,int _y):x(_x),y(_y){}
	friend Pnt operator+(cs Pnt &a,cs Pnt &b){return Pnt(a.x+b.x,a.y+b.y);}
	friend Pnt operator-(cs Pnt &a,cs Pnt &b){return Pnt(a.x-b.x,a.y-b.y);}
	friend int operator*(cs Pnt &a,cs Pnt &b){return a.x*b.y-b.x*a.y;}
	friend bool operator==(cs Pnt &a,cs Pnt &b){return a.x==b.x&&a.y==b.y;}
	friend int dot(cs Pnt &a,cs Pnt &b){return a.x*b.x+a.y*b.y;}
	int len()cs{return x*x+y*y;}
};

inline bool cmp(cs Pnt &a,cs Pnt &b){
	return a*b?a*b>0:a.len()<b.len();
}

cs int N=1e3+7;

int n,m;

int id[N],rk[N],vl[N];

Pnt p[N],q[N];

int sum[N][N];

inline bool cmp_q(int i,int j){return cmp(q[i],q[j]);}
inline bool cmp_p(int i,int j){return cmp(p[i],p[j]);}
void init_q(){
	for(int re i=1;i<=m;++i)id[i]=i;
	std::sort(id+1,id+m+1,cmp_q);
	for(int re i=1;i<=m;++i)rk[id[i]]=i;
	for(int re i=1;i<=m;++i)
		while(rk[i]!=i){
			int j=rk[i];
			std::swap(rk[i],rk[j]);
			std::swap(q[i],q[j]);
			std::swap(vl[i],vl[j]);
		}
}
void init_p(){
	for(int re i=1;i<=n;++i)id[i]=i;
	std::sort(id+1,id+n+1,cmp_p);
	for(int re i=1;i<=n;++i)rk[id[i]]=i;
	for(int re i=1;i<=n;++i)
		while(rk[i]!=i){
			int j=rk[i];
			std::swap(rk[i],rk[j]);
			std::swap(p[i],p[j]);
		}
	for(int re i=1;i<=n;++i)rk[id[i]]=i;
}

namespace SBT{
	cs int N=2e3+7;
	
	Pnt nd[N];int tot,rt;
	int sum[N],vl[N],siz[N],lc[N],rc[N];
	
	inline void clear(){tot=rt=0;}
	inline void push_up(int u){
		siz[u]=siz[lc[u]]+siz[rc[u]]+1;
		sum[u]=sum[lc[u]]+sum[rc[u]]+vl[u];
	}
	inline void Zig(int &u){
		int v=lc[u];lc[u]=rc[v],rc[v]=u;
		push_up(u),push_up(v);u=v;
	}
	inline void Zag(int &u){
		int v=rc[u];rc[u]=lc[v],lc[v]=u;
		push_up(u);push_up(v);u=v;
	}
	
	inline void ins(int &u,cs Pnt &p,int v){
		if(!u){
			u=++tot;nd[u]=p,sum[u]=vl[u]=v;
			siz[u]=1;lc[u]=rc[u]=0;return;
		}
		if(nd[u]*p==0){sum[u]+=v,vl[u]+=v;return;}
		if(nd[u]*p>0){
			ins(lc[u],p,v);push_up(u);
			if(siz[lc[u]]>siz[u]/2)Zig(u);
		}else {
			ins(rc[u],p,v);push_up(u);
			if(siz[rc[u]]>siz[u]/2)Zag(u);
		}
	}
	inline int query(int u,cs Pnt &p,bool f){
		if(!u)return 0;
		if(nd[u]*p<0||(nd[u]*p==0&&f))
			return sum[lc[u]]+vl[u]+query(rc[u],p,f);
		else return query(lc[u],p,f);
	}
}

void init_qy(int st){
	SBT::clear();int h=1;
	while(h<=m&&cmp(q[h],p[st]))++h;
	for(int re i=st+1;i<=n;++i){
		while(h<=m&&(cmp(q[h],p[i])||q[h]==p[i]))
			SBT::ins(SBT::rt,q[h]-p[st],vl[h]),++h;
		sum[st][i]=SBT::query(SBT::rt,p[i]-p[st],1);
		sum[i][st]=-SBT::query(SBT::rt,p[i]-p[st],0);
	}
}

int ps[N];
void Query(){int l=gi();ll area=0,ans=0;
	for(int re i=1;i<=l;++i)ps[i]=rk[gi()];ps[l+1]=ps[1];
	for(int re i=1;i<=l;++i)area+=p[ps[i]]*p[ps[i+1]];
	if(area<0)std::reverse(ps+1,ps+l+1),ps[l+1]=ps[1];
	for(int re i=1;i<=l;++i)ans+=sum[ps[i]][ps[i+1]];
	cout<<abs(ans)<<"\n";
}

void Main(){
	n=gi(),m=gi();
	for(int re i=1;i<=n;++i)p[i].x=gi()+10000,p[i].y=gi()+10001;
	for(int re i=1;i<=m;++i)q[i].x=gi()+10000,q[i].y=gi()+10001,vl[i]=gi();
	init_q();init_p();for(int re i=1;i<=n;++i)init_qy(i);
	int Q=gi();while(Q--)Query();
}

void file(){
#ifdef zxyoi
	freopen("freeze.in","r",stdin);
#endif
}
signed main(){file();Main();return 0;}