[bzoj4827][Hnoi2017]礼物_FFT
礼物 bzoj-4827 Hnoi-2017
题目大意:给定两个长度为$n$的手环,第一个手环上的$n$个权值为$x_i$,第二个为$y_i$。现在我可以同时将所有的$x_i$同时加上自然数$c$。我也可以将第一个手环任意旋转。旋转后每一个$x$对应一个$y$,那么代价为$\sum\limits_{i=0}^{n-1} (x_i-y_i)^2$。求最小代价。
注释:$1\le n\le 10^5$,$0\le maxval \le 100$。
想法:
水题啊.....
推推式子,我们假设就加了$c$,那么总代价为$\sum\limits_{i=0}^{n-1} (x_i-y_i+c)^2$
$\Leftrightarrow \sum\limits_{i=0}^{n-1} ((x_i+c)^2-2y_i\cdot (x_i+c)+y_i^2)$
$\Leftrightarrow \sum\limits_{i=0}^{n-1} (x_i^2+y_i^2)+nc^2+2c\cdot \sum\limits_{i=0}^{n-1}(x_i-y_i)-2\sum\limits_{i=0}^{n-1} x_i\cdot y_i$
我们发现,上述式子分成了三部分:
显然第一部分是定值。
第二部分的话我们设$sum_x=\sum\limits_{i=0}^{n-1} x_i$,$sum_y=\sum\limits_{i=0}^{n-1} y_i$,那么如果$sum_x>sum_y$,显然$c$为0的时候取得最小值,反之用二次函数的对称轴来判断最小值即可。
剩下了第三部分:
我们要求对应位相乘的最大值
还没有旋转的时候我们发现就是对应位相乘的和。
如果进行了一次旋转,那么就相当于$x_1$和$y_{n-1}$配对,剩下的依次对应。
故此我们发现,x序列的开始一段对应乘上$y$的后面一段,$x$剩下的部分对应上$y$开头剩下的那部分。
我们将$y$翻转,发现其实每个答案就是两个卷积的和。
用$FFT$优化然后$O(n)$更新答案即可。
Code:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#define N 50010
using namespace std; typedef double db;
const db pi=acos(-1);
inline void Max(int &x,int y) {x=max(x,y);}
inline char nc() {static char *p1,*p2,buf[100000]; return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;}
int rd() {int x=0; char c=nc(); while(!isdigit(c)) c=nc(); while(isdigit(c)) x=(x<<3)+(x<<1)+(c^48),c=nc(); return x;}
struct cp
{db x,y;cp() {x=y=0;}cp(db x_,db y_) {x=x_,y=y_;}cp operator + (const cp &a) const {return cp(x+a.x,y+a.y);}cp operator - (const cp &a) const {return cp(x-a.x,y-a.y);}cp operator * (const cp &a) const {return cp(x*a.x-y*a.y,x*a.y+y*a.x);}
}a[N<<2],b[N<<2];
void fft(cp *a,int len,int flg)
{int i,j,k,t;cp tmp,w,wn;for(i=k=0;i<len;i++){if(i>k) swap(a[i],a[k]);for(j=len>>1;(k^=j)<j;j>>=1);}for(k=2;k<=len;k<<=1){wn=cp(cos(2*pi*flg/k),sin(2*pi*flg/k));t=k>>1;for(i=0;i<len;i+=k){w=cp(1,0);for(j=i;j<i+t;j++){tmp=a[j+t]*w;a[j+t]=a[j]-tmp;a[j]=a[j]+tmp;w=w*wn;}}}if(flg==-1) for(i=0;i<len;i++) a[i].x/=len;
}
inline int sqr(int x) {return x*x;}
int main()
{int ans=0;int sum1=0,sum2=0;int all1=0,all2=0;int n=rd(),m=rd();for(int i=0;i<n;i++){a[i].x=rd();sum1+=a[i].x; all1+=sqr(a[i].x);}for(int i=n-1;~i;i--){b[i].x=rd();sum2+=b[i].x; all2+=sqr(b[i].x);}// printf("%d %d %d %d\n",sum1,sum2,all1,all2);int len=1; while(len<(n<<1)) len<<=1;fft(a,len,1); fft(b,len,1);for(int i=0;i<len;i++) a[i]=a[i]*b[i];fft(a,len,-1);ans=(int)(a[n-1].x+0.1);for(int i=1;i<n;i++){Max(ans,(int)(a[i-1].x+0.1)+(int)(a[n+i-1].x+0.1));}if(sum1>=sum2) printf("%d\n",all1+all2-2*ans);else{int ANS=all1+all2-2*ans;int dlt=(sum1-sum2)*2;int x=(-1*dlt)/(2*n);ANS+=min(n*sqr(x)+dlt*x,n*sqr(x+1)+dlt*(x+1));cout << ANS << endl ;}return 0;
}
小结:冷静下来还是非常水的。
转载于:https://www.cnblogs.com/ShuraK/p/10012396.html
[bzoj4827][Hnoi2017]礼物_FFT相关推荐
- bzoj4827: [Hnoi2017]礼物 FFT
bzoj4827: [Hnoi2017]礼物 Description 我的室友最近喜欢上了一个可爱的小女生.马上就要到她的生日了,他决定买一对情侣手 环,一个留给自己,一 个送给她.每个手环上各有 n ...
- [BZOJ4827][Hnoi2017]礼物(FFT)
4827: [Hnoi2017]礼物 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 1315 Solved: 915 [Submit][Statu ...
- bzoj千题计划303:bzoj4827: [Hnoi2017]礼物
https://www.lydsy.com/JudgeOnline/problem.php?id=4827 式子化简一下,发现最后只跟 Σ xi*yi 有关 第二个序列反转,就可以用FFT优化 注意: ...
- bzoj4827[Hnoi2017]礼物 (FFT,生成函数)
Description 我的室友最近喜欢上了一个可爱的小女生.马上就要到她的生日了,他决定买一对情侣手 环,一个留给自己,一 个送给她.每个手环上各有 n 个装饰物,并且每个装饰物都有一定的亮度.但是 ...
- (每日一题)P3723 [AH2017/HNOI2017]礼物(经典FFT)
整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 每日一题(莫反 / 多项式 / 母函数 / 群论) 2021.4.15 多项式 FFT Problem ...
- [BZOJ 4827][Hnoi2017]礼物
4827: [Hnoi2017]礼物 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 1091 Solved: 748 [Submit][Statu ...
- P3723 [AH2017/HNOI2017]礼物(FFT)
P3723 [AH2017/HNOI2017]礼物 式子化简 ∑i=1n(xi−yj)2\sum_{i = 1} ^{n} (x_i- y_j) ^2\\ i=1∑n(xi−yj)2 我们对第一 ...
- P3723 [AH2017/HNOI2017]礼物 FFT + 式子化简
传送门 文章目录 题意: 思路: 题意: 思路: 首先可以知道,我们对某个数组加上一个正数数的操作可以转换成对一个数组加上一个任意数,所以我们设变化量为xxx. 对于∑i=1n(ai−bi)2\sum ...
- AH2017/HNOI2017 礼物
P3723 [AH2017/HNOI2017]礼物 题目大意 给两个数列 a a a和 b b b, a a a和 b b b中的数都小于等于 m m m, b b b可以首尾相接地旋转,要选择一个整 ...
最新文章
- 日志服务Flink Connector《支持Exactly Once》
- linux配置ssh免密码,Linux下配置SSH免密通信 - “ssh-keygen”的基本用法
- SimpleUrlHandlerMapping 处理器映射的配置--转
- Docker:Docker的简介、安装、使用方法之详细攻略
- C# == equals 本质理解
- CPU核数跟多线程的关系
- 面向对象的代码生成方法
- [react] 在React中如何引入图片?哪种方式更好?
- 【动态规划】P1018 线性dp:乘积最大
- 陕西西咸新区沣西新城大数据产业发展思考
- STM32F103C8T6之SG90舵机控制
- 跑跑卡丁车手游怎么用电脑玩 跑跑卡丁车模拟器玩法教程
- ngx_http_core_module模块提供的变量
- 2023AB斗篷玩法全攻略
- 淘宝线上线下“出淘”欲打造零售业航母
- JMeter BeanShell 应用
- java控制台打印输出 中文乱码 解决办法
- 初中数学分几个模块_初中数学有几部分
- 为什么HashMap链表长度超过8会转成树结构
- 深职院入职培训心得2021
热门文章
- 【抖音拍摄小技巧】如何短期内获取快速播放流量?
- 计算机网络组建与维护,计算机网络组建与维护大作业.docx
- 从零开始做App 系列之项目立项+预估时间篇
- J2EE开发 考试总结
- python 数据分析(pandas)
- SyncFusion 的免费电子书
- php实战开发管理系统,深入浅出PHP框架Thinkphp实战开发(权限管理、CMS内容管理系统)...
- Java双亲委派机制
- 三星S10系列国内价格确定:5300元起步 最高1万2
- 软工个人项目作业——论文查重系统