F.鸡数题！(2024牛客寒假训练营)

思路解析：

条件四中若将和(i!=j)看作二进制的话，那么他们相&等于0意味着相应的二进制位不能全为一

条件三中转化为二进制后为111111...(n个)，结合条件四及条件一二来考虑的话题目的要求就是将n个不同的1分配给m个格子里，问有多少种分法，形象一点来说就是把n个不同的球分配到m个相同的盒子里，不能有空盒子（即条件一），有多少种分法？

举个简单的例子：

假设n=3,m=2，将3个1分到2个格子里，注意每个1都是不同的，假设三个1的代号分别为a,b,c，则可以有{a,{b,c}}，{b,{a,c}}，{c,{a,b}}三种分法，即a1，a2可以分别为{001，110}，{010，101}，{011，100}三种情况

题目本质来说就是求第二类Stirling数

第二类Stirling数的通项公式为：

观察通项公式最右边部分，分子假设为a，分母假设为b，计算过程需要求a除以b然后对1e9+7(假设为m)取模，这里a和b都是很大的数，容易溢出，导致取模错误。用逆可以避免除法计算，设b的逆元为，有：

(a/b)mod m=( (a/b) mod m )( (b) mod m )=(a/b  b) mod m = (a) mod m

=(a mod m)( mod m) mod m

可以通过费马小定理来求。

废话不多说，代码如下：

#include<bits/stdc++.h>
using namespace std;
using ll = long long;
using T = pair<int, int>;
//set<int>S;
//unordered_map<int, int>mp;
const int N = 1e5 + 10;
const ll mod = 1e9 + 7;
ll f[N];
//快速幂求a^n
ll fastPow(ll a, ll n, ll mod) {
	ll ans = 1;
	a %= mod;
	while (n) {
		if (n & 1)ans = (ans * a) % mod;
		a = (a * a) % mod;
		n >>= 1;
	}
	return ans;
}
//费马小定理求a模mod的逆
ll mod_inverse(ll a, ll mod) {
	return fastPow(a, mod - 2, mod);
}
//阶乘预处理
void init_fac() {
	f[0] = 1;
	f[1] = 1;
	ll i=2;
	for (; i <= 100000; ++i) {
		f[i] = (f[i - 1] * i) % mod;
	}
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	ll n, m; cin >> n >> m;
	if (n < m) {
		cout << 0 << endl;
		return 0;
	}
	ll i;
	init_fac();
	ll ans{};
	for (i = 0; i <= m; ++i) {
		if ((m - i) & 1) {
			ans += (-1 * fastPow(i, n,mod))* mod_inverse(f[i]*f[m-i],mod);
			ans %= mod;
		}
		else {
			ans += (fastPow(i, n,mod)) * mod_inverse(f[i] * f[m - i], mod);
			ans %= mod;
		}
	}
	cout << ans;
	return 0;
}