AωD!

摘要： 这里给出一种基于有限微积分的做法，由于求斯特林数时需要不定模数FNT，于是不太跑的过 …… 先转下降幂 $$\begin{align*}\sum_{t=0}^{n}{t^qq^t}&=\sum_{t=0}^{n}\sum_{p=0}^{q}{q\brace p}t^{\underline{p}}q^ 阅读全文

摘要： ？？？？？ 怎么找一个数列的规律（线性递推）呢？当然就用BM啦！ 估计这个东西我以后也遇不到几次…… 为什么这个东西会出现在模拟赛里？？？ 这个算法有什么用呢？比如说有一道题，在 $m * n$ 的网格上搞一些事情，$m$ 非常小，$n$ 非常大。显然是一个状压dp套矩阵快速幂~~裸题~~。不过呢， 阅读全文

摘要： 要省选啦！码码板子…… const int N = 200010; char st[N]; int rk[N], sa[N], nr[N], ns[N], tab[N], hi[N]; void calc(int n) { for (int i = 0; i <= 26; ++ i) tab[i] 阅读全文

摘要： Bluestein's algorithm 算法可以在\(O (n \log n) \)的时间内完成任意长度的 DFT 考虑DFT，有： \(\begin{align*} y_k &= \sum_{i = 0}^{n - 1} a_i \omega_n^{ki}\\ &= \sum_{i = 0}^ 阅读全文

摘要： 支配树模版题～ 这东西似乎非常好写…… 但是证明好长好长啊～ $$sdom(u) = min\{\{v|(v, u) \in E, dfn(v) < dfn(u)\}\cup\{sdom(q)|dfn(q)>dfn(u), \exists (p, x) \in E, q \in path(S, p) 阅读全文

摘要： 给定平面上的\(n\)个点，求一个\(n - 1\)阶多项式经过这些点…… 我以前应该是只会\(O(n^3)\)的高斯消元的……就是直接把方程列出来直接解的那种…… 考虑拉格朗日插值法： $$F=\sum_{1 \leq i \leq n}y_i\prod_{j\neq i \& 1 \leq j 阅读全文

摘要： Day-2 听说我能来NOI（捂脸熊 xy说我肩负重任（奶的漂亮），要MAKE XJ GREAT AGAIN （大雾 感觉我已经沦为暴力隔膜选手了啊…… Day-1 报到日 早上开super meat boy，把第5张图推的差不多了…… 然而我并没有把它拷出来，绍兴一中网络太慢……惨惨啊…… 没办法 阅读全文

摘要： 我们需要求\(g(P)\)，其中\(g\)是一个只有一项的多项式 暴力是\(n ^ 3 \log k\)的，显然过不去 怎么办？ 特征方程优化矩阵快速幂～ 考虑方程\(|xI - P| = 0\)，把det展开后可以得到一个方程\(f(x) = 0\)，这个方程称为\(P\)的特征方程，\(f\)称 阅读全文

摘要： 维护一个序列，初始全为\(1\) 支持两种操作： 1.对于所有的位置\(i\)，将它的值乘上\(i + a\) 2.询问\(a\)处的值 \(q=120000\) 20s 512M —————— 如果把第一个操作看成乘上一个\(x + a_i\)，第二个操作看成询问\(x = a_i\)处多项式的值 阅读全文

摘要： 怎么样快速提高姿势水平呢？自信说应该做TopCoder！ 于是我就看到了很多XJOI题，BZOJ题……好像可以搬些题呢！ 简单整理一下，方便未来搬题（雾 为什么好多题都没有官方题解啊QAQ 6.23:对比了下jiry的博客，感觉我写的好简略啊QAQ，不过这篇博客只是用来搬题的对吧（捂脸 301C：用 阅读全文