题面

传送门:洛咕


Solution

调得我头大,我好菜啊

好吧,我们来颓柿子吧:

我们可以只旋转其中一个手环。对于亮度的问题,因为可以在两个串上增加亮度,我们也可以看做是可以为负数的。

所以说,我们可以假设我们旋转\(B\)串,上下要加上的亮度差为\(p\),可以直接拍出一个最暴力的柿子:

设\(f(x)\)表示\(B\)串以\(x\)为开头的差异值,有:

\(f(x)=\sum_{i=0}^{x-1}(B[i]-A[i+n-x]+p)^2+\sum_{i=x}^{n-1}(B[i]-A[i-x]+p)^2\)

大力展开化简后有:

\(f(x)=\sum_{i=0}^{n-1}A[i]^2+\sum_{i=0}^{n-1}B[i]^2+n*p^2-2p\sum_{i=0}^{n-1}(A[i]-B[i])-2\sum_{i=0}^{x-1}(B[i]*A[i+n-x])-2\sum_{i=x}^{n-1}(B[i]*A[i-x])\)

前两项\(\sum_{i=0}^{n-1}A[i]^2+\sum_{i=0}^{n-1}B[i]^2\)显然\(O(n)预处理出来\)

中间两项\(n*p^2-2p\sum_{i=0}^{n-1}(A[i]-B[i])\)是一个关于\(p\)的二次函数,我们找最小值就好。(因为这题\(m\)非常小,我们也可以暴力枚举),复杂度\(O(1)\)或\(O(m)\)。

最后两项\(-2\sum_{i=0}^{x-1}(B[i]*A[i+n-x])-2\sum_{i=x}^{n-1}(B[i]*A[i-x])\)看起来非常像卷积,但是并不是,因此我们得做点处♂理。

蒟蒻本人是这样处理的:

首先,后面那个循环范围是肯定没法卷的,因此我们先把后面的循环处理一下得:

\(-2\sum_{i=0}^{x-1}(B[i]*A[i+n-x])-2\sum_{i=0}^{n-x-1}(A[i]*B[i+x])\)

然后我们可以考虑把前面那项的\(A\)反转(这样可以处理掉\(n\)来方便卷积),把后面那项的\(B\)反转(这样可以制造\(n\)与\(\sum\)对应)

\(-2\sum_{i=0}^{x-1}(B[i]*A'[x-1-i])-2\sum_{i=0}^{n-x-1}(A[i]*B'[n-1-i-x])\)

哦豁,卷积,搞定。

时间复杂度\(O(n*logn)\)


Code

我什么时候才能一次性写对FFT啊

