【BJOI2019】勘破神机（下降幂转自然幂）（第一类斯特林数）（特征方程）

最新推荐文章于 2023-03-31 18:27:09 发布
最新推荐文章于 2023-03-31 18:27:09 发布
阅读量345 收藏 1
点赞数
分类专栏： 组合数学
版权声明：本文为博主原创文章，遵循 CC 4.0 BY 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/zxyoi_dreamer/article/details/102630734
版权

传送门

题解：

完全自己推出来的第一道数学神题。

首先我们知道宽度为 2 2 2的部分方案数是斐波那契数列。

设 f n f_n fn​表示长度为 n n n的时候方案数，题目要求的实际上是这个东西：

∑ n = l r ( f n k ) \sum_{n=l}^r{f_n\choose k} n=l∑r​(kfn​​)

看上去非常地不可做，考虑表示成下降幂的形式，然后利用第一类斯特林数直接转成自然幂，直接上通项公式，然后利用二项式定理展开，后面就是一个等比数列求和：

∑ n = l r ( f n k ) = ∑ n = l r f n k ‾ k ! = 1 k ! ∑ n = l r s k , i ( − 1 ) k − i ∑ n = l r ( A x 1 n + B x 2 n ) i = 1 k ! ∑ n = l r s k , i ( − 1 ) k − i ∑ j = 0 i ( i j ) A j B i − j ∑ n = l r ( x 1 j x 2 i − j ) n \begin{aligned} &\sum_{n=l}^r{f_n\choose k}\\ =&\sum_{n=l}^r\frac{f_n^{\underline{k}}}{k!}\\ =&\frac{1}{k!}\sum_{n=l}^rs_{k,i}(-1)^{k-i}\sum_{n=l}^r(Ax_1^n+Bx_2^n)^i\\ =&\frac{1}{k!}\sum_{n=l}^rs_{k,i}(-1)^{k-i}\sum_{j=0}^i{i\choose j}A^jB^{i-j}\sum_{n=l}^{r}(x_1^jx_2^{i-j})^n \end{aligned} ===​n=l∑r​(kfn​​)n=l∑r​k!fnk​​​k!1​n=l∑r​sk,i​(−1)k−in=l∑r​(Ax1n​+Bx2n​)ik!1​n=l∑r​sk,i​(−1)k−ij=0∑i​(ji​)AjBi−jn=l∑r​(x1j​x2i−j​)n​

由于 k k k很小，直接做就行了。

所以现在的问题就是求通项公式。

宽度为 2 2 2的情况已经玩烂了，就是斐波那契数列。

考虑 3 3 3的情况，显然长度为奇数的情况方案为 0 0 0，设 g n g_n gn​表示长度为 2 n 2n 2n的时候的方案数。定义一个块不可分割当且仅当我们无法竖着切一刀且不切到任何一张牌。

注意到宽度为 3 3 3的时候，任意偶数长度的不可切割块都存在，长度为 2 2 2的有 3 3 3种构造方式，长度更长的有 2 2 2种构造方式。

我们有初始条件 g 0 = 1 , g 1 = 3 g_0=1,g_1=3 g0​=1,g1​=3，以及递推关系：

g n = 3 ⋅ g n − 1 + 2 ⋅ ∑ i = 0 n − 2 g i g_n=3\cdot g_{n-1}+2\cdot \sum_{i=0}^{n-2}g_i gn​=3⋅gn−1​+2⋅i=0∑n−2​gi​

写出 g n − 1 g_{n-1} gn−1​的表达式，差分得到 g n = 4 g n − 1 − g n − 2 g_n=4g_{n-1}-g_{n-2} gn​=4gn−1​−gn−2​

这个形式非常简洁，直接解特征方程可以得到 g n = 3 + 3 6 ( 2 + 3 ) n + 3 − 3 6 ( 2 − 3 ) n g_n=\frac{3+\sqrt 3}{6}(2+\sqrt 3)^n+\frac{3-\sqrt 3}{6}(2-\sqrt 3)^n gn​=63+3 ​​(2+3 ​)n+63−3 ​​(2−3 ​)n

代码（两个namespace 是直接复制粘贴的）：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

using std::cerr;
using std::cout;

