bzoj 4827 [Hnoi2017]礼物—

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=4827

式子就是 \sum_{i=0}^{n-1}(a[ i ] - b[ i+k ] + c)^2 。把 b 翻成两倍后卷积即可。关于 c 的部分是一个二次函数，注意 c 只能是整数！

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define db double
using namespace std;
const int N=5e4+5,M=N*6;
const db pi=acos(-1);
int n,m,sm,len,ca,cb,r[M];
struct cpl{db x,y;}a[M],b[M],I;
cpl operator+ (cpl a,cpl b){return (cpl){a.x+b.x,a.y+b.y};}
cpl operator- (cpl a,cpl b){return (cpl){a.x-b.x,a.y-b.y};}
cpl operator* (cpl a,cpl b){return (cpl){a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x};}
int rdn()
{int ret=0;bool fx=1;char ch=getchar();while(ch>'9'||ch<'0'){if(ch=='-')fx=0;ch=getchar();}while(ch>='0'&&ch<='9') ret=(ret<<3)+(ret<<1)+ch-'0',ch=getchar();return fx?ret:-ret;
}
void fft(cpl *a,bool fx)
{for(int i=0;i<len;i++)if(i<r[i])swap(a[i],a[r[i]]);for(int R=2;R<=len;R<<=1){int m=R>>1;cpl Wn=(cpl){ cos(pi/m),fx?-sin(pi/m):sin(pi/m) };for(int i=0;i<len;i+=R){cpl w=I;for(int j=0;j<m;j++,w=w*Wn){cpl tmp=w*a[i+m+j];a[i+m+j]=a[i+j]-tmp;a[i+j]=a[i+j]+tmp;}}}
}
int main()
{n=rdn();m=rdn(); I.x=1;for(int i=n-1;i>=0;i--){a[i].x=rdn();sm+=a[i].x*a[i].x;ca+=a[i].x;}for(int i=0;i<n;i++){b[i].x=b[i+n].x=rdn();sm+=b[i].x*b[i].x;cb+=b[i].x;}int c=floor((db)(cb-ca)/n),tmp=n*c*c+2*(ca-cb)*c;c++; tmp=min(tmp,n*c*c+2*(ca-cb)*c);sm+=tmp;len=1;for(;len<=n*3;len<<=1);for(int i=0;i<len;i++)r[i]=(r[i>>1]>>1)+((i&1)?len>>1:0);fft(a,0); fft(b,0);for(int i=0;i<len;i++)a[i]=a[i]*b[i];fft(a,1); tmp=0;for(int i=(n<<1)-1;i>=n-1;i--)tmp=max(tmp,int(a[i].x/len+0.5));sm-=tmp<<1;printf("%d\n",sm);return 0;
}

转载于:https://www.cnblogs.com/Narh/p/10024331.html

bzoj 4827 [Hnoi2017]礼物——FFT相关推荐

BZOJ 4827 [Hnoi2017]礼物 ——FFT
题目上要求一个循环卷积的最小值,直接破环成链然后FFT就可以了. 然后考虑计算的式子,可以分成两个部分分开计算. 前半部分FFT,后半部分扫一遍. #include <map> #incl ...
[BZOJ 4827][Hnoi2017]礼物
4827: [Hnoi2017]礼物 Time Limit: 20 Sec Memory Limit: 512 MB Submit: 1091 Solved: 748 [Submit][Statu ...
[4827][Hnoi2017]礼物——FFT
题目大意: 我的室友最近喜欢上了一个可爱的小女生.马上就要到她的生日了,他决定买一对情侣手环,一个留给自己,一个送给她.每个手环上各有 n 个装饰物,并且每个装饰物都有一定的亮度.但是在她生日的前一 ...
【刷题】BZOJ 4827 [Hnoi2017]礼物
Description 我的室友最近喜欢上了一个可爱的小女生.马上就要到她的生日了,他决定买一对情侣手环,一个留给自己,一个送给她.每个手环上各有 n 个装饰物,并且每个装饰物都有一定的亮度.但是在 ...
P3723 [AH2017/HNOI2017]礼物 FFT + 式子化简
传送门文章目录题意: 思路: 题意: 思路: 首先可以知道,我们对某个数组加上一个正数数的操作可以转换成对一个数组加上一个任意数,所以我们设变化量为xxx. 对于∑i=1n(ai−bi)2\sum ...
bzoj4827: [Hnoi2017]礼物 FFT
bzoj4827: [Hnoi2017]礼物 Description 我的室友最近喜欢上了一个可爱的小女生.马上就要到她的生日了,他决定买一对情侣手环,一个留给自己,一个送给她.每个手环上各有 n ...
[AH2017/HNOI2017]礼物(FFT)
题目描述我的室友最近喜欢上了一个可爱的小女生.马上就要到她的生日了,他决定买一对情侣手环,一个留给自己,一个送给她.每个手环上各有 n 个装饰物,并且每个装饰物都有一定的亮度.但是在她生日的前一 ...
4827: [Hnoi2017]礼物
首先知道旋转的一般套路就是开两倍. 然后此时发现可以把 ∑ni=1(ai−bi)2∑i=1n(ai−bi)2\sum_{i=1}^{n}(a_i-b_i)^2化作∑ni=1a2i+∑ni=1b2i−2 ...
bzoj4827[Hnoi2017]礼物 (FFT，生成函数)
Description 我的室友最近喜欢上了一个可爱的小女生.马上就要到她的生日了,他决定买一对情侣手环,一个留给自己,一个送给她.每个手环上各有 n 个装饰物,并且每个装饰物都有一定的亮度.但是 ...

bzoj 4827 [Hnoi2017]礼物——FFT

bzoj 4827 [Hnoi2017]礼物——FFT相关推荐

最新文章

热门文章