礼物 HYSBZ - 4827 (fft + 构造 )
作者：互联网
我的室友最近喜欢上了一个可爱的小女生。马上就要到她的生日了，他决定买一对情侣手 环，一个留给自己，一
个送给她。每个手环上各有 n 个装饰物，并且每个装饰物都有一定的亮度。但是在她生日的前一天，我的室友突
然发现他好像拿错了一个手环，而且已经没时间去更换它了！他只能使用一种特殊的方法，将其中一个手环中所有
装饰物的亮度增加一个相同的自然数 c（即非负整数）。并且由于这个手环是一个圆，可以以任意的角度旋转它，
但是由于上面 装饰物的方向是固定的，所以手环不能翻转。需要在经过亮度改造和旋转之后，使得两个手环的差
异值最小。在将两个手环旋转且装饰物对齐了之后，从对齐的某个位置开始逆时针方向对装饰物编号 1,2,…,n，
其中 n 为每个手环的装饰物个数，第 1 个手环的 i 号位置装饰物亮度为 xi，第 2 个手 环的 i 号位置装饰物
亮度为 yi，两个手环之间的差异值为(参见输入输出样例和样例解释)： \sum_{i=1}^{n}(x_i-y_i)^2麻烦你帮他
计算一下，进行调整（亮度改造和旋转），使得两个手环之间的差异值最小， 这个最小值是多少呢？
Input
输入数据的第一行有两个数n, m，代表每条手环的装饰物的数量为n，每个装饰物的初始 亮度小于等于m。
接下来两行，每行各有n个数，分别代表第一条手环和第二条手环上从某个位置开始逆时 针方向上各装饰物的亮度。
1≤n≤50000, 1≤m≤100, 1≤ai≤m
Output
输出一个数，表示两个手环能产生的最小差异值。
注意在将手环改造之后，装饰物的亮度 可以大于 m。
Sample Input5 6 1 2 3 4 5 6 3 3 4 5
Sample Output1 【样例解释】 需要将第一个手环的亮度增加1，第一个手环的亮度变为： 2 3 4 5 6 旋转一下第二个手环。对于该样例，是将第 二个手环的亮度6 3 3 4 5向左循环移动 2017-04-15 第 6 页，共 6 页 一个位置，使得第二手环的最终的亮度为 ：3 3 4 5 6。 此时两个手环的亮度差异值为1。
要点
把一个数组扩大一倍，相当于右移
这个题求的是sigma A[ i ] * B [ i + k]
把a翻转，就成了 sigma A[ n-i ] * B [ i + k]
这样之后 和就成了 n + k 也就是 n+k项的系数
存在n+k项的数组里面
1 #include"bits/stdc++.h"
2 #define sd(x) scanf("%d",&(x));
3 #define sf(x) scanf("%lf",&(x));
4 #define sld(x) scanf("%lld",&(x));
5 using namespace std;
6
7 const int maxn = 2e6+10;
8 const double Pi = acos(-1.0);
9 struct cp
10 {
11 double x,y;
12 cp (double xx=0,double yy=0)
13 {x=xx,y=yy;}
14 }a[maxn],g[maxn],b[maxn];
15
16 cp operator + (cp a,cp b){ return cp(a.x+b.x , a.y+b.y);}
17 cp operator - (cp a,cp b){ return cp(a.x-b.x , a.y-b.y);}
18 cp operator * (cp a,cp b){ return cp(a.x*b.x-a.y*b.y , a.x*b.y+a.y*b.x);}//不懂的看复数的运算那部分
19
20 int n,m;
21 int l,r[maxn];
22 int limit = 1;
23 int s1,s2,f1,f2;
24
25
26 inline void fft(cp *a,int ff)
27 {
28 for(int i=0;i<limit;i++)
29 if(i<r[i])swap(a[i],a[r[i]]);
30 for(int mid=1;mid<limit;mid<<=1)
31 {
32 cp wn(cos(Pi/mid) , ff*sin(Pi/mid));
33 for(int R=mid<<1,j=0;j<limit;j+=R)
34 {
35 cp w(1,0);
36 for(int k=0;k<mid;k++,w=w*wn)
37 {
38 cp x=a[j+k],y=w*a[j+mid+k];
39 a[j+k]=x+y;
40 a[j+mid+k]=x-y;
41 }
42 }
43
44 }
45
46 }
47
48 int main()
49 {
50
51 sd(n); sd(m);
52 for(int i=0;i<n;i++)sf(a[n-1-i].x);
53 for(int i=0;i<n;i++)sf(b[i].x);
54 for(int i=0;i<n;i++)b[n+i]=b[i];
55
56 for(int i=0;i<n;i++)
57 f1+=(int)a[i].x*(int)a[i].x,
58 f2+=(int)b[i].x*(int)b[i].x,
59 s1+=(int)a[i].x,
60 s2+=(int)b[i].x;
61
62 while(limit<=3*n) limit<<=1,l++;//最高次项为 3*n
63 // cout<<"LIMIT: "<<limit<<endl;
64 for(int i=0;i<limit;i++)
65 r[i]=(r[i>>1]>>1)|( (i&1)<<(l-1));
66
67
68 fft(a,1); fft(b,1);
69 for(int i=0;i<=limit;i++)a[i]=a[i]*b[i];
70 fft(a,-1);
71 int ans = 2e9;
72
73
74
75 for(int i=0;i<limit;i++)a[i].x=(int)(a[i].x/limit+0.5);
76 // cout<<f1<<" "<<f2<<" "<<s1<<" "<<s2<<endl;
77 for(int i=0;i<n;i++)
78 {
79 for(int c=0;c<m;c++)
80 {
81 ans=min(ans,f1+f2-2*(int)a[n+i].x+2*c*(s1-s2)+n*c*c);
82 }
83 }
84 cout<<ans;
85
86
87 return 0;
88
89
90
91
92 }
标签：int,HYSBZ,fft,手环,4827,亮度,maxn,装饰物,cp 来源： https://www.cnblogs.com/zhangbuang/p/11018878.html