//Luogu P3723 [AH2017/HNOI2017]礼物
//Jan,20th,2019
//颓柿子+FFT加速计算
#include<iostream>
#include<cstdio>
#include<cmath>
#include<complex>
#include<algorithm>
using namespace std;
long long read()
{
long long x=0,f=1; char c=getchar();
while(!isdigit(c)){if(c=='-') f=-1;c=getchar();}
while(isdigit(c)){x=x*10+c-'0';c=getchar();}
return x*f;
}
const int M=50000+100;
const int N=M*4;
const double PI=acos(-1);
typedef complex <double> cp;
inline cp omega(int K,int n)
{
return cp(cos(2*PI*K/n),sin(2*PI*K/n));
}
void FFT(cp a[],int n,bool type)
{
static int tmp[N],num=n-1,len=0;
while(num!=0) num/=2,len++;
for(int i=0,j;i<=n;i++)
{
for(j=0,num=i;j<len;j++)
tmp[j]=num%2,num/=2;
reverse(tmp,tmp+len);
for(j=0,num=0;j<len;j++)
num+=tmp[j]*(1<<j);
if(i<num) swap(a[i],a[num]);
}
for(int l=2;l<=n;l*=2)
{
cp x0=omega(1,l);
if(type==true) x0=conj(x0);
int m=l/2;
for(int j=0;j<n;j+=l)
{
cp x(1,0);
for(int k=0;k<m;k++,x*=x0)
{
cp temp=x*a[j+k+m];
a[j+k+m]=a[j+k]-temp;
a[j+k]=a[j+k]+temp;
}
}
}
}
int n,m,a[N],b[N];
cp B1[N],A1[N],B2[N],A2[N];
long long ans;
int main()
{
freopen("3723.in","r",stdin); n=read(),m=read();
for(int i=0;i<n;i++)
a[i]=read();
for(int i=0;i<n;i++)
b[i]=read(); long long dif=0;
for(int i=0;i<n;i++)
ans+=a[i]*a[i]+b[i]*b[i],dif+=a[i]-b[i];
long long t_ans=0x3f3f3f3f*0x3f3f3f3f;
for(int i=-m;i<=m;i++)
t_ans=min(t_ans,n*i*i-2*i*dif);
ans+=t_ans; int t=1;
while(t<2*n) t*=2;
reverse(a,a+n);
for(int i=0;i<n;i++)
A1[i]=a[i],B1[i]=b[i];
FFT(A1,t,false);
FFT(B1,t,false);
for(int i=0;i<t;i++)
A1[i]*=B1[i];
FFT(A1,t,true);
for(int i=0;i<t;i++)
A1[i].real()/=t; reverse(a,a+n);
reverse(b,b+n);
for(int i=0;i<n;i++)
A2[i]=a[i],B2[i]=b[i];
FFT(A2,t,false);
FFT(B2,t,false);
for(int i=0;i<t;i++)
A2[i]*=B2[i];
FFT(A2,t,true);
for(int i=0;i<t;i++)
A2[i].real()/=t; t_ans=(long long)(2*floor(A2[n-1].real()+0.5));
for(int i=1;i<n;i++)
t_ans=max(t_ans,(long long)(2*floor(A1[i-1].real()+A2[n-i-1].real()+0.5)));
ans-=t_ans;
printf("%lld",ans);
return 0;
}

