blank_monster-CSDN博客

原创 ICPC2021 济南站总结

大一省赛打星，所以这场济南算是我队第一次比赛了。 PKU出题真是离谱。一题手快铜，两题手快银。我队属于两题手慢的那个，最后铜首没拿银。虽从结果上来看，第一场就有铜，还是不错的。但还是落下了许多遗憾，包括C可以更快写出来， J也是可以做的。这锅我背了，数学场把队友给演了，数学选手简直崩溃。吸取上海vp时的教训，我们开场三人三线开题，我从后往前看，连着看了L、M、K三题都没什么想法，就去跟榜了。这时C、D、K都有人过，而K过的人更多。我就再回去看K了， K的大致意思是，一个人B

2021-11-16 21:12:30 1417

原创【模板】 Splay树

Messenger(板子题) #include <bits/stdc++.h> using namespace std; const int N = 2e5 + 10; const int inf = 0x3f3f3f3f; char str[N]; namespace Splay{ struct Node{ int son[2]; int fa, val, size, cnt; int rev; }t[N]; int root, tot; void push_up(

2021-11-07 19:00:39 146

原创【模板】点分治

点分治有两种统计答案的方案方案1：分别处理子树在处理完一颗子树的后，再将节点统一加进备选答案中，可以避免同颗子树上的两点相互干扰洛谷P3806 【模板】点分治1 #include<bits/stdc++.h> using namespace std; const int N = 1e4 + 10; const int LIMIT = 1e8 + 10; const int inf = 0x3f3f3f3f; int m; bitset<LIMIT> exist; int

2021-10-19 22:50:44 202

原创【模板】高斯消元

高斯消元

2021-10-16 18:11:11 139

原创【模板】可持久化trie

最大异或和 #include <bits/stdc++.h> using namespace std; #define rep(i, a, b) for(int (i) = (a); (i) <= (b); (i)++) typedef long long ll; const int N = 6e5 + 10; const int inf = 0x3f3f3f3f; const int LEN_OF_BIT = 24; int s[N]; namespace trie_persiste

2021-09-16 19:48:00 113

原创【模板】Treap

普通平衡树(模板题) #include <bits/stdc++.h> using namespace std; typedef long long ll; #define rep(i, a, b ) for(int (i)=(a);(i)<=(b);(i)++) const int N = 1e5 + 10; const int inf = 0x3f3f3f3f; namespace Treap{ struct Node{ int l, r; i

2021-08-26 16:55:18 71

原创【模板】莫队算法

小Z的礼物 #include <bits/stdc++.h> using namespace std; typedef long long ll; #define rep(i, a, b ) for(int (i)=(a);(i)<=(b);(i)++) const int N = 5e4 + 10; int n, m, a[N], col[N], fz[N], fm[N]; int gcd(int a, int b){ return b == 0 ? a : gcd(b, a

2021-08-25 21:53:12 103

原创【模板】线段树，权值线段树

P3372 【模板】线段树 1 #include <bits/stdc++.h> using namespace std; typedef long long ll; #define rep(i, a, b ) for(int (i)=(a);(i)<=(b);(i)++) #define lson pos << 1 #define rson pos << 1 | 1 const int N = 1e5 + 10; ll a[N]; //struct Segmen

2021-08-14 20:32:36 91

原创拉格朗日插值

简易版（只要求f(x)的值，不要求求出系数还原多项式） P4781 【模板】拉格朗日插值 #include <bits/stdc++.h> using namespace std; typedef long long ll; #define rep(i, a, b ) for(int (i) = (a); (i) <= (b); (i)++) const int N = 2e3 + 10; const int mod = 998244353; ll x[N], y[N]; ll fas

2021-08-14 09:47:13 76

原创【模板】FWT

还没弄懂FMT， FWT间得区别待填坑 P4717 【模板】快速莫比乌斯/沃尔什变换 (FMT/FWT) #include <bits/stdc++.h> using namespace std; typedef long long ll; #define rep(i, a, b ) for(int (i) = (a); (i) <= (b); (i)++) const int N = (1 << 17) + 10; const int mod = 998244353; ll

2021-08-14 09:31:06 114

原创【模板】树链剖分

O(nlog2n)O(nlog^2n)O(nlog2n) 【模板】树链剖分 #include <bits/stdc++.h> using namespace std; typedef long long ll; #define rep(i, a, b ) for(int (i) = (a); (i) <= (b); (i)++) #define lson pos << 1 #define rson pos << 1 | 1 const int N = 1e6 + 1

2021-08-13 08:56:36 103

原创【模板】网络流

Dinic算法(O(n2m)O(n^2m)O(n2m)，二分图O(mn)O(m\sqrt{n})O(mn)),可跑边数1e4~1e5的数据小M的作物 #include <bits/stdc++.h> using namespace std; typedef long long ll; #define rep(i, a, b) for(int i = a; i <= b; i++) const int N = 3e3 + 10; const int M = 2e6 + 2000 + 1

2021-08-11 23:55:47 88

原创 Min_25筛解题报告

