GODI2018基站

背景：

考场上30分，后来听说是FFT，觉得好高级，感觉这种问题好神仙。
改用$Markdown$。

题意：

在平面直角坐标系给$n$个点，每个点有一个对应的类型($012$)，求两个相同类型夹中点为不同类型的点的方案数。

思路：

由数学知识可以知道，对于三个点$(x_1,y_1)$，$(x_2,y_2)$，$(x_3,y_3)$，当$(x_2,y_2)$为另外两点中点时，必然有$x_2=\frac{x_1+x_3}{2},y_2=\frac{y_1+y_3}{2}$，同理，知道后者也可推出前者。

考虑将二维问题转化为一维问题。

对于每一个点的坐标$(x_i,y_i)$变为$1000\ast x_i+y_i$。那么对于三个点$1000\ast x_1+y_1$，$1000\ast x_2+y_2$，$1000\ast x_3+y_3$，若$1000\ast x_2+y_2$为其它两点的中点，必然有$1000\ast x_2+y_2=\frac{1000\ast x_1+y_1+1000\ast x_3+y_3}{2}$，即$(1000\ast x_2+y_2)\ast 2=1000\ast x_1+y_1+1000\ast x_3+y_3$。也就是说当前仅当满足$y\ast 2=x+z$时，$y$为$x,z$的中点。同理，知道后者也可推出前者。

问题就变得简单了。
考虑设$y$的类型为$0$，其它类型同理。
构造$A_i=[typ{e}_i=0]$，然后另$S_i=A_i\ast A_i$。我们有上面可以知道，若$y$点为$x,z$中点，其结果必然是偶数，而$x,z$，$z,x$为一种情况，因此结果需要除以$2$。所以由$S_i$来计算结果即可。

时间复杂度$\Theta (nlogn)$。

代码：

// luogu-judger-enable-o2
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const double pai=acos(-1.0);
struct comp
{
	double x,y;
	comp(double xx=0,double yy=0):x(xx),y(yy) {}
	friend comp operator+(const comp &x,const comp &y) {return comp(x.x+y.x,x.y+y.y);}
	friend comp operator-(const comp &x,const comp &y) {return comp(x.x-y.x,x.y-y.y);}
	friend comp operator*(const comp &a,const comp &b) {return comp(a.x*b.x-a.y*b.y,a.x*b.y+b.x*a.y);}
}a[2500000];
	int limit=1,n,l=0,ma=0;
	int r[2500000],ans[2500000],id[2500000];
	int type[2500000];
void init(int ma)
{
	while(limit<=(ma<<1))
		limit<<=1,l++;
	for(int i=1;i<=limit;i++)
		r[i]=((r[i>>1]>>1)|((i&1)<<(l-1)));
}
void FFT(comp *now,int ty)
{
	for(int i=0;i<limit;i++)
		if(i<r[i]) swap(now[i],now[r[i]]);
	for(int mid=1;mid<limit;mid<<=1)
	{
		comp wn(cos(pai/mid),ty*sin(pai/mid));
		for(int j=0,R=(mid<<1);j<limit;j+=R)
		{
			comp w(1,0);
			for(int k=0;k<mid;k++,w=w*wn)
			{
				comp x=now[j+k],y=w*now[j+k+mid];
				now[j+k]=x+y;
				now[j+k+mid]=x-y;
			}
		}
	}
}
void work(int x)
{
	for(int i=0;i<=limit;i++)
		a[i]=(comp){0,0};
	for(int i=0;i<=ma;i++)
		if(type[i]==x) a[i].x=1; else a[i].x=0;
	FFT(a,1);
	for(int i=0;i<=limit;i++)
		a[i]=a[i]*a[i];
	FFT(a,-1);
	for(int i=0;i<=(ma<<1);i+=2)
		if(type[i/2]&&type[i/2]!=x) ans[i/2]+=(int)floor(a[i].x/limit+0.001)/2;
}
int main()
{
	int x,y,z;
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	{
		scanf("%d %d %d",&x,&y,&z);
		type[x*1000+y]=z+1;
		id[i]=x*1000+y;
		ma=max(ma,x*1000+y);
	}
	init(ma);
	work(1),work(2),work(3);
	for(int i=1;i<=n;i++)
		printf("%d\n",ans[id[i]]);
}