cs int mod=998244353;
inline int add(int a,int b){a+=b-mod;return a+(a>>31&mod);}
inline int dec(int a,int b){a-=b;return a+(a>>31&mod);}
inline int mul(int a,int b){ll r=(ll)a*b;return r>=mod?r%mod:r;}
inline int power(int a,int b,int res=1){
	for(;b;b>>=1,a=mul(a,a))(b&1)&&(res=mul(res,a));
	return res;
}
inline void Inc(int &a,int b){a+=b-mod;a+=a>>31&mod;}
inline void Dec(int &a,int b){a-=b;a+=a>>31&mod;}
inline void Mul(int &a,int b){a=mul(a,b);}
inline void ex_gcd(int a,int b,int &x,int &y){
	if(!b){x=1,y=0;return ;}ex_gcd(b,a%b,y,x);y-=a/b*x;
}
inline int inv(int a){
	int x,y;ex_gcd(a,mod,x,y);
	return dec(x,0);
}

cs int N=5e2+7;
int fac[N];
int C[N][N],s[N][N];
inline void init(){
	fac[0]=fac[1]=1;s[0][0]=1;
	for(int re i=2;i<N;++i)fac[i]=mul(fac[i-1],i);
	for(int re i=0;i<N;++i){C[i][0]=1;
		for(int re j=1;j<=i;++j){
			C[i][j]=add(C[i-1][j],C[i-1][j-1]);
			s[i][j]=dec(s[i-1][j-1],mul(s[i-1][j],i-1));
		}
	}
}
//x_n = A * x_1^n + B * x_2^n
namespace Solve2{
	struct cp{
		int x,y;//x + y\sqrt 5
		cp(){}cp(int _x,int _y):x(_x),y(_y){}
		friend cp operator+(cs cp &a,cs cp &b){return cp(add(a.x,b.x),add(a.y,b.y));}
		friend cp operator-(cs cp &a,cs cp &b){return cp(dec(a.x,b.x),dec(a.y,b.y));}
		friend cp operator*(cs cp &a,cs cp &b){
			return cp(add(mul(a.x,b.x),mul(5,mul(a.y,b.y))),add(mul(a.x,b.y),mul(a.y,b.x)));
		}
		friend cp operator+(cs cp &a,int b){return cp(add(a.x,b),a.y);}
		friend cp operator-(cs cp &a,int b){return cp(dec(a.x,b),a.y);}
		friend cp operator*(cs cp &a,int b){return cp(mul(a.x,b),mul(a.y,b));}
		void operator+=(cs cp &b){Inc(x,b.x),Inc(y,b.y);}
	};
	inline cp inv(cs cp &a){return cp(a.x,dec(0,a.y))*::inv(dec(mul(a.x,a.x),mul(5,mul(a.y,a.y))));}
	inline cp power(cp a,ll b){
		cp res(1,0);
		for(;b;b>>=1,a=a*a)if(b&1)res=res*a;
		return res;
	}
	cs cp A=cp(0,::inv(5)),B=cp(0,mod-::inv(5)),x_1=cp(1,1)*::inv(2),x_2=cp(1,mod-1)*::inv(2);
	cp pw_A[N],pw_B[N],pw_1[N],pw_2[N];
	inline void init_pw(){
		pw_A[0]=pw_B[0]=pw_1[0]=pw_2[0]=cp(1,0);
		for(int re i=1;i<N;++i){
			pw_A[i]=pw_A[i-1]*A;
			pw_B[i]=pw_B[i-1]*B;
			pw_1[i]=pw_1[i-1]*x_1;
			pw_2[i]=pw_2[i-1]*x_2;
		}
	}
	inline cp calc(cp a,ll n){
		if(a.x==1&&a.y==0)return cp((n+1)%mod,0);
		return (power(a,n+1)-1)*inv(a-1);
	}
	inline int solve(ll l,ll r,int k){
		cp res(0,0);
		for(int re i=0;i<=k;++i){
			for(int re j=0;j<=i;++j){
				res+=pw_A[j]*pw_B[i-j]*mul(C[i][j],s[k][i])*
				(calc(pw_1[j]*pw_2[i-j],r)-calc(pw_1[j]*pw_2[i-j],l-1));
			}
		}
		return mul(res.x,::inv(fac[k]));
	}
	inline void main(int T){
		init_pw();while(T--){
			ll l,r;int k;scanf("%lld%lld%d",&l,&r,&k);
			cout<<mul(::inv((r-l+1)%mod),solve(l+1,r+1,k))<<"\n";
		}
	}
}

