BZOJ3257 [Zjoi2014]力 多项式 FFT
原文链接http://www.cnblogs.com/zhouzhendong/p/8762639.html
题目传送门 - BZOJ3527
题意
给出长度为$m$的序列$q_{1..m}$,让你输出长度为$m$的序列$E_{1..m}$。
其中:
$$E_i=\sum_{j=1}^{i-1}\frac{q_j}{(i-j)^2}-\sum_{j=i+1}^{m}\frac{q_j}{(i-j)^2}$$
题解
我们设
$$f_i=q_i,g_i=\frac 1{i^2}(g_0=0,且对于i>m,g_i=0),f_i^r=f_{m-i}$$
于是:
$$E_i=\sum_{j=0}^i f_jg_{i-j}-\sum_{j=i}^m f_jg_{j-i}$$
这个式子的前一半是个裸的卷积,可以直接$FFT$。后一半要稍微变变。(有dalao说是两个裸的卷积!Orz)
$$\sum_{j=i}^m f_{j}g_{j-i}\\=\sum_{j=i}^{m}f_{m-j}^rg_{i-j}\\=\sum_{j=0}^{m-i}f_{m-i-j}^rg_{j}$$
令$i'=m-i$,则:
$$\sum_{j=0}^{m-i}f_{m-i-j}^rg_{j}\\=\sum_{j=0}^{i'}f_{i'-j}^rg_{j}\\=\sum_{j=0}^{i'}f_{j}^rg_{i'-j}$$
于是也是一个裸的卷积形式了。
但是这题我做了很久。QAQ。
一开始把$m$和$n$搞错,而且还很自信的认为是等价的。($n$见代码)
然后发现$IDFT$之后忘记把数字除以$n$了。
然后还是挂了。
然后看(%)了看(%)网上AC的代码,发现我和他唯一的区别就是他把$i>m$的$g_i$都搞成$0$,而我没搞,但是我仍然很自信的认为是对的QAQ。
最后把我的自信全部打翻,然后就突然A掉了QAQ。
代码
#include <bits/stdc++.h>
using namespace std;
const int N=1<<19;
int m,n,L,R[N];
double q[N],PI=acos(-1.0);
struct C{
double r,i;
C(){}
C(double a,double b){r=a,i=b;}
C operator + (C x){return C(r+x.r,i+x.i);}
C operator - (C x){return C(r-x.r,i-x.i);}
C operator * (C x){return C(r*x.r-i*x.i,i*x.r+r*x.i);}
}A[N],B[N],ans1[N],ans2[N],w[N];
void FFT(C a[]){
for (int i=0;i<n;i++)
if (i<R[i])
swap(a[i],a[R[i]]);
for (int d=1,t=n>>1;d<n;d<<=1,t>>=1)
for (int i=0;i<n;i+=(d<<1))
for (int j=0;j<d;j++){
C tmp=w[t*j]*a[i+j+d];
a[i+j+d]=a[i+j]-tmp;
a[i+j]=a[i+j]+tmp;
}
}
void FFT_times(C A[],C B[],C C[]){
FFT(A),FFT(B);
for (int i=0;i<n;i++)
w[i].i*=-1.0,C[i]=A[i]*B[i];
FFT(C);
for (int i=0;i<n;i++)
w[i].i*=-1.0,C[i].r/=n;
}
int main(){
scanf("%d",&m);
for (int i=1;i<=m;i++)
scanf("%lf",&q[i]);
for (n=1,L=0;n<m*2;n<<=1,L++);
for (int i=0;i<n;i++){
w[i]=C(cos(2*i*PI/n),sin(2*i*PI/n));
R[i]=(R[i>>1]>>1)|((i&1)<<(L-1));
A[i]=C(q[i],0);
B[i]=C((i==0||i>m)?0:1.0/i/i,0);
}
FFT_times(A,B,ans1);
for (int i=0;i<n;i++){
A[i]=C(i<=m?q[m-i]:0,0);
B[i]=C((i==0||i>m)?0:1.0/i/i,0);
}
FFT_times(A,B,ans2);
for (int i=1;i<=m;i++)
printf("%.3lf\n",ans1[i].r-ans2[m-i].r);
return 0;
}
BZOJ3257 [Zjoi2014]力 多项式 FFT的更多相关文章
- 【BZOJ 3527】 3527: [Zjoi2014]力 (FFT)
3527: [Zjoi2014]力 Time Limit: 30 Sec Memory Limit: 256 MBSec Special JudgeSubmit: 2003 Solved: 11 ...
- [Luogu P3338] [ZJOI2014]力 (数论 FFT 卷积)
题面 传送门: 洛咕 BZOJ Solution 写到脑壳疼,我好菜啊 我们来颓柿子吧 \(F_j=\sum_{i<j}\frac{q_i*q_j}{(i-j)^2}-\sum_{i>j} ...
- [ZJOI2014] 力 - 多项式乘法 FFT
题意:给定 \({q_i}\),求 \[E_i = \sum_{i<j}{\frac{q_j}{(j-i)^2}} - \sum_{i>j}{\frac{q_j}{(j-i)^2}}\] ...
- BZOJ 3527: [ZJOI2014]力(FFT)
BZOJ 3527: [ZJOI2014]力(FFT) 题意: 给出\(n\)个数\(q_i\),给出\(Fj\)的定义如下: \[F_j=\sum \limits _ {i < j} \fra ...
- BZOJ3527 [Zjoi2014]力 【fft】
题目 给出n个数qi,给出Fj的定义如下: 令Ei=Fi/qi,求Ei. 输入格式 第一行一个整数n. 接下来n行每行输入一个数,第i行表示qi. 输出格式 n行,第i行输出Ei.与标准答案误差不超过 ...
- 洛谷P3338 [ZJOI2014]力(FFT)
传送门 题目要求$$E_i=\frac{F_i}{q_i}=\sum_{j=1}^{i-1}\frac{q_j}{(i-j)^2}-\sum_{j=i+1}^n\frac{q_j}{(j-i)^2}$ ...
- 【BZOJ3527】[ZJOI2014] 力(FFT)
题目: BZOJ3527 分析: FFT应用第一题-- 首先很明显能把\(F_j\)约掉,变成: \[E_j=\sum _{i<j} \frac{q_i}{(i-j)^2}-\sum_{i> ...
- [BZOJ3527][ZJOI2014]力:FFT
分析 整理得下式: \[E_i=\sum_{j<i}{\frac{q_i}{(i-j)^2}}-\sum_{j>i}{\frac{q_i}{(i-j)^2}}\] 假设\(n=5\),考虑 ...
- 【BZOJ】3527: [Zjoi2014]力(fft+卷积)
http://www.lydsy.com/JudgeOnline/problem.php?id=3527 好好的一道模板题,我自己被自己坑了好久.. 首先题目看错.......什么玩意.......首 ...
随机推荐
- 【APIO2016】【UOJ205】【LOJ2568】烟花表演 可合并堆
题目大意 有一棵树,每条边都有一个边权,现在你要修改边权,使得修改后根到所有叶子的距离相等. 要求所有边权非负. 修改的代价为\(\lvert\)每条边修改前的边权\(-\)修改后的边权\(\rver ...
- Hdoj 1847.Good Luck in CET-4 Everybody! 题解
Problem Description 大学英语四级考试就要来临了,你是不是在紧张的复习?也许紧张得连短学期的ACM都没工夫练习了,反正我知道的Kiki和Cici都是如此.当然,作为在考场浸润了十几载 ...
- docker_监控
目录 Docker 自带的监控子命令 sysdig Weave Scope cAdvisor Prometheus Prometheus 的架构 多维数据模型 Docker 自带的监控子命令 ps d ...
- 指数型生成函数 及 多项式求ln
指数型生成函数 我们知道普通型生成函数解决的是组合问题,而指数型生成函数解决的是排列问题 对于数列\(\{a_n\}\),我们定义其指数型生成函数为 \[G(x) = a_0 + a_1x + a_2 ...
- 使用vlc打开usb摄像头
打开vlc播放器 可以打开网络串流的方式打开摄像头,但只支持第一个摄像头 这一串地址拼凑方法看下面,下面可以选择摄像头 为什么只支持第一个摄像头可以参考下一篇使用Vlc.DotNet打开摄像头并截图 ...
- qt: 系统默认程序打开文件或者软件;
Qt提供了QDesktopServices类, 可以利用openUrl函数调用默认程序打开文件: 源码参考: #ifdef Q_OS_WIN32 m_szHelpDoc = QString(" ...
- SpringBoot项目@RestController使用 redirect 重定向无效
Spring MVC项目中页面重定向一般使用return "redirect:/other/controller/";即可. 而Spring Boot当我们使用了@RestCont ...
- java代码调用第三方接口
一.利用httpclient来字符串参数(url是第三方接口,不带参数,如:http://192.168.16.200:8081/faceInfo/list,param是url后面所要带的参数) pu ...
- [物理学与PDEs]第3章第5节 一维磁流体力学方程组 5.2 一维磁流体力学方程组的 Lagrange 形式
由 $$\bex \cfrac{\p \rho}{\p t}&+u_1\cfrac{\p \rho}{\p x}+\rho\cfrac{\p u_1}{\p x}=0, \eex$$ 我们可以 ...
- 关于JS中的常用表单验证+正则表达式
一.非空验证 trim:去空格(去掉前后的空格),任何字符串都可以用这个方法.写法为:if(v.trim().length==0),表示如果去掉空格后的字符串的长度为0. <body> & ...