题目大意：给定n个物品，可以用一个/两个/三个不同的物品凑出不同的价值，求每种价值有多少种拼凑方案(顺序不同算一种)
首先搞出这n个物品的母函数a
将a的每项的平方求和得到多项式b
将a的每项的立方求和得到多项式c
那么如果不考虑顺序和重复 那么方案数就是a+b+c
现在考虑顺序和重复后
三个物品的方案数为(a^3-3*a*b+2*c)/6
两个物品的方案数为(a^2-b)/2
一...
分类：
其他好文 时间：
2015-02-04 13:07:41
阅读次数：
308
最近经常怀念大学的数学课程，越来越觉得数学知识是那么的重要~~数项级数与函数项级数快速傅里叶变换概率论与数理统计测度论代数与几何分形与混沌收敛与发散我经常想我的核心竞争力在哪里，可能就是大学时候多花了一些时间在数学上，仅此而已~
分类：
其他好文 时间：
2015-01-27 12:57:34
阅读次数：
165
快速傅里叶变换模板题算法理解请看《算法导论》第30章《多项式与快速傅里叶变换》，至于证明插值唯一性什么的看不懂也没关系啦～只要明白这个过程是怎么算的就ok。递归版：(4252ms 23468kb) 1 //UOJ 34 递归版 2 #include 3 #include 4 #include 5 #...
分类：
其他好文 时间：
2015-01-18 19:42:15
阅读次数：
189
题目大意：给定一个由'a'和'b'构成的字符串，求不连续回文子序列的个数
首先回文一定是将字符串倍增 由于求的是不连续回文子序列的个数 因此我们可以求出总回文子序列的个数，然后减掉连续的
连续的就是回文子串 用Manacher算法可以O(n)求解
不连续的就有些难搞了
首先我们令f[i]表示以i为中心的对称字符对个数
比如s[]=$#a#b#a 那么s[4]='b' f[4]=2
那么...
分类：
编程语言 时间：
2014-12-27 17:40:35
阅读次数：
243
题目大意：给定n个点，第i个点和第j个点之间的库仑力为(qi*qj)/(i-j)^2，定义左侧为正方向，求每个点受的合力与电荷量的比值
题解详见 http://eolv.farbox.com/post/shui-yu-zheng-feng/2014-12-07 我懒得打了- -
#include
#include
#include
#include
#include
#define ...
分类：
其他好文 时间：
2014-12-25 18:29:14
阅读次数：
156
题目大意：给定两个长度为n的序列a和b，求c[k]=Σa[i]*b[i-k]
这东西不是卷积的形式，因此我们将b数组反转，之后就是卷积的形式了
然后就愉♂悦地上FFT吧
#include
#include
#include
#include
#include
#define M 263000
#define PI 3.14159265358979323846264338327950...
分类：
其他好文 时间：
2014-12-25 16:28:41
阅读次数：
148
题目大意：给定两个高精度整数，求两个数的乘积
FFT大法好
系统的complex比手写慢了2.5倍 简直吓死人- -
#include
#include
#include
#include
#include
#define M 131080
#define PI 3.1415926535897932384626433832795028841971
using namespace s...
分类：
其他好文 时间：
2014-12-25 14:40:21
阅读次数：
178
题意：给出n个数qi，给出Fj的定义如下：令Ei=Fi/qi，求Ei. fft的那一堆东西还是背不到啊。。。这次写虽说完全自己写的，但是还是在参见了以前fft程序的情况下调了很久，主要在如下几点写错：1、非递归中内层数组调用中下表忘掉加k 2、每次转换乘的那个数是cos(...)+isin(.......
分类：
其他好文 时间：
2014-12-20 02:06:41
阅读次数：
212
二维傅里叶变换，图像傅里叶变换，图像傅里叶谱...
分类：
其他好文 时间：
2014-12-16 01:21:46
阅读次数：
4888
埋了一天的算导就当我看懂了？。。。目前仅限于学到FFT计算多项式的卷积，什么频域什么东西的那些我都不懂。。。。最让我感到神奇的是复数这个概念，，，好强大。。复数的话我大概懂得这点？：复数有实部和虚部，其中虚部的单位是$i=$，定义为$e=a+bi$，带$i$的是虚部然后当虚部为0时，这个复数就是实数...
分类：
其他好文 时间：
2014-11-26 20:32:19
阅读次数：
285
