题意:给定两个 \(n\) 元环,环上每个点有权值,分别为 \(x_i, y_i\)。定义两个环的差值为

\[\sum_{i=0}^{n-1}{(x_i-y_i)^2}
\]

可以旋转其中的一个环,或者将其中一个环的每种权值加上一个数。求最小化的差值。

Solution: 加数只需要加在一个上面即可(假设可以为负),那么差值可以写成

\[\sum_{i=0}^{n-1}{(x_i-y_{i+k}+c)^2}
\]

我们可以将差值定义为旋转位数\(k\)与加数\(c\)的函数,即 \(f(k,c)\) 。我们现在要找的就是二元函数的峰值。展开得

\[f(k,c)=\sum{x_i^2} + \sum{y_i^2} + nc^2 + 2c \sum{x_i} - 2c\sum{y_i} - 2 \sum{x_i y_{(i+k)\%n}}
\]

我们惊喜地发现没有交叉项,即

\[f(k,c) = g(k) + h(c)
\]

那么我们只需要分别最小化两部分即可。

考虑\(g(k)\),翻转序列\(x\),这是一个循环卷积的形式。我们可以将序列\(y\)扩增一倍来转化为线性卷积。那么此时

\[g(k) = \sum_{i=0}^{n-1}{x_{n-1-i} y_{i+k}}
\]

对应到多项式乘法上,找\(g(k)\)的最小值,即在幂次为 \([n-1,2n-1)\) 的项中找最大系数即可。

考虑\(h(c)\),由于\(m \leq 100\),暴力枚举即可。

#include<bits/stdc++.h>
#define pi acos(-1)
using namespace std;
struct poly {
typedef complex<double> E;
int n,m;
vector <double> c;
void read(int deg) {
c.resize(deg+1);
for(int i=0;i<=deg;i++) scanf("%lf",&c[i]);
}
void write() {
for(int i=0;i<c.size();i++) printf("%lf ",c[i]);
printf("\n");
}
void fft(E *a,int f,int *R){
for(int i=0;i<n;i++)if(i<R[i])swap(a[i],a[R[i]]);
for(int i=1;i<n;i<<=1){
E wn(cos(pi/i),f*sin(pi/i));
for(int p=i<<1,j=0;j<n;j+=p){
E w(1,0);
for(int k=0;k<i;k++,w*=wn){
E x=a[j+k],y=w*a[j+k+i];
a[j+k]=x+y;a[j+k+i]=x-y;
}
}
}
}
void mul(vector<double> A, vector<double> B){
int *R; E *a,*b; int L=0;
int _s=(A.size()+B.size()+2)<<5;
a=(E*)malloc(_s); b=(E*)malloc(_s); R=(int*)malloc(_s);
memset(a,0,_s);memset(b,0,_s);memset(R,0,_s);
n=A.size()-1; for(int i=0;i<=n;i++) a[i]=A[i];
m=B.size()-1; for(int i=0;i<=m;i++) b[i]=B[i];
m=n+m;for(n=1;n<=m;n<<=1)L++;
for(int i=0;i<n;i++)R[i]=(R[i>>1]>>1)|((i&1)<<(L-1));
fft(a,1,R);fft(b,1,R);
for(int i=0;i<=n;i++)a[i]=a[i]*b[i];
fft(a,-1,R);
c.resize(m+1);
for(int i=0;i<=m;i++)c[i]=a[i].real()/n;
free(a); free(b); free(R);
}
poly mul(poly pa,poly pb) {
poly pc; pc.mul(pa.c,pb.c); return pc;
}
poly operator * (const poly &pb) {
poly ret = mul(*this,pb); return ret;
}
poly operator *= (const poly &pb) {
mul(this->c,pb.c); return *this;
}
}; int n,m,a[100005],b[100005]; int main() {
ios::sync_with_stdio(false);
cin>>n>>m;
for(int i=0;i<n;i++) cin>>a[i];
for(int i=0;i<n;i++) cin>>b[i];
long long sx=0,sy=0;
for(int i=0;i<n;i++) sx+=a[i],sy+=b[i];
long long px=0,py=0;
for(int i=0;i<n;i++) px+=a[i]*a[i],py+=b[i]*b[i];
long long ans = 1e+9;
for(int i=-100;i<=100;i++) {
ans = min(ans, n*i*i+2*i*sx-2*i*sy);
}
ans += px + py;
poly x,y;
x.c.resize(n);
for(int i=0;i<n;i++) x.c[i]=a[n-i-1];
y.c.resize(2*n);
for(int i=0;i<n;i++) y.c[i]=y.c[i+n]=b[i];
poly z=x*y;
long long mx=0;
for(int i=0;i<n;i++) mx=max(mx,(long long)(z.c[n+i-1]+0.5));
cout<<ans-2*mx<<endl;
}

