結果

問題 No.2356 Back Door Tour in Four Seasons
コンテスト
ユーザー AngrySadEight
提出日時 2023-04-18 14:37:50
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 584 ms / 2,000 ms
コード長 2,791 bytes
コンパイル時間 849 ms
コンパイル使用メモリ 77,064 KB
最終ジャッジ日時 2025-02-12 09:55:52
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
using namespace std;
typedef long long ll;

ll mod2 = 998244353;

ll my_pow(ll x, ll n, ll mod){
    //　繰り返し二乗法．x^nをmodで割った余り．
    ll ret;
    if (n == 0){
        ret = 1;
    }
    else if (n % 2 == 1){
        ret = (x * my_pow((x * x) % mod, n / 2, mod)) % mod;
    }
    else{
        ret = my_pow((x * x) % mod, n / 2, mod);
    }
    return ret;
}

ll inv(ll x, ll mod){
    return my_pow(x, mod - 2, mod);
}

int main(){
    ll N;
    cin >> N;
    vector<char> S(N);
    vector<ll> A(N);
    ll all_sum = 0;
    ll spring = 0;
    for (ll i = 0; i < N; i++){
        cin >> S[i] >> A[i];
        all_sum += A[i];
        if (S[i] == 'P') spring++;
    }
    ll summer = 0;
    ll fall = 0;
    ll winter = 0;
    for (ll i = 0; i < N; i++){
        ll prd = my_pow(N - 1, A[i], mod2);
        prd = (prd + mod2 - my_pow(N - 2, A[i], mod2)) % mod2;
        prd = (prd * inv(my_pow(N - 1, A[i], mod2) , mod2)) % mod2;
        if (S[i] == 'U') summer = (summer + prd) % mod2;
        if (S[i] == 'F') fall = (fall + prd) % mod2;
        if (S[i] == 'W') winter = (winter + prd) % mod2;
    }
    /*
    vector<ll> summer(201, 0);
    vector<ll> fall(201, 0);
    vector<ll> winter(201, 0);
    for (ll i = 0; i < N; i++){
        if (S[i] == 'U') summer[A[i]]++;
        if (S[i] == 'F') fall[A[i]]++;
        if (S[i] == 'W') winter[A[i]]++;
    }
    vector<ll> pow_1(202, 0);
    vector<ll> pow_2(202, 0);
    vector<ll> pow_3(606, 0);
    for (ll i = 0; i <= 200; i++){
        pow_1[i] = my_pow(N - 1, i, mod2);
        pow_2[i] = my_pow(N - 2, i, mod2);
    }
    for (ll i = 0; i <= 600; i++){
        if (i <= all_sum) pow_3[i] = my_pow(N - 1, all_sum - i, mod2);
    }
    ll ans = 0;
    for (ll i = 1; i <= 200; i++){
        for (ll j = 1; j <= 200; j++){
            for (ll k = 1; k <= 200; k++){
                if (summer[i] == 0 || fall[j] == 0 || winter[k] == 0) continue;
                ll prod = 1;
                prod = (prod * ((pow_1[i] - pow_2[i] + mod2) % mod2)) % mod2;
                prod = (prod * ((pow_1[j] - pow_2[j] + mod2) % mod2)) % mod2;
                prod = (prod * ((pow_1[k] - pow_2[k] + mod2) % mod2)) % mod2;
                prod = (prod * pow_3[i + j + k]) % mod2;
                prod = (prod * summer[i]) % mod2;
                prod = (prod * fall[j]) % mod2;
                prod = (prod * winter[k]) % mod2;
                prod = (prod * spring) % mod2;
                ans = (ans + prod) % mod2;
            }
        }
    }*/
    ll ans = 1;
    ans = (ans * summer) % mod2;
    ans = (ans * fall) % mod2;
    ans = (ans * winter) % mod2;
    ans = (ans * spring) % mod2;
    ans = (ans * my_pow(N - 1, all_sum, mod2)) % mod2;
    cout << ans << endl;
    
}
0