其他分享
首页 > 其他分享> > 礼物 HYSBZ - 4827 (fft + 构造 )

礼物 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