- B1: 前缀和与差分
- B2: 快速组合
- B3: 二分查找
- B4: 高精度加法
- B4: 高精度除法
- B4: 高精度乘法
- B4: 高精度减法
- D1: 单调栈
- D2: 单调队列
- D3: Spare Table
- D4: 字典树
- D5: 并查集
- D6: 树状数组-限定区间计数
- D6: 树状数组-单点更新区间求和
- D7: 线段树-基础
- G1: 二分图-判定
- G2: 二分图-最大匹配
- G3: Astar k短路 on matrix
- G4: 拓扑排序
- G5: 连通无环无向图的重心
- G6: MST-Kruskal
- G7: MST-Prim
- G8: ShortestPath-BellmanFord
- G9: ShortestPath-Dijkstra
- G9: ShortestPath-Dijkstra Heap OPT
- G10: ShortestPath-SPFA
- G11: ShortestPath-Floyd
- M1: 快速幂与龟速乘
- M2: 欧几里得-最大公约数
- M2: EX欧几里得-翡蜀定理
- M2: EX欧几里得-线性同余方程
- M2: EX欧几里得-乘法逆元
- M3: 费马小定理-乘法逆元
- M4: 分解质因数
- M4: 分解质因数-欧拉函数
- M4: 分解质因数-多数乘积约数计数
- M4: 分解质因数-多数乘积约数求和
- M5: 欧拉筛
- M5: 欧拉筛-质数筛
- M6: 埃氏筛-质数筛
- M7: Stein算法-最大公约数
- M8: 矩阵乘法与快速幂
- O1: 快速排序
- O2: 并归排序
- S1: KMP
- S2: 字符串哈希
typedef long long llong;
static unordered_map< int, int > pfs;
const int mod = 1e9 + 7;
template< class type >
inline int lowbit(const type x) { return -x & x; }
template< class type >
void pfactor(type x) {
static int k;
static type cur;
k = lowbit(x), cur = x / k;
while (k >>= 1u) ++pfs[2];
for (int i = 3; i <= x / i; i += 2)
while (cur % i == 0) ++pfs[i], cur /= i;
if (cur > 1) ++pfs[cur];
}
llong sumdivs(const int x) {
llong ret = 1;
for (const auto &pf : pfs) {
static llong t;
t = 1;
for (int i = 0; i < pf.second; i++) t = (t * pf.first + 1) % mod;
(ret *= t) %= mod;
}
return ret;
}