题面

思路

这题很像bzoj4827礼物

还是一样的思路,我们把$y$倍长,$y[i+k]=y[i]+n$

然后令$f(s,c)$表示从$y$的第$s$个开始匹配,位置偏移量为$c$的答案

可以得到$f(s,c)=\sum_{i=0}{n-1}(x_i-y_{i+s}+c)2=\sum_{i=0}{n-1}(x_i2+y_{i+s}+c^2+2x_ic-2y_{i+s}x-2x_iy_{i+s})$

我们可以把右边这个式子视为关于$c$的二次函数

用$FFT$可以快速得到不同的$s$下,$2x_iy_{i+s}$这一项的值,又因为其他系数都是确定的,所以我们可以用$O(k\log k)$确定不同的s对应的c的函数

然后可以算出来这个函数的最佳取值,再从所有最佳取值中选出最优解即可

Code

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define ll long long
using namespace std;
inline int read(){
int re=0,flag=1;char ch=getchar();
while(!isdigit(ch)){
if(ch=='-') flag=-1;
ch=getchar();
}
while(isdigit(ch)) re=(re<<1)+(re<<3)+ch-'0',ch=getchar();
return re*flag;
}
namespace FFT{
struct complex{
double x,y;
complex(double xx=0,double yy=0){x=xx;y=yy;}
inline complex operator +(complex &b){return complex(x+b.x,y+b.y);}
inline complex operator -(complex &b){return complex(x-b.x,y-b.y);}
inline complex operator *(complex &b){return complex(x*b.x-y*b.y,x*b.y+y*b.x);}
}A[100010],B[100010];
int cnt,lim,r[100010];
const double pi=acos(-1.0);
void fft(complex *a,double type){
int i,j,k,mid;complex x,y,w,wn;
for(i=0;i<lim;i++) if(i<r[i]) swap(a[i],a[r[i]]);
for(mid=1;mid<lim;mid<<=1){
wn=complex(cos(pi/mid),type*sin(pi/mid));
for(j=0;j<lim;j+=(mid<<1)){
w=complex(1,0);
for(k=0;k<mid;k++,w=w*wn){
x=a[j+k];y=a[j+k+mid]*w;
a[j+k]=x+y;
a[j+k+mid]=x-y;
}
}
}
if(type==-1) for(i=0;i<lim;i++) a[i].x=(ll)(a[i].x/lim+0.5);
}
void init(int len){
memset(A,0,sizeof(A));memset(B,0,sizeof(B));
cnt=0;lim=1;
while(lim<=len) lim<<=1,cnt++;
for(int i=0;i<lim;i++) r[i]=((r[i>>1]>>1)|((i&1)<<(cnt-1)));
}
}
void mul(){
using namespace FFT;
fft(A,1);
fft(B,1);
for(int i=0;i<lim;i++) A[i]=A[i]*B[i];
fft(A,-1);
}
double sqr(double x){
return x*x;
}
ll n,k,x[20010],y[20010],pre1[20010],pre2[20010];
int main(){
int T=read();
while(T--){
n=read();k=read();
ll t1=0,t2=0,i,j,ans=1e15,tl,tm; for(i=0;i<k;i++){
x[i]=read();
t1+=(x[i]*x[i]);
t2+=2*x[i];
}
for(i=0;i<k;i++){
y[i]=read();
y[i+k]=y[i]+n;
}
pre1[0]=y[0]*y[0];
pre2[0]=y[0];
for(i=1;i<(k<<1);i++){
pre1[i]=pre1[i-1]+y[i]*y[i];
pre2[i]=pre2[i-1]+y[i];
} FFT::init(k*3);
for(i=0;i<k;i++) FFT::A[i].x=x[k-1-i];
for(i=0;i<(k<<1);i++) FFT::B[i].x=y[i];
mul(); for(i=k-1;i<(k<<1)-1;i++){
j=i-k+1;
tl=t1+pre1[j+k-1]-pre1[j-1];
tm=t2-2*(pre2[j+k-1]-pre2[j-1]);//这里处理的是y的前缀和,以及y的平方的前缀和
ll tmpc=(ll)(-(double)tm/(2.0*(double)k)); ans=min(ans,(ll)(k*sqr(tmpc)+tmpc*tm+tl-2*FFT::A[i].x));//这里需要上下都摸♂索一下
tmpc--;
ans=min(ans,(ll)(k*sqr(tmpc)+tmpc*tm+tl-2*FFT::A[i].x));
tmpc+=2;
ans=min(ans,(ll)(k*sqr(tmpc)+tmpc*tm+tl-2*FFT::A[i].x));
}
printf("%lld\n",ans);
}
}