namespace Solve3{
	struct cp{
		int x,y;//x + y\sqrt 3
		cp(){}cp(int _x,int _y):x(_x),y(_y){}
		friend cp operator+(cs cp &a,cs cp &b){return cp(add(a.x,b.x),add(a.y,b.y));}
		friend cp operator-(cs cp &a,cs cp &b){return cp(dec(a.x,b.x),dec(a.y,b.y));}
		friend cp operator*(cs cp &a,cs cp &b){
			return cp(add(mul(a.x,b.x),mul(3,mul(a.y,b.y))),add(mul(a.x,b.y),mul(a.y,b.x)));
		}
		friend cp operator+(cs cp &a,int b){return cp(add(a.x,b),a.y);}
		friend cp operator-(cs cp &a,int b){return cp(dec(a.x,b),a.y);}
		friend cp operator*(cs cp &a,int b){return cp(mul(a.x,b),mul(a.y,b));}
		void operator+=(cs cp &b){Inc(x,b.x),Inc(y,b.y);}
	};
	inline cp inv(cs cp &a){return cp(a.x,dec(0,a.y))*::inv(dec(mul(a.x,a.x),mul(3,mul(a.y,a.y))));}
	inline cp power(cp a,ll b){
		cp res(1,0);
		for(;b;b>>=1,a=a*a)if(b&1)res=res*a;
		return res;
	}
	cs cp A=cp(3,1)*::inv(6),B=cp(3,mod-1)*::inv(6),x_1=cp(2,1),x_2=cp(2,mod-1);
	cp pw_A[N],pw_B[N],pw_1[N],pw_2[N];
	inline void init_pw(){
		pw_A[0]=pw_B[0]=pw_1[0]=pw_2[0]=cp(1,0);
		for(int re i=1;i<N;++i){
			pw_A[i]=pw_A[i-1]*A;
			pw_B[i]=pw_B[i-1]*B;
			pw_1[i]=pw_1[i-1]*x_1;
			pw_2[i]=pw_2[i-1]*x_2;
		}
	}
	inline cp calc(cp a,ll n){
		if(a.x==1&&a.y==0)return cp((n+1)%mod,0);
		return (power(a,n+1)-1)*inv(a-1);
	}
	inline int solve(ll l,ll r,int k){
		cp res(0,0);
		for(int re i=0;i<=k;++i){
			for(int re j=0;j<=i;++j){
				res+=pw_A[j]*pw_B[i-j]*mul(C[i][j],s[k][i])*
				(calc(pw_1[j]*pw_2[i-j],r)-calc(pw_1[j]*pw_2[i-j],l-1));
			}
		}
		return mul(res.x,::inv(fac[k]));
	}
	inline void main(int T){
		init_pw();while(T--){
			ll l,r;int k;scanf("%lld%lld%d",&l,&r,&k);
			cout<<mul(::inv((r-l+1)%mod),solve(l+1>>1,r>>1,k))<<"\n";
		}
	}
}

