#include using namespace std; const long long mod = 998244353; template< typename T > vector< pair< pair< T, T >, T > > quotient_range(T N) { T M; vector< pair< pair< T, T >, T > > ret; for(M = 1; M * M <= N; M++) { ret.emplace_back(make_pair(M, M), N / M); } for(T i = M; i >= 1; i--) { T L = N / (i + 1) + 1; T R = N / i; if(L <= R && ret.back().first.second < L) ret.emplace_back(make_pair(L, R), N / L); } return ret; } using ll = long long; using mat = vector>; /// 行列積 mat mat_mul(mat &a, mat &b) { mat res(a.size(), vector(b[0].size())); for (int i = 0; i < a.size(); i++) { for (int j = 0; j < b[0].size(); j++) { for (int k = 0; k < b.size(); k++) { (res[i][j] += a[i][k] * b[k][j]) %= mod; } } } return res; } /// 行列累乗 mat mat_pow(mat a, long long n) { mat res(a.size(), vector(a.size())); // 単位行列で初期化 for (int i = 0; i < a.size(); i++) res[i][i] = 1; // 繰り返し二乗法 while (n > 0) { if (n & 1) res = mat_mul(a, res); a = mat_mul(a, a); n >>= 1; } return res; } int main() { int N, M, K; cin >> N >> M >> K; auto range = quotient_range(M); int k(size(range)); mat dp(1, vector(k)); for(int i = 0; i < k; i++) { dp[0][i] = (range[i].first.second - range[i].first.first + 1); } mat P(k, vector(k)); for(int i = 0; i < k; i++) { for(int j = 0; j < k; j++) { if(K < abs(range[i].second-range[j].second)) continue; P[i][j] = (range[j].first.second - range[j].first.first + 1); } } P = mat_pow(P, N-1); dp = mat_mul(dp, P); long long ans(0LL); for(int i = 0; i < k; i++) { ans = (ans + dp[0][i]) % mod; } cout << ans << endl; return 0; }