赛时看懂这题就是把n个有区别球放入m个无区别箱子中的问题,但是那时候并不知道它是 第二类斯特林数,于是赛时也没去推过式子。赛后看了利用容斥原理对其的相关证明,然后就发现原来是这么回事,于是赛后进行了补题。
第二类斯特林数:
F题的题面如下:
所有数按位或后等于2^n - 1说明n位二进制每一位都需要有1存在,而任意两个数按位与等于0说明不存在相同位存在1的数,3、4两个条件可以类比为将n个有区别球放到m个有区别箱子中,而第二个条件使得任意组合的a数组都会按照从小到大排序,等于排除了顺序因素,将使得原题类比为
将n个有区别球放到m个无区别箱子中。
然后就是我补题的过程,
刚开始就是纯把这个式子实现了一遍,但是没有考虑到m的1e5范围将会使得求组和数的时间非常庞大,然后交上去就T了,只过了13%的样例。
然后我考虑能不能将求组合数变成线性的时间复杂度,于是想到了
这样就可以通过记录上一步算出来的组合数求当前步的组合数达到线性的复杂度。
调整完代码后又交了一次,结果居然Wa在了81%,当时就很不理解了为什么会wa??
后面手写了几个很大的样例,结果发现出现了负数,原因是因为算出的结果可能会到1e-18,只加一个MOD再MOD并不能让它变成正数,改成循环就AC了。下面是我的代码:
#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int N = 2e5 + 10;
const int p = 1e9 + 7;
ll o[N], k[N], m, n, gg = 1, res = 1, y = 1, u = 0;
ll qmi(ll a, ll b)
{
a %= p;
b %= p;
ll res = 1;
while (b)
{
if (b & 1)
res = res * a % p;
a = a * a % p;
b /= 2;
}
return res % p;
}
void solve()
{
cin >> n >> m;
for (int i = 1; i <= m; i++)
k[i] = qmi(i, p - 2) % p;//预处理1~m的逆元
for (int i = 0; i <= m; i++)
o[i] = qmi(m - i, n) % p;//预处理0~m每个数的n次方
for (int i = 1; i <= m; i++)
y = y * i % p;//处理m!
res = qmi(y, p - 2) % p;
u = o[0] % p;
for (int i = 1; i <= m; i++)
{
gg = gg * (m - i + 1) % p;
gg = gg * k[i] % p;
gg %= p;
if (i & 1)
u -= gg * o[i] % p;
else
u += gg * o[i] % p;
}
while (u < 0)
u = (u + p) % p;
cout << u % p * res % p << endl;
}
int main()
{
ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
solve();
}