第二类斯特林数总结

最主要的两个式子:

套路1:

$$\begin{array}{l}x^k=\sum_{i=0}^x\begin{pmatrix}x\\i\end{pmatrix}\begin{Bmatrix}k\\i\end{Bmatrix}i!\\\end{array}$$

左边的式子可以看成将k个球放到x个有标号的盒子中,总共的方案数。

右边的式子可以看成先选出i个非空的盒子,将k个球放到这i个盒子中且不空的方案数。但由于第二类斯特林数是无序的,还要乘上i阶乘。

套路2:

$$\begin{array}{l}\begin{Bmatrix}n\\m\end{Bmatrix}=\frac1{m!}\;\sum_{i=0}^m\;(-1)^i\;\;\begin{pmatrix}m\\i\end{pmatrix}(m-i)^n\\\end{array}$$

这是第二类斯特林数的容斥形式。

右边相当于枚举有多少个空的盒子,然后剩下的盒子随便放。由于这是有标号的结果,最后要除以m阶乘。

有了套路2,就能NTT在O(nlogn)时间内算出$$\begin{array}{l}\\\begin{Bmatrix}n\\0\end{Bmatrix}\sim\begin{Bmatrix}n\\1\end{Bmatrix}\end{array}$$>。

1.求第一类斯特林数S(n,0)~S(n,m),n,m<=5E5,模998244353。

朴素NTT即可。

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 typedef long long int ll;
 4 const ll mod=998244353;
 5 const ll G=3;
 6 const ll Gi=332748118;
 7 const ll maxn=1E6+5;
 8 ll a[maxn],b[maxn],fac[maxn],n,m,len,limit,sum[maxn],ans;
 9 int r[maxn];
10 int re(int x)
11 {
12     int ans=0;
13     for(int i=0;i<len;++i)ans=ans*2+((x&(1<<i))>0);
14     return ans;
15 }
16 ll qpow(ll x,ll y)
17 {
18     ll ans=1,base=x;
19     while(y)
20     {
21         if(y&1)ans=ans*base%mod;
22         base=base*base%mod;
23         y>>=1;
24     }
25     return ans;
26 }
27 void init()
28 {
29     fac[0]=1;
30     for(int i=1;i<=n;++i)fac[i]=fac[i-1]*i%mod;
31     for(int i=0;i<=n;++i)b[i]=qpow(fac[i],mod-2)*qpow(i,n)%mod;
32     for(int i=0;i<=n;++i)a[i]=(qpow(-1,i)*qpow(fac[i],mod-2)%mod+mod)%mod;
33 }
34 void NTT(ll*A,int g)
35 {
36     for(int i=0;i<limit;++i)
37         if(r[i]<i)swap(A[i],A[r[i]]);
38     for(int i=2;i<=limit;i*=2)
39     {
40         ll w;
41         if(g==1)w=qpow(G,(mod-1)/i);
42         else w=qpow(Gi,(mod-1)/i);
43         for(int j=0;j<limit/i;++j)
44         {
45             ll d=1;
46             for(int k=0;k<i/2;++k)
47             {
48                 ll a=A[i*j+k],b=A[i*j+i/2+k]*d%mod;
49                 A[i*j+k]=(a+b)%mod;
50                 A[i*j+i/2+k]=(a-b+mod)%mod;
51                 d=d*w%mod;
52             }
53         }
54     }
55 }
56 int main()
57 {
58     ios::sync_with_stdio(false);
59     cin>>n;
60     init();
61     len=0,limit=1;
62     while(limit<=n*2)
63     {
64         ++len;
65         limit<<=1;
66     }
67     for(int i=n+1;i<=limit;++i)a[i]=0;
68     for(int i=n+1;i<=limit;++i)b[i]=0;
69     for(int i=0;i<limit;++i)r[i]=re(i);
70     NTT(a,1);
71     NTT(b,1);
72     for(int i=0;i<limit;++i)a[i]=a[i]*b[i]%mod;
73     NTT(a,-1);
74     ll g=qpow(limit,mod-2);
75     for(int i=0;i<limit;++i)a[i]=a[i]*g%mod;
76     for(int i=0;i<=n;++i)cout<<a[i]<<" ";cout<<endl;
77     return 0;
78 }

2.求$$\sum_{i=0}^n\;\sum_{j=0}^i\;\begin{Bmatrix}i\\j\end{Bmatrix}\times2^j\times(j!)\;$$

原文地址:https://www.cnblogs.com/GreenDuck/p/10758629.html

时间: 2024-11-08 04:18:49

第二类斯特林数总结的相关文章

Gym 101147G 第二类斯特林数

