牛客练习赛76 C-CG的通关秘籍 组合+数学推导
传送门: https://ac.nowcoder.com/acm/contest/10845/C
题意
有
n
个
位
置
,
[
1..
m
]
中
的
数
无
穷
多
个
,
构
造
所
有
可
能
的
序
列
。
有n个位置,[1..m]中的数无穷多个,构造所有可能的序列。
有n个位置,[1..m]中的数无穷多个,构造所有可能的序列。
每
个
序
列
的
贡
献
取
决
于
,
若
a
i
>
a
i
−
1
,
贡
献
+
1
,
若
a
i
<
a
i
−
1
,
贡
献
−
1.
相
等
没
有
贡
献
。
每个序列的贡献取决于,若a_i>a_{i-1},贡献+1,若a_i<a_{i-1},贡献-1.相等没有贡献。
每个序列的贡献取决于,若ai>ai−1,贡献+1,若ai<ai−1,贡献−1.相等没有贡献。
问 : 构 造 完 所 有 可 能 的 序 列 之 后 的 贡 献 和 , 请 m o d 1 e 9 + 7 问:构造完所有可能的序列之后的贡献和,请mod\;1e9+7 问:构造完所有可能的序列之后的贡献和,请mod1e9+7
思路
因为贡献只取决于第i位和第i-1位,然后其他位置随便取,则为该贡献 ∗ m n − 2 . *m^{n-2}. ∗mn−2.
然后从第二个位置开始枚举,假设枚举到了第i位。
第i位上可以取[1,m]之间的所有数,假设取一个j,则对于前一位来说有j-1个比它小,m-j个比它大。
所以贡献为 ( j − 1 ) + 2 ∗ ( m − j ) 。 (j-1)+2*(m-j)。 (j−1)+2∗(m−j)。
所以该位置上的总贡献为 [ ∑ j = 1 m ( j − 1 ) + ( m − j ) ] ∗ m n − 2 [\sum_{j=1}^m(j-1)+(m-j)]*m^{n-2} [∑j=1m(j−1)+(m−j)]∗mn−2
最后结合一下所有位置(注意是从第二位开始取),贡献和为 ∑ i = 2 n [ ∑ j = 1 m ( j − 1 ) + 2 ∗ ( m − j ) ] ∗ m n − 2 \sum_{i=2}^n[\sum_{j=1}^m(j-1)+2*(m-j)]*m^{n-2} i=2∑n[j=1∑m(j−1)+2∗(m−j)]∗mn−2
整
理
一
下
得
:
整理一下得:
整理一下得:
3
m
2
−
3
m
2
∗
m
n
−
2
∗
(
n
−
1
)
\frac{3m^2-3m}{2}*m^{n-2}*(n-1)
23m2−3m∗mn−2∗(n−1)
Code(702MS)
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<int, int> pdd;
#define INF 0x3f3f3f3f
#define lowbit(x) x & (-x)
#define mem(a, b) memset(a , b , sizeof(a))
#define FOR(i, x, n) for(int i = x;i <= n; i++)
// const ll mod = 998244353;
const ll mod = 1e9 + 7;
// const double eps = 1e-6;
// const double PI = acos(-1);
// const double R = 0.57721566490153286060651209;
ll quick_pow(ll a, ll b) {
ll ans = 1;
while(b) {
if(b & 1) ans = ans * a % mod;
a = a * a % mod;
b >>= 1;
}
return ans % mod;
}
void solve() {
int _; scanf("%d",&_);
ll inv2 = quick_pow(2, mod - 2);
while(_--) {
ll n, m; scanf("%lld%lld",&n,&m);
if(m == 1) {
printf("0\n");
continue;
}
ll ans = m * (3 * m - 3) % mod * inv2 % mod;
ans = ans * quick_pow(m, n - 2) % mod * (n - 1) % mod;
printf("%lld\n",(ans + mod) % mod);
}
}
signed main() {
solve();
}