P3338 [ZJOI2014]力(FFT)

P3338 [ZJOI2014]力

Fj=∑i=1j−1qi×qj(i−j)2−∑i=j+1nqi×qj(i−j)2Ej=∑i=1j−1qi(i−j)2−∑i=j+1nqi(i−j)2f(i)=qi,g(i)=1i2,f(0)=0,g(0)=0Ej=∑i=0jf(i)g(j−i)−∑i=jnf(i)g(i−j)F_j = \sum_{i = 1} ^{j - 1} \frac{q_i \times q_j}{(i - j) ^ 2} - \sum_{i = j + 1} ^{n} \frac{q_i \times q_j}{(i - j) ^ 2}\\ E_j = \sum_{i = 1} ^{j - 1} \frac{q_i}{(i - j) ^ 2} - \sum_{i = j + 1} ^{n} \frac{q_i}{(i - j) ^ 2}\\ f(i) = q_i, g(i) = \frac{1}{i ^ 2}, f(0) = 0, g(0) = 0\\ E_j = \sum_{i = 0} ^{j} f(i) g(j - i) - \sum_{i = j} ^{n} f(i) g(i - j)\\ Fj=i=1j1(ij)2qi×qji=j+1n(ij)2qi×qjEj=i=1j1(ij)2qii=j+1n(ij)2qif(i)=qi,g(i)=i21,f(0)=0,g(0)=0Ej=i=0jf(i)g(ji)i=jnf(i)g(ij)

前项是标准的多项式卷积,后项我们另考虑后项。

∑i=jnf(i)g(i−j)\sum\limits_{i = j} ^{n} f(i) g(i - j)i=jnf(i)g(ij)T=i−jT = i - jT=iji=T+ji = T + ji=T+j,原式子∑i=0n−jf(i+j)g(i)\sum\limits_{i = 0} ^{n - j} f(i + j)g(i)i=0njf(i+j)g(i)

我们翻转f(i)f(i)f(i)得到f′(n−i)=f(i)f'(n - i) = f(i)f(ni)=f(i),代入原式子∑i=0n−jf′(n−(i+j))g(i)\sum\limits_{i= 0} ^{n - j} f'(n - (i + j)) g(i)i=0njf(n(i+j))g(i)

n−j=Tn - j= Tnj=T∑i=0Tf′(T−i)g(i)\sum\limits_{i = 0} ^{T} f'(T - i)g(i)i=0Tf(Ti)g(i),也就是一个多项式卷积了。

#include <bits/stdc++.h>using namespace std;struct Complex {double r, i;Complex(double _r = 0, double _i = 0) : r(_r), i(_i) {}
};Complex operator + (const Complex &a, const Complex &b) {return Complex(a.r + b.r, a.i + b.i);
}Complex operator - (const Complex &a, const Complex &b) {return Complex(a.r - b.r, a.i - b.i);
}Complex operator * (const Complex &a, const Complex &b) {return Complex(a.r * b.r - a.i * b.i, a.r * b.i + a.i * b.r);
}Complex operator / (const Complex &a, const Complex &b) {return Complex((a.r * b.r + a.i * b.i) / (b.r * b.r + b.i * b.i), (a.i * b.r - a.r * b.i) / (b.r * b.r + b.i * b.i));
}typedef long long ll;const int N = 1e6 + 10;int r[N];Complex a[N], b[N], c[N];void get_r(int lim) {for (int i = 0; i < lim; i++) {r[i] = (i & 1) * (lim >> 1) + (r[i >> 1] >> 1);}
}void FFT(Complex *f, int lim, int rev) {for (int i = 0; i < lim; i++) {if (i < r[i]) {swap(f[i], f[r[i]]);}}const double pi = acos(-1.0);for (int mid = 1; mid < lim; mid <<= 1) {Complex wn = Complex(cos(pi / mid), rev * sin(pi / mid));for (int len = mid << 1, cur = 0; cur < lim; cur += len) {Complex w = Complex(1, 0);for (int k = 0; k < mid; k++, w = w * wn) {Complex x = f[cur + k], y = w * f[cur + mid + k];f[cur + k] = x + y, f[cur + mid + k] = x - y;}}}if (rev == -1) {for (int i = 0; i < lim; i++) {f[i].r /= lim;}}
}int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);// ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);int n, lim = 1;scanf("%d", &n);for (int i = 1; i <= n; i++) {scanf("%lf", &a[i].r);b[n - i].r = a[i].r;c[i].r = 1.0 / i / i;}n <<= 1;while (lim <= n) {lim <<= 1;}get_r(lim);FFT(a, lim, 1);FFT(b, lim, 1);FFT(c, lim, 1);for (int i = 0; i < lim; i++) {a[i] = a[i] * c[i];b[i] = b[i] * c[i];}FFT(a, lim, -1);FFT(b, lim, -1);n >>= 1;for (int i = 1; i <= n; i++) {printf("%.3f\n", a[i].r - b[n - i].r);}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314123.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

