数学
文章平均质量分 84
FSYo
我不去想是否能够成功,
既然选择了远方,
便只顾风雨兼程
展开
-
牛顿插值多项式
感觉是很好玩的东西,而且好像比拉格朗日插值有用,可以用来动态更新插值多项式(维护每个区间的均差。但我只理解了均差的定义和它的巧妙,还没有完全 get 背后的原理 …然后又可以继续往前和,最后发现剩下的就是。然后再拆,发现可以变成。原创 2022-11-21 13:53:27 · 788 阅读 · 0 评论 -
Van der Waerden 定理
https://zhuanlan.zhihu.com/p/150790806https://www.bilibili.com/video/av86801472/题外话:我发现 Van der Waerden 很帅 /se定理内容大概为,将自然数划分成 mmm 个集合,存在任意长度的等差数列。这里给出小学二年级 都能懂的证明(我乱想的,有错请指出)(参考了上面知乎的回答和 B 站的那个视频,不过我感觉知乎的太深奥,B 站的讲得不是很清楚)我们设 F(m,k)F(m,k)F(m,k) 表示若将 [1原创 2021-11-05 00:14:09 · 926 阅读 · 2 评论 -
正态分布?
全是个人理解正态分布是一种广泛出现的连续概率分布,比如升高,分数二项分布是离散情况下的概率分布比如仍硬币,正面的可能性是 ppp,那么仍 nnn 次,xxx 次正面的概率为 (nx)px(1−p)n−x\binom nxp^x(1-p)^{n-x}(xn)px(1−p)n−x容易得到均值 μ=np\mu=npμ=np,方差 σ2=np(1−p)\sigma^2=np(1-p)σ2=np(1−p)并且画柱状图画出来就是钟形,而且和正态分布的概率密度函数特别像对于均值为 μ\muμ 方差为 σ2\原创 2021-09-05 23:28:19 · 490 阅读 · 0 评论 -
杨表与钩子公式
搬运 2019 袁方舟论文我们先来简单认识一下杨表:一个标准杨表是说,每一行严格上升,每一列严格上升设杨表 λ=(λ1,λ2,…,λm)\lambda=(\lambda_1,\lambda_2,\dots,\lambda_m)λ=(λ1,λ2,…,λm) 表示有 mmm 行,每行有 λi\lambda_iλi 个数,并且 λ1≥λ2≥⋯≥λm\lambda_1\ge \lambda_2\ge\dots\ge \lambda_mλ1≥λ2≥⋯≥λm杨表的边角说的是 {(i,j)}\{(原创 2021-03-05 20:59:32 · 825 阅读 · 0 评论 -
「ZJOI2018」树
设大小为 nnn 的不同的树为 Tn,1,…,Tn,mT_{n,1},\dots,T_{n,m}Tn,1,…,Tn,m其中每个 TTT 表示这类树的一个集合那么求的就是 1(n−1)!k∑i∣Tn,i∣k\frac{1}{(n-1)!^k}\sum_{i}|T_{n,i}|^k(n−1)!k1∑i∣Tn,i∣k参考:传送门设 置换 fff 下的不动点为 fix(f)fix(f)fix(f)所有 fff 使得状态 XXX 在该置换下为不动点的集合为 SG(X)SG(X)SG(X)我们设原创 2021-02-02 21:44:56 · 1118 阅读 · 0 评论 -
数学题
第三波数学题独孤玖 orz设 f(n)=lcm(1,2,…,n)f(n)=\text{lcm}(1,2,\dots,n)f(n)=lcm(1,2,…,n)证明:f(2n)∣f(n)×(2nn)f(2n)\mid f(n)\times \binom{2n}{n}f(2n)∣f(n)×(n2n)我们设 f(n)=n!p,f(2n)=2nn‾qf(n)=\frac{n!}{p},f(2n)=\frac{2n^{\underline n}}{q}f(n)=pn!,f(2n)=q2nn我们用 2nn原创 2020-11-08 16:00:51 · 270 阅读 · 0 评论 -
数树
即统计交集的连通块个数,并有 ycy^cyc 的贡献对于问题 111:考虑点 n−cn-cn−c 条边作为交集,即硬点出 ccc 个联通块,设每个联通块大小为 aia_iai,那么方案数就是 ∏ai×nc−2\prod a_i\times n^{c-2}∏ai×nc−2,求出 fc=∏ai×nn−c−2f_c=\prod a_i \times n^{n-c-2}fc=∏ai×nn−c−2 表示点 ccc 条边断掉的方案数, 然后二项式反演(注意我们直接点联通块的话不好容斥)Ans=∑cyn∑原创 2020-10-20 19:47:31 · 729 阅读 · 0 评论 -
小 H 爱染色
写出答案的式子∑if(i)[(n−im)2−(n−i−1m)2]\sum_{i}^{}f(i)\Big[\binom{n-i}{m}^2-\binom{n-i-1}{m}^2\Big]i∑f(i)[(mn−i)2−(mn−i−1)2]是关于 nnn 的 3m+13m+13m+1 次多项式,我们需要知道各个点值利用 f(k)=∑if(i)∏j≠ii−jk−j=∑if(i)(m−i)!i!(−1)m−ikm+1‾1k−if(k)=\sum_if(i)\prod_{j\neq i}\frac{i-j}原创 2020-10-20 19:00:09 · 278 阅读 · 0 评论 -
【清华集训2017】生成树计数
∑∑di=n−2(∏aidi+1)(∏(di+1)m)(∑(di+1)m)=∏ai∑∑di=n−2(∏aidi)(∑i(di+1)2m∏j≠i(dj+1)m)(n−2)!∏di!Fi(z)=∑jaij(j+1)mzjj!,Gi(z)=∑jaij(j+1)2mzjj!=∏ai∑i[zn−2]Gi(z)∏j≠iFi(z)=∏ai∏Fi(z)∑iGi(z)Fi(z)=∏aiexp[∑ilnFi(z)]×∑iGi(z)Fi(z)\sum_{\sum d_i=n-2}(\prod a_i^{d_i+1})\Bi原创 2020-10-18 12:31:33 · 195 阅读 · 0 评论 -
PKUWC D1T3
月考考得稀撇 /kk然后这题改了一年Loser\mathcal{Loser}Loser∑ilengthi×vi∑j[gcd(pi,j)=1][gcd(t,j)=1]=∑ilengthi×vi∑j∑l∣j,l∣piμ(l)∑r∣j,r∣tμ(r)=∑ilengthi×vi∑l∣pi∑r∣tμ(l)μ(r)⌊nlcm(l,r)⌋\sum_ilength_i\times v_i\sum_{j}[gcd(p_i,j)=1][gcd(t,j)=1]\\=\sum_i length_i\times v_i\sum原创 2020-10-15 00:08:14 · 178 阅读 · 0 评论 -
F2 重做
CF1349F2从小到大,从后往前放,即写出一个放东西的顺序序列 {pi}\{p_i\}{pi},将其划分成若干段极长的下降序列,每一段从小到大对应一个数字。假设要求出 ttt 的答案,我们可以枚举每个处于第 ttt 个下降序列的位置并统计它的贡献 ∑iwaysi,t(ni)(n−i)!\sum_{i}\text{ways}_{i,t}\binom{n}{i}(n-i)!∑iwaysi,t(in)(n−i)!其中 waysi,j\text{ways}_{i,j}waysi,j 表示长为 ii原创 2020-10-08 21:27:31 · 222 阅读 · 0 评论 -
如何优雅地求和?
给定 m,n,p,f(1),…,f(m)m,n,p,f(1),\dots,f(m)m,n,p,f(1),…,f(m),求出∑i=0nf(i)(ni)pi(1−p)n−i\sum_{i=0}^nf(i)\binom{n}{i}p^i(1-p)^{n-i}i=0∑nf(i)(in)pi(1−p)n−iGF\mathcal{GF}GF 告诉我们一种很优雅的方式:∑ct[zt](pez+1−p)n\sum c_t[z^t](pe^z+1-p)^n∑ct[zt](pez+1−p)n是很优美,但是给点值很原创 2020-10-04 21:56:47 · 1666 阅读 · 0 评论 -
【CTSC2017】游戏(贝叶斯公式)(线段树)
传送门考虑 iii 为 1 的概率更前后第一个确定的游戏有关,即 iii 为 1 的概率为P(xi=1∣xl=a,xr=b)=P(xi=1,xl=a,xr=b)P(xl=a,xr=b)=P(xi=1,xr=b∣xl=a)P(xl=a)P(xr=b∣xl=a)P(xl=a)=P(xi=1,xr=b∣xl=a)P(xr=b∣xl=a)P(x_i=1|x_l=a,x_r=b)=\frac{P(x_i=1,x_l=a,x_r=b)}{P(x_l=a,x_r=b)}\\=\frac{P(x_i=1,x_r=b|原创 2020-06-21 17:29:11 · 267 阅读 · 0 评论 -
【CF643F】Bears and Juice
传送门成都七中外地生考试加强版考虑喝过一桶酒或果汁的熊的集合,集合中的元素是一个二元组 (x,tx)(x,t_x)(x,tx) 表示熊的编号以及第几天和的,注意到如果一只熊和了这桶果汁没倒那么不会再喝,如果倒了那么也不会再和,那么每桶果汁或酒的集合就形如 {(x,tx)}\{(x,t_x)\}{(x,tx)} 其中 xxx 只出现一次,若两个集合不同则可以区分它们,床的限制是对集合大小的限制,那么容易发现答案即为anst=∑i=0min(n−1,p)(ni)tians_t=\sum_{i=0}原创 2020-06-10 11:04:18 · 245 阅读 · 0 评论 -
二次剩余
二次剩余:对于整数 q,nq,nq,n,若存在 xxx 使得x2≡qmod nx^2\equiv q\mod nx2≡qmodn则称 qqq 是 nnn 的二次剩余,否则 qqq 是 nnn 的非二次剩余对于奇质数 ppp,引入LegendreLegendreLegendre 符号定义:(ap)={0p∣a−1p∤a,a不是模p意义下的二次剩余1p∤a,a是模p意义下的二次剩余 \left(\frac{a}{p}\right)=\left\{\begin{aligned}& 0原创 2020-06-10 09:55:48 · 376 阅读 · 0 评论 -
斐波那契循环节
Source: here斐波那契数列:fi=fi−1+fi−2,i≥2f0=0,f1=1fi=ϕi−φi5f_i=f_{i-1}+f_{i-2}, i\ge 2\\f_0=0,f_1=1\\f_i=\frac{\phi^i-\varphi^i}{\sqrt 5}fi=fi−1+fi−2,i≥2f0=0,f1=1fi=5ϕi−φi其中 ϕ,φ\phi,\varphiϕ,φ 为两个特征根,ϕ=5+12,φ=5−12\phi =\frac{\sqrt 5+1}{2},\varphi = \f原创 2020-06-10 08:27:19 · 729 阅读 · 0 评论 -
「2019 集训队互测 Day 1」整点计数(min25)
传送门以下整点皆在第一象限讨论令 f(x)=∑i≥1∑j≥1[i2+j2=x2]f(x)=\sum_{i\ge 1}\sum_{j\ge 1}[i^2+j^2=x^2]f(x)=∑i≥1∑j≥1[i2+j2=x2],要求的就是 ∑i=1Nf(i)\sum_{i=1}^Nf(i)∑i=1Nf(i)考虑把整点 (x,y)(x,y)(x,y) 用复平面的高斯整数 x+yix+yix+yi 表示,而 (x+yi)(x−yi)=x2+y2(x+yi)(x-yi)=x^2+y^2(x+yi)(x−yi原创 2020-05-26 16:38:13 · 500 阅读 · 0 评论 -
P5591 小猪佩奇学数学 (单位根反演)
传送门思路:题目提示得很明显是单位根反演,我们将 ⌊ik⌋\lfloor\frac{i}{k}\rfloor⌊ki⌋ 拆成 i−i%kk\frac{i-i\%k}{k}ki−i%k,前半部分用二项式定理,用单位根反演处理后半部分,即∑t=1k−1t∑i=0n(ni)pi∑j=0k−1ωkj(i−t)k=1k∑tt∑j=0k−1ωk−jt(pωkj+1)n=1k∑j(pωkj+1)n∑tt...原创 2020-04-12 12:08:51 · 240 阅读 · 0 评论 -
【省选模拟】20/04/07 考试 (组合意义)(生成函数)(多项式)
有 nnn 个独立的变量(实数)xi∈[0,ai]x_i\in[0,a_i]xi∈[0,ai],求 (∑xi)m(\sum x_i)^m(∑xi)m 的期望直接考虑组合意义:就是 nnn 个里面任选 mmm 个乘起来,先不考虑顺序,那么我们可以枚举每个 xix_ixi 被选的次数,基于此构造 EGFEGFEGFFi(x)=∑j≥0∫0aitjdtj!xt=∑j≥0aij+1(j...原创 2020-04-07 18:48:17 · 295 阅读 · 0 评论 -
【省选模拟】20/04/07 亲 (推式子)(生成函数)(多项式)(拉格朗日插值)
题意:求 ∑i=1n(ni)Qi(∑j=1ijk)\sum_{i=1}^n\binom{n}{i}Q^{i}(\sum_{j=1}^ij^k)∑i=1n(in)Qi(∑j=1ijk),n≤1e9,k≤3e3,Q≤1e9n\le 1e9,k\le3e3,Q\le 1e9n≤1e9,k≤3e3,Q≤1e9发现 3 个大神都有不同的神仙方法,只有我什么都不会无敌的 ldxldxldx...原创 2020-04-07 17:59:58 · 312 阅读 · 0 评论 -
数学题
第二波数学题,独孤玖 orz对于给定的 a>1a>1a>1,证明:∀n∈N∗\forall n\in N^{*}∀n∈N∗,总存在整系数 P(x)P(x)P(x) 使得 P(0),P(1),…,P(n)P(0),P(1),\dots,P(n)P(0),P(1),…,P(n) 互不相同且均形如 ak(k∈Z,ak∈N∗)a^k(k\in Z,a^{k}\in N^{*})ak(...原创 2020-04-05 19:00:00 · 290 阅读 · 0 评论 -
Tutte 矩阵解决一般图图的匹配问题
详见 《2017国家集训队论文》“基于线性代数的一般图匹配”前言:这个东西比带花树好记多了符号规定:Ai,jA_{i,j}Ai,j 表示第 iii 行第 jjj 列,AI,JA_{I,J}AI,J 表示暴力集合 III 中的行和集合 JJJ 中的列剩下的矩阵,AI,JA^{I,J}AI,J 表示除去 III 中的行和 JJJ 中的列的矩阵,rank Arank\ Arank...原创 2020-04-02 16:46:20 · 2020 阅读 · 0 评论 -
【NOI2012】骑行川藏(拉格朗日乘数法)
传送门limits:φ(v1,v2,…,vn)=∑ki(vi−vi′)2si=Emaximize{f(v1,v2,…,vn)=∑sivi}L(v1,v2,…,vn)=f+λφlimits:\varphi(v_1,v_2,\dots,v_n)=\sum k_i(v_i-v_i')^2s_i=E\\ maximize\{f(v_1,v_2,\dots,v_n)=\sum \frac{s_i}{v_...原创 2020-03-26 13:55:39 · 285 阅读 · 3 评论 -
拉格朗日乘数法
求 z=f(x,y)z=f(x,y)z=f(x,y) 在 φ(x,y)=0\varphi(x,y)=0φ(x,y)=0 条件下的最值我们知道,若在 (x0,y0)(x_0,y_0)(x0,y0) 取得最值,在 (x0,y0)(x_0,y_0)(x0,y0) 的某个邻域内 φ(x,y)=0\varphi(x,y)=0φ(x,y)=0 可以确定一个连续且具有连续导数的函数 y=ψ(x)y...原创 2020-03-26 13:34:30 · 316 阅读 · 0 评论 -
数学题
对于 ttt 个集合 A1,A2,…,At∈ZA_1,A_2,\dots,A_t\in ZA1,A2,…,At∈Z,若存在 b1,b2,…bk∈Zb_1,b_2,\dots b_k\in Zb1,b2,…bk∈Z使得 ∀i≠j,(bi+A)∩(bj+A)=∅\forall i\neq j,(b_i+A)\cap(b_j+A)=\empty∀i=j,(bi+A)∩(bj+A)=...原创 2020-03-22 19:57:22 · 270 阅读 · 0 评论 -
常见函数积分
∫kdx=kx+C\int k\text{d}x=kx+C∫kdx=kx+C∫xμdx=xμ+1μ+1+C(μ≠1)\int x^{\mu}\text{d}x=\frac{x^{\mu+1}}{\mu+1}+C(\mu\neq 1)∫xμdx=μ+1xμ+1+C(μ=1)∫dxx=ln∣x∣+C\int \frac{\text{d}x}{x}=\text{ln}|x|+C∫xdx=l...原创 2020-03-12 10:52:21 · 8708 阅读 · 0 评论 -
「LibreOJ NOI Round #1」动态几何问题(整除分块)
传送门转换一下题意就是 ∑i≤n,j≤mij∈{x∣x=k2}\sum_{i\le n,j\le m}ij\in \{x|x=k^2\}∑i≤n,j≤mij∈{x∣x=k2},除去 i,ji,ji,j 中的平方因子,那么剩下的两个数相等,枚举这个数(不含平方因子)可以得到ans=∑i=1nμ(i)2⌊ni⌋⌊mi⌋ans=\sum_{i=1}^{n}\mu(i)^2\sqrt {\lflo...原创 2020-03-10 22:51:30 · 225 阅读 · 0 评论 -
Taylor Formula
Taylor Formula:Taylor\ Formula:Taylor Formula:目的:我们试图用一个多项式去拟合一个复杂的函数,这样在处理问题上提供极大的便利提出问题:设 f(x)f(x)f(x) 在 x0x_0x0 具有 nnn 阶导数,试图找出一个关于 x−x0x-x_0x−x0 的 nnn 次多项式g(x)=a0+a1(x−x0)+a2(x−...原创 2020-03-05 23:44:25 · 280 阅读 · 0 评论 -
L'hospital Rule
L′hospital Rule:L'hospital\ Rule:L′hospital Rule:若 x→a,x→∞x\to a,x\to \inftyx→a,x→∞ 时两个函数 f(x),F(x)f(x),F(x)f(x),F(x) 都趋近于 000 或无穷大,那么极限 limx→af(x)F(x)\lim_{x\to a}\frac{f(x)}{F(x)}limx→a...原创 2020-03-05 22:12:41 · 962 阅读 · 2 评论 -
常见函数求导及求导法则
函数连续:若 f(x)f(x)f(x) 满足, limx→ 0[f(x0+Δx)−f(x0)]=0\lim_{x\to \ 0}[f(x_0+\Delta x)-f(x_0)]=0x→ 0lim[f(x0+Δx)−f(x0)]=0则成 f(x)f(x)f(x) 在 x0x_0x0 连续定义:设函数 f(x)f(x)f(x) 在点 x0x_0x0 的...原创 2020-03-02 21:15:35 · 5737 阅读 · 2 评论 -
线性递推与 Cayley - Hamilton 定理的简要证明
引言:本篇对 OIOIOI 领域 - “线性递推” 中用的重要定理 Cayley−HamiltonCayley -HamiltonCayley−Hamilton 定理做出简要证明符号规定,∣A∣|A|∣A∣ 表示矩阵 AAA 的行列式,AijA_{ij}Aij 表示 AAA (i,j)(i,j)(i,j) 元的代数余子式,A∗A^{*}A∗ 表示 AAA 的伴随矩阵, 根据定义我们知道...原创 2020-02-20 18:03:29 · 839 阅读 · 0 评论 -
【模板】矩阵求逆(矩阵初等变换)
O(n5)O(n^5)O(n5) 做法:先求出 AAA 的伴随矩阵 A∗A^{*}A∗,后利用 A∗A∗=∣A∣∗E⇒A−1=A∗∣A∣A*A^{*}=|A|*E\Rightarrow A^{-1}=\frac{A^{*}}{|A|}A∗A∗=∣A∣∗E⇒A−1=∣A∣A∗ 求解需要求 O(n2)O(n^2)O(n2) 次行列式O(n4)O(n^4)O(n4) 做法:对每一行来一波...原创 2020-02-19 13:20:28 · 1638 阅读 · 0 评论 -
线性代数初探(行列式,矩阵初等变换,矩阵的秩)
引言:边听网课边看线性代数,爽!线代太好玩了,鉴于博主的老年记忆,赶紧记录下来本文主要介绍行列式的一些性质与应用,还有矩阵的一些运算大概是《线性代数》的精简版外加一些自己的理解行列式的定义:令 ppp 为 1,2,...,n1,2,...,n1,2,...,n 的一个排列,排列中的逆序对个数为 ttt,那么行列式为∑(−1)tai,p1a2,p2...an,pn\sum(-1)...原创 2020-02-18 21:54:45 · 472 阅读 · 0 评论