結果

問題 No.1304 あなたは基本が何か知っていますか？私は知っています．
ユーザー mfbgjsczmfbgjscz
提出日時 2021-02-05 18:12:01
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 261 ms / 2,000 ms
コード長 2,105 bytes
コンパイル時間 1,803 ms
コンパイル使用メモリ 175,504 KB
実行使用メモリ 21,592 KB
最終ジャッジ日時 2023-09-03 04:13:15
合計ジャッジ時間 37,052 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 113 ms
21,456 KB
testcase_01 AC 152 ms
21,356 KB
testcase_02 AC 114 ms
21,464 KB
testcase_03 AC 261 ms
21,388 KB
testcase_04 AC 112 ms
21,476 KB
testcase_05 AC 113 ms
21,460 KB
testcase_06 AC 113 ms
21,336 KB
testcase_07 AC 112 ms
21,460 KB
testcase_08 AC 187 ms
21,476 KB
testcase_09 AC 112 ms
21,444 KB
testcase_10 AC 114 ms
21,524 KB
testcase_11 AC 112 ms
21,460 KB
testcase_12 AC 113 ms
21,468 KB
testcase_13 AC 114 ms
21,404 KB
testcase_14 AC 113 ms
21,332 KB
testcase_15 AC 113 ms
21,524 KB
testcase_16 AC 113 ms
21,464 KB
testcase_17 AC 153 ms
21,456 KB
testcase_18 AC 113 ms
21,340 KB
testcase_19 AC 113 ms
21,336 KB
testcase_20 AC 114 ms
21,508 KB
testcase_21 AC 153 ms
21,572 KB
testcase_22 AC 113 ms
21,456 KB
testcase_23 AC 112 ms
21,436 KB
testcase_24 AC 114 ms
21,352 KB
testcase_25 AC 226 ms
21,456 KB
testcase_26 AC 191 ms
21,340 KB
testcase_27 AC 188 ms
21,392 KB
testcase_28 AC 197 ms
21,348 KB
testcase_29 AC 189 ms
21,348 KB
testcase_30 AC 218 ms
21,348 KB
testcase_31 AC 229 ms
21,592 KB
testcase_32 AC 188 ms
21,460 KB
testcase_33 AC 227 ms
21,456 KB
testcase_34 AC 188 ms
21,460 KB
testcase_35 AC 217 ms
21,500 KB
testcase_36 AC 226 ms
21,392 KB
testcase_37 AC 221 ms
21,444 KB
testcase_38 AC 219 ms
21,460 KB
testcase_39 AC 248 ms
21,340 KB
testcase_40 AC 220 ms
21,356 KB
testcase_41 AC 247 ms
21,444 KB
testcase_42 AC 217 ms
21,400 KB
testcase_43 AC 218 ms
21,436 KB
testcase_44 AC 225 ms
21,464 KB
04_evil_A_01 AC 265 ms
21,456 KB
04_evil_A_02 AC 229 ms
21,336 KB
04_evil_A_03 AC 250 ms
21,340 KB
04_evil_A_04 AC 224 ms
21,432 KB
04_evil_A_05 AC 226 ms
21,456 KB
04_evil_A_06 AC 117 ms
21,436 KB
04_evil_A_07 AC 114 ms
21,460 KB
04_evil_A_08 AC 117 ms
21,460 KB
04_evil_A_09 AC 115 ms
21,460 KB
04_evil_A_10 AC 116 ms
21,436 KB
05_evil_B_01 AC 292 ms
21,464 KB
05_evil_B_02 AC 283 ms
21,456 KB
05_evil_B_03 AC 288 ms
21,400 KB
05_evil_B_04 AC 290 ms
21,336 KB
05_evil_B_05 AC 285 ms
21,460 KB
05_evil_B_06 AC 284 ms
21,400 KB
05_evil_B_07 AC 284 ms
21,340 KB
05_evil_B_08 AC 283 ms
21,456 KB
05_evil_B_09 AC 296 ms
21,464 KB
05_evil_B_10 AC 298 ms
21,344 KB
06_evil_C_01 TLE -
06_evil_C_02 TLE -
06_evil_C_03 TLE -
06_evil_C_04 TLE -
06_evil_C_05 TLE -
06_evil_C_06 TLE -
06_evil_C_07 TLE -
06_evil_C_08 TLE -
06_evil_C_09 TLE -
06_evil_C_10 AC 76 ms
21,548 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#pragma GCC target("avx")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
using namespace std;
#define ll long long
#define rep(i, m, n) for (ll i = m; i < n; i++)
#define mod 998244353

void fwt(vector<ll>& f) {
    int n = f.size();
    for(int i = 1; i < n; i <<= 1) {
        for(int j = 0; j < n; j++) {
            if((j & i) == 0) {
                ll x = f[j], y = f[j | i];
                f[j] = x + y, f[j | i] = x - y;
            }
        }
    }
}

void ifwt(vector<ll>& f) {
    int n = f.size();
    for(int i = 1; i < n; i <<= 1) {
        for(int j = 0; j < n; j++) {
            if((j & i) == 0) {
                ll x = f[j], y = f[j | i];
                f[j] = (x + y) * 499122177 % mod, f[j | i] = (x - y + mod) % mod * 499122177 % mod;
            }
        }
    }
}

vector<ll> mul(vector<ll> &mt1, vector<vector<ll>> &mt2) {
    vector<ll> res(2);
    res[0] = (mt1[0] * mt2[0][0] % mod + mt1[1] * mt2[0][1] % mod) % mod;
    res[1] = (mt1[0] * mt2[1][0] % mod + mt1[1] * mt2[1][1] % mod) % mod;
    return res;
}

vector<vector<ll>> update(vector<vector<ll>> &mt) {
    vector<vector<ll>> res(2, vector<ll>(2, 0));
    rep(i, 0, 2) rep(j, 0, 2) rep(k, 0, 2) res[i][j] = (res[i][j] + mt[i][k] * mt[k][j] % mod) % mod;
    return res; 
}

int main(){
	cin.tie(nullptr);
	ios::sync_with_stdio(false);
	ll N, K, X, Y; 
    cin >> N >> K >> X >> Y;
	ll M = 262144;
	vector<vector<ll>> dp(M, vector<ll>(2));
    vector<ll> H(M, 0);
    rep(i, 0, K) {
		ll A;
		cin >> A;
		H[A] = 1;
	}
	fwt(H);
    rep(v, 0, M) dp[v][1] = H[v];
    rep(v, 0, M) dp[v][0] = (dp[v][1] * H[v] - K + mod) % mod;
	rep(v, 0, M) {
        vector<vector<ll>> mt(2, vector<ll>(2));
        mt[0] = {H[v], (1 - K + mod) % mod};
        mt[1] = {1, 0};
        ll P = N - 1;
        while(P) {
            if(P & 1) dp[v] = mul(dp[v], mt);
            mt = update(mt);
            P /= 2;
        }
    }
    vector<ll> ep(M);
    rep(v, 0, M) ep[v] = dp[v][1];
	ifwt(ep);
	ll res = 0;
	rep(v, X, Y + 1) res = (res + ep[v]) % mod;
	cout << (res + mod) % mod << endl;
}
0