結果

問題 No.2527 H and W
ユーザー syndro_6syndro_6
提出日時 2023-11-03 21:47:07
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 29 ms / 2,000 ms
コード長 4,410 bytes
コンパイル時間 2,114 ms
コンパイル使用メモリ 205,940 KB
実行使用メモリ 19,340 KB
最終ジャッジ日時 2024-09-25 19:54:57
合計ジャッジ時間 3,333 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,940 KB
testcase_02 AC 18 ms
18,828 KB
testcase_03 AC 2 ms
6,944 KB
testcase_04 AC 29 ms
18,896 KB
testcase_05 AC 12 ms
6,944 KB
testcase_06 AC 28 ms
18,856 KB
testcase_07 AC 2 ms
6,940 KB
testcase_08 AC 19 ms
18,856 KB
testcase_09 AC 18 ms
18,840 KB
testcase_10 AC 19 ms
18,960 KB
testcase_11 AC 12 ms
6,940 KB
testcase_12 AC 28 ms
18,832 KB
testcase_13 AC 26 ms
18,792 KB
testcase_14 AC 20 ms
18,788 KB
testcase_15 AC 12 ms
6,940 KB
testcase_16 AC 18 ms
18,864 KB
testcase_17 AC 19 ms
18,784 KB
testcase_18 AC 18 ms
18,832 KB
testcase_19 AC 17 ms
18,996 KB
testcase_20 AC 18 ms
19,340 KB
testcase_21 AC 18 ms
18,796 KB
testcase_22 AC 18 ms
18,668 KB
testcase_23 AC 18 ms
18,868 KB
testcase_24 AC 14 ms
13,480 KB
testcase_25 AC 16 ms
16,492 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

// (⁠◕⁠ᴗ⁠◕⁠✿⁠)
 
#include <bits/stdc++.h>
#define rep(i, n) for (ll i = 0; i < (n); i++)
#define srep(i, s, n) for (ll i = s; i < (n); i++)
#define len(x) ((int)(x).size())
#define all(x) (x).begin(), (x).end()
using namespace std;
template<typename T> using vc = vector<T>;
template<typename T> using vv = vc<vc<T>>;
template<typename T> using vvv = vv<vc<T>>;
using vi = vc<int>;using vvi = vv<int>; using vvvi = vv<vi>;
using ll = long long;using vl = vc<ll>;using vvl = vv<ll>; using vvvl = vv<vl>;
using ld = long double; using vld = vc<ld>; using vvld = vc<vld>; using vvvld = vc<vvld>;
using uint = unsigned int;
using ull = unsigned long long;
const ld pi = 3.141592653589793;
const int inf = 0x3f3f3f3f;
const ll INF = 0x3f3f3f3f3f3f3f3f;
// const ll mod = 1000000007;
const ll mod = 998244353;
inline bool inside(long long y, long long x, long long H, long long W) {return 0 <= (y) and (y) < (H) and 0 <= (x) and (x) < (W); }

#define debug(var)  do{std::cout << #var << " : \n";view(var);}while(0)
template<typename T> void view(T e){cout << e << endl;}
template<typename T> void view(const vc<T>& v){for(const auto& e : v){ cout << e << " "; } cout << endl;}
template<typename T> void view(const vv<T>& vv){ for(const auto& v : vv){ view(v); } }

struct modint {
    ll val;
    using mint = modint;
    modint() {};
    modint(ll v) {
        if (0 <= v && v < mod) val = v;
        else{
            v %= mod;
            if (v < 0) v += mod;
            val = v;
        }
    }
    mint operator++(int){
        val++;
        if (val == mod) val = 0;
        return *this;
    }
    mint operator--(int){
        if (val == 0) val = mod;
        val--;
        return *this;
    }
    mint& operator+=(const mint& v){
        val += v.val;
        if (val > mod) val -= mod;
        return *this;
    }
    friend mint operator+(const mint& u, const mint& v){
        return mint(u) += v;
    }
    mint& operator-=(const mint& v){
        val -= v.val;
        if (val < 0) val += mod;
        return *this;
    }
    friend mint operator-(const mint& u, const mint& v){
        return mint(u) -= v;
    }
    mint& operator*=(const mint& v){
        val *= v.val;
        if (val > mod) val %= mod;
        return *this;
    }
    friend mint operator*(const mint& u, const mint& v){
        return mint(u) *= v;
    }
    mint pow(ll n){
        mint a = *this, ret = 1;
        while (n > 0){
            if (n & 1) ret *= a;
            a *= a;
            n >>= 1;
        }
        return ret;
    }
    mint operator/=(const mint& v){
        val *= mint(v).pow(mod - 2).val;
        if (val > mod) val %= mod;
        return *this;
    }
    friend mint operator/(const mint& u, const mint& v){
        return mint(u) /= v;
    }
    friend bool operator==(const mint& u, const mint& v){
        return u.val == v.val;
    }
    friend bool operator!=(const mint& u, const mint& v){
        return u.val != v.val;
    }
    friend bool operator>(const mint& u, const mint& v){
        return u.val > v.val;
    }
    friend bool operator<(const mint& u, const mint& v){
        return u.val < v.val;
    }
};

using mint = modint;

struct Binoms {
    public : 
        Binoms(int n){
            fact = vc<mint>(1, 1);
            inv_fact = vc<mint>(1, 1);
            update(n);
        };

        void update(int n){
            if (sz >= n) return;
            fact.resize(n + 1, 1);
            inv_fact.resize(n + 1, 1);
            for (int i = sz; i < n; i++) fact[i + 1] = fact[i] * (i + 1);
            inv_fact[n] = 1 / fact[n];
            for (int i = n - 1; i >= sz + 1; i--) inv_fact[i] = inv_fact[i + 1] * (i + 1);
            sz = n;
            return;
        }

        mint C(int n, int k){
            if (sz < n) update(n);
            return fact[n] * inv_fact[n - k] * inv_fact[k];
        }

        mint fac(int n){
            if (sz < n) update(n);
            return fact[n];
        }

    private :
        int sz = 0;
        vc<mint> fact, inv_fact;
};

int main(){
    ll H, W, K; cin >> H >> W >> K;
    mint ret = 0;
    ll i = 1;
    Binoms B(1);
    while (i * i <= K){
        if (K % i == 0){
            ll p = i, q = K / i;
            if (p <= H && q <= W) ret += B.C(H, p) * B.C(W, q);
            if (q <= H && p <= W && p != q) ret += B.C(H, q) * B.C(W, p);
        }
        i++;
    }
    cout << ret.val << endl;
}
0