(2019 CCPC 秦皇岛) A - Angle Beats (计算几何)

解：没做过几何题，但还是看了这个题，思想还是很简单的，分两类，一个就是以询问点为直角点，再选出一个点求出一个斜率，然后求在另一条边上的点；第二类就是询问点不是直角点，那就其他的点作为直角点，来更新这q次询问的点。

然后就是怎么写不T呢，看到一种map维护的写法很妙啊，键值为向量，就是斜率，然后值就是个数，但是我们可以将所有的向量都变到一二象限，因为对于一条边，顺时针的直角和逆时针的直角都是一样的，我们在这将其归为一类，所以可以通过（-x,-y）转化即可。

#include<bits/stdc++.h>
#define il inline
#define pb push_back
#define ms(_data,v) memset(_data,v,sizeof(_data))
#define SZ(a) int((a).size())
using namespace std;
typedef long long ll;
const ll inf=0x3f3f3f3f;
const int N=2e3+5; 
//il int Add(int &x,ll y) {return x=x+y>=mod?x+y-mod:x+y;}
//il int Mul(int &x,ll y) {return x=x*y>=mod?x*y%mod:x*y;}
struct node{
	ll x,y;
	il node change()const{
		if(x<0 || (x==0 && y<0)) return {-x,-y};
		return {x,y};
	}
	node operator -(const node &tp) const{ //保证在一二象限，角度减小180，但求得直角，所以不影响 
		return {x-tp.x,y-tp.y};
	}
	bool operator <(const node &tp)const{//角度从小到大 
		node na=change(),nb=tp.change();
		return na.x*nb.y < na.y*nb.x;
	}
}a[N],qu[N];
map<node,int> mp; 
int n,q,ans[N];
int main(){
	std::ios::sync_with_stdio(0);cin.tie(0);
	cin>>n>>q;
	for(int i=1;i<=n;++i){
		cin>>a[i].x>>a[i].y;
	}
	for(int i=1;i<=q;++i){
		cin>>qu[i].x>>qu[i].y;
	}
	node k;
	for(int i=1;i<=q;++i){//询问点作为直角点 
		mp.clear(); 
		for(int j=1;j<=n;++j){
			mp[a[j]-qu[i]]++;
		}
		for(int j=1;j<=n;++j){
			k=a[j]-qu[i];
			k={-k.y,k.x};
			ans[i]+=mp[k];
		}
		ans[i]/=2;//因为一个直角被算了两遍 
	}
	for(int i=1;i<=n;++i){//询问点不作为直角点 
		mp.clear(); 
		for(int j=1;j<=n;++j){
			if(i!=j) mp[a[j]-a[i]]++;
		}
		for(int j=1;j<=q;++j){
			k=qu[j]-a[i];
			k={-k.y,k.x};
			ans[j]+=mp[k];
		}
	} 
	for(int i=1;i<=q;++i) cout<<ans[i]<<endl;
	return 0;
}