Sankkl1的博客

摘自：https://www.cnblogs.com/LiuRunky/p/Cartesian_Tree.html利用单调栈O(n)O(n)O(n)：#include <cstdio>#include <cstring>#include <algorithm>#include <vector>using namespace std;const int N=100005;int n;int a[N];int root;int ls[N

题目链接：https://codeforces.com/problemset/problem/1437/E分析简单地说，就是求被几个确定的数分成的每个区间内的最长上升子序列，我们可以对每个区间求一遍最长上升子序列。题目中会有约束：区间的两边是已经确定的数。可以先去除已经在这两个数范围之外的数，因为这样的数是一定要进行操作来修改的。那么如何求最长上升子序列呢，可以参考文章，时间复杂度为O(nlogn)O(nlogn)O(nlogn)。但是要将一个数列修改成严格单调是不能仅仅看最长上升子序列的，例如[

摘自：https://blog.youkuaiyun.com/qq_37774171/article/details/81203890假设存在一个序列d[1…9] ={ 2，1 ，5 ，3 ，6，4， 8 ，9， 7}，可以看出来它的LIS长度为5。下面一步一步试着找出它。我们定义一个序列B，然后令 i = 1 to 9 逐个考察这个序列。此外，我们用一个变量Len来记录现在最长算到多少了首先，把d[1]有序地放到B里，令B[1] = 2，就是说当只有1一个数字2的时候，长度为1的LIS的最小末尾是2。这时L

摘自百度百科：https://baike.baidu.com/item/lucas/4326261?fr=aladdinint Lucas (ll n , ll m , int p) { return m == 0 ? 1 : 1ll*comb (n%p , m%p , p) * Lucas (n/p,m/p,p)%p ;}//comb()函数中，因为q ， r < p , 所以这部分暴力完成即可。//C++求C(n, m) mod 10007 版本二 要求p z在100000左右

题目链接：https://codeforces.com/problemset/problem/1430/E分析先将字符串reverse一下，然后从后往前遍历，每次遍历一个字符，找到这个字符在原串中的最后一个位置，将其移动过来即可。简化一下，相当于每次将某个字符的最右边一个移动到最右边（用过以后要去掉），先记录每个字符的位置pospospos，然后记录这个字符前面有xxx个字符已经被移走了，答案就加上n−pos−xn-pos-xn−pos−x。每个位置的xxx值可以用树状数组来维护，维护区间修改。代

欧拉函数及其性质：https://blog.youkuaiyun.com/qq_37493070/article/details/81988725

题目链接：https://acm.hdu.edu.cn/showproblem.php?pid=7084代码#include<bits/stdc++.h>using namespace std;const int N = 3e5 + 10;int q, len;char s[N];int ne[N];vector<int> g1[N], g2[N];int dfn1[N], dfn2[N], cnt1, cnt2, rr1[N], rr2[N], re[N], h

#include<bits/stdc++.h>using namespace std;const int N = 1e5;char s[N], t[N];int ne[N], len;void getnext(){ int j = 0, k = -1; ne[0] = -1; len = strlen(t); while(j < len) { if(k == -1 || t[j] == t[k]) { j++; k++; ne[j] =

关系：https://www.cnblogs.com/Milkor/p/4464515.html相关应用：https://www.cnblogs.com/Milkor/p/4474835.html

题目链接：https://acm.hdu.edu.cn/showproblem.php?pid=7073主要是为了学习一下mt随机的用法代码#include<bits/stdc++.h>using namespace std;typedef long long ll;const int N=2000000;const int maxn=2e6+10;const int K=30;mt19937_64 gen(time(0));int tot,prime[maxn],_,n,a

维护可持久化数组题目链接：https://www.luogu.com.cn/problem/P3919#include <bits/stdc++.h>using namespace std;typedef long long ll;const int N = 1000005;int n,m;int a[N];int val[N*20], ls[N * 20], rs[N * 20], cnt, t[N * 20];void build(int &rt, int l,

题目链接：https://ac.nowcoder.com/acm/contest/11260/E分析先跑一边dfs处理出每个点的dfs序就能得到以每个点为根节点的子树所拥有的点的dfs序的区间，如图由于是dfs序，每个点的子树所包含的所有点的dfs序是连续的。我们先根据每个查询的lll和rrr，向上找到距离开始点最远的在区间内的点，答案就是以该点为根的子树大小减去该点的dfs序区间内温度小于等于x−1x-1x−1的点的数量。由于查询次数较多，我们可以记录每次要查询的dfs序区间以及x−1x-1x