项链 [FFT]的更多相关文章

  1. [2016北京集训试题15]项链-[FFT]

    Description Solution 设y[i+k]=y[i]+n. 由于我们要最优解,则假如将x[i]和y[σ[i]]连线的话,线是一定不会交叉的. 所以,$ans=\sum (x_{i}-y_ ...

  2. HDU5730 Shell Necklace(DP + CDQ分治 + FFT)

    题目 Source http://acm.hdu.edu.cn/showproblem.php?pid=5730 Description Perhaps the sea‘s definition of ...

  3. HDU 5730 Shell Necklace cdq分治+FFT

    题意:一段长为 i 的项链有 a[i] 种装饰方式,问长度为n的相连共有多少种装饰方式 分析:采用dp做法,dp[i]=∑dp[j]*a[i-j]+a[i],(1<=j<=i-1) 然后对 ...

  4. #8 //HDU 5730 Shell Necklace(CDQ分治+FFT)

    Description 给出长度分别为1~n的珠子,长度为i的珠子有a[i]种,每种珠子有无限个,问用这些珠子串成长度为n的链有多少种方案 题解: dp[i]表示组合成包含i个贝壳的项链的总方案数 转 ...

  5. HDU - 5730 :Shell Necklace(CDQ分治+FFT)

    Perhaps the sea‘s definition of a shell is the pearl. However, in my view, a shell necklace with n b ...

  6. 并行计算提升32K*32K点(32位浮点数) FFT计算速度(4核八线程E3处理器)

    对32K*32K的随机数矩阵进行FFT变换,数的格式是32位浮点数.将产生的数据存放在堆上,对每一行数据进行N=32K的FFT,记录32K次fft的时间. 比较串行for循环和并行for循环的运行时间 ...

  7. BZOJ1878[SDOI2009]HH的项链

    Description HH有一串由各种漂亮的贝壳组成的项链.HH相信不同的贝壳会带来好运,所以每次散步 完后,他都会随意取出一段贝壳,思考它们所表达的含义.HH不断地收集新的贝壳,因此, 他的项链变 ...

  8. 【BZOJ-2179&2194】FFT快速傅里叶&快速傅里叶之二 FFT

    2179: FFT快速傅立叶 Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 2978  Solved: 1523[Submit][Status][Di ...

  9. 为什么FFT时域补0后,经FFT变换就是频域进行内插?

    应该这样来理解这个问题: 补0后的DFT(FFT是DFT的快速算法),实际上公式并没变,变化的只是频域项(如:补0前FFT计算得到的是m*2*pi/M处的频域值, 而补0后得到的是n*2*pi/N处的 ...

随机推荐

  1. 使用nsis开发自定义安装包使用心得,以及遇到坑

    因为新公司需要开发pc应用的自定义安装包,开始时候计划使用nsis开发,论坛上面有很多不错的例子,而且完成度很强, 随便拿来修改使用,但是后续的开发过程中遇到的问题就逐个出现. 首先说一下nsis的优 ...

  2. Python 文件访问模式

    f = open('xxx文件', '访问模式') r    以只读方式打开文件(read).文件的指针将会放在文件的开头.默认模式. w   打开一个文件只用于写入(write).如果该文件已存在则 ...

  3. 【读书笔记】你不知道的JavaScript(上卷)--作用域是什么

    第一章 作用域 1.理解作用域 几个名词的介绍 引擎:从头到尾负责整个JavaScript程序的编译及执行过程 编译器:负责语法分析及代码生成器等脏活累活 作用域:负责收集并维护由所有声明的标识符(变 ...

  4. 日期格式兼容iOS

    iOS不支持2016-02-11 12:21:12格式的日期 目前Safari可以支持的标准格式: MM-dd-yyyy yyyy/MM/dd MM/dd/yyyy MMMM dd, yyyy MMM ...

  5. HTML语义化的重要性

    语义化标签就是尽量使用有相对应的结构的含义的Html的标签 1.结构更好,更利于搜索引擎的抓取(SEO的优化)和开发人员的维护(可维护性更高,因为结构清晰,so易于阅读). 2.更有利于特殊终端的阅读 ...

  6. 【jQeury】input输入框状态,input事件,blur事件,focus事件

    //输入框正在输入时 $("#test1").on('input',function(){ alert('正在输入'); }) //输入框得到焦点时 $("#test2& ...

  7. 包围轮廓的矩形边界 opencv

    #include<opencv2/opencv.hpp> #include<iostream> using namespace std; using namespace cv; ...

  8. Hacker Cups and Balls Gym - 101234A 二分+线段树

    题目:题目链接 题意:有编号从1到n的n个球和n个杯子. 每一个杯子里有一个球, 进行m次排序操作,每次操作给出l,r. 如果l<r,将[l,r]范围内的球按升序排序, 否则降序排, 问中间位置 ...

  9. [Link-Cut-Tree][BZOJ2002]弹飞绵羊

    题面 Description 某天,Lostmonkey发明了一种超级弹力装置,为了在他的绵羊朋友面前显摆,他邀请小绵羊一起玩个游戏.游戏一开始,Lostmonkey在地上沿着一条直线摆上\(n\)个 ...

  10. 零基础学css

    选择器:标签选择器.id选择器.类选择器 ---------------------------------------------------------------------------- 标签 ...