- 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: 字符串哈希
This doc has been marked needs-review. Following information are in upcoming update plan.
int ar[maxn];
int sum[maxn];
class BIT {
static inline int lowbit(int x) { return x & -x; }
inline int prefix(int x) {
static int ret;
ret = 0;
while (x) {
ret += cnt[x];
x -= lowbit(x);
}
return ret;
}
public:
int cnt[maxh + 1];
void modify(int x, int c) {
while (x <= maxh) {
cnt[x] += c;
x += lowbit(x);
}
}
int solve(int l, int r) { return prefix(r) - prefix(l - 1); }
} bit;
Usage
// typedef unsigned long long ullong;
int n;
cin >> n;
for (int i = 0; i < n; i++) {
cin >> ar[i];
++ar[i];
}
for (int i = 0; i < n; i++) {
sum[i] = bit.solve(ar[i] + 1, maxh);
bit.modify(ar[i], 1);
}
fill(bit.cnt + 1, bit.cnt + maxh + 1, 0);
for (int i = n - 1; i >= 0; --i) {
sum[i] += bit.solve(1, ar[i] - 1);
bit.modify(ar[i], 1);
}
ullong ret = 0;
for (int i = 0; i < n; i++)
ret += ((ullong)(1 + sum[i]) * sum[i]) >> 1u;
printf("%llu\n", ret);