设斐波那契数列f1=1,f2=1,f3=2...f1=1,f2=1,f3=2...f1=1,f2=1,f3=2... 求Sn mod mS_n \ mod\ mSn​ mod m的值#include<bits/stdc++.h>using namespace std;typedef long long ll;const int N = 3;int n,m;void mul(int c[], int a[], int b[][N]){

摘自：https://blog.youkuaiyun.com/weixin_41162823/article/details/83990051#include<iostream>#include<cstdio>#include<cstring>using namespace std;const int maxn = 105;struct Matrix{ int n,m; int v[maxn][maxn]; Matrix(int n,int m):n(n),m(m)

题目链接：https://acm.hdu.edu.cn/showproblem.php?pid=7063#include<iostream>#include<cstdio>#include<cstring>#include<cstdio>#include<cmath>#include<algorithm>#define esp 1e-8const double pi=acos(-1.0);using namespac

摘自：https://www.luogu.com.cn/problem/solution/P3808#include<cstdio>#include<cstring>#include<algorithm> #include<queue>#define maxn 5000000+10using namespace std;char str[maxn*2];struct node{ int fail;//失配指针； int cnt;//单词出现

设数列y=(an+b)qn−1y=(an+b)q^{n-1}y=(an+b)qn−1那么∑1ny=(An+B)qn−B\sum_1^ny=(An+B)q^n-B∑1n​y=(An+B)qn−B其中A=aq−1A=\frac{a}{q-1}A=q−1a​B=b−Aq−1B=\frac{b-A}{q-1}B=q−1b−A​

题目链接：https://acm.hdu.edu.cn/showproblem.php?pid=7020分析如果用线段树会超时，可以用树状数组，参考：https://www.luogu.com.cn/problem/solution/P4062线段树代码（超时）：#include<bits/stdc++.h>using namespace std;const int N=1e6+10;vector<int>g[N];int a[N],ans[N];struct no

题目链接：https://ac.nowcoder.com/acm/contest/11257/H分析每一行的全部扫描线都加入后，判断线段树总点的sum是否是小于d的，如果是则答案就在这一行中，递归去找答案即可。代码#include<bits/stdc++.h>using namespace std;//#define debug 1const int N=1e6+10;typedef long long LL;int lazy[N];struct node1{ i

摘至 OI Wiki#include <stdio.h>#include <stdlib.h>int main() { // For Windows //对拍时不开文件输入输出 //当然，这段程序也可以改写成批处理的形式 while (true) { system("gen > test.in"); //数据生成器将生成数据写入输入文件 system("test1.exe < test.in > a.out"); //获取程

转自OI Wiki

定义如果图G（有向图或者无向图）中所有边一次仅且一次行遍所有顶点的通路称作欧拉通路。如果图G中所有边一次仅且一次行遍所有顶点的回路称作欧拉回路。具有欧拉回路的图称为欧拉图（简称E图）。具有欧拉通路但不具有欧拉回路的图称为半欧拉图。存在的判断判断欧拉通路有向图：图连通，有一个顶点出度大入度1，有一个顶点入度大出度1，其余都是出度=入度。无向图：图连通，只有两个顶点是奇数度，其余都是偶数度的。判断欧拉回路有向图：图连通，所有的顶点出度=入度。无向图：图连通，所有顶点都是偶数度。

题目链接：https://acm.hdu.edu.cn/showproblem.php?pid=6988分析二分答案，再check答案：所有后缀的所有前缀即所有子串，利用前缀和二分找到小于等于当前答案的数量，去除重复的，即height[i]，注意要取min代码#include <iostream>#include <cstring>#include <algorithm>using namespace std;typedef long long ll;

题目链接：https://www.luogu.com.cn/problem/P4173代码#include<bits/stdc++.h>#define MAXN 2000005#define reg register#define inl inline#define db double#define eps 1e-6using namespace std;const int Mod=998244353;const db Pi=acos(-1.0);struct Comple

求多项式乘法#include<bits/stdc++.h>using namespace std;const int N=5e5+10,M=1<<20;const double PI=acos(-1);struct Complex{ double x,y; Complex operator+ (const Complex &t) const { return {x+t.x,y+t.y}; } Complex

hash//POJ 2417//baby_step giant_step// a^x = b (mod n) n为素数，a,b < n// 求解上式 0 <= x < n的解#include <cmath>#include <cstdio>#include <cstring>#define MOD 76543using namespace std;int hs[MOD], head[MOD], next[MOD], id[MOD],

