結果

問題 No.3532 Non-Fourth-Power Sets
コンテスト
ユーザー GaLLium
提出日時 2026-03-12 09:48:47
言語 C++23
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 1,489 ms / 6,000 ms
コード長 1,069 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 3,778 ms
コンパイル使用メモリ 340,196 KB
実行使用メモリ 275,928 KB
最終ジャッジ日時 2026-05-04 20:53:04
合計ジャッジ時間 40,046 ms
ジャッジサーバーID
（参考情報） 		judge1_1 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
#include <atcoder/modint>
using mint = atcoder::modint1000000007;

int main(){
    ios_base::sync_with_stdio(false); 
    cin.tie(nullptr);
    
    vector<int> mpf(20000005), primes;
    for(int i = 2;i < ssize(mpf);++i){
        if (mpf[i] == 0){
            mpf[i] = i;
            primes.push_back(i);
        }
        for (int p:primes){
            if (p > mpf[i] || p*i >= ssize(mpf)) break;
            mpf[p*i] = p;
        }
    }
    
    vector<mint> f(48000005);
    f[0] = 1; f[1] = 2; f[2] = 5; f[3] = 15;
    for (int n = 4;n < ssize(f);++n){
        mint n1 = (n+3)/4, n2 = (n+2)/4, n3 = (n+1)/4;
        f[n] = f[n-1] + (n1+n2+n3)*f[n-2] + (n2*n2+n2*n3+n3*n3)*f[n-3] + n3*n3*n3*f[n-4];
    }
    
    int T; cin >> T;
    while(T--){
        int A,B; cin >> A >> B;
        mint ans = 1;
        while (A > 1){
            int p = mpf[A], vp = 0;
            while (mpf[A] == p){
                A /= p; ++vp;
            }
            ans *= f[B*vp];
        }
        cout << ans.val() << '\n';
    }
}
0