题目:https://loj.ac/problem/3090 题解:https://www.luogu.org/blog/rqy/solution-p5320 1.用斯特林数把下降幂化为普通的幂次求和 2.找出通项公式,使得幂次变成二项式,进而将 [ l , r ] 的部分变成等比数列求和 3.模 998244353 下没有 \( \sqrt{5} \) ,所以“扩域”,就是把数表示成 \( a+b*\sqrt{5} \) :\( \sqrt{3} \) 也同理 注意扩域之后,不满足费马小定理,…

[BJOI2019]勘破神机(斯特林数,数论) 题面 洛谷 题解 先考虑\(m=2\)的情况. 显然方案数就是\(f_i=f_{i-1}+f_{i-2}\),即斐波那契数,虽然这里求出来是斐波那契的第\(n+1\)项,但是本质上没什么区别,就默认是斐波那契数列了. 斐波那契数列的特征根是\(\alpha=\frac{1+\sqrt 5}{2},\beta=\frac{1-\sqrt 5}{2}\),然后大力设一下通项是\(f_n=A\alpha^n+B\beta^n\),可以解出\(f_n=\f…

[BJOI2019]勘破神机 推式子好题 m=2,斐波那契数列,$f_{n+1}$项 不妨$++l,++r$,直接求$f_n$ 求$\sum C(f_n,k)$,下降幂转化成阶乘幂,这样都是多项式了,方便交换求和号 最后面的斐波那契数列用通项公式求.二项式展开. 交换求和号之后,枚举i,j 最后一项是等比数列求和. %rqy m=3, n为奇数是0 n是偶数时,令n=n/2 递推公式:$g_n=4\times g_{n-1}+g_{n-2}$ 证明:枚举从后往前第一个完全分出的块,除了块长为2的…

题意:f[i],g[i]分别表示用1*2的骨牌铺2*n和3*n网格的方案数,求ΣC(f(i),k)和ΣC(g(i),k),对998244353取模,其中l<=i<=r,1<=l<=r<=1e18 题解:显然打表发现f[i]为斐波那契数列,g[2i+1]=0,g[2i]=4g[2i-2]-g[2i-4]. 然后考虑m=2的斐波那契部分:k是给定的,仅需求斐波那契数列的下降幂,然后可以用第一类斯特林数去转换,然后求斐波那契数列的幂之和,假设斐波那契数列的两个特征根为a,b,则f(…

洛谷题面传送门 神仙题(为什么就没能自己想出来呢/zk/zk) 这是我 AC 的第 \(2\times 10^3\) 道题哦 首先考虑 \(m=2\) 的情况,我们首先可以想到一个非常 trivial 的 DP:\(dp_i\) 表示填好前 \(i\) 列的方案数,那么第 \(i\) 列显然有横着放和竖着放两种可能,方案数分别是 \(dp_{i-2}\) 和 \(dp_{i-1}\),因此我们有 \(dp_i=dp_{i-2}+dp_{i-1}\),边界条件 \(dp_0=1\),显然这个递推式…

真的是好题,只不过强行多合一有点过分了…… 题目大意: $T$ 组数据.每个测试点中 $m$ 相同. 对于每组数据,给定 $l,r,k$,请求出 $\dfrac{1}{r-l+1}\sum\limits_{n=l}^r\dbinom{f(n,m)}{k}\bmod 998244353$. 其中 $f(n,m)$ 表示用 $1\times 2$ 的骨牌(可以变成 $2\times 1$)填满 $n\times m$ 的网格的方案数. $1\le T\le 5,1\le l\le r\le 10^{…

传送门 首先我们要知道要求什么.显然每次放方块要放一大段不能从中间分开的部分.设\(m=2\)方案为\(f\),\(m=3\)方案为\(g\),\(m=2\)可以放一个竖的,或者两个横的,所以\(f_i=f_{i-1}+f_{i-2}\);\(m=3\),因为只有\(i\)为偶数有值,所以后面的\(i\)其实是\(2i\),然后可以发现要么放三个横的,要么像下面这样放长度为偶数的块 |--| ------ |--| |----| ---- |----| 可以注意到每种长度都有两种方案(长度为\(…

简单线性代数练习题 首先翻开具体数学生成函数一章,可以发现\(F(n),G(n)\)满足以下递推式 \[F(n)=F(n-1)+F(n-2),F(0)=1,F(1)=1\] \[G(n)=4G(n-2)-G(n-4),G(2)=3,G(0)=1\] 我们发现\(G\)只有偶数项有值(证明的话可以把网格黑白染色什么的来证明) 那么我们可以设\(T(n)=G(2n)\) 那么\(T\)服从以下递推式 \[T(n)=4T(n-1)-T(n-2),T(1)=1,T(0)=1\] 那么我们发现\(F\)和…

LOJ#3090. 「BJOI2019」勘破神机 为了这题我去学习了一下BM算法.. 很容易发现这2的地方是\(F_{1} = 1,F_{2} = 2\)的斐波那契数列 3的地方是\(G_{1} = 3,G_{2} = 11\)其中下标表示长度的\(\frac{1}{2}\),可以得到\(G_{3} = 4G_{2} - G_{1}\) 然后我们列一波特征根方程,可以得到 \(m = 2\)时 $$ \left{\begin{matrix} x_{1} = \frac{1 + \sqrt{5}}…

https://www.cnblogs.com/cjyyb/p/10747543.html 特征方程+斯特林反演化简式子,要注意在模998244353意义下5没有二次剩余,所以每个数都要用$a+b\sqrt{5}$的形式表示,运算类似复数. 斯特林反演的几个用法: 1.下降幂转幂:连续求和时可以通过等比数列求和公式加速. 2.幂转下降幂:类似自然数幂和地用有限微积分加速,或帮助设计DP状态. #include<cstdio> #include<algorithm> #define…

题目传送门 传送门 题目大意 设$F_{n}$表示用$1\times 2$的骨牌填$2\times n$的网格的方案数,设$G_{n}$$表示用$1\times 2$的骨牌填$3\times n$的网格的方案数. 给定$l, r, k$,求$\frac{1}{r - l + 1}\sum_{i = l}^{r} \binom{F_{i}}{k}$. 给定$l, r, k$,求$\frac{1}{r - l + 1}\sum_{i = l}^{r} \binom{G_{i}}{k}$. 之前好像在…

题目链接:洛谷 我一开始不知道$N,M$有什么用处,懵逼了一会儿,结果才发现是输入数据范围... $$\begin{aligned}\binom{n}{k}Ans&=\sum_{i=0}^k\binom{m}{i}\binom{n-m}{k-i}i^L \\&=\sum_{i=0}^k\binom{m}{i}\binom{n-m}{k-i}\sum_{j=0}^Lj!\binom{i}{j}\begin{Bmatrix}L \\ j\end{Bmatrix} \\&=\sum_{j…

题目:http://codeforces.com/contest/717/problem/A 是 BJOI2019 勘破神机 的弱化版. 令 \( g[i] \) 表示长为 i .以 1 结尾的方案数,有 \( g[i]=g[i-1]+g[i-2] , g[0]=g[1]=1 \) : 令 \( f[i] \) 表示长为 i 的方案数,有 \( f[i]=g[i]+g[i-1] \) 发现 \( f[i]=f[i-1]+f[i-2] , f[0]=1 , f[1]=2 \) 那么令 l+=2 ,…

据国内媒体报道,一些苹果上游供应商已经接到通知,iPhone6系列将会在5月底彻底停产,一时间,竟在网络上引发汹涌的怀念之情.iPhone6的特别之处在于它是苹果第一款大屏幕的智能手机,标志着库克彻底推翻乔布斯"3.5英寸黄金比例"的学说,屏幕扩大之后,不但增加了显示元素,更重要的是扩大了电池的容量,虽然掉色问题.屁股一压就弯的问题以及凸出的摄像头等等,都曾遭遇诟病,但"More than bigger"的设计非常受亚洲消费者欢迎,当年销量大涨,仅2014年Q4就卖…

目录 参考资料 前言 暴力 nlog^2n的做法 nlogn的做法 代码 参考资料 百度百科 斯特林数 学习笔记-by zhouzhendong 前言 首先是因为这道题,才去研究了这个玩意:[2019雅礼集训][第一类斯特林数][NTT&多项式]permutation 感觉这个东西非常的...巧妙. 暴力 第一类斯特林树S(n,k)就是将n个数字划分为k个不相区分的圆排列的方案数(即忽略顺序). 首先,第一类斯特林数有一个人尽皆知的\(O(n^2)\)递推式: \[S(n,k)=S(n-1,k-…

目录 题意 输入格式 输出格式 思路 代码 题意 找有多少个长度为n的排列,使得从左往右数,有a个元素比之前的所有数字都大,从右往左数,有b个元素比之后的所有数字都大. n<=2*10^5,a,b<=n 输入格式 输入三个整数n,a,b. 输出格式 输出一个整数,表示答案. 思路 这道题是真的神啊... 首先,根据官方题解的思路,首先有一个n^2的DP: 定义dp[i][j]表示一个长度为i的排列,从前往后数一共有j个数字大于所有排在它前面的数字. 首先有转移式: \[dp[i][j]=dp[…

题目描述 有一个\(n\)个元素的随机置换\(P\),求\(P\)分解出的轮换个数的\(m\)次方的期望\(\times n!\) \(n\leq 100000,m\leq 30\) 题解 解法一 有一种暴力的做法:设\(f_{i,j}\)为\(i\)个元素的随机置换\(P\),分解出的轮换个数的\(j\)次方的期望\(\times i!\) 考虑第\(P_i\)是什么. 如果是\(i\),那么就多了一个轮换,用二项式定理展开得到\(\sum_{k=0}^jf_{i-1,k}\binom{j}{…

BZOJ 洛谷 挺套路但并不难的一道题 \(Description\) 给定一棵\(n\)个点的树和\(K\),边权为\(1\).对于每个点\(x\),求\(S(x)=\sum_{i=1}^ndis(x,i)^K\). \(n\leq50000,\ k\leq150\). \(Solution\) 和其它求\(x^k\)的题一样,依旧用第二类斯特林数展开.(二项式定理依旧可以得到部分分,依旧不想看=-=) \[\begin{aligned}S(x)&=\sum_{i=1}^ndis(x,i)^K…

[题目]G. Partitions [题意]n个数$w_i$,每个非空子集S的价值是$W(S)=|S|\sum_{i\in S}w_i$,一种划分方案的价值是所有非空子集的价值和,求所有划分成k个非空子集的方案的价值和.1<=k<=n<=2*10^5,1<=wi<=10^9. [算法]斯特林数 [题解]首先价值与具体数字没有关系,即: $$ans=num*\sum_{i=1}^{n}w_i$$ 其中num表示1在每个k划分方案中所在集合的大小的和. 考虑一种角度,所在集合的大…

题目链接 洛谷P4609 题解 感性理解一下: 一神带\(n\)坑 所以我们只需将除了\(n\)外的\(n - 1\)个元素分成\(A + B - 2\)个集合,每个集合选出最大的在一端,剩余进行排列,然后选出\(A - 1\)个集合放左边,剩余放右边 容易发现分割集合并内部排列实质对应第一类斯特林数\[\begin{bmatrix} n - 1 \\ A + B - 2 \end{bmatrix}\] 所以答案就是 \[\begin{bmatrix} n - 1 \\ A + B - 2 \e…

数学杂烩总结(多项式/形式幂级数+FWT+特征多项式+生成函数+斯特林数+二次剩余+单位根反演+置换群) 因为不会做目录所以请善用ctrl+F 本来想的是笔记之类的,写着写着就变成了资源整理 一些有的没的的前置 导数 \(f'(x)=\lim\limits_{\triangle x\rightarrow 0}\frac{f(x+\triangle x)-f(x)}{\triangle x}\) \(\sin x:\cos x\) \(\cos x:-\sin x\) \(\ln x:\frac{…

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3625 题意: n个房间,房间里面放着钥匙,允许破门而入k个,拿到房间里面的钥匙后可以打开对应的门,但是1号门不能破门而入,求这样检查完所有房间,概率是多少? 分析: 钥匙随机放到房间,全排列有n!: n个房间,破k个门进入,就是第一类斯特林数S(n,k): 但是,第一个门不能破门而入,就是要减去S(n-1,k-1): 然后求和SUM = S(n,i)  {1<=i<=k} 概率就是 SUM / N…

Description “简单无向图”是指无重边.无自环的无向图(不一定连通). 一个带标号的图的价值定义为每个点度数的k次方的和. 给定n和k,请计算所有n个点的带标号的简单无向图的价值之和. 因为答案很大,请对998244353取模输出.$n \le 10^9,k \le 200000$ 化学学考时含义推式子+手动打表找规律得到了一个$O(nlogn)$的式子开心的很我以为我要AC了回来看数据范围就升天了. 问NC大神这题用到了什么:斯特林数/伯努利数.然后就自闭了学了一天的知识点还去做了点…

题意 ​ 题目链接:https://www.luogu.org/problem/P4827 ​ 给定一棵 \(n\) 个节点的树和一个常数 \(k\) ,对于树上的每一个节点 \(i\) ,求出 \(\displaystyle \sum_{i=1}^n{\rm dist}^k(i,j)\),其中 \(\rm{dist}\) 函数表示树上两点距离. ​ \(1 \leq n \leq 50000\) ​ \(1\leq k \leq 150\) 思路 ​ 看到求答案 \(k\) 次方的问题,应该联…

传送门 解题思路 比较有意思的一道数学题.首先\(n*k^2\)的做法比较好想,就是维护一个\(x^i\)这种东西,然后转移的时候用二项式定理拆开转移.然后有一个比较有意思的结论就是把求\(x^i\)这种东西变成组合数去求,具体来说就是\(n^k=\sum\limits_{i=1}^k\dbinom{n}{i}*S[k][i]*i!\),\(S\)表示第二类斯特林数,第二类斯特林数可以表示为有\(n\)个盒子要装\(m\)个小球,然后在给盒子和求加上编号就可以得出上面的式子.这样的话在根据帕斯卡…

正题 题目链接:https://www.luogu.com.cn/problem/P6620 题目大意 给出\(n,x,p,m\)和一个\(m\)次多项式\(f\)求 \[\sum_{k=0}^nf(k)\times x^k\times \binom{n}{k} \] 答案对\(p\)取模. \(1\leq n\leq 10^9,1\leq m\leq 1000\) 解题思路 什么混凝土数学题 首先我们发现这个组合数\(\binom{n}{k}\)处理的十分难受,有一个下降幂的结论正好可以把这个…

题目 题目里要求的是: \[\sum_{k=0}^n f(k) \times X^k \times \binom nk \] 这里面出现了给定的多项式,还有组合数,这种题目的套路就是先把给定的普通多项式转成下降幂多项式.这一步可以做到\(O(mlogm)\),(模板)但是这题不需要,这个后面再说.假设现在已经得出了f的下降幂多项式的系数\(b_i\),则: \[\begin{align} f(k)&=\sum_{i=0}^m b_ik^{\underline i}\\ ans&=\sum_…

插板法基础知识 斯特林数见百科 #include<iostream> #include<cmath> #include<cstdio> #include<cstring> #define LL long long #define eps 1e-7 #define MOD 1000000007 using namespace std; ][]={},stir2[][]={}; int main(){ ;i<=;i++){ c[i][]=c[i][i]=;…

Machine scheduling Time Limit: 5000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 1000    Accepted Submission(s): 363 Problem Description A Baidu's engineer needs to analyze and process large amount of data o…

Count the Buildings Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)Total Submission(s): 1249    Accepted Submission(s): 408 Problem Description There are N buildings standing in a straight line in the City, numbere…