[HNOI2017] 礼物 - 多项式乘法FFT的更多相关文章

  1. 多项式乘法(FFT)学习笔记

    ------------------------------------------本文只探讨多项式乘法(FFT)在信息学中的应用如有错误或不明欢迎指出或提问,在此不胜感激 多项式 1.系数表示法  ...

  2. 【learning】多项式乘法&fft

    [吐槽] 以前一直觉得这个东西十分高端完全不会qwq 但是向lyy.yxq.yww.dtz等dalao们学习之后发现这个东西的代码实现其实极其简洁 于是趁着还没有忘记赶紧来写一篇博 (说起来这篇东西的 ...

  3. 洛谷.3803.[模板]多项式乘法(FFT)

    题目链接:洛谷.LOJ. FFT相关:快速傅里叶变换(FFT)详解.FFT总结.从多项式乘法到快速傅里叶变换. 5.4 又看了一遍,这个也不错. 2019.3.7 叕看了一遍,推荐这个. #inclu ...

  4. @总结 - 1@ 多项式乘法 —— FFT

    目录 @0 - 参考资料@ @1 - 一些概念@ @2 - 傅里叶正变换@ @3 - 傅里叶逆变换@ @4 - 迭代实现 FFT@ @5 - 参考代码实现@ @6 - 快速数论变换 NTT@ @7 - ...

  5. [uoj#34] [洛谷P3803] 多项式乘法(FFT)

    新技能--FFT. 可在 \(O(nlogn)\) 时间内完成多项式在系数表达与点值表达之间的转换. 其中最关键的一点便为单位复数根,有神奇的折半性质. 多项式乘法(即为卷积)的常见形式: \[ C_ ...

  6. BZOJ4827 [Hnoi2017]礼物 多项式 FFT

    原文链接http://www.cnblogs.com/zhouzhendong/p/8823962.html 题目传送门 - BZOJ4827 题意 有两个长为$n$的序列$x$和$y$,序列$x,y ...

  7. UOJ 34 多项式乘法 FFT 模板

    这是一道模板题. 给你两个多项式,请输出乘起来后的多项式. 输入格式 第一行两个整数 nn 和 mm,分别表示两个多项式的次数. 第二行 n+1n+1 个整数,表示第一个多项式的 00 到 nn 次项 ...

  8. 【Luogu3808】多项式乘法FFT(FFT)

    题目戳我 一道模板题 自己尝试证明了大部分... 剩下的还是没太证出来... 所以就是一个模板放在这里 以后再来补东西吧.... #include<iostream> #include&l ...

  9. 【模板】多项式乘法(FFT)

    题目描述 给定一个n次多项式F(x),和一个m次多项式G(x). 请求出F(x)和G(x)的卷积. 输入输出格式 输入格式: 第一行2个正整数n,m. 接下来一行n+1个数字,从低到高表示F(x)的系 ...

随机推荐

  1. spring cloud微服务快速教程之(九) Spring Cloud Alibaba--sentinel-限流、熔断降级

    0.前言 sentinel的限流.降级功能强大,可以在控制面板中任意制定规则,然后推送到微服务中: 可以根据URL单独制定规则,也可以根据资源名批量制定规则: 需要注意的地方是:1.GITHUB文件在 ...

  2. Dart中类的getter和setter

    Dart类Getters和Setter Getters和Setter(也称为访问器和更改器)允许程序分别初始化和检索类字段的值. 使用get关键字定义getter或访问器.Setter或存取器是使用s ...

  3. mysql第八课

    开启事务: START TRANSACTION; 提交事务: COMMIT; 回滚事务: ROLLBACK; 事务的概念:原子性,一致性,隔离性,持久性 READ UNCOMMITTED(读了未提交) ...

  4. vector内存分配简单介绍

    众所周知,vector的size()其实并不代表它占用的空间,它实际占用空间可以用capacity()查看 众所周知,push_back()时,如果size==capacity则会使capacity从 ...

  5. 吴裕雄--天生自然HADOOP操作实验学习笔记:使用hive操作hbase

    实验目的 熟悉hive和hbase的操作 熟悉hadoop.hbase.hive.zookeeper的关系 熟练大数据环境的搭建 学会分析日志排除问题 实验原理 1.hive整合hbase原理 前面大 ...

  6. form表单中给input 添加 数量可以增减的按钮

    只需给input表单增加type=number即可

  7. JAVA8对象属性的计算

    Men men = new Men(); men.setName("UU"); men.setAge("56"); Men men1 = new Men(); ...

  8. LeetCode 237. 删除链表中的节点 (单链表遍历)

    题目链接:https://leetcode-cn.com/problems/delete-node-in-a-linked-list/ 请编写一个函数,使其可以删除某个链表中给定的(非末尾)节点,你将 ...

  9. 三分钟学会.NET Core Jwt 策略授权认证

    一.前言 大家好我又回来了,前几天讲过一个关于Jwt的身份验证最简单的案例,但是功能还是不够强大,不适用于真正的项目,是的,在真正面对复杂而又苛刻的客户中,我们会不知所措,就现在需要将认证授权这一块也 ...

  10. 自定义结构体 Map

    struct Hashmap{ ,maxe=; ],nxt[maxe+],w[maxe+]; ]; ;;} void Add(int x,int y){son[++E]=y;nxt[E]=lnk[x] ...