小鸟游六花

vanishment this world

一、算法实现

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
#define AHO_CORASICK_ENABLE_WORDS_COUNT
// https://dl.acm.org/doi/pdf/10.1145/360825.360855
// https://oi.wiki/string/ac-automaton/
struct aho_corasick {

static const int ALPHABET_SIZE = 27;

// 字典树(在build的过程中优化为下一个匹配位置)
vector<array<int, ALPHABET_SIZE>> trie_nodes;
// fail的指针
vector<int> fail;
// 启动后缀链接 用于快速匹配所有出现的单词
// https://leetcode.cn/problems/construct-string-with-minimum-cost/description/
#ifdef AHO_CORASICK_ENABLE_SUFFIX_LINK
// 后缀链接(Suffix Link),用于快速寻找fail指针指向的单词尾,0表示nullptr
vector<int> last;
// 字典树节点是否为单词结尾
vector<bool> trie_exists;
#endif

// 快速匹配某个单词是否出现
// https://leetcode.cn/problems/stream-of-characters/description/
#ifdef AHO_CORASICK_FAST_MATCHED
// 是否匹配到了任意一个单词(用于快速匹配)
vector<bool> trie_matched;
#endif

// 匹配每个单词在文章中出现的次数
// https://www.luogu.com.cn/problem/P5357
#ifdef AHO_CORASICK_ENABLE_WORDS_COUNT
// fail的指针构建反向树
vector<vector<int>> fail_tree;
// 每一个前缀被匹配的数量
vector<int> trie_count;
// dfs统计匹配节点 fail_tree时是否访问 用于优化递归
vector<bool> trie_visited;
#endif

aho_corasick() {
clear();
}
void clear() {
// 必须
trie_nodes.clear();
fail.clear();
#ifdef AHO_CORASICK_ENABLE_SUFFIX_LINK
trie_exists.clear();
last.clear();
#endif
#ifdef AHO_CORASICK_FAST_MATCHED
trie_matched.clear();
#endif
#ifdef AHO_CORASICK_ENABLE_WORDS_COUNT
fail_tree.clear();
trie_count.clear();
trie_visited.clear();
#endif
new_node();
}
int new_node() {
trie_nodes.push_back({});
fail.push_back(0);
#ifdef AHO_CORASICK_ENABLE_SUFFIX_LINK
trie_exists.push_back(false);
last.push_back(0);
#endif
#ifdef AHO_CORASICK_FAST_MATCHED
trie_matched.push_back(false);
#endif
#ifdef AHO_CORASICK_ENABLE_WORDS_COUNT
fail_tree.emplace_back();
trie_visited.push_back(false);
trie_count.push_back(0);
#endif
return int(trie_nodes.size() - 1);
}
int insert(const string& s) {
int p = 0;
for (char c : s) {
if (trie_nodes[p][c & 31] == 0) {
trie_nodes[p][c & 31] = new_node();
}
p = trie_nodes[p][c & 31];
}
#ifdef AHO_CORASICK_FAST_MATCHED
trie_matched[p] = true;
#endif
#ifdef AHO_CORASICK_ENABLE_SUFFIX_LINK
trie_exists[p] = true;
#endif
return p;
}

void build() {
queue<int> q;
for (int i = 1; i < ALPHABET_SIZE; ++i) {
int u = trie_nodes[0][i];
if (u != 0) {
q.push(u);
}
fail[u] = 0;
}
while (!q.empty()) {
int u = q.front();
q.pop();
#ifdef AHO_CORASICK_FAST_MATCHED
trie_matched[u] = trie_matched[fail[u]] || trie_matched[u];
#endif
for (int i = 1; i < ALPHABET_SIZE; ++i) {
int v = trie_nodes[u][i];
if (v == 0) {
trie_nodes[u][i] = trie_nodes[fail[u]][i];
} else {
fail[v] = trie_nodes[fail[u]][i];
#ifdef AHO_CORASICK_ENABLE_SUFFIX_LINK
last[v] = trie_exists[fail[v]] ? fail[v] : last[fail[v]];
#endif
q.push(v);
}
}
}
#ifdef AHO_CORASICK_ENABLE_WORDS_COUNT
int n = int(trie_nodes.size());
for (int i = 1; i < n; ++i) {
fail_tree[fail[i]].push_back(i);
}
#endif
}

#ifdef AHO_CORASICK_ENABLE_WORDS_COUNT
void search_count(const string& text) {
fill(trie_count.begin(), trie_count.end(), 0);
fill(trie_visited.begin(), trie_visited.end(), false);
++trie_count[0];
int p = 0;
for (char c : text) {
p = trie_nodes[p][c & 31];
++trie_count[p];
}
}
// 先调用search_count 然后通过这个方法统计词频
int calc_count(int u) {
if (trie_visited[u]) {
return trie_count[u];
}
trie_visited[u] = true;
for (int v : fail_tree[u]) {
trie_count[u] += calc_count(v);
}
return trie_count[u];
}
#endif
};

