题目链接:

https://www.luogu.org/problemnew/show/P3723

题目:

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

但是在她生日的前一天,我的室友突然发现他好像拿错了一个手环,而且已经没时间去更换它了!他只能使用一种特殊的方法,将其中一个手环中所有装饰物的亮度增加一个相同的自然数 c(即非负整数)。并且由于这个手环是一个圆,可以以任意的角度旋转它,但是由于上面装饰物的方向是固定的,所以手环不能翻转。需要在经过亮度改造和旋转之后,使得两个手环的差异值最小。

在将两个手环旋转且装饰物对齐了之后,从对齐的某个位置开始逆时针方向对装饰物编号1,2,…,n,其中 n 为每个手环的装饰物个数, 第 1 个手环的 i 号位置装饰物亮度为 xi,第 2 个手环的 i 号位置装饰物亮度为 yi,两个手环之间的差异值为(参见输入输出样例和样例解释):

$\sum_{i=1}^{n} (x_i-y_i)^2$

麻烦你帮他计算一下,进行调整(亮度改造和旋转),使得两个手环之间的差异值最小,这个最小值是多少呢?

题解:

差异值=$\sum_{i=1}^{n}(a_i+x-b_i)^2$

$\sum_{i=1}^na_i^2+\sum_{i=1}^{n}b_i^2+nx^2+2x(\sum_{i=1}{n}a_i-\sum_{i=1}{n}b_i)-2\sum_{i=1}^{n}a_ib_i$

我们枚举$x(-m<=x<=m)$,发现除了最后一项都是定值

那么我们另最后一项最大即可

由于$a$,$b$其实都是可以旋转的,那么我们把$a$倍长

末项$=\sum_{i=x}^{n+x-1}a_ib_{i-x+1}$

再按照套路把$b$反向,$b_i=b_{n-i+1}$

末项$=\sum_{i=x}^{n-x+1}a_ib_{n-i+x}$

$=\sum_{i=1}^{n}a_{i-1+x}b_{n-i+1}$

这显然是一个卷积的形式,即把$a$与$b$卷起来的第$n+x$项

代码:

#include<algorithm>
#include<cstring>
#include<cstdio>
#include<iostream>
#include<cmath>
using namespace std;
typedef double db;
typedef long long ll; const int N=1e6+;
const ll inf=1e18;
const db pi=acos(-1.0);
int r[N];
struct complex
{
db x,y;
complex (db xx=,db yy=) {x=xx;y=yy;}
}A[N],B[N];
complex operator + (complex a,complex b) {return complex(a.x+b.x,a.y+b.y);}
complex operator - (complex a,complex b) {return complex(a.x-b.x,a.y-b.y);}
complex operator * (complex a,complex b) {return complex(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x);}
inline int read()
{
char ch=getchar();int s=,f=;
while (ch<''||ch>'') {if (ch=='-') f=-;ch=getchar();}
while (ch>=''&&ch<='') {s=(s<<)+(s<<)+ch-'';ch=getchar();}
return s*f;
}
void fft(int limit,complex *a,int type)
{
for (int i=;i<limit;i++) if (i<r[i]) swap(a[i],a[r[i]]);
for (int len=;len<limit;len<<=)
{
complex wn=complex(cos(pi/len),type*sin(pi/len));
for (int k=;k<limit;k+=(len<<))
{
complex w=complex(,);
for (int l=;l<len;l++,w=w*wn)
{
complex Nx=a[k+l],Ny=w*a[k+len+l];
a[k+l]=Nx+Ny;
a[k+len+l]=Nx-Ny;
}
}
}
}
int n,m;
int a[N],b[N];
int main()
{
ll a1=,b1=,a2=,b2=;
n=read();m=read();
for (int i=;i<=n;i++) a[i]=read(),a1+=a[i],a2+=a[i]*a[i];
for (int i=;i<=n;i++) b[i]=read(),b1+=b[i],b2+=b[i]*b[i]; for (int i=;i<=n;i++)
{
A[i].x=A[i+n].x=a[i];
B[i].x=b[n-i+];
} int limit=,l=;
while (limit<n+n+n) limit<<=,++l;
for (int i=;i<limit;i++) r[i]=(r[i>>]>>)|((i&)<<(l-)); fft(limit,A,);fft(limit,B,);
for (int i=;i<=limit;i++) A[i]=A[i]*B[i];
fft(limit,A,-);
for (int i=;i<=limit;i++) A[i].x=(ll)(A[i].x/limit+0.5); ll ans=inf;
for (int x=;x<=n;x++)
for (int z=-m;z<=m;z++)
ans=min(ans,a2+b2+n*z*z+2ll*z*(a1-b1)-2ll*(ll)A[x+n].x);
printf("%lld\n",ans);
return ;
}

