- 博客(22)
- 收藏
- 关注
RSS订阅
原创 「置顶」NOI To do list
未学习、待重学及待刷题的知识点Treap / Splay / 替罪羊树最大流2-SAT匈牙利算法二分图最大匹配拉格朗日插值FFT 快速傅里叶变换NTT 数论变换Miller-RabinPollard RhoBSGSexBSGS *exLucas *生成函数Burnside引理Polya定理crtexcrt容斥 *概率期望 *高斯消元法 *Matrix...
2018-08-29 21:01:40
495
原创 「SPOJ 5971」LCMSUM
题目传送门题意求∑ni=1lcm(i,n)∑i=1nlcm(i,n)\sum_{i=1}^{n} lcm(i,n)，TTT组询问1≤T≤3×105,1≤n≤1061≤T≤3×105,1≤n≤1061\leq T \leq 3\times 10^5, 1 \leq n \leq 10^6 题解∑i=1nlcm(i,n)∑i=1nlcm(i,n)\sum_{i=1}^{n}...
2018-08-31 21:33:49
335
原创 「BZOJ 2154」Crash的数字表格「莫比乌斯反演」
题目传送门题意求∑ni=1∑mj=1lcm(i,j)∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^{n} \sum_{j=1}^{m} lcm(i,j)，对201010092010100920101009取模题解∑ni=1∑mj=1lcm(i,j)∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^{n} \sum_{j=1}^{m} lcm(i,j)=∑...
2018-08-30 16:13:41
259
原创 「POI2007」ZAP-Queries「莫比乌斯反演」
题目传送门题意T组询问求，∑i=1a∑j=1b[gcd(i,j)=d]∑i=1a∑j=1b[gcd(i,j)=d]\sum_{i=1}^{a} \sum_{j=1}^{b} [gcd(i, j)=d]，T,a,b,d≤5×104T,a,b,d≤5×104T,a,b,d\leq 5 \times 10^4题解不妨设a≤ba≤ba \leq b.∑i=1a∑j=1b[gcd(i...
2018-08-29 16:10:11
250
原创 「学习笔记」矩阵乘法与矩阵快速幂
矩阵乘法定义一个n×mn×mn \times m矩阵AAA乘上一个m×pm×pm\times p的矩阵BBB，得到一个n×pn×pn \times p的矩阵CCC：Ci,j=∑mk=1Ai,k×Bk,j(1≤i≤n,1≤j≤p)Ci,j=∑k=1mAi,k×Bk,j(1≤i≤n,1≤j≤p)C_{i,j}=\sum_{k=1}^{m} A_{i,k} \times B_{k,j}\...
2018-08-28 22:36:13
238
原创 「Luogu P2257」YY的GCD「莫比乌斯反演」
题目传送门题意求 ∑ni=1∑mj=1[gcd(i,j)isprime]∑i=1n∑j=1m[gcd(i,j)isprime]\sum_{i=1}^{n} \sum_{j=1}^{m} [gcd(i, j)\;is\;prime]题解不妨设n≤mn≤mn \leq m，推一下式子∑i=1n∑j=1m[gcd(i,j)isprime]∑i=1n∑j=1m[gcd(i,j)ispr...
2018-08-28 16:34:09
259
原创 「NOI 2018」屠龙勇士「扩展中国剩余定理」
题目传送门题解首先，每条龙使用哪把剑是确定的，我们可以根据题意使用multisetmultisetmultiset（可重集）求出.#include <iostream>#include <cstdio>#include <cmath>#include <set>using namespace std;#
2018-08-28 15:49:52
333
原创 「ZJOI 2008」骑士「基环树DP」
BZOJ 题目传送门 Luogu 题目传送门题意求由基环树组成的森林的最大权独立集，n<=106n<=106nfi,0fi,0f_{i,0}和fi,1fi,1f_{i,1}分别表示在以结点iii为根的子树内，iii选/不选的最大独立集的权值转移：fi,0=∑j∈son(i)max(fj,0,fj,1)fi,0=∑j∈son(i)max(fj,0,fj,1)f_{i,0}=...
2018-08-26 18:00:48
397
原创 「学习笔记」BSGS及扩展
离散对数问题：给定求y,z,p,y,z,p,y,z,p,求 yx≡zyx≡zy^x \equiv z (mod(mod(mod p)p)p)的最小整数解.BSGSShank的大步小步算法(Shank’s Baby-Step-Giant-Step Algorithm)这里介绍一种避开求逆元的BSGS（常数小令m=⌈p–√⌉m=⌈p⌉ m = \lceil\sqrt p\rceil，...
2018-08-26 15:49:48
260
原创 「学习笔记」泰勒级数
多项式函数是长这样的函数：f(x)=a0+a1x+a2x2+…+anxnf(x)=a0+a1x+a2x2+…+anxnf(x)=a_0+a_1x+a_2x^2+\ldots+a_nx^n它有一个很NiceNiceNice的特点：代人xxx，在O(n)O(n)O(n)的时间内就可以求出f(x)f(x)f(x)，没有任何障碍.但是这样的函数：g(x)=exg(x)=exg(x)...
2018-08-26 15:15:34
794
原创 「学习笔记」级数与检验法
Warning：这篇文章都是定理公式若级数∑∞n=1an=S∑n=1∞an=S\sum_{n=1}^{\infty} a_n=S，则说这个级数收敛到SSS。否则说该级数发散。级数1+1+1+…1+1+1+…1+1+1+\ldots的发散的.级数12+14+18+…12+14+18+…\frac{1}{2}+\frac{1}{4}+\frac{1}{8}+\ldots是收敛的，收...
2018-08-26 15:09:47
1252
原创 「学习笔记」自适应辛普森法
引例计算积分：∫RLcx+dax+bdx∫LRcx+dax+bdx\int_{L}^{R} \frac{cx+d}{ax+b} dx题解其实可以直接求导然后做，然而太弱根本推不出来.于是就可以使用自适应辛普森积分。辛普森（Simpson）公式是牛顿-科特斯公式当n=2时的情形。∫RLf(x)dx≈(R−L)[f(L)+f(R)+4f(mid)]6∫LRf(x)dx≈...
2018-08-26 15:07:24
830
原创 「学习笔记」牛顿迭代法
牛顿迭代法（Newton′smethodNewton′smethodNewton's method）一般用于求函数的一个零点。牛顿迭代法牛顿迭代法三个步骤：随机 / 猜一个ppp值求x=px=px=p时的切线，即求导数令p=p=p=切线零点，返回步骤222，重复若干次.即：xn+1=xn−f(xn)f′(xn)xn+1=xn−f(xn)f′(xn)x_{n+1} = x...
2018-08-26 14:51:34
923
原创 「学习笔记」基环树
基环树是一种图，它由一个环组成，环上每个点都是一棵树点树根，所以称为基环树。当然，一棵树上连一条边也会变成基环树。下图是一个基环树.基环树一般分成环和树来分别处理（显然环的处理较为麻烦），那首先得找到环.找环大概就是dfs一下，找到一个在此结点之前走过的相邻结点就开始记录环.vector<int> G[MAXN]; //基环树int fa[MAXN]; ...
2018-08-26 14:34:33
9749
3
原创 「CF535D」Tavas and Malekas「KMP」
题意一个有小写字母组成的长度为nnn文本串，给出匹配串及匹配的mmm个位置（按顺序给出），求原串有多少种可能.题解可以转化为判定这些匹配点有没有冲突。有冲突就是000，没有冲突就是26k26k26^k，kkk是没有被匹配覆盖的位置个数.显然两个匹配有冲突他们一定有重叠。重叠的一定是匹配串的后缀和前缀。所以假设两个匹配位置他们重叠的部分长度为xxx，我们只需要判断匹配串长度为xxx的...
2018-08-23 16:39:40
255
原创 「CF126B」 Password「KMP」
题意在字符串中找一个子串，满足既是前缀也是后缀，还在中间出现过.题解首先预处理出KMPKMPKMP的nextnextnext数组，然后从nextnnextnnext_n开始跳，不停让i=nextii=nextii=next_i。考虑nextnnextnnext_n的含义：前缀等于后缀。每跳一次相当于这个前缀等于后缀的长度减小了。每次枚举的时候使用KMPKMPKMP匹配中间位置，...
2018-08-23 15:07:04
332
原创 「Luogu P2398」GCD SUM「莫比乌斯反演」
题目传送门题意求∑ni=1∑nj=1gcd(i,j)∑i=1n∑j=1ngcd(i,j)\sum_{i=1}^{n}\sum_{j=1}^{n} gcd(i,j)，n<=105n<=105n∑i=1n∑j=1ngcd(i,j)∑i=1n∑j=1ngcd(i,j)\sum_{i=1}^{n}\sum_{j=1}^{n} gcd(i,j)=∑d=1n∑i=1n∑j=1n[g...
2018-08-23 14:09:11
226
原创 「学习笔记」爬山算法与模拟退火
爬山算法（Hill Climbing ）爬山算法（Hill Climbing ）是一种局部择优的方法，采用启发式方法，是对深度优先搜索的一种改进，它利用反馈信息帮助生成解的决策。 属于人工智能算法的一种。有一些问题，是找全局最大值的。题目所述出自变量与值之间的关系是函数，而这函数图像往往像山的形状，可能有许多局部最大、局部最小。于是“爬山”、“高地”、“山脊”等词形象地表示了算法的作用与函...
2018-08-22 09:50:45
965
原创 「学习笔记」中国剩余定理及扩展
中国剩余定理总是不太熟练，写一篇BlogBlogBlog复习一下好了。中国剩余定理有nnn个线性同余方程:x≡ai(modmi)x≡ai(modmi)x \equiv a_i \;(mod\; m_i)，mimim_i都互质，求解这个线性同余方程组使用中国剩余定理（ChineseremaindertheoremChineseremaindertheoremChinese\; r...
2018-08-19 23:10:18
298
原创 「学习笔记」莫比乌斯反演
前言（今天入坑莫比乌斯反演，感觉全都是不可做题..）莫比乌斯反演属于数论中较难的部分吧，做这种题一般长片的推导，化简，最后用枚举或者整除分块等求出答案（我做题少只见过这些）.下面要引入一些东西做铺垫.数论函数与积性函数首先定义数论函数f(x)f(x)f(x)：定义域和值域都是整数的函数.定义积性函数：∀a,b,gcd(a,b)=1∀a,b,gcd(a,b)=1\forall...
2018-08-19 19:04:53
483
原创 「HDU 6314」Matrix「容斥」
题目传送门题意对一个n∗mn∗mn*m的二维方格涂上黑白两色，求至少aaa行全黑，bbb列全黑的方案数.题解可以先考虑行，因为行列是独立的.我们用g(k)g(k)g(k)表示至少有kkk行全黑的方案数，容斥求一下：g(k)=∑i=knf(i)∗Cin∗2n−ig(k)=∑i=knf(i)∗Cni∗2n−ig(k) = \sum_{i=k}^{n}f(i)*C_{n}^{...
2018-08-19 16:57:29
355
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人