[Luogu P3723] [AH2017/HNOI2017]礼物 (FFT 卷积)的更多相关文章

  1. LUOGU P3723 [AH2017/HNOI2017]礼物 (fft)

    传送门 解题思路 首先我们设变化量为\(r\),那么最终的答案就可以写成 : \[ ans=min(\sum\limits_{i=1}^n(a_i-b_i+r)^2) \] \[ ans=min(\s ...

  2. 洛谷 P3723 [AH2017/HNOI2017]礼物 解题报告

    P3723 [AH2017/HNOI2017]礼物 题目描述 我的室友最近喜欢上了一个可爱的小女生.马上就要到她的生日了,他决定买一对情侣手环,一个留给自己,一个送给她.每个手环上各有 \(n\) 个 ...

  3. P3723 [AH2017/HNOI2017]礼物

    题目链接:[AH2017/HNOI2017]礼物 题意: 两个环x, y 长度都为n k可取 0 ~ n - 1      c可取任意值 求 ∑ ( x[i] - y[(i + k) % n + 1] ...

  4. 洛谷P3723 [AH2017/HNOI2017]礼物(FFT)

    传送门 首先,两个数同时增加自然数值相当于只有其中一个数增加(此增加量可以小于0) 我们令$x$为当前的增加量,${a},{b}$分别为旋转后的两个数列,那么$$ans=\sum_{i=1}^n(a_ ...

  5. [AH2017/HNOI2017]礼物(FFT)

    题目描述 我的室友最近喜欢上了一个可爱的小女生.马上就要到她的生日了,他决定买一对情侣手 环,一个留给自己,一 个送给她.每个手环上各有 n 个装饰物,并且每个装饰物都有一定的亮度.但是在她生日的前一 ...

  6. Luogu 3723 [AH2017/HNOI2017]礼物

    BZOJ 4827 $$\sum_{i = 1}^{n}(x_i - y_i + c)^2 = \sum_{i = 1}^{n}(x_i^2 + y_i^2 + c^2 - 2 * x_iy_i + ...

  7. 洛谷P3723 [AH2017/HNOI2017]礼物

    吴迪说他化学会考上十分钟就想出来了,太神了%%%不过我也十分钟 但是调了一个多小时啊大草 懒得人话翻译了,自己康吧: 我的室友(真的是室友吗?)最近喜欢上了一个可爱的小女生.马上就要到她的生日了,他决 ...

  8. 笔记-[AH2017/HNOI2017]礼物

    笔记-[AH2017/HNOI2017]礼物 [AH2017/HNOI2017]礼物 \[\begin{split} ans_i=&\sum_{j=1}^n(a_j-b_j+i)^2\\ =& ...

  9. bzoj 4827: [Hnoi2017]礼物 [fft]

    4827: [Hnoi2017]礼物 题意:略 以前做的了 化一化式子就是一个卷积和一些常数项 我记着确定调整值还要求一下导... #include <iostream> #include ...

随机推荐

  1. 中秋国庆8天挑战赛 之 挑战8天掌握微信小程序

    中秋国庆8天挑战赛 挑战8天掌握微信小程序 当前学习进度: // 10.1​// 学习内容:​// 10.2​// 学习内容:​// 10.3​// 学习内容:​// 10.4​// 学习内容:​// ...

  2. Java知识系统回顾整理01基础07类和对象01引用

    一.引用的定义 引用的概念,如果一个变量的类型是 类类型,而非基本类型,那么该变量又叫做引用. 二.引用和指向 new Hero(); 代表创建了一个Hero对象 但是也仅仅是创建了一个对象,没有办法 ...

  3. 【题解】CF375D Tree and Queries

    Link \(\text{Solution:}\) 讲实话这题有点烦,不知道为啥改了下\(\text{dfs}\)就过了--原版本\(dfs\)好像没啥错啊-- 其实对于子树问题,我们求出原来树的\( ...

  4. C# Socket TCP发送图片与接收图片

    如果需要查看更多文章,请微信搜索公众号 csharp编程大全,需要进C#交流群群请加微信z438679770,备注进群, 我邀请你进群! ! ! --------------------------- ...

  5. Python基础笔记1-Python读写yaml文件(使用PyYAML库)

    最近在搭建自动化测试项目过程中经常遇到yaml文件的读写,为了方便后续使用,决定记下笔记. 一,YAML 简介 YAML,Yet Another Markup Language的简写,通常用来编写项目 ...

  6. 前端传递的json格式与SpringMVC接收实体类的对应关系

    这篇文章主要是帮助刚刚入行的猿猿尽快适应Restful风格的搬砖生活 @RequestBody注解 基本介绍:@RequestBody主要用来接收前端传递给后端的json字符串中的数据的(请求体中的数 ...

  7. 多测师讲解 _接口自动化框架设计分层思想(001)_高级讲师肖sir

    第一层: 第二层:调用接口层 VOQGWBZYNBOAVZGE

  8. 修改LiveChart的提示显示位置

    问题:修改LiveChart的提示显示位置 摘要:相信WPF开发者在用LiveCharts的时候会有遇到这个需求.就是产品要求折线图的提示要显示的正常点. 需求:如下图所示.原本显示是在下方并且没有小 ...

  9. java 画 哆啦A梦

    package Demo;import javax.imageio.ImageIO;import java.awt.*;import java.awt.image.BufferedImage;impo ...

  10. 编程语言那么多,为什么偏偏是C语言成了大学的必修课?

    谁叫你不幸生在中国了?--何祚庥(中国科学院院士) 这是一本给非计算机专业的大学生的C语言的书."我不是学计算机的,为啥要学C语言?"这个问题每年在中华大地都会被问上几百万次. 被 ...