摘要：
const int MAXN = 2e5 + 10; vector<int> G[MAXN]; int dis[MAXN], mxdis; void dfs(int u, int p) { dis[u] = dis[p] + 1; if(dis[u] > dis[mxdis]) mxdis = u; 阅读全文
摘要：
不打开O2的话效率不够高，但是这种程度上的优化确实应该留给编译器去做了，学个计算机假如连编译器能做的优化都要手动去做并且牺牲程序的简洁性，那还是早点改行去工地搬砖吧。 const int MAXLOGN = 20; const int MAXN = 1 << MAXLOGN; const int M 阅读全文
摘要：
把正整数n分成若干个正整数之和的方法数。 http://www.51nod.com/Challenge/Problem.html#problemId=1259 动态规划 答案是g[n] const int MAXF = 1e3 + 10; const int MAXG = 2e5 + 10; int 阅读全文
摘要：
斯特林子集数/第二类斯特林数 从把n个不同的小球放到k个相同的盒子里，且每个盒子至少要有一个小球的选法。 快速计算一行斯特林数： 使用这个式子 \({n \brace m} = \frac{1}{m!} \sum\limits_{i=0}^m (-1)^i \binom {m}{i} (m-i)^n 阅读全文
摘要：
边 \((u,v)\) 表示 u 的拓扑序在 v 之前。 因为是先从入度为0的点开始消。 struct Graph { static const int MAXN = 2e5 + 10; int n; vector<int> eout[MAXN]; int din[MAXN]; int topo[M 阅读全文
摘要：
多项式加减乘法 \[ h1(x)=\sum\limits_{i=0}^{n1-1}h_{1i}\cdot x^i; \forall i \ge n1, h1_{i} = 0.\\ h2(x)=\sum\limits_{i=0}^{n2-1}h_{2i}\cdot x^i; \forall i \ge 阅读全文
摘要：
序列 \(a=[a_0,a_1,a_2,a_3,a_4,...]\) 的普通型生成函数，定义为形式幂级数： \(F(x)=\sum\limits_{i=0}^{n}a_ix^i\) 。 若序列 \(a\) 拥有通项公式，那么 \(a_i\) 就等于通项公式。 再设 \(G(x)=\sum\limit 阅读全文
摘要：
最长上升子序列 const int MAXN = 2e5 + 10; int n, a[MAXN], dp[MAXN]; int LIS1() { memset(dp, INF, sizeof(dp[0]) * (n + 1)); for(int i = 1; i <= n; ++i) *lower 阅读全文
摘要：
质因数分解法求Mobius函数 const int MAXN = 1e6 + 10; int p[MAXN], ptop; int pm[MAXN], pk[MAXN], mu[MAXN]; void sieve(int n) { memset(pm, 0, sizeof(pm[0]) * (n + 阅读全文
摘要：
$s_1(n)=\sum\limits_i=\frac{1}{2}n(n+1) \ s_2(n)=\sum\limits_i2=\frac{1}{6}n(n+1)(2n+1)\ s_3(n)=\sum\limits_i3=\frac{1}{4}n2(n+1)^2$ 注意 \(\sum\limits_ 阅读全文