[HNOI2011] 卡农 题解

众所周知卡农是一种复调音乐的写作技法，小余在听卡农音乐时灵感大发，发明了一种新的音乐谱写规则。

他将声音分成 nn 个音阶，并将音乐分成若干个片段。音乐的每个片段都是由 11 到 nn 个音阶构成的和声，即从 nn 个音阶中挑选若干个音阶同时演奏出来。

为了强调与卡农的不同，他规定任意两个片段所包含的音阶集合都不同。同时为了保持音乐的规律性，他还规定在一段音乐中每个音阶被奏响的次数为偶数。

现在的问题是：小余想知道包含 mm 个片段的音乐一共有多少种。
两段音乐 aa 和 bb 同种当且仅当将 aa 的片段重新排列后可以得到 bb。例如：假设 aa 为 \{\{1,2\},\{2,3\}\}{{1,2},{2,3}}，bb 为 \{\{2,3\},\{1,2\}\}{{2,3},{1,2}}，那么 aa 与 bb 就是同种音乐。

答案对 10^8+7108+7 取模。

输入格式

仅一行两个正整数 n,mn,m

输出格式

输出一行一个整数表示答案。

输入输出样例

输入 #1复制

2 3

输出 #1复制

1

说明/提示

【数据范围】
对于 20\%20% 的数据，1\le n,m \le 51≤n,m≤5；
对于 50\%50% 的数据，1\le n,m \le 30001≤n,m≤3000；
对于 100\%100% 的数据，1\le n,m \le 10^61≤n,m≤106。

【样例解释】
音乐为 \{\{1\},\{2\},\{1,2\}\}{{1},{2},{1,2}}

发现题目中定义的同种音乐其实是假的，只要除以m!m!就可以了（利用逆元）。

这样就是在集合S=\{1,2,...,n\}S={1,2,...,n}中选出mm个子集，满足三点性质：

（1）所有选出的mm个子集都不能为空。

（2）所有选出的mm个子集中，不能存在两个完全一样的集合。

（3）所有选出的mm个子集中，11到nn每个元素出现的次数必须是偶数。

DP。定义状态f[i]f[i]表示到第ii个子集，满足所有三点性质的方案数。

考虑从容斥的角度分析，得出转移。由性质（3）得到，确定了前i-1i−1个子集后，第ii个子集也随之确定。这时方案数为A_{2^n-1}^{i-1}A2n−1i−1​。

然后去掉不满足性质（1）的方案数。可以得出，如果第ii个子集为空，那么前i-1i−1个子集可以凑成一个合法的方案。因此不满足性质（1）的方案数为f[i-1]f[i−1]。

最后去掉不满足性质（2）的方案数。假设第jj个子集与第ii个子集重复，那么如果把第jj个子集和第ii个子集去掉，剩下的i-2i−2个子集可以凑成一个合法的方案，即f[i-2]f[i−2]。而jj有i-1i−1种取值，子集ii的取法为2^n-1-(i-2)2n−1−(i−2)。

因此不满足性质（2）的方案数为f[i-2]\times(i-1)\times(2^n-i+1)f[i−2]×(i−1)×(2n−i+1)。 得出转移：

f[i]=A_{2^n-1}^{i-1}-f[i-1]-f[i-2]\times(i-1)\times(2^n-i+1)f[i]=A2n−1i−1​−f[i−1]−f[i−2]×(i−1)×(2n−i+1)

代码：

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int N = 1e6 + 5, MX = 1e8 + 7; int n, m, f[N], orz = 1, mx = 1, A[N];
int qpow(int a, int b) {
	int res = 1; while (b) b & 1 ? res = 1ll * res * a % MX : 0,
		a = 1ll * a * a % MX, b >>= 1; return res;
}
int main() {
	int i; cin >> n >> m; for (i = 1; i <= n; i++) orz = orz * 2 % MX;
	orz = (orz - 1 + MX) % MX; A[0] = 1; for (i = 1; i <= m; i++)
		A[i] = 1ll * A[i - 1] * ((orz - i + 1 + MX) % MX) % MX,
	mx = 1ll * mx * i % MX; f[f[1] = 0] = 1; for (i = 2; i <= m; i++)
		f[i] = (A[i - 1] - f[i - 1] + MX - 1ll * f[i - 2] *
		(i - 1) % MX * (orz - i + 2 + MX) % MX + MX) % MX;
	cout << 1ll * f[m] * qpow(mx, MX - 2) % MX << endl; return 0;
}