一、算法实现

1、平衡二叉树

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
// https://www.luogu.com.cn/problem/P3369
template<typename T, int N>
struct AVL {
int root, pos;
int lc[N], rc[N], dep[N], cnt[N], size[N];
T val[N];
void clear() {
// 简化空值处理,0位置作为空哨兵
pos = root = 0;
memset(lc, 0, sizeof(lc));
memset(rc, 0, sizeof(rc));
memset(dep, 0, sizeof(dep));
memset(cnt, 0, sizeof(cnt));
memset(val, 0, sizeof(val));
memset(size, 0, sizeof(size));
}
// 添加元素 允许相同
void add(T x) {
root = add(root, x);
}
// 删除元素 存在多个只删除一个
void remove(T x) {
root = remove(root, x);
}
// 获得元素的排名,定义排名为比当前数小的数的个数+1
int rank(T x) {
return less_cnt(root, x) + 1;
}
// 查询数据结构中排名为i的数(i不应该超过size[root])
T index(int i) {
if (size[root] < i) return -1;
return index(root, i);
}
// 求x的前驱,前驱定义为小于x且最大的数
T predecessor(int x) {
return predecessor(root, x);
}
// 求x的后继,后继定义为大于x且最小的数
T successor(int x) {
return successor(root, x);
}

private:
T predecessor(int r, T x) {
if (r == 0) return INT_MIN;
if (val[r] >= x) return predecessor(lc[r], x);
return max(val[r], predecessor(rc[r], x));
}
T successor(int r, T x) {
if (r == 0) return INT_MAX;
if (val[r] <= x) return successor(rc[r], x);
return min(val[r], successor(lc[r], x));
}
T index(int r, int i) {
if (size[lc[r]] >= i)
return index(lc[r], i);
else if (size[lc[r]] + cnt[r] < i)
return index(rc[r], i - size[lc[r]] - cnt[r]);
return val[r];
}
// 查询有多少个元素比x小
int less_cnt(int r, T x) {
if (r == 0) return 0;
if (val[r] >= x) return less_cnt(lc[r], x);
return cnt[r] + size[lc[r]] + less_cnt(rc[r], x);
}


// --------- Basic Function -------------
// 调整节点信息
void push_up(int r) {
dep[r] = max(dep[lc[r]], dep[rc[r]]) + 1;
size[r] = size[lc[r]] + size[rc[r]] + cnt[r];
}
// 左旋
int left_rotate(int p) {
int r = rc[p];
rc[p] = lc[r];
lc[r] = p;
push_up(p);
push_up(r);
return r;
}
// 右旋
int right_rotate(int p) {
int l = lc[p];
lc[p] = rc[l];
rc[l] = p;
push_up(p);
push_up(l);
return l;
}
// 维护平衡
int maintain(int p) {
int lh = dep[lc[p]];
int rh = dep[rc[p]];
// LL LR RR RL
if (lh - rh > 1) {
if (dep[lc[lc[p]]] >= dep[rc[lc[p]]]) {
p = right_rotate(p);
} else {
lc[p] = left_rotate(lc[p]);
p = right_rotate(p);
}
} else if (rh - lh > 1) {
if (dep[rc[rc[p]]] >= dep[lc[rc[p]]]) {
p = left_rotate(p);
} else {
rc[p] = right_rotate(rc[p]);
p = left_rotate(p);
}
}
return p;
}
// 添加元素
int add(int r, T x) {
if (r == 0) {
cnt[++pos] = 1;
val[pos] = x;
push_up(pos);
return pos;
}
if (val[r] > x) lc[r] = add(lc[r], x);
else if (val[r] < x) rc[r] = add(rc[r], x);
else ++cnt[r];
push_up(r);
return maintain(r);
}
// 删除元素
int remove(int r, int x) {
if (r == 0) {
return r;
}
if (val[r] > x) {
lc[r] = remove(lc[r], x);
} else if (val[r] < x) {
rc[r] = remove(rc[r], x);
} else {
if (cnt[r] > 1) {
--cnt[r];
push_up(r);
return maintain(r);
}
if (lc[r] == 0 && rc[r] == 0) {
return 0;
}
if (lc[r] == 0 && rc[r] != 0) {
// 用右孩子(后继)替换当前节点
r = rc[r];
} else if (lc[r] != 0 && rc[r] == 0) {
// 用左孩子(前驱)替换当前节点
r = lc[r];
} else {
// 用后继替换当前节点
int s = rc[r];
while (lc[s] != 0) s = lc[s];
rc[r] = remove_most_left(rc[r], s);
lc[s] = lc[r];
rc[s] = rc[r];
r = s;
}
}

push_up(r);
return maintain(r);
}
// 删除r的后继s
int remove_most_left(int r, int s) {
if (r == s) return rc[r];
lc[r] = remove_most_left(lc[r], s);
push_up(r);
return maintain(r);
}
};
AVL<int, 100001> tree;
int main(){
tree.clear();
int n, op, x;
scanf("%d", &n);
while (n--) {
scanf("%d%d", &op, &x);
if (op == 1) tree.add(x);
else if (op == 2) tree.remove(x);
else if (op == 3) printf("%d\n", tree.rank(x));
else if (op == 4) printf("%d\n", tree.index(x));
else if (op == 5) printf("%d\n", tree.predecessor(x));
else printf("%d\n", tree.successor(x));
}
return 0;
}