大致题意: n个孩子,k场比赛,每个孩子至少参加一场比赛,且每场比赛只能由一个孩子参加.问有多少种分配方式. 分析: k>n,就无法分配了. k<=n.把n分成k堆的方案数乘以n的阶乘.N分成k堆得方案数即第二类斯特林数 http://blog.csdn.net/acdreamers/article/details/8521134 #include <bits/stdc++.h> using namespace std; typedef long long ll; const ll

Light OJ 1236 Race 第二类斯特林数

第二类斯特林数 n 匹马 分成1 2 3... n组 每一组就是相同排名 没有先后 然后组与组之间是有顺序的 在乘以组数的阶乘 #include <cstdio> #include <cstring> using namespace std; int dp[1010][1010]; int a[1010]; int main() { a[0] = 1; dp[0][0] = 1; for(int i = 1; i <= 1000; i++) { dp[i][0] = 0; d

swjtu oj Paint Box 第二类斯特林数

http://swjtuoj.cn/problem/2382/ 题目的难点在于,用k种颜色,去染n个盒子,并且一定要用完这k种颜色,并且相邻的格子不能有相同的颜色, 打了个表发现,这个数是s(n, k) * k! s(n, k)表示求第二类斯特林数. 那么关键是怎么快速求第二类斯特林数. 这里提供一种O(k)的算法. 第二类斯特林数: #include <cstdio> #include <cstdlib> #include <cstring> #include <

hdu 2512 一卡通大冒险(第二类斯特林数)

递推思路如下,i张卡片分成j堆,那么分为两种情况:第i张卡片自成一堆或没有自成一堆. 那么自成一堆的话就是dp[i-1][j-1]种情况 不自成一堆的话就是就能在j堆种任意挑一堆放入,所以有dp[i-1][j]*j种情况 综上,如下: dp[i][j]=dp[i-1][j]*j+dp[i-1][j-1]. 关于第二类斯特林数,百度就好. 具体代码 #include <iostream> using namespace std; int dp[2005][2005]; int main() {

poj 3088 组合计数,第二类斯特林数

题意:给出n个数字[1,n],问你可以组成多少种不同的序列,其中每一个序列由几个部分组成,每个部分包含几个数字,每部分内数字无序,部分之间数字有序.每个数字最多使用一次,可以不用. 思路:枚举从n个数字中选出i个数字(组合数),再枚举将这i个数字分成j个部分(第二类斯特林数),然后乘上j的全排列. 1 #include <iostream> 2 #include <cstring> 3 #include <cstdio> 4 using namespace std; 5

LightOJ 1326 - Race(第二类斯特林数啊 )

题目链接:http://lightoj.com/volume_showproblem.php?problem=1326 百度百科:斯特林数 ACdreamer:第一类Stirling数和第二类Stirling Disky and Sooma, two of the biggest mega minds of Bangladesh went to a far country. They ate, coded and wandered around, even in their holidays.

Gym Gym 101147G 第二类斯特林数

题目链接:http://codeforces.com/gym/101147/problem/G 题意:n个人,去参加k个游戏,k个游戏必须非空,有多少种放法? 分析: 第二类斯特林数,划分好k个集合后乘以阶乘: 1 #include <bits/stdc++.h> 2 3 using namespace std; 4 5 const int maxn = 1010; 6 const long long MOD = 1000000000 + 7L; 7 long long stir[maxn][

poj 1430 第二类斯特林数

1 #include <iostream> 2 #include <cmath> 3 #include <algorithm> 4 using namespace std; 5 6 int get2(long long n){ 7 if(n==0) 8 return 0; 9 int cnt =0; 10 while(n){ 11 cnt += n/2; 12 n = n/2; 13 } 14 return cnt; 15 } 16 int main(){ 17 18

HDU 4045 Machine scheduling (第二类斯特林数+DP)

A Baidu's engineer needs to analyze and process large amount of data on machines every day. The machines are labeled from 1 to n. On each day, the engineer chooses r machines to process data. He allocates the r machines to no more than m groups ,and

关于第二类斯特林数的一丢丢东西

关于第二类斯特林数的一丢丢东西 第二类斯特林数 S(n,m)表示有\(n\)个有区别小球,要放进\(m\)个相同盒子里,且每个盒子非空的方案数 考虑一个很容易的递推: \[S(n,m)=S(n-1,m-1)+m*S(n-1,m)\] 考虑组合意义: 假设前面的\(n-1\)个球丢进了\(m-1\)个组,因为每个组非空,所以这个球只有一种选择--自己一组 如果前面的球已经分成了\(m\)组,那么,这个球就有\(m\)种放法 所以这个递推式就是这样来的 那么,只考虑组合意义可不可以算? 当然是可以的