結果

問題 No.3776 Magic Square of Product
コンテスト
ユーザー cho435
提出日時 2026-10-23 23:13:08
言語 C++23
(gcc 15.3.0 + boost 1.92.0 + ACL)
コンパイル:
g++-15 -O2 -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 68 ms / 2,500 ms
+ 125µs
コード長 1,901 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 6,316 ms
コンパイル使用メモリ 385,484 KB
実行使用メモリ 42,560 KB
最終ジャッジ日時 2026-10-23 23:13:32
合計ジャッジ時間 8,514 ms
ジャッジサーバーID
（参考情報） 		judge4_0 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
#include "atcoder/all"

using namespace std;
using ll = long long;
#define rep(i, s, t) for (ll i = s; i < (ll)(t); i++)
#define all(x) begin(x), end(x)
template <class T> bool chmin(T& x, T y) {
	return x > y ? (x = y, true) : false;
}
template <class T> bool chmax(T& x, T y) {
	return x < y ? (x = y, true) : false;
}


std::vector<int> prime_table(int n) {
	std::vector<int> p(n + 1, -1);
	for (int i = 2; i * i <= n; i++) {
		if (p[i] != -1) continue;
		for (int j = i * i; j <= n; j += i) p[j] = i;
	}
	return p;
}

std::map<int, int> table_factor(int n) {
	const static std::vector<int> table = prime_table(1e7);
	assert(n <= 1e7);
	std::map<int, int> res;
	while (n > 1) {
		if (table[n] == -1) {
			res[n]++;
			break;
		}
		res[table[n]]++;
		n /= table[n];
	}
	return res;
}

using mint = atcoder::modint998244353;
vector<mint> C = {1, 6, 21, 55, 120, 231, 406, 666, 1035, 1540, 2211, 3081, 4186, 5565, 7260, 9316, 11781, 14706, 18145, 22155, 26796, 32131, 38226, 45150, 52975, 61776, 71631, 82621, 94830, 108345, 123256};

void solve() {
	// vector<vector<int>> add = {
	// 	{1, 0, 0, 0, 1, 0, 0, 0, 1},
	// 	{1, 0, 0, 0, 0, 1, 0, 1, 0},
	// 	{0, 1, 0, 1, 0, 0, 0, 0, 1},
	// 	{0, 1, 0, 0, 0, 1, 1, 0, 0},
	// 	{0, 0, 1, 1, 0, 0, 0, 1, 0},
	// 	{0, 0, 1, 0, 1, 0, 1, 0, 0},
	// };
	// set<vector<int>> nw;
	// nw.insert(vector<int>(9, 0));
	// rep(lp, 0, 30) {
	// 	set<vector<int>> nx;
	// 	for(auto v : nw) {
	// 		for(auto w : add) {
	// 			rep(i, 0, 9) v[i] += w[i];
	// 			nx.insert(v);
	// 			rep(i, 0, 9) v[i] -= w[i];
	// 		}
	// 	}
	// 	swap(nw, nx);
	// 	cout << nw.size() << endl;
	// }
	int n;
	cin >> n;
	mint ans = 1;
	for(auto [x, c] : table_factor(n)) {
		ans *= C[c];
	}
	cout << ans.val() << '\n';
}

int main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	cout << fixed << setprecision(15);
	int t = 1;
	cin >> t;
	while (t--) solve();
}
0