1、Treap(有旋)

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
// https://www.luogu.com.cn/problem/P3369

// 维护小根堆
template<typename T, int N>
struct Treap {
int root, pos;
int lc[N], rc[N], pri[N], cnt[N], size[N];
T val[N];

void clear() {
// 简化空值处理,0位置作为空哨兵
pos = root = 0;
memset(lc, 0, sizeof(lc));
memset(rc, 0, sizeof(rc));
memset(pri, 0, sizeof(pri));
memset(cnt, 0, sizeof(cnt));
memset(val, 0, sizeof(val));
memset(size, 0, sizeof(size));
pri[0] = INT_MAX;
}
// 添加元素 允许相同
void add(T x) {
root = add(root, x);
}
// 删除元素 存在多个只删除一个
void remove(T x) {
root = remove(root, x);
}
// 获得元素的排名,定义排名为比当前数小的数的个数+1
int rank(T x) {
return less_cnt(root, x) + 1;
}
// 查询数据结构中排名为i的数(i不应该超过size[root])
T index(int i) {
if (size[root] < i) return -1;
return index(root, i);
}
// 求x的前驱,前驱定义为小于x且最大的数
T predecessor(int x) {
return predecessor(root, x);
}
// 求x的后继,后继定义为大于x且最小的数
T successor(int x) {
return successor(root, x);
}

private:
T predecessor(int r, T x) {
if (r == 0) return INT_MIN;
if (val[r] >= x) return predecessor(lc[r], x);
return max(val[r], predecessor(rc[r], x));
}
T successor(int r, T x) {
if (r == 0) return INT_MAX;
if (val[r] <= x) return successor(rc[r], x);
return min(val[r], successor(lc[r], x));
}
T index(int r, int i) {
if (size[lc[r]] >= i) return index(lc[r], i);
else if (size[lc[r]] + cnt[r] < i) return index(rc[r], i - size[lc[r]] - cnt[r]);
return val[r];
}
// 查询有多少个元素比x小
int less_cnt(int r, T x) {
if (r == 0) return 0;
if (val[r] >= x) return less_cnt(lc[r], x);
return cnt[r] + size[lc[r]] + less_cnt(rc[r], x);
}


// --------- Basic Function -------------
// 调整节点信息
void push_up(int r) {
size[r] = size[lc[r]] + size[rc[r]] + cnt[r];
}
// 左旋
int left_rotate(int p) {
int r = rc[p];
rc[p] = lc[r];
lc[r] = p;
push_up(p);
push_up(r);
return r;
}
// 右旋
int right_rotate(int p) {
int l = lc[p];
lc[p] = rc[l];
rc[l] = p;
push_up(p);
push_up(l);
return l;
}
// 添加元素
int add(int r, T x) {
if (r == 0) {
cnt[++pos] = 1;
val[pos] = x;
pri[pos] = rand();
size[pos] = 1;
return pos;
}
if (val[r] > x) {
lc[r] = add(lc[r], x);
if (pri[r] > pri[lc[r]]) r = right_rotate(r);
} else if (val[r] < x) {
rc[r] = add(rc[r], x);
if (pri[r] > pri[rc[r]]) r = left_rotate(r);
} else {
++cnt[r];
}
push_up(r);
return r;
}
// 删除元素
int remove(int r, int x) {
if (r == 0) {
return r;
}
if (val[r] > x) {
lc[r] = remove(lc[r], x);
} else if (val[r] < x) {
rc[r] = remove(rc[r], x);
} else {
if (cnt[r] > 1) {
--cnt[r];
push_up(r);
return r;
}
if (lc[r] == 0 && rc[r] == 0) {
return 0;
}
if (lc[r] == 0 && rc[r] != 0) {
r = rc[r];
} else if (lc[r] != 0 && rc[r] == 0) {
r = lc[r];
} else {
if (pri[lc[r]] < pri[rc[r]]) {
r = right_rotate(r);
rc[r] = remove(rc[r], x);
}
else {
r = left_rotate(r);
lc[r] = remove(lc[r], x);
}
}
}
push_up(r);
return r;
}

};
Treap<int, 100001> tree;
int main() {

tree.clear();
int n, op, x;
scanf("%d", &n);
while (n--) {
scanf("%d%d", &op, &x);
if (op == 1) tree.add(x);
else if (op == 2) tree.remove(x);
else if (op == 3) printf("%d\n", tree.rank(x));
else if (op == 4) printf("%d\n", tree.index(x));
else if (op == 5) printf("%d\n", tree.predecessor(x));
else printf("%d\n", tree.successor(x));
}
return 0;
}

