結果

問題 No.2818 A Game I Play to Pass the Time
ユーザー zeta7532zeta7532
提出日時 2024-07-19 13:42:23
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 1,773 bytes
コンパイル時間 2,037 ms
コンパイル使用メモリ 210,516 KB
実行使用メモリ 10,624 KB
最終ジャッジ日時 2024-07-19 13:42:33
合計ジャッジ時間 9,815 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
10,624 KB
testcase_01 TLE -
testcase_02 -- -
testcase_03 -- -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
using namespace std;
using ll = long long;
const ll mod = 998244353;
#define fi first
#define se second
#define rep(i,n) for(ll i=0;i<n;i++)
#define all(x) x.begin(),x.end()
#define faster ios::sync_with_stdio(false);cin.tie(nullptr)

ll modpow(ll A,ll N){
    ll res=1;
    while(N){
        if(N&1) res=res*A%mod;
        A=A*A%mod;
        N>>=1;
    }
    return res;
}

ll modinv(ll A){return modpow(A,mod-2);}

int main() {
    vector<ll> P;
    for(ll i=2;i<=100;i++){
        bool ok=true;
        for(ll j=2;j<i;j++){
            if(i%j==0) ok=false;
        }
        if(ok) P.push_back(i);
    }
    vector<vector<ll>> Q(25,vector<ll>(100,1));
    rep(i,25){
        rep(j,100-1){
            Q[i][j+1]=Q[i][j];
            Q[i][j+1]*=P[i];
            Q[i][j+1]%=mod;
        }
    }
    ll T;
    cin >> T;
    while(T--){
        ll X,N;
        cin >> X >> N;
        vector<ll> R(100,1);
        rep(i,100-1){
            R[i+1]=R[i];
            R[i+1]*=((N+i-1)%mod);
            R[i+1]%=mod;
            R[i+1]*=modinv(i+1);
            R[i+1]%=mod;
        }
        ll loop=1;
        vector<ll> A(25,0);
        rep(i,25){
            while(X%P[i]==0){
                A[i]++;
                X/=P[i];
            }
            loop*=(A[i]+1);
        }
        ll ans=0;
        rep(LOOP,loop){
            ll Y=LOOP;
            ll Z=N;
            rep(i,25){
                Z*=Q[i][(Y%(A[i]+1))];
                Z%=mod;
                Z*=R[A[i]-(Y%(A[i]+1))];
                Z%=mod;
                Y/=(A[i]+1);
            }
            ans+=Z;
            ans%=mod;
        }
        cout << ans << endl;
    }
    return 0;
}
0