signed main(){
#ifdef zxyoi
	freopen("calc.in","r",stdin);
#endif
	int T,m;scanf("%d%d",&T,&m);init();
	if(m==2)Solve2::main(T);else Solve3::main(T);
	return 0;
}
洛谷P5393 【模板】下降幂多项式转普通多项式（多项式快速插值）
苟为蒟蒻又何妨
08-03 525
传送门 已知下降幂多项式A(x)=∑i=0naixi‾A(x)=\sum_{i=0}^na_ix^{\underline i}A(x)=∑i=0n​ai​xi​ 现在求普通多项式B(x)=∑i=0nbixiB(x)=\sum_{i=0}^nb_ix^iB(x)=∑i=0n​bi​xi使得A(x)=B(x)A(x)=B(x)A(x)=B(x) 思路： 考虑先ffpffpffp转点值，然后快速插值即可...
P5324 [BJOI2019]删数 （线段树 求区间 0 的个数）
Top_xiao的博客
05-16 938
题目描述 对于任意一个数列，如果能在有限次进行下列删数操作后将其删为空数列，则称这个数列可以删空。一次删数操作定义如下： 记当前数列长度为kk，则删掉数列中所有等于kk的数。 现有一个长度为nn的数列aa，有mm次修改操作，第ii次修改后你要回答： 经过ii次修改后的数列aa，至少还需要修改几个数才可删空？ 每次修改操作为单点修改或数列整体加一或数列整体减一。 输入输出格式 输入格...
【校内模拟】矩阵求和（组合数转下降幂转自然幂）（斯特林数）（树状数组）
zxyoi_dreamer的博客（不定期诈尸）
02-05 487
简要题意： 一个 n×mn\times mn×m 矩阵，第 iii 行第 jjj 列的权值为 (i−1)⋅m+j(i-1)\cdot m + j(i−1)⋅m+j，需要你支持一下三种操作： R，交换两行 C，交换两列 Q，询问对一个子矩阵求 kkk 次二维前缀和后矩阵中元素之和。 数据范围n,m,Q≤1e5,k≤10n,m,Q\leq 1e5,k\leq 10n,m,Q≤1e5,k≤10 题...
斯特林数与上升幂与下降幂
00ffcc的博客
09-01 2344
重新审视 Stirling Number 转载自http://www.mamicode.com/info-detail-1982164.html 原文地址http://www.cnblogs.com/Sdchr/p/7462045.html 原文markdown挂掉没法看了，转载顺便修一下 无标号计数, 带标号计数 我们经常处理的有两种计数类型: 组合计数, 排列计数 但是, 我更喜欢这样称呼: 组合计数为无标号计数,排列计数为带标号计数 因为它能反应我们直观看到的 无标号计数有一个等价的概念:定序计数
【模板】【洛谷P5393】下降幂多项式转普通多项式（拉格朗日插值）（NTT）
zxyoi_dreamer的博客（不定期诈尸）
08-15 270
传送门 题解： 首先对于下降幂多项式，其系数表达的点值EGF只需要乘上一个exe^xex。 于是我们得到了若干个点值，直接份治拉格朗日插值即可。 代码： #include<bits/stdc++.h> #define ll long long #define re register #define cs const using std::cerr; using std::cout...
【第一类斯特林数】自然数幂和
07-12 753
用带符号第一类斯特林数求自然数幂和， 先推一发第一类斯特林数，设F(x,n)=x∗(x−1)∗.....∗(x−n+1)F(x,n)=x*(x-1)*.....*(x-n+1) =∑i=1nxi∗Sn,i=\sum_{i=1}^nx^i*S_{n,i} 又因为F(x,n−1)∗(x−n+1)=F(x,n)F(x,n-1)*(x-n+1)=F(x,n) 把两个式子拆开，搞一下即可得出： Si
洛谷 P5322 [BJOI2019]排兵布阵
Red_Sky_ta的博客
12-04 454
P5322 [BJOI2019]排兵布阵 算法思路：分组背包
【背包DP练习】洛谷 P5020货币系统 P1757通天之分组背包 P1064[NOIP2006 提高组]金明的预算方案 P5322 [BJOI2019]排兵布阵
wldcmzy的博客
03-17 259
洛谷 P5020货币系统 https://www.luogu.com.cn/problem/P5020 思路是把货币从小到大排序，然后按顺序依次完全背包dp，每次dp检查i-1种面值的货币能不能凑出第i种货币 . 最后四个样例数据范围比较大，但应该都不是需要时间特别长的，而是卡一些技巧。 80分代码：修改一（不是亿）点点就可以满分 #include <cstdio> #include <iostream> #include <cstring> #include <a
BJOI2018 双人猜数游戏 题解
最新发布
09-14
状态转移方程可以这样设置：如果范围只有一个数，那么猜测次数就是0；如果范围内的数已经被排除，那么猜测次数就是上一个状态的值；否则，我们需要考虑两种情况，一是猜测范围的左边界，二是猜测右边界，取这两个...
BJOI 2013. 压力1
08-03
【BJOI 2013. 压力1】是一道典型的图论问题，主要涉及无向图、路径操作和树剖等算法。题目要求处理一个包含N个顶点和M条边的无向图，以及Q个点对，对于每个点对(p, q)，需要将从p到q的所有路径上的所有顶点的权值...
下降幂+第二类斯特林数 推柿子大题
sophilex的博客
03-31 382
然后还有一个东西极大阻碍了我们的推柿子进程，就是f(k),它跟外面的组合数并不能进行什么神奇操作，所以我们可以考虑将其转化一下形式。可以稍微想一想，不难发现，有了上述的性质，我们可以把一个跟k有关的下降幂变成与k无关的一个下降幂，就多了很多操作的空间。,没错它其实就是组合数，只是换了一个名字，但正是因为它与幂函数有许多可以类比的地方，所以也叫它下降幂。不难发现，前面已经是一个二项式定理的样子了，这也就意味着我们成功把n的复杂度降下来了！考虑这样一个组合意义：有k种颜色，每一个人选择任意一种的方案数。
P6620 [省选联考 2020 A 卷] 组合数问题（斯特林数、下降幂）
BUG_Creater_jie的博客
03-16 514
好吃的柿子
斯特林数
yuedongdehuoyan的博客
07-26 234
斯特林数有两大类, 第一类:将n个不同元素构成m个圆排列的数目. 1.对无符号整数:dp[i][j]=dp[i-1][j-1]+(i-1)*dp[i-1][j]; 2.对有符号整数:dp[i][j]=dp[i-1][j-1]-(i-1)dp[i-1][j]; 第二类:表示将n个不同的元素分成m个集合的方案数 dp[i][j]=dp[i-1][j-1]+j*dp[i-1][j]; 具体原...
题解 [联合省选 2020 A] 组合数问题（LOJ #3300 / 洛谷 P6620）【第二类斯特林数 下降幂多项式】
破壁人五号的博客
06-24 423
题目链接：洛谷 P6620 / LOJ #3300 题意 求：∑k=0nf(k)×xk×(nk) mod p\sum_{k=0}^n f(k)\times x^k\times{n\choose k}\bmod p∑k=0n​f(k)×xk×(kn​)modp，其中 fff 为一个 mmm 次多项式，n,xn,xn,x 给定，ppp 为给定的数（不一定是质数）。n,x,p≤109n,x,p\leq 10^9n,x,p≤109，m≤103m\leq 10^3m≤103。 题解 多项式与组合数不太搭，我们假装已经
下降幂与多项式乘法
Freopen的博客
07-11 2764
下降幂指数生成函数：f(x)=∑i=0aii!xi‾f(x) = \sum_{i=0}\frac {a_i}{i!} x^{\underline i}f(x)=∑i=0​i!ai​​xi​ 指数生成函数：f2(x)=∑i=0aii!xif_2(x) = \sum_{i=0}\frac {a_i}{i!} x^{i}f2​(x)=∑i=0​i!ai​​xi 对于f(x)f(x)f(x)的点值指数生成...
普通多项式转下降幂多项式小结
a_forever_dream的博客
07-02 801
前置芝士 斯特林数 正题 对于给定的多项式 F(x)=∑i=0naixiF(x)=\sum_{i=0}^n a_i x^iF(x)=∑i=0n​ai​xi，求出一个下降幂多项式 G(x)=∑i=0nbixi‾G(x)=\sum_{i=0}^n b_i x^{\underline i}G(x)=∑i=0n​bi​xi​，满足 F(x)=G(x)F(x)=G(x)F(x)=G(x)。 主要用到的技巧就是利用第二类斯特林数展开 xnx^nxn，即 xn=∑j=0nCxjj!S(n,j)x^n=\sum_{j=
省选模拟赛20200212 算算算（巧用斯特林数+维护下降幂）
cqbzcsq的博客
02-13 303
题解 以后看到幂和并且幂次比较小的时候尽量往斯特林数方面想 关于题解，有一个巧妙的巧妙的实现就是把 i！乘进组合数 接下来就只需要维护每个值+1之后的下降k次幂之和（把k取0~100的下降幂之和都要动态维护） 想了我好久。。。 后来看到了标程，就一句话f[i][k]=f[i-1][k-1]*i+f[i-1][k]，妙啊！！！！ 展开一下，发现每个数的下降k-1...
具体数学-第6课（下降阶乘幂）
算法码上来
04-18 988
原文链接： 具体数学-第6课 - WeiYang Blog 上节课讲到下降阶乘幂和差分运算，这节课继续讲它和差分的各种性质。 性质1 首先在后面章节会证明， 的二项展开形式和普通的 是一样的，这里提一下，暂时用不到。 性质2 接下来给出下降阶乘幂为负数的定义： 性质3 和普通幂 不同，下降阶乘幂有如下性质： ...
自然数幂和 斯特林数
LYD729
07-12 1573
2018 UPD： 其实第二类斯特林数做自然数幂和更简单，这里简单写一下： 由一个基本式子出发 nk=∑i=0k{k i}[n]in^k=\sum_{i=0}^k\big\{^k_{\ i}\big\}[n]_i 考虑对nn求和 Ans=∑i=0nik=∑i=0n∑j=0k{kj}[i]jAns=\sum_{i=0}^ni^k=\sum_{i=0}^n\sum_{j=0}^k\big\{
BJOI2018双人猜数游戏题解深度解析
在解题的过程中，理解游戏规则是第一步。通常，游戏会有一个明确的数域范围，玩家必须在规定次数内猜出目标数字。每次猜测后，系统会提示猜测结果是太高、太低还是正确。玩家必须根据提示调整其下一次的猜测值。 接...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值