找素数要求求出1e11以内的素数个数这题有趣的是代码只求了g函数不用求S函数，试分析原因 g函数的意义是：前n个数中所有素数的贡献Σi=1n[i∈Prime]f(i)\Sigma_{i = 1}^{n}[i \in Prime]f(i)Σi=1n[i∈Prime]f(i) 而S函数是：前n个数的贡献Σi=1nf(i)\Sigma_{i = 1}^{n}f(i)Σi=1nf(i)（即积性函数的前缀和）提要求素数个数，可构造函数f(x)=[x∈Prime]f(x) = [x \in Prime]f(

2021-08-10 13:20:24 167

原创单调队列的做题思路（自看）

挺难理解单调队列的，就写些自己平时思考到的小结论，可能对解题有点帮助不妨设队列是个单调递增队列,若当前遍历到第i个元素，队列尾是第j个元素（可能已进行部分出栈操作） 1.1 a[k] > a[j], j < k < i 1.2 用que[tail]表示队尾元素编号,若a[i] < a[que[tail]], 则a[i] < a[que[tail - 1] + 1]. 也可单开一个left数组, 存大于a[j]的连续左区间内的最大值 Feel good ...

2021-07-20 13:48:43 129

原创线段树的应用

已知线段树最基本的操作是对于数组的单点修改和区间查询，因此有了以下推广：之一：对一棵树进行单点修改和子树查询（不论这棵树是否为二叉树，不论节点是否顺序编号）：时间复杂度O(logn)即线段树的操作时间，退化成链也一样原因，线段树查询需要知道l 和r，当节点顺序编号时， l就是该节点的编号，所以我们只需知道r即可完成查询操作。将整棵树进行一遍dfs，就有r = l + size[l] - 1，其中size[l]为以l为根节点的子树大小(包括l) 而当节点不是顺序编号时，另开两个数组idx[], 和

2021-04-14 11:27:01 233

原创算术分解定理取根号的注意点

很多时候，对于一个数据范围很大的数，如n∈\in∈(1, 2e9)分解质因数，若循环n次就会超时，所以通常做法是循环n\sqrt{n}n次，再有因数的性质，算出大于n\sqrt{n}n的情况，然而此方法有些部分需要引起注意，以下代码为例 void find_prime_factor(int n) { for (int i = 2; i <= sqrt(n); i++) { if (n % i == 0) { prime[++cnt] = i; while (n % i == 0)

2021-03-12 22:55:23 402

原创最小生成树

Kruskal的优化：以Highways 这题为例思路，点之间两两建边,利用kruskal和并查集，先连上已有的边(将边的两个端点放入一个并查集中)，再将不同集合的边连接直到边数等于n 做题过程1:结构体2个（Point记录点坐标，Edge记录每条边的情况），优先队列(至今为止Kruskal用的好像都是这个?，每两个点之间建边，但是跳过在同意集合的两个点。其中建边的时间复杂度是O(n^2),KruskalO(nlogn)，怎么看都不会超时的情况下，还是因为常数过大TLE了(黄色部分仍可优化) 每两个点之

2021-02-15 08:34:27 147

原创背包

01背包: 一维数组dp方程 for(int j = n; j >= w[i];j--）{ dp[j] = dp[j - w[i] ] +v[i]; } 多重背包: 1：将一个有x份物品的拆成x一份的物品，循环数倍扩大，时间复杂度受不了 2：二进制优化:原理：将一个数x用二进制表示后，每一位的0或1可表示1-x的所有数 int num = 0; for (int i = 1; i <= m; i++) { int a, b, c; scanf("%d%d%d", &a

2021-02-02 22:49:40 105

原创最短路学习笔记

三种算法 floyd 复杂度：O(n^3) 应用场景：多源最短路,非负边邻接矩阵存图最短路 #include <cstdio> #include <cstring> #include <algorithm> using namespace std; const int INF = 0x3f3f3f3f; int n,m; const int MAX_N = 105; int d[MAX_N][MAX_N]; int main() { while(scanf("%d

2021-01-30 22:37:18 89

原创 Hash

二维哈希公式： sum[i][j]=sum[i-1][j]*p+h[i][j-1]*p-h[i-1]j-1]*p *p+a[i][j] 先来看二维矩阵前缀和：如图每个矩阵右下方的数字代表这个矩阵的前缀和下标；我们由数学公式sum4 = sum2 + sum3- sum1 +a[i][j]可得sum4 即sum[i][j] = sum[i-1][j] + sum[i][j-1] -sum[i-1][j-1] +a[i][j] 模仿一维哈希函数有: h[i] = h[i-1] *p1 h[j] = h[j-

2021-01-29 23:17:43 104

原创 dp

不知道是今天的题目有挑战性，还是自己太菜。很多题思路都不清晰。学长也是太看得起我们了，直接放一道EC-final赛的银牌题，成为第一道没人AC的题今天学的迷迷糊糊吧，总结也写不来，状态转移方程也不知道从哪下手。把感想记下来，以后好好补dp吧 Catching Cheaters 这题不难(前提是要找到状态转移方程) dp[i][j]表示在a中取前i个组成一个子串a‘，b中取前j个组成一个子串b’，求a‘和b’的最大score 有1)(a[i]==b[j]时)dp[i][j]=（dp[i-1][j-1] +

2021-01-27 23:07:27 135