- 博客(26)
- 资源 (1)
- 收藏
- 关注
RSS订阅
原创 存板子
主席树#include<bits/stdc++.h>using namespace std;#define SIZE 200010struct SegmentTree { int lc, rc; int sum; #define lc(x) tree[x].lc #define rc(x) tree[x].rc #define sum(...
2019-07-19 12:18:25
185
原创 古代猪文
题意就是给你n,g,求g∑d∣n Cndmod 999911659g^{\sum d|n\:C^{d}_{n}}mod\:999911659g∑d∣nCndmod999911659
2019-06-27 19:58:33
333
1
原创 序列终结者
原题地址：序列终结者(做了这道题就不能再作序列了qaq)标记下传的进阶之题和会有负数，初始化要-INF（敲重点!）fhq treap大法好啊！贴代码#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<cstdlib>#i...
2019-06-18 20:03:25
174
原创 fhq treap
旋转！合并！easy！quickly！上代码// luogu-judger-enable-o2#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<cstdlib>#include<ctime>#include&l...
2019-06-16 17:04:18
188
原创 宝石管理系统
fhq treap的练手题原题地址：宝石管理系统没什么好讲的，直接上代码#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<cstdlib>#include<ctime>#include<algorithm...
2019-06-11 21:05:35
177
原创 CF118D Caesar's Legions
As（这个神秘的男人)给我随机分到的题一开题看数据范围就知道不用优化，于是直接给Cy神仙说是红题（LuoGu紫题肯定是被恶意评分了），于是又被Zxh大神仙嘲讽了真的不难，计数Dp，方程易推f[a][b][c][d]f[a][b][c][d]f[a][b][c][d]表示对于前a个士兵，第a个是兵种（0为步兵，1为骑兵），末尾有连续c个该兵种，总共有d个该兵种然后可知如果第a个兵种的方案可以...
2019-06-02 14:42:35
218
1
原创 Simpson&自适应Simpson
Simpson公式：∫lrf(x)dx≈(r−l)(f(l)+f(r)+4f(l+r2))6\int^{r}_{l}f(x)dx\approx\frac{(r-l)(f(l)+f(r)+4f(\frac{l+r}{2}))}{6}∫lrf(x)dx≈6(r−l)(f(l)+f(r)+4f(2l+r))inline double simpson(double l, double r) {...
2019-05-31 21:43:22
315
原创 矩阵求逆
#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<algorithm>using namespace std;typedef long long LL;const int N=405;const LL mod=1e9+7;...
2019-05-29 20:46:48
164
2
原创 [FFT]快速傅里叶变换（分治版）
FFT真好恰！#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<algorithm>#include<complex>using namespace std;#define cp complex<doubl...
2019-05-29 20:30:10
1138
7
原创 快读模板
inline int read() { int s = 0, w = 1; char c = getchar(); for (; !isdigit(c); c = getchar()) if (c == '-') w = -1; for (; isdigit(c); c = getchar()) s = (s << 1) + (s << 3) + (c ^ 48);...
2019-05-27 19:36:53
171
原创 LuoGu4735 最大异或和
不要万能头！！！不要scanf！！！#include<cstdio>#include<cstring>#include<algorithm>#include<iostream>using namespace std;#define SIZE 480010int n, m, tot;int trie[24 * SIZE][2], la...
2019-05-21 21:36:05
120
原创 欧拉函数
欧拉函数欧拉函数是小于xxx的整数中与xxx互质的数的个数，一般用φ(x)\varphi(x)φ(x)表示。特殊的，φ(1)=1\varphi(1)=1φ(1)=1。欧拉函数计算公式φ(x)={x∏i=1n(1−1pi)，若n>11，若n=1\varphi(x)=\left\{\begin{array}{cc}x\prod^{n}_{i=1}(1- \frac{1}{p_i})...
2019-05-12 20:13:44
591
原创 整除分块（数论分块）
问：求∑i=1N⌊Ni⌋,N≤1012问：求\sum^{N}_{i=1}⌊\frac{N}{i}⌋,N\leq10^12问：求∑i=1N⌊iN⌋,N≤1012思考性质:思考性质:思考性质:(1).⌊Ni⌋最多只有2N种取值(1).⌊\frac{N}{i}⌋最多只有2N种取值(1).⌊iN⌋最多只有2N种取值(2).设⌊Ni⌋⌊Ni′⌋(2).设⌊\frac{N}{i}⌋⌊\frac{N...
2019-05-10 21:57:54
480
原创 莫比乌斯函数&莫比乌斯反演
莫比乌斯反演：莫比乌斯反演：莫比乌斯反演：f(n)=∑d∣ng(d)⇒g(n)=∑d∣nf(nd)∗μ(d)f(n)=\sum^{}_{d|n}g(d)\Rightarrow g(n)=\sum^{}_{d|n}f(\frac{n}{d})*\mu(d)f(n)=d∣n∑g(d)⇒g(n)=d∣n∑f(dn)∗μ(d)莫比乌斯函数：莫比乌斯函数：莫比乌斯函数：(1)若d=1,则μ(d)...
2019-05-09 19:50:30
309
原创 积性函数和完全积性函数
若a,b互质且f(ab)=f(a)*f(b),则f(x)是积性函数若a,b不互质且也有f(ab)=f(a)*f(b),则f(x)是完全积性函数
2019-05-09 18:00:55
2393
原创 拉格朗日插值
∑i=1nyi∗(∏i=1nj≠ix−xjxi−xj)\sum^{n}_{i=1}yi*(\prod^{n}_{i=1}j\neq i\frac{x-xj}{xi-xj})i=1∑nyi∗(i=1∏nj̸=ixi−xjx−xj)
2019-05-04 20:51:44
158
转载 c++模拟质点在万有引力下的运动轨迹
FxkkksNBFxkkksNBFxkkksNB女装大佬不同寻常女装大佬不同寻常女装大佬不同寻常她的blog她的原地址PhysicsIStheMostBeautiful\mathcal{Physics IStheMostBeautiful}PhysicsIStheMostBeautiful通过c++模拟质点在万有引力下的运动轨迹0.ForwardItroduction\mathcal{F...
2019-04-23 16:26:34
804
转载 球体积公式推导（积分）
转载自吊打我的全机房最强的女装大佬 ——Fxkkks她的博客她的原地址球的体积刚刚学了定积分的一点皮毛…来玩一玩废话不多说，Let’s start设球的半径为RRR我们把球（我们通过半球来考虑）切成好多好多（nnn片）薄片，就是一个一个的圆，设圆的半径分别为rrr面积S(r)=πr2S(r)=\pi r^2S(r)=πr2到圆心距离为xxx的圆的半径f(x)=R2−x2f(x...
2019-04-23 16:20:59
6689
原创 LuoGu P2503 [HAOI2006]均分数据
原题地址比平衡点还裸的模拟退火模板题，但是窝就是不会打，太菜了n个数分m组，让方差最小，就是对于ai把ai放到最小的xi中（1<=x<=n）对于n个数随机重组退火很迷的我对着大佬的参数一个一个试下来（真开心嘤嘤嘤）走好// luogu-judger-enable-o2#include<bits/stdc++.h>using namespace std;...
2019-04-22 21:09:56
168
转载 卡常火车头
#pragma GCC diagnostic error “-std=c++11”#pragma GCC target(“avx”)#pragma GCC optimize(3)#pragma GCC optimize(“Ofast”)#pragma GCC optimize(“inline”)#pragma GCC optimize("-fgcse")#pragma GCC opti...
2019-04-12 21:12:20
432
原创 状压DP仓促
状压DP学习笔记【Part1 前言&闲谈】作为一个<全机房 && 全浙江 && 全中国>最菜的蒟蒻，自己对状压这个算法真的不怎么熟悉，所以决定花清明的时间好好补一补状压和状压DP这个（想必您们各位大佬都会，一看题目就轻松秒掉的）算法。对于状压，就是通过二进制位运算使一些状态被压缩成一个数，从而化简算法的时间复杂度（毕竟不是人人都可以像xk一...
2019-04-08 14:27:14
188
5
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人