The Preliminary Contest for ICPC Asia Shanghai 2019 C Triple（FFT+暴力） 传送门：https://nanti.jisuanke.com/t/41400 题意： 给你三个数组a，b，c，要你求有多少个三元组（i，j，k），使得 $$ \b ...
分类：
其他好文 时间：
2019-09-16 10:01:48
阅读次数：
100
https://www.luogu.org/problem/P1919 ...
分类：
其他好文 时间：
2019-09-13 01:06:43
阅读次数：
83
信噪比(SNR) 有用信号功率与噪声功率的比（此处功率为平均功率），也等于幅度比的平方 $$SNR(dB)=10\log_{10}\frac{\sum_{n=0}^{N-1}s^2(n)}{\sum_{n=0}^{N-1}d^2(n)}=10*\log_{10}(\frac{P_{signal}}{ ...
分类：
其他好文 时间：
2019-09-11 17:53:29
阅读次数：
261
自己整理出来的模板 存在的问题： 1.多项式求逆常数过大(尤其是浮点数FFT) 2.log只支持f[0]=1的情况，exp只支持f[0]=0的情况 有待进一步修改和完善 FFT： 1 #include<bits/stdc++.h> 2 using namespace std; 3 typedef l ...
分类：
其他好文 时间：
2019-09-10 18:00:11
阅读次数：
95
除了之前说的从空间滤波器中获得频率域滤波器，还可以从频率域中直接生成滤波器，这些滤波器被规定为距滤波器中心点的距离不同的函数。可以创建一个用于实现频率滤波器的网格数组，最主要的是需要计算任何点到频率矩形中一个指定点的距离函数，FFT（快速傅里叶）算法是假设变换的原点位于频率矩形的左上角，因此需要将原 ...
分类：
其他好文 时间：
2019-09-01 01:30:09
阅读次数：
121
题目大意：给定长度为 $n 1$ 的序列 $g$，求 $f$ 序列，其中 $f$ 为 $$ f[i]=\sum_{j=1}^{i} f[i j] g[j] $$ 学会了分治 $fft$。 发现这个式子中也含有卷积，但是这是一个递推式，即：$f$ 数组是未知的。 考虑分治策略，即：假设已经算出区间 $ ...
分类：
其他好文 时间：
2019-08-31 00:47:54
阅读次数：
102
FFT Introduction $FFT$ 是一种利用神奇操作，在 $nlog$ 的时间内代替 $n^2$ 的朴素多项式乘法的算法。由 $DFT$ 和 $IDFT$ 两个部分组成。 原理来自于，对于一个 $n$ 次多项式 $f(x)$ ，它一般会被表达为： $$ f(x)=\sum_{i=0}^n ...
分类：
其他好文 时间：
2019-08-30 13:12:51
阅读次数：
73
任意模数FFT 这是一个神奇的魔法，但是和往常一样，在这之前，先 $\texttt{orz}\ \color{orange}{\texttt{matthew99}}$ 问题描述 给定 2 个多项式 $F(x), G(x)$ ，请求出 $F(x) G(x)$。 系数对 p 取模 ，$2 \le p \ ...
分类：
其他好文 时间：
2019-08-29 13:28:35
阅读次数：
67
在教练的要求下开始学习多项式算法了，不过因为不太会积分和求导先把多项式牛顿迭代，多项式指数函数，多项式幂函数，多项式快速幂等内容咕掉了，于是这一篇博客就是其他基础多项式内容的总结。 Fast Fourier Transform $FFT$，快速傅里叶变换，可以在$O(n\log_2n)$的时间内计算 ...
分类：
编程语言 时间：
2019-08-27 23:03:35
阅读次数：
118
"题目" 题目大意 也懒得解释题目大意了…… 正解 正解居然是$FFT$？ 不要看题目的那个式子这么长，也不要在那个式子上下手。 其实我们会发现，不同的$(x_i x_j,y_i y_j,z_i z_j)$并不多。 如果我们求出每个三元组的出现次数，后面的就好做了。 那怎么求呢？ 祭出我们的大杀器— ...
分类：
其他好文 时间：
2019-08-27 22:39:31
阅读次数：
89
