結果

問題 No.8123 Calculated N !
ユーザー KumaTachiRen
提出日時 2025-03-29 11:01:21
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,024 ms / 2,000 ms
コード長 996 bytes
コンパイル時間 3,498 ms
コンパイル使用メモリ 287,984 KB
実行使用メモリ 18,884 KB
最終ジャッジ日時 2025-04-01 20:50:45
合計ジャッジ時間 14,262 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 6
other AC * 16
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

#include <atcoder/modint>
using mint = atcoder::modint1000000007;

using ll = long long;
using vl = vector<ll>;

ll floor_sqrt(ll n) {
	ll l = 0, r = 1e9;
	while (r - l > 1) {
		ll x = (l + r) / 2;
		(x * x <= n ? l : r) = x;
	}
	return l;
}

pair<vl, vl> prime_count(ll n) {
	ll sq = floor_sqrt(n);
	vl x, c;
	for (ll i = 1; i <= sq; i++) x.push_back(n / i);
	for (ll i = x.back() - 1; i; i--) x.push_back(i);
	unordered_map<ll, ll> s;
	for (auto i : x) s[i] = i - 1;
	for (ll p = 2; p <= sq; p++) {
		if (s[p] <= s[p - 1]) continue;
		ll t = s[p - 1];
		for (auto i : x) {
			if (i < p * p) break;
			s[i] -= s[i / p] - t;
		}
	}
	reverse(x.begin(), x.end());
	for (auto i : x) c.push_back(s[i]);
	return {x, c};
}

int main() {
	ll n;
	cin >> n;
	auto [x, c] = prime_count(n);
	mint ans = 1;
	for (int i = 1; i < x.size(); i++) {
		mint s = 1;
		for (ll v = n; v;) s += v /= x[i];
		ans *= s.pow(c[i] - c[i - 1]);
	}
	cout << ans.val() << "\n";
}
0