- 博客(59)
- 收藏
- 关注
RSS订阅原创 线性基模板
线性基：struct LB { ll d[100]; int tot; LB() { tot = 0; memset(d, 0, sizeof(d)); } bool add(ll x) { for (int i = 63; ~i; i--) { if (x & (1ll << i)) { //注意，如果i大于31，前面的1的后面一定要加ll...
2019-11-09 14:26:08
140
原创 向天杰的动态凸包
还记得普通凸包里的那个while循环吗？我们现在只要把这个while循环给改一下，让每次新加入的点 按照极角序 往前 / 后遍历，我们就可以在平均复杂度 O(?) 的情况下求出来新的凸包。 我也不知道具体多少，但大多数情况下很快，至少不用遍历所有点如果题目有额外要求，我们可以直接在插入函数里面进行维护，如本模板就是维护了新凸包的周长。本模板其实并非插入点，而是删除点，大家可以自己...
2019-11-03 11:01:22
188
原创 两圆交点模板 2018icpc沈阳L
链接 ：https://codeforces.com/gym/101955/problem/L本来是过了的…但看了大佬的代码之后开始怀疑自己是怎么卡过去的大佬的代码：#include <bits/stdc++.h>using namespace std;const double Pi=acos(-1);const double eps=1e-12;const i...
2019-10-31 21:59:24
235
原创 HDU 2202 旋转卡壳
HDU 2202#include<bits/stdc++.h>#define ll long longusing namespace std;const int N=5e4+100;const double eps=0.00;struct point { double x,y; point() {} point(double _x,double _y): x(...
2019-10-26 15:50:12
150
原创 幂塔函数 欧拉降幂模板
链接 ：https://codeforces.com/contest/906/problem/D#include<bits/stdc++.h>#define ll long long#define Mod(a,mod) a>=mod?a%mod+mod:ausing namespace std;const int N=1e5+10;ll qpow(ll a,ll...
2019-10-19 14:01:41
269
原创 螺旋矩阵规律
#include<bits/stdc++.h>#define ll long longusing namespace std;const int mod=1e9+7;const int N=3e2+10;int work(int n, int m, int i, int j) {//n行m列的螺旋矩阵，i行j列的序号 if (i == 1) ret...
2019-10-18 16:52:04
949
原创 拉格朗日插值法
当 X 的取值不一定是 1-n 连续的时候，我们只能 O(n²) 地插出结果洛谷 P4781#include<bits/stdc++.h>#define ll long longusing namespace std;const int N=2e3+10;const int mod=998244353;int x[N],y[N];ll quick_pow(ll a...
2019-10-17 14:26:51
258
原创 HDU6608 大素数检测模板
链接 ：http://acm.hdu.edu.cn/showproblem.php?pid=6608具体分析看这里 ： 大佬的博客#include <cstdio>#include <cstdlib>#include <map>using namespace std;long long gcd(long long a,long long b)...
2019-10-17 13:13:29
139
原创 抄的 FFT 和 NTT 模板
洛谷P1919 ：https://www.luogu.org/problem/P1919FFT ：#include<cmath>#include<cstdio>#include<cstring>#include<algorithm>const int N=1e6+10;const double Pi=acos(-1.0);st...
2019-10-10 16:40:16
191
原创 高斯消元 模板
题目链接：https://ac.nowcoder.com/acm/contest/1109/A#include<cstdio>#include<cstring>#include<cstdlib>#include<cmath>#include<ctime>#include<iostream>#include&...
2019-10-04 17:45:26
144
原创 最小球覆盖
2018 南京现场赛 D ：https://nanti.jisuanke.com/t/33673三分：#include <iostream>#include <math.h>#include <stdio.h>#include <algorithm>#include <string.h>using namespace...
2019-09-19 15:34:46
135
原创 判断两直线的位置关系
POJ 1269 ：http://poj.org/problem?id=1269#include <algorithm>#include <iostream>#include <cstring>#include <cstdlib>#include <vector>#include <queue>#includ...
2019-09-17 22:24:07
851
原创 凸包问题
POJ 3348 ：http://poj.org/problem?id=3348#include<stdio.h>#include<iostream>#include<algorithm>#include<string.h>#include<math.h>using namespace std;typedef long ...
2019-09-17 21:43:19
242
原创 计算几何常用公式+模板
const double pi=acos(-1);const double eps=1e-6;//精度int sgn(double a, double b) { if(fabs(a-b)<eps) return 0; if(a<b) return -1; return 1;}struct point { double x,y,z; point() {} point...
2019-09-12 16:42:09
268
原创 HDU6559 几何基础 解方程
题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6559过程太过于麻烦，直接贴大佬博客：https://www.cnblogs.com/Dillonh/p/11196418.html核心方程：#include<bits/stdc++.h>using namespace std;const int maxn=1e5+1...
2019-09-12 14:45:20
292
原创 字符串的最小表示法
对于一个字符串S，我们认为他是循环的，及最后一位的下一位是第一位，其他以此类推字符串的起点是任意的，我们需要找到一个起点，使得从该起点开始的字符串最小，即为一个字符串的最小表示题目链接：https://ac.nowcoder.com/acm/contest/887/A#include <bits/stdc++.h>#define ll long longusing...
2019-09-03 20:15:55
191
原创 蔡勒公式 模板
直接看大佬博客https://blog.csdn.net/SongBai1997/article/details/85227008 https://blog.csdn.net/qq_42804678/article/details/82793438没有证明，背会就完事了...
2019-08-28 11:31:15
209
原创 2019 牛客 第十场G 计算几何 枚举
又是一个看半天看不懂，忽然茅塞顿开的题目，就是一个枚举的思想https://www.cnblogs.com/ccsu-kid/p/11380666.html#include<bits/stdc++.h>#define ll long longusing namespace std;const int maxn=1e6+10;double x[maxn],y[maxn],...
2019-08-27 13:54:40
105
原创 STL max_element min_element nth_element
nth_elemet作为stl类的一个求第k大的数，据说平均复杂度是o(n)nth_element(a+1,a+k,a+1+n),从a[1]到a[n]中找第k小的数同时nth_element会对原数组排序，第k个数刚好为第k大的数，它的左边都比他小，右边都比他大，但左右两边都是无序的max_element min_element 用法一样max_element(a+1,a...
2019-08-27 13:27:44
125
转载 2019 牛客多校 第九场C 神仙思维
大佬博客：https://blog.csdn.net/qq_35950004/article/details/99683772#include<bits/stdc++.h>#define maxn 100005#define mod 1000000007using namespace std;int n,b;int pw[maxn],pws[maxn],pwsm[m...
2019-08-27 12:26:21
90
原创 洛谷 P1447 能量采集
新发现一种巨简洁的做法，利用容斥原理，复杂度为 O(nlogn)#include <bits/stdc++.h>#define ll long longusing namespace std;const int maxn=1e5+10;ll a[maxn];int main(){ ll ans=0,n,m; scanf("%lld%lld",&n,&...
2019-08-26 16:33:41
164
原创 欧拉函数 线性筛
//写在前面//对于欧拉函数有如下三条性质 p为素数 //1、phi(p)=p-1//2、当p与i互质时有: phi(p*i)=phi(p)*phi(i)//3、当i%p==0时有：phi(p*i)=p*phi(i)#include <iostream>#include <string.h>#define maxn 1005using namespace ...
2019-08-26 16:32:22
97
原创 杜教BM 模板
#include <cstdio>#include <cstring>#include <cmath>#include <algorithm>#include <vector>#include <string>#include <map>#include <set>#include &...
2019-08-24 16:52:36
141
原创 欧拉筛 模板 and 2018 ICPC 南京网络赛 J
欧拉筛模板：const int maxn=1e7+10;int cnt=0,prime[maxn];bool isvisited[maxn];void solve(int n){//筛 1-n 的数 for(int i = 2; i < n; i ++) { if(!isvisited[i]) prime[cnt ++] = i; for...
2019-08-24 16:19:47
133
原创 高斯消元 模板
求线性方程的整数解，思路就是线性代数的矩阵化简直接上kuangbin大神的代码 #include<stdio.h> #include<algorithm> #include<iostream> #include<string.h> #include<math.h> using namesp...
2019-08-21 17:26:50
91
转载 LGV 算法 非降路径
主要就是这个定理的内容路径的起点和终点根据题意自行选择通常情况下，这个定理会与非降路径结合起来附上一道例题帮助理解https://blog.csdn.net/weixin_42165981/article/details/81154921其实是我写的太垃圾，不看大佬博客无法理解...
2019-08-21 14:25:47
368
原创 补个SG模板
以前的丢掉了，在这里新补一个#define MAXN 1000 + 10#define N 20int f[N],SG[MAXN],S[MAXN];void getSG(int n){ int i,j; memset(SG,0,sizeof(SG)); for(i = 1; i <= n; i++){ memset(S,0,sizeof(S...
2019-08-21 12:48:45
101
原创 类欧几里得算法 模板
这个算法和欧几里得并无直接联系，它是一个完完全全的新算法，我们省略推导过程，直接记录结论。首先，我们有如下定义： 注意：接下来出现的出发均为向下取整！对于上面这三个式子，我们分别有当a>=c or b>=c的时候，我们有当a<c and b<c的时候，我们有当a>=c or b&...
2019-08-20 16:47:25
341
原创 BSGS算法 扩展BSGS算法 模板
BSGS算法 可以快速的求出，满足a^x ≡ b(mod p)的最小的非负整数x （p必须为质数，即 gcd(a,p) == 1）我们先将x拆分成i*m-j的形式（其中m为sqrt(p)向上取整的值，则原式化为 ≡ b (mod p) ) 显然 0<=j<m 1<=i<=m移向后得 ≡ b*(mod p)我们从 0-m 枚举 j，并将 ...
2019-08-16 15:43:49
430
1
原创 二次剩余定理 模板
某些坑已经挖好…随时回来填补勒让德符号：https://www.bk.gugeso.site/wiki/%E5%8B%92%E8%AE%A9%E5%BE%B7%E7%AC%A6%E5%8F%B7关于雅可比符号的介绍较少，大致意思就是，当p不为质数的时候，我们需要根据唯一分解定理把p分解成数个质数相乘，分别计算答案接下来 参考了这位大佬的博客（基本全抄）https...
2019-08-16 10:43:01
361
原创 卢卡斯定理 扩展卢卡斯定理 模板
卢卡斯定理 模板根据此题所写：https://www.luogu.org/problem/P3807由于此题求的是所以根据具体题意要作出相应的改变 （这个组合数的格式并不标准，但我看着比较顺眼）#include<iostream>#include<algorithm>#include<cstdio>#define ed 100005us...
2019-08-15 10:58:41
267
原创 求逆元的几种方法
Plus ： 当a为正数时，无论n的正负，a%n结果一样 当a为负数时，a%n 的符合人的习惯的结果= a%n+n扩展欧几里得#include<bits/stdc++.h>using namespace std;int exgcd(int a,int b,int &x,int &y)//扩展欧几里得算法{ ...
2019-08-14 13:52:42
682
原创 中国剩余定理 扩展中国剩余定理 模板
中国剩余定理解线性同余线性方程/*long long gcd(LL a,LL b){ return b==0?a:gcd(b,a%b);}*/#include<bits/stdc++.h>#define ll long longconst int maxn=1e5+10;//扩展欧几里得算法 void gcd(ll a,ll b,ll &d,ll ...
2019-08-13 16:50:08
247
原创 组合数 模板
当组合数要求范围不是很大时，我们可以直接数组打表算出void init(){ for(int i=0; i<=3000; i++){ c[i][0]=1; c[i][i]=1; for(int j=1; j<i; j++){ c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod; ...
2019-08-13 12:47:41
162
原创 矩阵快速幂 模板
根据牛客多校所写：https://ac.nowcoder.com/acm/contest/885/B#include<bits/stdc++.h>using namespace std;typedef long long ll;const int maxn=1000100;const int NUM=2;//当前矩阵是几行几列的ll MOD;//取模char s[m...
2019-08-07 11:50:33
108
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人