数学
SC.ldxcaicai
我很菜=_=
展开
-
bzoj4487: [Jsoi2015]染色问题(容斥原理)
传送门题意简述:用ccc中颜色给一个n∗mn*mn∗m的方格染色,每个格子可涂可不涂,问最后每行每列都涂过色且ccc中颜色都出现过的方案数。思路:令fi,j,kf_{i,j,k}fi,j,k表示至少有iii行没涂色,至少有jjj列没涂色,至少有ccc种颜色没涂色的方案数。于是fi,j,k=CniCmjCck(c−k+1)(n−i)(m−j)f_{i,j,k}=C_n^iC_m^jC_...原创 2019-02-09 11:37:53 · 405 阅读 · 0 评论 -
常系数齐次线性递推算法学习
文章目录简介求法矩阵快速幂特征多项式一些定义Cayley-Hamilton定理递推优化简介定义:设有数列{an}满足递推关系an=∑i=1kan−ifi\{a_n\}满足递推关系a_n=\sum\limits_{i=1}^{k}a_{n-i}f_i{an}满足递推关系an=i=1∑kan−ifi,则称该数列满足kkk阶齐次线性递推关系。求法现在我们从最基础的矩阵快速幂开始一步一步...原创 2019-02-09 16:32:02 · 1324 阅读 · 0 评论 -
bzoj3944: Sum(杜教筛)
传送门题意:思路:直接上杜教筛。知道怎么推导就很简单了,注意预处理的范围。然后我因为预处理范围不对被zxyoi教育了(ldx你这个傻×两倍常数活该被卡TLE) 喜闻乐见代码:#include<bits/stdc++.h>#define ri register intusing namespace std;const int N=7500005,lim=750000...原创 2019-02-12 21:25:34 · 233 阅读 · 0 评论 -
bzoj2655: calc(生成函数+拉格朗日插值)
传送门题意简述:问有多少数列满足如下条件:所有数在[1,A][1,A][1,A]之间。没有相同的数数列长度为nnn一个数列的贡献是所有数之积,问所有满足条件的数列的贡献之和。A≤1e9,n≤500A\le1e9,n\le500A≤1e9,n≤500思路:肯定不能枚举所有情况。我们先规定这个数列满足a1&lt;a2&lt;⋅˙⋅⋅&lt;ana_1&a...原创 2019-02-19 14:16:29 · 311 阅读 · 0 评论 -
bzoj4568: [Scoi2016]幸运数字(倍增+线性基)
传送门题意:给你一棵带点权的树,多次询问路径的最大异或和。思路:线性基上树??倍增维护一下就完了。时间复杂度O(nlog3n)O(nlog^3n)O(nlog3n)代码:#include<bits/stdc++.h>#define ri register int#define fi first#define se secondusing namespace st...原创 2019-03-25 21:56:27 · 153 阅读 · 0 评论 -
第一类斯特林数小结
第一类斯特林数s1nms1_n^ms1nm表示将nnn个数放进mmm个圆排列的方案数。有一个显然的递推式:s1nm=s1n−1m−1+(n−1)s1n−1ms1_n^m=s1_{n-1}^{m-1}+(n-1)s1_{n-1}^ms1nm=s1n−1m−1+(n−1)s1n−1m,对应的意义:要么第nnn个单独构成一个新的圆排列,要么放在之前某个数的后面。还有一种组合意义:一共进行...原创 2019-04-02 21:33:18 · 842 阅读 · 0 评论 -
第二类斯特林数小结
第二类斯特林数s2nms2_n^ms2nm表示将nnn个数放进mmm个集合的方案数。有一个显然的递推式:s2nm=s2n−1m−1+m∗s2n−1ms2_n^m=s2_{n-1}^{m-1}+m*s2_{n-1}^ms2nm=s2n−1m−1+m∗s2n−1m,对应的意义:要么第nnn个单独构成一个新的集合,要么放在之前某个集合中。然而第二类斯特林数并没有什么生成函数。但它可以花...原创 2019-04-02 22:02:18 · 904 阅读 · 0 评论 -
bzoj3157&&3516国王奇遇记(dp)
传送门传送门随便推一波式子然后O(m^2)dp即可。令f(i)=∑k=1nkimkf(i)=\sum_{k=1}^nk^im^kf(i)=∑k=1nkimk然后用扰乱法化简:(m−1)f(i)=∑k=1nkimk+1−∑k=1nkimk(m-1)f(i)=\sum_{k=1}^nk^im^{k+1}-\sum_{k=1}^nk^im^k(m−1)f(i)=∑k=1nkimk+1−∑k...原创 2019-06-26 17:34:00 · 190 阅读 · 0 评论 -
SRM 548 div1 KingdomAndCities(补集转化+分类讨论)
传送门m=0m=0m=0的情况就是一个补集转化,用个O(n2k2)O(n^2k^2)O(n2k2)的dpdpdp即可。m=1,2m=1,2m=1,2的情况分类讨论即可。代码:#include<bits/stdc++.h>#define ri register intusing namespace std;const int N=105;typedef long long...原创 2019-07-15 23:53:04 · 150 阅读 · 0 评论 -
圆的反演数学证明
Part0点的反演:本质上是点与点的映射。是一个点(x,y)(x,y)(x,y)根据另一个定点(a,b)(a,b)(a,b)(一般取原点)和一个常数kkk,进行变换得到一个新点(x′,y′)(x',y')(x′,y′),满足(x,y),(x′,y′),(a,b)(x,y),(x',y'),(a,b)(x,y),(x′,y...原创 2019-08-03 14:27:37 · 811 阅读 · 0 评论 -
LOJ#3157. 「NOI2019」机器人(dp+下降幂多项式)
传送门神题OrzOrzOrz首先50pts50pts50pts大家应该都会。。。直接设fl,r,xf_{l,r,x}fl,r,x表示区间[l,r][l,r][l,r]中最大值为xxx的方案数,然后枚举最大值出现位置转移即可。复杂度?由于对于每个区间[l,r][l,r][l,r]只有2/32/32/3个点可以用于转移,所以前两维总状态数大概是O(n)O(n)O(n)级别的。而对于最后一...原创 2019-08-06 00:38:02 · 715 阅读 · 2 评论 -
hdu5628 Clarke and math(狄利克雷卷积快速幂)
传送门一眼题。。。注意不能有行末空格233分析一下显然有g=Ik∗fg=I^k*fg=Ik∗f直接上狄利克雷卷积+快速幂即可。。。复杂度O(nlogn2)O(n\log n^2)O(nlogn2)代码:#include<bits/stdc++.h>#define ri register intusing namespace std;const int rlen=1...原创 2019-08-02 13:27:39 · 382 阅读 · 0 评论 -
bzoj4710: [Jsoi2011]分特产(容斥原理)
传送门题意简述:有nnn个人,mmm种物品,给出每种物品的数量aia_iai,问每个人至少分得一个物品的方案数(n,m,每种物品数≤1000n,m,每种物品数\le1000n,m,每种物品数≤1000)。思路:我们算出fif_ifi表示至少有iii个人没有分到物品的方案数容斥一下即可。于是fi=Cni∏j=1mCn−i−1+ajn−i−1f_i=C_n^i\prod_{j=1}^mC...原创 2019-02-09 11:03:04 · 169 阅读 · 0 评论 -
bzoj2839: 集合计数(容斥原理)
传送门题意简述:对于一个有N个元素的集合在其2^N个子集中取出若干集合(至少一个),使得它们的交集的元素个数为K,求取法的方案数。思路:考虑枚举相交的是哪kkk个,有CnkC_n^kCnk种方案,然后考虑剩下的可选可不选一共有22n−k2^{2^{n-k}}22n−k种选法,但是这样选出来的集合可能有其余的数相交,因此我们容斥掉多余的:ans=Cnk∗∑i=0n−kCni22n−k−ia...原创 2019-02-09 10:36:51 · 274 阅读 · 0 评论 -
2018.10.31 bzoj4737: 组合数问题(lucas定理+容斥原理+数位dp)
传送门这是一道让我重新认识lucaslucaslucas的题。考虑到lucaslucaslucas定理:(nm)≡(n%pm%p)∗(npmp)\binom n m \equiv \binom {n\%p} {m\%p}*\binom{\frac n p}{\frac m p}(mn)≡(m%pn%p)∗(pmpn) (mod(mod(mod p)p)p)所以可以看成(nm)\b...原创 2018-10-31 21:41:54 · 289 阅读 · 0 评论 -
2018.11.06 NOIP训练 最大获利(profit)(01分数规划+最大权闭合子图)
传送门好题啊。∑i&lt;jpi,jK∗(200−K)&gt;X\frac{\sum_{i&lt;j}p_{i,j}}{K*(200-K)}&gt;XK∗(200−K)∑i<jpi,j>X=>∑i&lt;jpi,j−XK(200−K)&gt;0\sum_{i&lt;j}p_{i,j}-XK(200-K)&...原创 2018-11-06 12:02:29 · 228 阅读 · 0 评论 -
2018.11.02 NOIP模拟 优美的序列(数论+单调栈/链表)
传送门考虑如果一个区间满足最小值等于最大公约数那么这个区间是合法的。因此我们对于每一个点维护可以延展到的最左/右端点保证这一段区间的gcdgcdgcd等于这个点的值。这个可以用之前同类的链表或者用单调栈维护。代码...原创 2018-11-02 16:40:58 · 248 阅读 · 0 评论 -
2018.11.02 洛谷P2312 解方程(数论)
传送门直接做肯定会TLETLETLE.于是考验乱搞能力的时候到了。我们随便选几个质数来checkcheckcheck合法解,如果一个数无论怎么checkcheckcheck都是合法的那么就有很大概率是正确答案了。事实证明这个做法是对的。因此对于某一个质数pri[i]pri[i]pri[i]我们把所有系数模一个pri[i]pri[i]pri[i]之后带入1&nbsp;pri[i−1]1~p...原创 2018-11-02 21:02:59 · 162 阅读 · 0 评论 -
poj3734 Blocks(生成函数)
传送门生成函数入门题。按照题意构造函数:对于限定必须是出现偶数次的颜色:1+x22!+x44!+...=ex+e−x21+\frac {x^2}{2!}+\frac {x^4}{4!}+...=\frac{e^x+e^{-x}}21+2!x2+4!x4+...=2ex+e−x对于无限定的颜色:1+x1!+x22!+...=ex1+\frac x{1!}+\frac{x^2}{2!}...原创 2018-12-30 19:44:29 · 211 阅读 · 0 评论 -
NOIP训练 偶数个5(简单数论)
传送门对于出题人zxyoizxyoizxyoi先%\%%为敬题目需要龟速乘差评。题意简述:5e55e55e5组数据,给出n,请你求出所有n位数中有偶数个5的有多少,n≤1e18n\le1e18n≤1e18思路:一眼数位dpdpdp,哎哟这nnn怎么这么大绝望.jpg既然是zxyoizxyoizxyoi大毒瘤的题自然要推一推式子了无奈.jpg考虑对每一位构造生成函数:首位:F(x)=...原创 2018-12-31 16:08:52 · 466 阅读 · 0 评论 -
codeforces 528D. Fuzzy Search(fft)
传送门fftfftfft好题。题意简述:给两个字符串s,ts,ts,t,问ttt在sss中出现了几次,字符串只由A,T,C,GA,T,C,GA,T,C,G构成。两个字符匹配的定义:当si−k,si−k+1,...,si+k−1,si+ks_{i-k},s_{i-k+1},...,s_{i+k-1},s_{i+k}si−k,si−k+1,...,si+k−1,si+k中存在至少一个字...原创 2019-01-26 09:49:51 · 416 阅读 · 0 评论 -
codechef Chef at the Food Fair(线段树+泰勒展开)
传送门题意:现在有nnn个位置,每个位置上有一个值aia_iai.要求支持如下两种操作:区间乘vvv求区间的(1−ai)(1-a_i)(1−ai)之积思路:考虑转换式子:Ans=∏i=lr(1−ai)=e∑i=lrln(1−ai)Ans=\prod_{i=l}^r(1-a_i)=e^{\sum_{i=l}^rln(1-a_i)}Ans=∏i=lr(1−ai)=e∑i=l...原创 2019-02-14 18:13:23 · 313 阅读 · 0 评论 -
类欧几里得算法学习
文章目录几个常用等式类欧几里得介绍求f求g求h几个常用等式a≤⌊bc⌋⇔ac≤ba\le\left\lfloor\frac bc\right\rfloor\Leftrightarrow ac\le ba≤⌊cb⌋⇔ac≤ba&amp;lt;⌈bc⌉⇔ac&amp;lt;ba&amp;lt;\left\lceil\frac bc\right\rceil\Leftrightarrow ac&...原创 2019-02-06 19:01:17 · 3283 阅读 · 4 评论 -
bzoj2987: Earthquake(类欧几里得)
传送门题意简述:求满足ax+by+c≤0ax+by+c\le0ax+by+c≤0的二元组(x,y)(x,y)(x,y)对数。思路:类欧几里得算法模板题。把式子变化一下变成:求满足0≤y≤−ax+cb0\le y\le\frac{-ax+c}b0≤y≤b−ax+c的二元组(x,y)(x,y)(x,y)对数。然后就变成求∑i=0⌊ca⌋⌊−ax+cb⌋+1\sum_{i=0}^{\lef...原创 2019-02-06 20:24:24 · 297 阅读 · 0 评论 -
bzoj2187: fraction(类欧几里得)
传送门题意简述:多组询问,每次给出a,b,c,da,b,c,da,b,c,d,求满足ab&lt;pq&lt;cd\frac ab&lt;\frac pq&lt;\frac cdba<qp<dc的所有二元组(p,q)(p,q)(p,q)中qqq为第一关键字,ppp为第二关键字排出来字典序最小的那一对。思路:设计函数f(a,b,p,q,c,d)...原创 2019-02-06 22:24:57 · 210 阅读 · 0 评论 -
数学小知识点整理(TBC)
文章目录前言素数与同余线性筛部分素数线性递推逆元指数循环节降幂当求逆元时模数与求逆元的数有可能不互质时的处理方法一个神奇的结论拓展欧拉定理杂乱的一些性质/技巧二进制枚举子集异或前缀和O(n)预处理popcount多项式一类组合数学卡特兰数通项斯特林数错排公式二项式反演前言感觉距离scoi2019scoi2019scoi2019的时间不多了博主因为太弱所以现在慌得一批,现在尝试梳理一些小知识点顺...原创 2019-02-08 21:57:59 · 1322 阅读 · 0 评论 -
2018.10.09 NOIP训练 蛋糕(数学)
描述今天是鲍勃的生日,爱丽丝打算做一个蛋糕送给他。这是鲍勃的 n 岁生日,所以爱丽丝的蛋糕必须是正 n 边形。而且,鲍勃很喜 欢数字 m,所以这个蛋糕必须放在一个正 m 边形的盒子里。为了让气氛更加浪漫, 爱丽丝将在蛋糕的中心插上一根蜡烛,显然,蜡烛既在蛋糕的中心,又在盒子的 中心是最好的。换句话说,爱丽丝应该使正 n 边形的蛋糕能被容纳在正 m 边形的盒子里,且 使其中心重合。事实上,爱丽...原创 2018-10-09 21:44:59 · 603 阅读 · 0 评论