2021-08-22

Add or Multiply 1
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 131072/131072 K (Java/Others)
Total Submission(s): 150 Accepted Submission(s): 36

Problem Description
前两段和第五题相同，但你不需要阅读第五题就可以完成这个题目。

你有一个数字 x 和若干个操作，每个操作是 +ai 或者乘 ×ai 中的一种。你可以重新排列这些操作的顺序，然后对数字 x 执行这些操作。

比如说三个操作是 +a1,+a2,×a3。如果按顺序执行这三个操作，那么得到的结果是 ((x+a1)+a2)×a3。如果排列成 +a2,×a3,+a1，那么得到的结果是 ((x+a2)×a3)+a1。

我们会发现，有一些操作顺序计算出来的结果是本质相同的，比如说+a1,+a2,×a3和+a2,+a1,×a3这样运算下来结果是一样的。我们认为两个操作顺序计算的结果本质相同，当且仅当无论代入什么数，计算出来的结果都是一样的。

请问有多少种本质不同的操作序列。换句话说就是最多能找到多少个操作序列，使得这些操作序列任意两个都不是本质相同的。由于答案很大，输出对109+7取模的结果。

在这个题目中，我们只会给出加法操作和乘法操作的个数，分别是n,m，并不会给出具体的顺序和数字。容易发现，答案与具体的顺序和数字无关。

Input
本题有多组测试数据。

第一行一个整数 T(1≤T≤104) 表示数据组数。

对于每组数据，一行两个整数 n,m(1≤n,m≤3000) 描述一组询问，表示操作中加法个数和乘法个数。

Output
对于每组数据输出一行，表示答案。

Sample Input
3
2 1
5 5
100 100

Sample Output
4
329462
294770659

Source
2021 年百度之星·程序设计大赛 - 复赛

很显然，连续的加号和乘号之间交换顺序是本质相同的运算，所以题目可以转换成将加号和乘号进行分组，由于他们的组之间是相互分割的，所以组数之差不超过1，也就是说有三种情况：
1、加号比乘号组数多一个，此时只有一种情况也就是 +×+
2、乘号比加号组数多一个，此时也只有一种情况也就是 ×+×
3、加号和乘号组数一样多，此时有两种情况，也就是×+×+和+×+×
所以我们穷举其中一个的组数，然后分类讨论另一种组数即可。
那么如何求出分组的情况呢，我们用ans[i] [j]表示将i个物品分成j组（组内顺序无关但是组间顺序有关），考虑递推求解。将i个物品分成j组可以转换成将i个物品放入j个盒子中且每个盒子至少有一个物品。ans[i-1] [j]表示当前已经将i-1个物品放入了j个盒子中，需要将最后一个物品放入盒子中，那么现在有两种情况：
1、当前还剩一个空盒子，也就是说当前将i-1个物品放入了j-1个盒子中得到情况数为ans[i-1] [j-1]，此时最后一个物品必须放入这个空盒子中（否则情况不合法，若有两个及以上的空盒子也不合法），此时也就是ans[i-1] [j-1]*1。
2、没有空盒子，也就是当前的情况数为ans[i-1] [j]，那么最后一个物品可以放在任意一个盒子，也就是有j种选择，所以这种情况得到的方案数为
ans[i-1] [j]×j。
ans[0][0]显然是1，ans[i][1]也显然是1（只放入一个盒子当然只有一种情况）。于是我们就可以预处理得到所有情况的方案数。
接下来看代码吧

#include<cstdio>
#include<iostream>
#include<math.h>
#include<algorithm>
#include<cstring>
#include<unordered_map>
#include<string>
#include<queue>
#include<map>
#include<ctime>
#include<limits.h>
#include <stack>
#include<set>
#include<iomanip>
#define ll long long
#define int long long
#define um unordered_map
#pragma warning(disable : 4996)
using namespace std;
#define rep(l,a,b) for(ll l=a;l<=b;l++)
#define repp(l,a,b) for(ll l=a;l>=b;l--)
#define mem(x) memset(x,0,sizeof(x))
#define fr fast_IO::read()
#define INF 0X3f3f3f3f
namespace fast_IO {
	ll read() {
		ll num = 0;
		char c;
		bool tag = false;
		while ((c = getchar()) != '-' && c != '+' && (c < '0' || c > '9') && ~c);
		if (!~c)return EOF;
		if (c == '-')tag = true;
		else if (c == '+')tag = false;
		else num = c ^ 48;
		while ((c = getchar()) >= '0' && c <= '9')
			num = (num << 1) + (num << 3) + (c ^ 48);
		if (tag)return -num;
		return num;
	}
}
//--------------------------------------------------------------
const int mod = 1e9 + 7;
#define N 6005
int n, m, k;
ll mul[N];
ll cc[N][N];
int ans[N][N], jie[N];//把i个操作分成j组的方案

void init()
{
	jie[0] = 1;//预处理阶乘
	for (int i = 1; i <= 3000; ++i) ans[i][1] = 1, jie[i] = jie[i - 1] * i % mod;
	for (int i = 2; i <= 3000; ++i)//预处理方案数
		for (int j = 2; j <= i; ++j)
			ans[i][j] = (ans[i - 1][j - 1] + ans[i - 1][j] * j % mod) % mod;
	for (int i = 2; i <= 3000; ++i)//对组间进行全排序
		for (int j = 2; j <= i; ++j)//所以要乘上j!
			ans[i][j] = ans[i][j] * jie[j] % mod;
}
signed main()
{
	//freopen("in.txt", "r", stdin);
	init();
	int t;
	cin >> t;
	while (t--)
	{
		n = fr, m = fr;
		ll res = 0;
		rep(i, 1, n)//穷举组数
		{//对应上面所说的三种情况，因为ll最大有2e19，所以取两次模即可，取模运算可是很慢的
			res = (res + ans[n][i] * (ans[m][i - 1] + 2 * ans[m][i] + ans[m][i + 1]) % mod) % mod;
		}
		printf("%lld\n", res);
	}
}