[AH2017/HNOI2017] 礼物 解题报告 (FFT)的更多相关文章

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

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

  2. [Luogu P3723] [AH2017/HNOI2017]礼物 (FFT 卷积)

    题面 传送门:洛咕 Solution 调得我头大,我好菜啊 好吧,我们来颓柿子吧: 我们可以只旋转其中一个手环.对于亮度的问题,因为可以在两个串上增加亮度,我们也可以看做是可以为负数的. 所以说,我们 ...

  3. P3723 [AH2017/HNOI2017]礼物

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

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

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

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

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

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

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

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

    [Luogu3723] [DarkBZOJ4827] 题解 首先,有一个结论:两个手环增加非负整数亮度,等于其中一个增加一个整数亮度(可以为负) 设增加亮度为x.求\(\sum_{i=1}^{n}(a ...

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

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

  9. [AH2017/HNOI2017]礼物

    题解: 水题 化简一波式子会发现就是个二次函数再加上一个常数 而只有常数中的-2sigma(xiyi)是随移动而变化的 所以只要o(1)求出二次函数最大值然后搞出sigma(xiyi)就可以了 这个东 ...

随机推荐

  1. 利用JDBC或者事物或者调用存储过程实现往MySQL插入百万级数据

    转自:http://www.cnblogs.com/fnz0/p/5713102.html 想往某个表中插入几百万条数据做下测试, 原先的想法,直接写个循环10W次随便插入点数据试试吧,好吧,我真的很 ...

  2. Python学习历程之面对对象浅识

    # ===============================封装====================================# class Bar:# def __init__(se ...

  3. Makefile经典教程(转)

    转自:http://blog.csdn.net/ruglcc/article/details/7814546/ makefile很重要 什么是makefile?或许很多Winodws的程序员都不知道这 ...

  4. HD-ACM算法专攻系列(23)——Crixalis's Equipment

    题目描述: AC源码:此次考察贪心算法,解题思路:贪心的原则是使留下的空间最大,优先选择Bi与Ai差值最大的,至于为什么?这里用只有2个设备为例,(A1,B1)与(A2,B2),假设先搬运A1,搬运的 ...

  5. 获取类似QQ似的时间,昨天或具体日期

    最近在做一个聊天功能,并且要在用户列表上显示最后聊天时间,类似QQ的日期显示. 问群里和百度后,群里没人鸟我,网上也没搜到,最后苦于无奈只能自己封装了. 不过话说回来了,大哥与小弟的区别就是大哥写好封 ...

  6. 解决问题方法:没有设置对象,app.Config没有配置

  7. 达夫设备之js

    最近阅读<高性能JavaScript>时,在书中的“达夫设备“ . 对此,有些感悟,同时有些疑问,希望看到的朋友,能帮忙解释下,在此先提前感谢了. 1. 先说自己的理解吧: ”达夫设备“的 ...

  8. Broadcast Receiver广播接收器

    1.概述 广播接收器不仅能接受来自系统的内容,也可以接受来自其他app的内容.广播分为标准广播和有序广播. 2.标准广播 一种完全异步执行的广播,在广播发出之后几乎所有的广播接收器都在同一时刻接受到广 ...

  9. SQL Server-简单查询语句,疑惑篇

      前言 对于一些原理性文章园中已有大量的文章尤其是关于索引这一块,我也是花费大量时间去学习,对于了解索引原理对于后续理解查询计划和性能调优有很大的帮助,而我们只是一些内容进行概括和总结,这一节我们开 ...

  10. zabbix监控自身为监控机(server)

    Zabbix 监控主机 添加自身为被监控机 这里我的环境已经部署完毕(没有经验的小伙伴可以看我上一篇文章) 等待一会刷新如下页面 字体出现以下亮度就成功了