群论基础

一、群的定义及性质

1、定义

1)若$V=<S, ·>$ 是代数系统,如果$·$是可结合的,称$V$为半群。
2)若$V=<S, ·>$ 存在单位元$e$,则称$V$是幺半群或独异点。
3)若$V=<S, ·>$ 是独异点, 若$\forall a \in S$ 有 $a^{-1} \in S$, 则称$V$为群, 通常记为$G$。
4)若群$G$是有穷集, 则称$G$为有限群,否则称为无限群,群$G$的基数成为群$G$阶。
5)只含有单位元的群称作平凡群。
6)若群$G$中的二元运算是可交换的,则称$G$为交换群或阿贝尔(Abel)群。
7)设$G$是群,$a \in G$, 使得等式$a^{k} = e$成立的最小正整数$k$称为$a$的阶(或者周期),记作$|a|=k$, 这时也称$a$为$k$阶元,若不存在这样的正整数$k$,则称$a$为无限阶元。

2、定理

1)$\forall a \in G, a^{-1}$ 唯一。
2)$\forall a \in G, (a^{-1})^{-1} = a$。
3)$\forall a,b \in G, (ab)^{-1} = b^{-1}a^{-1}$。
4)$\forall a \in G, n,m \in Z, (a)^{n}(a)^{m} =a^{n+m}$。
5)$\forall a \in G, n,m \in Z,(a^{n})^{m} =a^{nm}$。
6)若$G$为交换群,则$(ab)^{n} = a^{n}b^{n}$。
7)$\forall a,b,c \in G$,若$ab = ac$ 则 $b=c$。
8)$\forall a,b,c \in G$,若$ba = ca$ 则 $b=c$。
9)若$a \in G$,且$|a|=r$, 设$k$是整数、则$a^{k} = e$ 当且仅当 $r|k$。
10)若$a \in G$,且$|a|=r$, 则$|a^{-1}| = |a|$。

1-8根据定义容易证明,下面证明9:
充分性:因为$r|k, a^{r}=e$, 所以$k = mr$, 则 $a^{k} = a^{mr} = (a^{r})^{m} = e^{m} = e$。
必要性:因为$a^{k}=e, a^{r}=e,r \le k$,所以有$k=mr+i(0 \le i \lt r), e = a^{k} = a^{mr+i} = a^{mr}a^{i}$,$i$ 必须等于0。

下面证明10:

​ $(a^{-1})^{r} = e$, 设$|a^{-1}| = t$, 则 $t|r$, 这说明$a$的逆元的阶是$a$的阶的因子,又$a$的逆和$a$互逆,所以有 $r|t$, 所以$r=t$。

​ 也可以考虑对称展开(在例题1有所体现):

