結果

問題 No.2857 Div Array
ユーザー applejam
提出日時 2025-10-02 13:42:40
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 30 ms / 2,000 ms
コード長 1,825 bytes
コンパイル時間 7,288 ms
コンパイル使用メモリ 334,096 KB
実行使用メモリ 7,720 KB
最終ジャッジ日時 2025-10-02 13:42:52
合計ジャッジ時間 9,589 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <algorithm>
#include <cmath>
#include <bitset>
#include <set>
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
using ll = long long;
using mint = modint998244353;
using vi = vector<int>;
using vvi = vector<vi>;
using vvvi = vector<vvi>;
using vll = vector<ll>;
using vvll = vector<vll>;
using vvvll = vector<vvll>;
using vmi = vector<mint>;
using vvmi = vector<vmi>;
using vvvmi = vector<vvmi>;
#define all(a) (a).begin(), (a).end()
#define rep2(i, m, n) for (int i = (m); i < (n); ++i)
#define rep(i, n) rep2(i, 0, n)
#define drep2(i, m, n) for (int i = (m)-1; i >= (n); --i)
#define drep(i, n) drep2(i, n, 0)

vvmi mat(const vvmi &a, const vvmi &b, int len){
    vvmi res(len, vmi(len));
    rep(i, len)rep(j, len){
        mint tmp = mint(0);
        rep(k, len)tmp += a[i][k]*b[k][j];
        res[i][j] = tmp;
    }
    return res;
}


int main(){

    int n, m, k; cin >> n >> m >> k;
    set<int> st;
    rep2(i, 1, m+1)st.insert(m/i);
    vi id; for(auto i : st)id.push_back(i);
    vi cnt(m+1, 0); rep2(i, 1, m+1)cnt[m/i]++;
    int len = id.size();

    vvmi base(len, vmi(len, mint(0)));
    rep(i, len)rep(j, len){
        if(abs(id[i] - id[j]) <= k)base[i][j] = mint(cnt[id[i]]);
    }

    vvvmi matrix; matrix.push_back(base);
    rep(i, 31){
        int l = matrix.size();
        matrix.push_back(mat(matrix[l-1], matrix[l-1], len));
    }

    vvmi e(len, vmi(len, mint(0))); rep(i, len)e[i][i] = mint(1);
    bitset<32> bs(n-1);
    rep(i, 32){
        if(bs[i])e = mat(e, matrix[i], len);
    }
    
    vmi ft(len); rep(i, len)ft[i] = mint(cnt[id[i]]);
    mint ans = mint(0);
    rep(i, len){
        rep(j, len){
            ans += e[i][j] * ft[j];
        }
    }
    cout << ans.val() << endl;
    return 0;
}
0