#include <bits/stdc++.h>using namespace std;#define fi first#define se second#define MP make_pair#define pb push_backtypedef long long ll;int n,m;int a[100005];int res[100005],cnt[100005],ans,unit;void add(int x){ x = a[x]; if(!cnt[x])

题目链接：https://www.acwing.com/problem/content/description/3170/在二维平面上有 n 个点，第 i 个点的坐标为 (xi,yi)。请你找出一个点，使得该点到这 n 个点的距离之和最小。该点可以选择在平面中的任意位置，甚至与这 n 个点的位置重合。输入格式第一行包含一个整数 n。接下来 n 行，每行包含两个整数 xi,yi，表示其中一个点的位置坐标。输出格式输出最小距离和，答案四舍五入取整。数据范围1≤n≤100,0≤xi,yi≤1

#include<cstdio>#include<cmath>#include<algorithm>#include<cstring>#include<iostream>#define CLR(a,b) memset(a,b,sizeof(a));const int inf=0x3f3f3f3f;using namespace std;const double PI = acos(-1.0);typedef unsigned lo

#include <bits/stdc++.h>using namespace std;typedef long long ll;#define N 10000006#define fi first#define se second#define pb push_backint vis[N],prime[N],cnt;int phi[N];void Euler_phi(int n){ phi[1] = 1; for(int i=2;i<=n;i++) {

#include <bits/stdc++.h>using namespace std;typedef long long ll;#define N 10006#define fi first#define se second#define pb push_backint x,y;int exgcd(int a, int b){ if(b == 0) { x = 1; y = 0; return a; } int res =

文章目录大数乘法大数加法大数减法大数除法大数乘法#include<bits/stdc++.h>using namespace std;char a[10005],b[10005];int x[10005],y[10005],z[20005];int main(){ int l1,l2,i,j,k; while(~scanf(" %s %s",a,b)) { memset(z,0,sizeof(z)); l1=strlen(a)

∑i=1nni \sum_{i=1}^{n} \frac{n}{i} i=1∑n​in​当我们要求上面式子的时候，很自然想到 O(n) 复杂的解法，但是如果在多组输入的情况下，这样的复杂度可能不够优秀，这时我们就要用到整除分块的方法。for(int l=1,r;l<=n;l=r+1){ r=n/(n/l); ans+=(r-l+1)*(n/l);}...

单调队列的元素保证是单调递增或者是单调递减的，即队头元素是最大值或最小值。#include <vector>#include <deque>using namespace std;class MonotonicQueue{ private: deque<int>data; public: void push(int n) { while(!data.empty()&&

多重背包问题有 N 种物品和一个容量是 V 的背包。第 i 种物品最多有 si 件，每件体积是 vi，价值是 wi。求解将哪些物品装入背包，可使物品体积总和不超过背包容量，且价值总和最大。输出最大价值。基础解法include<iostream>#include<algorithm>using namespace std;const int array_size = 101;int f[array_size], N, V, v, w, s;struct Goods

#include<bits/stdc++.h>#define MAXN 200100using namespace std;int n,m,size_disc; int n_init[MAXN],n_disc[MAXN]; //n_init:原始数组; n_disc:离散化数组int rt[MAXN],lc[MAXN<<5],rc[MAXN<<5],sum[MAXN<<5]; //rt:不同版本的根节点; lc/rc:左儿子右儿子

复杂度：O(n^1/4)链接：https://blog.youkuaiyun.com/water_zero_saber/article/details/79901588#include<cstdio>#include<cstring>#include<cmath>#include<algorithm>#include<iostream>#include<queue>using namespace std; typedef lon

#include<bits/stdc++.h>using namespace std;#define re register int#define ll long long#define INF 0x3f3f3f3f#define maxn 50005#define maxm 500005inline ll read(){ ll x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')f=

https://www.cnblogs.com/nublity/p/11620792.htmlhttps://blog.youkuaiyun.com/wu_tongtong/article/details/79362339?utm_medium=distribute.pc_relevant.none-task-blog-baidujs_utm_term-0&spm=1001.2101.3001.4242#pragma GCC optimize(2)#pragma GCC optimize(3, "Ofas

详解链接：https://blog.youkuaiyun.com/viafcccy/article/details/87483567#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>#define PI 3.1415926535using namespace std;struct node{ int x,y;};