社区讨论
很疑惑为什么会错
AT_arc004_4 [ARC004D] 表現の自由 ( Freedom of expression )参与者 6已保存回复 10
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 10 条
- 当前快照
- 1 份
- 快照标识符
- @m0rwhyz8
- 此快照首次捕获于
- 2024/09/07 16:46 去年
- 此快照最后确认于
- 2024/09/07 18:12 去年
CPP
#include <bits/stdc++.h>
using namespace std;
#define int long long
int mod = 1e9 + 7;
int jie[200005], inv[200005];
int qpow(int a, int b) {
int ans = 1;
while (b) {
if (b % 2 == 1)
ans = ans * a % mod;
a = a * a % mod;
b /= 2;
}
return ans % mod;
}
int c(int n, int m) {
if (n < m)
return 0;
return jie[n] * inv[m] % mod * inv[n - m] % mod;
}
map<int, int>q;
signed main() {
jie[0] = 1;
inv[0] = qpow(jie[0], mod - 2);
for (int i = 1; i <= 200000; i++) {
jie[i] = jie[i - 1] * i % mod;
inv[i] = qpow(jie[i], mod - 2);
}
int n, m;
cin >> n >> m;
if (n < 0)
n *= -1;
for (int i = 2; i * i <= n; i++) {
while (n % i == 0) {
q[i]++;
n /= i;
}
}
if (n > 1)
q[n]++;
int ans = 1;
for (auto x : q) {
int y = x.second;
ans = ans * c(y + m - 1, y) % mod;
}
cout << ans * qpow(2, m - 1) % mod;
}
服了一直 。
回复
共 10 条回复,欢迎继续交流。
正在加载回复...