CF755G PolandBall and Many Other Balls(多项式/倍增fft)

CF755G PolandBall and Many Other Balls 倍增fft模版题 这种问题关键在于两个状态转移&#xff0c;一个是单点加1&#xff0c;还有一个是乘2&#xff0c;这样我们就相当于拥有了二进制下的左移操作和加1操作&#xff0c;那么可以在O(logn)的时间内表示出任何一个数。 然后对…

VS, VS Code, VS Online, VS xxx, 你都分清了吗?

首先说说部分童鞋容易混淆的 Visual Studio 和 Visual Studio Code 吧。其实&#xff0c;它们俩的关系&#xff0c;就相当于 Java 和 JavaScript&#xff0c;没啥关系。再说说 Visual Studio Online。这就复杂了。历史上&#xff0c;出现过两个 Visual Studio Online&#xff0…

Java修炼之路——基础篇——Java关键字

1&#xff1a;transient 当对象被序列化时&#xff0c;transient阻止其修饰的对象进行序列化&#xff1b;当反序列化时&#xff0c;此对象的值不会被恢复。 2&#xff1a;instanceof 判断引用指向的对象&#xff0c;是不是某个类及其子类的实例对象&#xff1b; class Person …

.NetCore从零开始使用Skywalking分布式追踪系统

将本文从0开始搭建两个webapi项目&#xff0c;使用Skywalking来追踪他们之间的调用关系及响应时间。开发环境为VisualStudio20191&#xff1a;安装Skywalking,可参考&#xff1a;https://www.cnblogs.com/sunyuliang/p/11422576.html&#xff0c;本列中搭建好后的Skywalking服务…

P3723 [AH2017/HNOI2017]礼物(FFT)

P3723 [AH2017/HNOI2017]礼物 式子化简 ∑i1n(xi−yj)2\sum_{i 1} ^{n} (x_i- y_j) ^2\\ i1∑n​(xi​−yj​)2 我们对第一个手环ccc&#xff0c;相当于(xic−yi)2(x_i c - y_i) ^ 2(xi​c−yi​)2&#xff0c;对第二个手环ccc相当于(xi−yi−c)2(x_i - y_i - c) ^2(xi​−…

CF623E Transforming Sequence(多项式/倍增fft/动态规划)

CF623E Transforming Sequence 经典的倍增NTT题目&#xff0c;但是由于万恶的模数导致这道题变成了倍增MTT 要求n个数前缀或严格递增的序列个数&#xff0c;一共有k位。 然后我们考虑进行dp&#xff0c;然后我的思路就是fi,jf_{i,j}fi,j​表示前i位在k位中有j位的方案数&…

2019-02-26-算法-进化(字符串转换成整数)

题目描述&#xff1a; 请你来实现一个 atoi 函数&#xff0c;使其能将字符串转换成整数。 首先&#xff0c;该函数会根据需要丢弃无用的开头空格字符&#xff0c;直到寻找到第一个非空格的字符为止。 当我们寻找到的第一个非空字符为正或者负号时&#xff0c;则将该符号与之后…

通过Service访问应用 (1)

目录通过Service访问应用 通过Pod IP访问应用 通过ClusterIP Service在集群内部访问 通过Service访问应用通过之前的操作&#xff0c;应用部署完成了&#xff0c;我们的Demo网站已经成功启动了&#xff0c;那么如何访问网站呢&#xff1f;通过Pod IP访问应用我们可以通过Pod IP…

P4239 任意模数多项式乘法逆(多项式/ MTT)

P4239 任意模数多项式乘法逆 这个题目简直就是毒瘤&#xff0c;不过还好我们可以使用vector封装要不然真的没法看&#xff0c;现在我们就会用vector封装MTT了&#xff0c;然后有一个代码细节就是这里的求逆还是在模意义下的&#xff0c;所以我们还是需要求逆。 #include<b…

多项式求逆模板(NTT + mod)

【模板】多项式乘法逆 /*Author : lifehappy */ #include <bits/stdc.h>using namespace std;typedef long long ll;const int N 1e6 10, mod 998244353, G 3;int r[N], n;ll a[N], b[N], c[N];ll quick_pow(ll a, int n) {ll ans 1;while(n) {if(n & 1) ans …

