博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
[BZOJ 3527][Zjoi2014]力(FFT)
阅读量:5872 次
发布时间:2019-06-19

本文共 791 字,大约阅读时间需要 2 分钟。

Description

给出n个数qi,给出Fj的定义如下:
令Ei=Fi/qi,求Ei.

Solution

设fi=qi gi=1/i/i(这里如果写成i*i可能会爆int)

那前半部分就是∑fi*gj-i 发现是一个卷积的形式

后半部分把数组反一下同理

#include
#include
#include
#include
#include
#define MAXN 100005#define pi acos(-1)using namespace std;int n;struct cp{ double r,i; cp(double r=0,double i=0):r(r),i(i){} cp operator + (const cp& x) {
return cp(r+x.r,i+x.i);} cp operator - (const cp& x) {
return cp(r-x.r,i-x.i);} cp operator * (const cp& x) {
return cp(r*x.r-i*x.i,r*x.i+i*x.r);}}a[MAXN*4],b[MAXN*4],c[MAXN*4];void brc(cp* x,int l){ int k=l/2; for(int i=1;i
=j) { k-=j; j>>=1; } if(k

 

转载于:https://www.cnblogs.com/Zars19/p/6833165.html

你可能感兴趣的文章
Do 32-bit build only with XCode 5.1
查看>>
-Java基础-方法
查看>>
IBM Bluemix云计算大会见闻
查看>>
业务安全漏洞挖掘归纳总结
查看>>
ansible批量安装服务器思路
查看>>
设计模式之行为模式(1)-状态、策略、责任链、访问者
查看>>
Spring Boot:在Spring Boot中使用定时任务
查看>>
如何在单例模式下禁止init
查看>>
JVM系列三:JVM参数设置、分析
查看>>
MySql Cluster 集成安装,Centos,坑点集锦
查看>>
Arrays.copyOfRange
查看>>
Java CyclicBarrier介绍
查看>>
Solaris 10 x86 Mono 三次折腾准备休战了
查看>>
PE文件感染
查看>>
网站目录爆破的扫描器的思路
查看>>
谷歌黑科技:gVisor轻量级容器运行时沙箱
查看>>
Tengine新增nginx upstream模块的使用
查看>>
第三节 整型和浮点型
查看>>
JDK1.8
查看>>
nohup命令简介
查看>>