​ $(a^{-1})^{r} = e$, 设$|a^{-1}| = t$, 得到 $t|r$, 又$((a^{-1})^{-1})^{t}=((a^{-1})^{t})^{-1}=e, |((a^{-1})^{-1})| = r$ 得到$r=t$。

​ 因此 $|a^{-1}| = |a|$。

3、例题

1)$\forall a,b \in G$ 是有限阶元,证明:$|b^{-1}ab|=|a|$。

证明:

​ 设$|a|=r, |b^{-1}ab|=t$,则有:

$$
\begin{alignedat}{1}
(b^{-1}ab)^r &= (b^{-1}ab)(b^{-1}ab)…(b^{-1}ab) \
&= (b^{-1}a^{r}b) \
&= e
\end{alignedat}
$$
因此$t|r$。

​ 另一方面,$a=b(b^{-1}ab)b^{-1}$,则:

$$
\begin{alignedat}{1}
(a)^t &= (b(b^{-1}ab)b^{-1})(b(b^{-1}ab)b^{-1})…(b(b^{-1}ab)b^{-1}) \
&= (b^{-1}(b^{-1}ab)^{t}b) \
&= e
\end{alignedat}
$$

因此$r|t$。

所以 $|b^{-1}ab|=|a|$。

2)$G$ 是有限阶元,证明:$G$ 中阶大于2的元素有偶数个。

证明:显然阶大于2的元素与它的逆不相同且成对出现。

二、子群与群的陪集分解

1、定义

1)设$G$是群,$H$是$G$的非空子集,如果$H$关于$G$中的运算构成群,则称$H$是$G$的子群,记作$H \le G$。
2)若$H$是$G$的子群,且$H \subset G$,则称$H$是$G$的真子群,记作$H \lt G$。
3)设$G$是群,$a \in G$,$a$所有的幂构成的集合称作由$a$生成的子群,记作$\lt a \gt$。
4)设$G$是群,$a \in G, \forall b \in G, ab=ba$,所有$a$构成的集合称作群$G$的中心。
5)若$H$是$G$的子群,且$a \in G$,令$Ha={ha|h \in H}$ 称$Ha$是子群$H$在$G$中的右陪集。称$a$为$Ha$的代表元素。
6)若$H$是$G$的子群,且$a \in G$,令$aH={ah|h \in H}$ 称$aH$是子群$H$在$G$中的左陪集。称$a$为$aH$的代表元素。

2、 定理

1)设$G$是群,$H$是$G$的非空子集,$H$是$G$的子群当且仅当:$\forall a,b \in H$ 有 $ab\in H$且$\forall a \in H$有$a^{-1} \in H$。
2)设$G$是群,$H$是$G$的非空子集,$H$是$G$的子群当且仅当:$\forall a,b \in H$ 有 $ab^{-1}\in H$。
3)设$G$是群,$H$是$G$的非空子集,$H$是有穷集,$H$是$G$的子群当且仅当:$\forall a,b \in H$ 有 $ab\in H$。
4)设$G$是群,$H$是$G$的非空子集,$H$是有穷集,$H$是$G$的子群当且仅当:$\forall a,b \in H$ 有 $ab\in H$。
5)设$ G $ 是一个有限群,$ H $ 是 $ G $ 的子群,则 $ |H| $(子群 $ H $ 的阶)是 $ |G| $(群 $ G $ 的阶)的约数,即: $|G| = [G : H] \cdot |H|$ 其中 $ [G : H] $ 是 $ H $ 在 $ G $ 中的 左陪集的个数(拉格朗日定理,不妨设$|H|=m$, 显然$|aH|=m$)。

3、例题

1)设$G$是群,$H$,$K$是$G$的子群, 证明$ H \cap G$也是$G$的子群。

易证,这说明交集部分只包含那些既符合 $H$的封闭性要求、又符合 $K$的封闭性要求的元素。

2)设$G$是群,$H$,$K$是$G$的子群, 证明$ H \cup G$也是$G$的子群,则$ H \subseteq K$ 或 $ K \subseteq H$ 。

证明:任取$h \in H, k \in K, k \notin H$, 若$ H \subseteq K$ 或 $ K \subseteq H$ 不成立,则$hk \notin H$,如若不然,则$hk \in H$ 则 $h^{-1}hk = k \in H$, 同理可得 $hk \notin K$, 矛盾,故结论成立。这说明了两个子群的并集能成为一个新的子群,只有当这两个子群之间有包含关系时。

0%