2019-02-26-算法-进化(回文数)

题目描述&#xff1a; 判断一个整数是否是回文数。回文数是指正序&#xff08;从左向右&#xff09;和倒序&#xff08;从右向左&#xff09;读都是一样的整数。 示例 1: 输入: 121 输出: true示例 2: 输入: -121 输出: false 解释: 从左向右读, 为 -121 。 从右向左读, 为…

POJ1742 Coins(DP)

Coins 思路 没分析复杂度写了个二进制拆分&#xff0c;然后做010101背包O(nlog(c)m)>10e7了O(nlog(c)m) > 10e7了O(nlog(c)m)>10e7了&#xff0c;所以还是想办法优化吧。 我们引入一个needneedneed数组&#xff0c;need[j]need[j]need[j]表示&#xff0c;在枚举到第…

13张PPT带你了解主动式消息队列处理集群

前言偷偷和你们说&#xff0c;我搞了一份内部资料&#xff0c;该内部资料共有13张PPT&#xff0c;据作者透露&#xff0c;该PPT至少花了整整1周时间才编写完成&#xff0c;其内容简洁明了&#xff0c;内容深度足够&#xff0c;易于初学者理解&#xff0c;也给深度开发人员分享了…

CF773F Test Data Generation(倍增FFT/动态规划)

CF773F Test Data Generation https://www.luogu.com.cn/problem/CF773F 这个题还是挺巧妙的&#xff0c;最后需要我们求解的实际上值域为a选择奇数个数最大数是奇数的方案数&#xff0c;然后这个东西显然包含了3个信息&#xff0c;值域、个数、最后一个数的奇偶性&#xff0c…

2019-02-27-算法-进化(寻找两个有序数组的中位数)

题目描述 给定两个大小为 m 和 n 的有序数组 nums1 和 nums2。 请你找出这两个有序数组的中位数&#xff0c;并且要求算法的时间复杂度为 O(log(m n))。 你可以假设 nums1 和 nums2 不会同时为空。 示例 1: nums1 [1, 3] nums2 [2]则中位数是 2.0示例 2: nums1 [1, 2]…

.Net Core 三大Redis客户端对比和使用心得

前言稍微复杂一点的互联网项目&#xff0c;技术选型都可能会涉及Redis&#xff0c;.NetCore的生态越发完善&#xff0c;支持.NetCore的Redis客户端越来越多&#xff0c;下面三款常见的Redis客户端&#xff0c;相信大家平时或多或少用到一些&#xff0c;结合平时对三款客户端的使…

2020第十一届蓝桥杯软件类省赛第二场C/C++ 大学 B 组(题解)

试题 A: 门牌制作 问题描述 小蓝要为一条街的住户制作门牌号。 这条街一共有 2020 位住户&#xff0c;门牌号从 1 到 2020 编号。 小蓝制作门牌的方法是先制作 0 到 9 这几个数字字符&#xff0c;最后根据需要将字 符粘贴到门牌上&#xff0c;例如门牌 1017 需要依次粘贴字符 …

队长开卖自家产“翠香”猕猴桃

猕猴桃品种有很多&#xff0c;但不是所有的果子都叫翠香。这两天我在公众号里卖了这个翠香猕猴桃&#xff0c;可能是有同学以为是做广告卖水果&#xff0c;其实是家里的亲戚猕猴桃成熟了&#xff0c;辛苦一年下来地里一共结了3000斤猕猴桃&#xff0c;遇到了一个难题就是如何把…

Java修炼之路——基础篇——Java集合类

集合类的全景图 常用集合类特性 1. Collection&#xff1a;每个位置对应一个元素1.1: List 存放有序元素&#xff0c;允许重复元素&#xff0c;允许元素为null1.1.1: ArrayList&#xff1a;内部结构为数组&#xff1b;初始容量为10&#xff1b;插入、删除的移动速度慢&#x…

1575 Gcd and Lcm

1575 Gcd and Lcm ∑i1n∑j1i∑k1ilcm(gcd(i,j),gcd(i,k))设f(n)∑i1n∑j1nlcm(gcd(i,n),gcd(j,n))f(p)3p2−3p1f(pk)(2k1)(p2k−p2k−1)pk−1\sum_{i 1} ^{n} \sum_{j 1} ^{i} \sum_{k 1} ^{i} lcm(gcd(i, j), gcd(i, k))\\ 设f(n) \sum_{i 1} ^{n} \sum_{j 1} ^{n} lcm…