社区讨论

很疑惑为什么会错

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;
}
服了一直 WAWA onon 11

回复

10 条回复,欢迎继续交流。

正在加载回复...