bzoj-3527 力

题意：

给出n个数qi，给出Fj的定义如下：

令Ei=Fi/qi。试求Ei。

题解：

这是一道难得的久违的自己推出来的数学题；

首先显然Ej和qj没啥关系。。。

然后观察分母实际上就是个i与j距离的平方；

展开一点的话可以发现

这样的形式(这个公式好模糊。。但是没什么关系嘛)；

那么显然这样稍微构造一下，可以从卷积得来；

然后画一画搞一搞就上FFT了；

FFT模板总敲错怎么办= =。。。

不要用long double ，因为cout会RE。。

代码：

#include<math.h>
#include<stdio.h>
#include<iomanip>
#include<string.h>
#include<iostream>
#include<algorithm>
#define N 262144
using namespace std;
typedef double ld;
const ld pi=acos(-1.0);
struct cp
{
	ld x,y;
	cp(ld a,ld b):x(a),y(b){}
	cp(){};
	cp operator +(const cp &a)
	{
		return cp(x+a.x,y+a.y);
	}
	cp operator -(const cp &a)
	{
		return cp(x-a.x,y-a.y);
	}
	cp operator *(const cp &a)
	{
		return cp(x*a.x-y*a.y,x*a.y+y*a.x);
	}
}a[N],b[N],c[N];
ld ans[N];
void FFT(cp *a,int len,int type)
{
	int i,j,h,t;
	for(i=0,t=0;i<len;i++)
	{
		if(i>t)	swap(a[i],a[t]);
		for(j=(len>>1);(t^=j)<j;j>>=1);
	}
	for(h=2;h<=len;h<<=1)
	{
		cp wn(cos(2*pi*type/h),sin(2*pi*type/h));
		for(i=0;i<len;i+=h)
		{
			cp temp,w(1,0);
			for(j=0;j<(h>>1);j++,w=w*wn)
			{
				temp=w*a[i+j+(h>>1)];
				a[i+j+(h>>1)]=a[i+j]-temp;
				a[i+j]=a[i+j]+temp;
			}
		}
	}
}
void slove(cp *a,cp *b,int len)
{
	FFT(a,len,1),FFT(b,len,1);
	for(int i=0;i<len;i++)	c[i]=a[i]*b[i];
	FFT(c,len,-1);
	for(int i=0;i<len;i++)	ans[i]=c[i].x/len;
}
int main()
{
	int n,m,i,j,k;
	scanf("%d",&n);
	for(i=0;i<n;i++)
	{
		scanf("%lf",&a[i].x);
		b[i].x=-1.0/(n-i)/(n-i);
	}
	for(i=1;i<=n;i++)
	{
		b[i+n].x=1.0/i/i;
	}
	for(i=(1<<30);i;i>>=1)
		if(n&i)
			{m=(i<<1);break;}
	m<<=1;
	slove(a,b,m);
	cout<<fixed<<setprecision(3);
	for(i=0;i<n;i++)
	{
		printf("%.3lf\n",ans[i+n]);
	}
	return 0;
}