bzoj 4827: [HNOI2017]礼物 (FFT)
一道FFT
然而据说暴力可以水70分
然而我省选的时候看到了直接吓傻了 连暴力都没打
太弱了啊QAQ
emmmm
详细的拆开就看其他题解吧233
最后那一步卷积其实我一直没明白
后来画画图终于懂了
只要把其中一个反过来
多项式乘法的结果中的每一项系数就对应某一个Σx[i] * y[j] 的结果
前面几项是不完全的结果
但是太小了就被忽略啦
代码如下
/**************************************************************
Problem: 4827
User: cminus
Language: C++
Result: Accepted
Time:5644 ms
Memory:24568 kb
****************************************************************/ #include <cstdio>
#include <cmath>
#include <complex>
using namespace std; const int N = ;
typedef long long ll;
typedef complex<double> cp;
const double pi = acos(-1.0);
cp A[N], B[N]; void FFT(cp *y, int n, int type) {
if (n == ) return ;
cp l[n >> ], r[n >> ];
for (int i = ; i <= n; i++)
if (i & ) r[i >> ] = y[i];
else l[i >> ] = y[i];
FFT(l, n >> , type); FFT(r, n >> , type);
cp omegan(cos( * pi / n), sin( * pi * type / n)), omega(, );
for (int i = ; i < n >> ; i++) {
y[i] = l[i] + r[i] * omega;
y[i + (n >> )] = l[i] - r[i] * omega;
omega *= omegan;
}
} int main() {
int n, m, ans = , y = ;
scanf("%d %d", &n, &m);
for (int i = ; i < n; i++) {
int x; scanf("%d", &x);
A[n - i - ] = x;
ans += x * x;
}
for (int i = ; i < n; i++) {
int x; scanf("%d", &x);
B[i] = x;
ans += x * x;
y += (int)B[i].real() - A[n - i - ].real();
}
int n1; for (n1 = ; n1 <= n * ; n1 <<= );
for (int i = ; i < n; i++)
B[i + n] = B[i];
FFT(A, n1, ); FFT(B, n1, );
for (int i = ; i <= n1; i++)
A[i] *= B[i];
FFT(A, n1, -);
int temp = , z = (-y) / n;
for (int i = ; i < n; i++) temp = max(temp, (int)(A[i + n - ].real() / n1 + 0.5));
ans -= temp * ;
temp = z * z * n + y * z * ;
z += ; temp = min(temp, z * z * n + y * z * );
z -= ; temp = min(temp, z * z * n + y * z * );
// 有理有据的精度优化
ans += temp;
printf("%d\n", ans);
return ;
}
bzoj 4827: [HNOI2017]礼物 (FFT)的更多相关文章
- bzoj 4827: [Hnoi2017]礼物 [fft]
4827: [Hnoi2017]礼物 题意:略 以前做的了 化一化式子就是一个卷积和一些常数项 我记着确定调整值还要求一下导... #include <iostream> #include ...
- bzoj 4827 [Hnoi2017]礼物——FFT
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4827 式子就是 \sum_{i=0}^{n-1}(a[ i ] - b[ i+k ] + c ...
- bzoj 4827 [Hnoi2017] 礼物 —— FFT
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4827 首先,旋转对应,可以把 b 序列扩展成2倍,则 a 序列对应到的还是一段区间: 再把 ...
- BZOJ 4827 [Hnoi2017]礼物 ——FFT
题目上要求一个循环卷积的最小值,直接破环成链然后FFT就可以了. 然后考虑计算的式子,可以分成两个部分分开计算. 前半部分FFT,后半部分扫一遍. #include <map> #incl ...
- bzoj 4827: [Hnoi2017]礼物【FFT】
记得FFT要开大数组!!开到快MLE的那种!!我这个就是例子TAT,5e5都RE了 在这题上花的时间太多了,还是FFT不太熟练. 首先看70分的n方做法:从0下标开始存,先n--,把a数组倍增,然后枚 ...
- 【刷题】BZOJ 4827 [Hnoi2017]礼物
Description 我的室友最近喜欢上了一个可爱的小女生.马上就要到她的生日了,他决定买一对情侣手 环,一个留给自己,一个送给她.每个手环上各有 n 个装饰物,并且每个装饰物都有一定的亮度.但是在 ...
- BZOJ 4827: [Hnoi2017]礼物 FFT_多项式_卷积
题解稍后在笔记本中更新 Code: #include <bits/stdc++.h> #define setIO(s) freopen(s".in","r&q ...
- BZOJ:4827: [Hnoi2017]礼物
[问题描述] 我的室友最近喜欢上了一个可爱的小女生.马上就要到她的生日了,他决定买一对情侣手 环,一个留给自己,一个送给她.每个手环上各有 n 个装饰物,并且每个装饰物都有一定的亮度. 但是在她生日的 ...
- 4827: [Hnoi2017]礼物
4827: [Hnoi2017]礼物 链接 分析: 求最小的$\sum_{i=1}^{n}(x_i-y_i)^2$ 设旋转了j位,每一位加上了c. $\sum\limits_{i=1}^{n}(x_{ ...
随机推荐
- CF1300E Water Balance
题目链接 problem 给出一个长度为n的序列,每次可以选择一个区间\([l,r]\)并将区间\([l,r]\)内的数字全部变为这些数字的平均数.该操作可以进行任意多次. 求出进行任意次操作后可以得 ...
- gulp常用插件之http-proxy-middleware使用
更多gulp常用插件使用请访问:gulp常用插件汇总 http-proxy-middleware这是一个用于后台将请求转发给其它服务器.其实这并不是转给gulp使用的,在其它构建工具也可以用. 更多使 ...
- xctf-i-got-id-200(perl网页文件+ARGV上传造成任意文件读取)
打开url发现有三个链接,点进去都是.pl文件,且只有files可以上传文件. .pl文件都是用perl编写的网页文件 这里上传了又将文件的内容全部打印出来,那么猜想后台应该用了param()函数. ...
- unity踩坑2020-01-21
这几天一直在测试一个类似于传奇的2d界面游戏,目前做的测试为: 人物动作响应,主要是8方向的判断和资源文件精灵的刷新. 学到的知识点: 1,Enum.GetHashCode() 可以得到这个枚举的索引 ...
- [Python]jieba切词 添加字典 去除停用词、单字 python 2020.2.10
源码如下: import jieba import io import re #jieba.load_userdict("E:/xinxi2.txt") patton=re.com ...
- Coxeter积分计算
\begin{align*}&\int_0^{\frac{\pi}{3}}{\arccos \left( \frac{1-\cos x}{\text{2}\cos x} \right) dx} ...
- javaWeb快速入门+——初体验-HelloWorld
文章转载自 https://www.cnblogs.com/1906859953Lucas/p/10821840.html 练习成品下载 https://www.lanzous.com/i9fljkj ...
- code码说明
https://www.cnblogs.com/wanglaowu/p/6229843.html
- vscode设置成中文
打开 VS Code Ctrl + Shift +p打开搜索框 搜索框内输入Configure Display Language 回车 修改代码中“locale”后面引号内内容为zh-CH 重新启动V ...
- vue 报错碰到的一些问题及其规范
报错信息:Expected error to be handled(需要处理的错误) 这是因为回调函数里面的参数error没有运用到,所以可以不设置参数,或者在回调函数内console.log(err ...