![](https://img-blog.csdnimg.cn/20201014180756928.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
---------mathematics----------
Iking123
这个作者很懒,什么都没留下…
展开
-
【JZOJ3657】【NOI2014模拟】随机游走(walk)(期望+高斯消元)
ProblemInputOutputHintSolution 乍一看,外加无穷级数的期望裸题,估计分分钟切掉。 细一想,好像有什么不对——我推不出到底怎么列方程。我试了好几个方法,但在手玩了几个数据之后,发现都是错的。 看了看那惜字如金的题解,绞尽脑汁地理解了一番,终于懂了。 首先,我们可以设点(i,j)到点(n,m)的期望距...原创 2018-01-30 19:59:43 · 657 阅读 · 0 评论 -
Pollard_Rho大数分解学习小记
前言该算法用于对大数进行质因子分解。它的核心思想是:对于一个数N,若N不为质数,设N=c*d(c、d均大于1),我们找到一个c,递归分解c和d。生日悖论一个班有k个同学,至少有两名学生的生日是同一天的概率P是多少?正难则反,考虑容斥。P=1-任意两名学生的生日都不同的概率。k=1,P=0k=1,P=0k=1,P=0k=2,P=1−1∗364365k=2,P=1−1∗3...原创 2018-08-24 16:50:02 · 769 阅读 · 1 评论 -
【JZOJ4732】【NOIP2016提高A组模拟8.23】函数(线筛求欧拉函数+Pollard_Rho大数分解)
Problem定义一个函数f(x),定义域和值域均为N∗N∗N^*,且满足∑d|nf(d)=n∑d|nf(d)=n\sum_{d|n} f(d)=n。现给定N个幸运数字A1,A2,...,AnA1,A2,...,AnA_1,A_2,...,A_n,求∑ni=1f(Ai)∑i=1nf(Ai)\sum_{i=1}^n f(A_i)。HintSolution解谜f(x)...原创 2018-08-25 16:56:28 · 269 阅读 · 0 评论 -
欧拉函数详解
前言关于这东西,我们都是几百年前就接触了。虽然说这东西现在看来很简单了,但是有些公式的证明也不简单,在此做个详解。定义在数论,对正整数n,欧拉函数(Euler’s totient function)是小于或等于n的正整数中与n互质的数的数目。严谨地,φ(n)=∑ni=1δgcd(i,n)1φ(n)=∑i=1nδgcd(i,n)1\varphi(n)=\sum_{i=1}^...原创 2018-08-25 17:33:00 · 1150 阅读 · 1 评论 -
【JZOJ5870】【NOIP2018模拟9.15】地图 (递推+DP+组合数学)
ProblemHintSolution首先,图中只会存在链和环。记图中有one个度数为1的点,two个度数为2的点。囿于每条链有两个度数为1的点(链的两端),链的数量是确定的:one2one2\frac{one}2。这时,我灵(nao)光(zi)一(wa)闪(te),想到了一个优(sha)美(bi)的方法。我的SB方法:递推+组合数学+容斥观察到链和环...原创 2018-09-15 17:25:15 · 366 阅读 · 0 评论 -
【JZOJ5894】【NOIP2018模拟10.5】同余方程(容斥+计数问题+类数位DP)
ProblemSolution首先,考虑类似差分的容斥。设fa,bf_{a,b}fa,b表示x在[0,a-1]范围内,y在[0,b-1]范围内的答案。则原Ans=fr1+1,r2+1−fl1,r2+1−fr1+1,l2+fl1,l2Ans=f_{r_1+1,r_2+1}-f_{l_1,r_2+1}-f_{r_1+1,l_2}+f_{l_1,l_2}Ans=fr1+1,r2+1−f...原创 2018-10-05 20:58:06 · 240 阅读 · 0 评论 -
【JZOJ5803】【2018.8.12省选模拟】girls(容斥+乱搞+三元环计数)
Description给定n个点,m对冲突(如点x与点y有冲突)。要求选不冲突的三个点i、j、k。三元组(i,j,k)(i,j,k)(i,j,k)贡献为A∗i+B∗j+C∗kA*i+B*j+C*kA∗i+B∗j+C∗k。求总贡献。n,m≤2∗105n,m≤2*10^5n,m≤2∗105。Solution正难则反,考虑容斥。记fif_ifi表示至少存在i对冲突的总贡献,则Ans...原创 2018-11-03 17:51:36 · 295 阅读 · 0 评论 -
【JZOJ2700】【GDKOI2012模拟02.01】数字
DescriptionT组询问。Data ConstraintSolution先说一个结论:D(n)=(n−1) mod 9+1D(n)=(n-1)\ mod\ 9+1D(n)=(n−1) mod 9+1。证明如下:首先,快速判断一个大数n是否为9的倍数有一个黑科技:那就是将其每个位数相加,判断得到的和是否为9的倍数。(即判断S(n)是否...原创 2018-12-26 21:03:25 · 199 阅读 · 0 评论 -
斯特林数小结(update gradually...)
概念斯特林数(Stirling number)是由18世纪数学家James Stirling提出的两类数。第一类:轮换数,将n个数排成k个非空的环(圆排列,有序)的方案数,用s(n,m)s(n,m)s(n,m)或[nk]\left[ \begin{matrix}n\\k\end{matrix} \right][nk]表示。还分为无符号第一类斯特林数su(n,m)s_u(n,m)su(n,...原创 2019-01-28 08:07:27 · 312 阅读 · 0 评论 -
【JZOJ6242】【NOI2019模拟2019.7.1】刀剑传奇(容斥+高维前/后缀和)
Description给定N(≤5×105)N(≤5\times10^5)N(≤5×105)个长度为20的01串,每次可选择其中恰好KKK个。Q(≤5×105)Q(≤5\times10^5)Q(≤5×105)个询问,每个询问表示为一个长度为20的01串,若其某一位为1则选择的KKK个串中至少要有一个的该位为1。求方案数。Solution考虑容斥。设f[i]f[i]f[i]表示选择了K...原创 2019-07-02 21:51:52 · 745 阅读 · 0 评论 -
Miller_Rabin素性测试学习小记
问题给出一个正整数n,判断它是不是质数。有一个简单暴力的方法:试除法,从2枚举到n−−√n\sqrt n,如果有一个数能整除n,则表明n为合数。但这么做的时间复杂度是O(n−−√)O(n)O(\sqrt n)的,当n的规模达到1016101610^{16}甚至更高时,基本上就会跑1s以上了。有什么更快的方法?费马小定理通过费马小定理我们知道,对于一个质数p,有:ap−1≡1...原创 2018-08-24 14:40:50 · 738 阅读 · 0 评论 -
【JZOJ4424】【HNOI2016模拟4.2】道路(DP+组合数学+二项式展开)
Problemn(≤50)个点,初始时没有边。要建造若干边,使得其联通。对于一个方案,若其有k条边,对答案贡献为k2k2k^2。Solution首先,假设对答案贡献为k0k0k^0怎么做(即求方案数)。设f(n)表示n个点,要求联通的方案数;设g(n)表示不要求联通的方案数。显然有:g(n)=2C2ng(n)=2Cn2g(n)=2^{C_n^2}。考虑容斥。f(n)=...原创 2018-08-20 22:37:28 · 179 阅读 · 0 评论 -
快速沃尔什变换——FWT学习小记
用途FWT主要解决这样的问题,给出A(x),B(x),求Ck=∑i⊗j=kAjBkCk=∑i⊗j=kAjBkC_k=\sum_{i\otimes j=k}A_jB_k,其中⊗⊗\otimes表示and(按位与)/or(按位或)/xor(按位异或)运算中的一种。定义首先,A(x)可以表示成一个n维向量(a0,a1,...,an−1)(a0,a1,...,an−1)(a_0,a_...原创 2018-08-17 19:50:27 · 612 阅读 · 0 评论 -
【JZOJ4021】【CF251D】Two Sets(异或方程组+高斯消元)
Problem 给出n(≤100000)个数字,要求将其分成两个集合,集合1的异或和为x1,集合2的异或和为x2,求出一种划分方案,使得在x1+x2最大的前提下x1尽量小。Solution 这道题我刚看感觉很不可做,于是打了个暴力。 正解竟然是神奇的高斯消元。 首先,我们瞄到“异或”这个词,马上联想到位运算,继而联想到从高到低逐位来搞。当然我在比赛时也想到了这一点。...原创 2018-03-24 17:29:56 · 258 阅读 · 0 评论 -
【JZOJ3661】【SHTSC2014】概率充电器(charger) (期望+树形DP)
Problem 给出一颗n(≤500000)个点的树,点i开始时有qi%的概率被直接充电,每条边有p%的概率可以导电,随后电能可以从直接充电的点经过通电的边使得其他点进行间接充电。求能充电的点数的期望,四舍五入到小数点后6位小数。Solution 这道期望题我乍一看感觉很简单,但是仔细想想,又不会了,所以一开始打了暴力。 正解是树形DP。 但这并不是普通的树形DP,...原创 2018-01-31 16:49:20 · 286 阅读 · 0 评论 -
【JZOJ5622】【NOI2018模拟4.2】table(组合数学)
ProblemInputOutputHintSolution40points 前8个点的n和m较小,可以直接暴力求出数表中的每个数。 时间复杂度:O(nm)O(nm)O(nm)。Code Code就不贴了,毕竟后面的60points方法也有这部分。60points 我们考虑第9~12个点,即p=1的情况...原创 2018-04-06 21:26:08 · 196 阅读 · 0 评论 -
莫比乌斯(Mobius)反演知识整合
前言 虽然我去年就已经学过这东西了,但是一直没有整理、归纳。而且我发现我对它的套路不够熟悉、对它的理解不够透彻,卒致GDOIDay2T1没切(详细事故记录戳这里)。所以,在此做一个知识整合。引子 先来看bzoj的一道题:Problem b。Problem 给出n组询问,每组询问给出a、b、c、d、k五个整数,求有多少个数对(x,y),满足a≤x≤b,c≤y≤d,且...原创 2018-05-01 16:21:57 · 942 阅读 · 0 评论 -
【JZOJ1161】【NOI2002第二试第三题】机器人M号(欧拉函数+递推)
Problemhttps://jzoj.net/senior/#main/show/1161Solution首先要喷一下OJ上的样例输出:明明要输出三行,样例却只输出了一行,害得我也跟着输出了一行,于是被坑掉100points!然后,分析题意。显然,对于x,若y(y>1且y≠x)为其约数,则y为其老师。而x的独立数则显然为欧拉函数ϕ(x)ϕ(x)\phi(x)...原创 2018-07-09 15:55:17 · 251 阅读 · 0 评论 -
【JZOJ4196】二分图计数(容斥+dfs)
ProblemHintSolution 看到这道题时我也想了一会容斥,只不过当时比较傻逼,想到的是2n2n2^n暴枚一个S,再2n2n2^n暴枚一个S1⊆S,然后就误以为会gg。。。竟然没有想到直接O(3n)O(3n)O(3^n)枚举。。。 下面是TJ: 这个容斥应该很好理解。当S1为空集时,S-S1=S,则S中所有点都可以任意匹配;这样可能会使某个点...原创 2018-07-06 17:23:26 · 444 阅读 · 1 评论 -
【JZOJ5796】【2018提高组】模拟A组&省选 划分(容斥+扩展中国剩余定理+扩展欧几里得)
Problem 有一个未知的序列x,长度为n。它的K-划分序列y指的是每连续K个数的和得到划分序列,y[1]=x[1]+x[2]+….+x[K],y[2]=x[K+1]+x[K+2]+….+x[K+K]….。若n不被K整除,则y[n/K+1]可以由少于K个数加起来。比如n=13,K=5,则y[1]=x[1]+…+x[5],y[2]=x[6]+….+x[10],y[3]=x[11]+x[12]...原创 2018-08-10 19:25:41 · 228 阅读 · 0 评论 -
【JZOJ5813】【NOIP提高A组模拟2018.8.14】 计算(质因数分解+DP+思维)
ProblemHintSolution这道题是妥妥的送了45points。因为100以内的数的约数个数均≤12,我们找出n的约数后,暴力dfs填数即可。时间复杂度O(σ(n)2m)O(σ(n)2m)O(\sigma(n)^{2m})。不过,满分做法还是需要一点思维的。假设我们现在的x数列满足条件I(∀i∈[1,2m],xi∈Z+,xi|n∀i∈[1,2m...原创 2018-08-14 15:59:04 · 259 阅读 · 0 评论 -
Codeforces 662C Binary Table(快速沃尔什变换)
Problem给定一个n(≤20)*m(≤100 000)的01矩阵,每次操作可以将一行或一列取反。求最终1的最少个数。Solution前置技能:快速沃尔什变换(FWT)。观察到n较小,考虑O(2n)O(2n)O(2^n)枚举每一行选或不选。不妨设f(x)表示行的操作状态为x时(我们可用一个二进制数表示状态),经过各种列操作后所得到的最少的1的个数。可以O(m)O(m...原创 2018-08-17 22:39:13 · 265 阅读 · 0 评论 -
【JZOJ6233】【NOI2019模拟2019.6.27】心的旋律(思维题)
Description对于一个两边各有nnn个点的二分图的左边一个点集AAA,定义F(A)F(A)F(A)为右边至少和AAA中的一个点相邻的点的集合。给定一个KKK,要求构造一个二分图,使得对于所有点集AAA,∣F(A)∣<∣A∣|F(A)|<|A|∣F(A)∣<∣A∣有恰好KKK次成立。Data Constraintn≤32,K∈[0,2n]n≤3...原创 2019-06-27 20:38:05 · 373 阅读 · 0 评论