結果

問題 No.3776 Magic Square of Product
コンテスト
ユーザー t98slider
提出日時 2026-10-23 22:12:59
言語 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  
実行時間 56 ms / 2,500 ms
+ 335µs
コード長 2,156 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 4,273 ms
コンパイル使用メモリ 382,688 KB
実行使用メモリ 46,480 KB
最終ジャッジ日時 2026-10-23 22:13:09
合計ジャッジ時間 8,538 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge4_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;
using mint = atcoder::modint998244353;

template <class T> struct csr {
    struct Node {
        csr* g;
        int u;
        template<class... Args>
        void emplace_back(Args&&... args){
            g->add_edge(u, T(std::forward<Args>(args)...));
        }
        auto begin(){ return g->E.begin() + g->start[u]; }
        auto end(){ return g->E.begin() + g->start[u + 1]; }
        int size(){ return g->start[u + 1] - g->start[u]; }
        T& operator[](int p){ return *(begin() + p); }
    };
    int N;
    std::vector<int> start;
    std::vector<T> E;
    std::vector<std::pair<int,T>> edge;
    csr(int n) : N(n), start(n + 1) {edge.reserve(n);}
    void add_edge(int u, T v){
		assert(0 <= u && u < N);
        start[u + 1]++;
        edge.emplace_back(u, v);
    }
    void build(){
        E.resize(edge.size());
        for(int i = 0; i < N; i++) start[i + 1] += start[i];
        auto cnt = start;
        for(auto [u, v] : edge) E[cnt[u]++] = v;
    }
	const int size() {return N;}
    Node operator[](int u) {return Node{this, u};}
};

int main(){
	ios::sync_with_stdio(false);
	cin.tie(0);
	constexpr int th = 1'000'000;
	vector<bool> pr(th + 1);
	csr<int> tb(th + 1);
	for(int i = 2; i <= th; i++){
		if(pr[i]) continue;
		for(int j = i; j <= th; j += i){
			pr[j] = true;
			int v = j, cnt = 0;
			do{
				v /= i;
				cnt++;
			}while(v % i == 0);
			tb[j].emplace_back(cnt);
		}
	}
	tb.build();
	vector<mint> dp(20);
	for(int i = 1; i < 20; i++){
		for(int a = 0; a <= i; a++){
			for(int b = 0; a + b <= i; b++){
				int c = i - (a + b);
				for(int a2 = 0; a2 <= i && a + a2 <= i; a2++){
					for(int b2 = 0; a2 + b2 <= i && b + b2 <= i; b2++){
						int c2 = i - (a2 + b2);
						
						int a3 = i - (a + a2);
						int b3 = i - (b + b2);
						int c3 = i - (c + c2);
						if(a3 + b3 + c3 == i && c3 >= 0) dp[i]++;
					}
				}
			}
		}

	}
	int T;
	cin >> T;
	while(T--){
		int v;
		cin >> v;
		if(v == 1){
			cout << "1\n";
			continue;
		}
		mint ans = 1;
		for(auto &&vl : tb[v]) ans *= dp[vl];
		cout << ans.val() << '\n';
	}
}
0