結果

問題 No.2758 RDQ
ユーザー InTheBloomInTheBloom
提出日時 2024-05-17 22:01:34
言語 C++23(gcc13)
(gcc 13.2.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 2,342 bytes
コンパイル時間 1,439 ms
コンパイル使用メモリ 126,696 KB
実行使用メモリ 13,760 KB
最終ジャッジ日時 2024-05-17 22:01:40
合計ジャッジ時間 5,250 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
13,760 KB
testcase_01 AC 2 ms
6,944 KB
testcase_02 AC 2 ms
6,944 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 1 ms
6,940 KB
testcase_05 TLE -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <algorithm>
#include <tuple>
#include <cmath>

using namespace std;
using ll = long long;

int main () {
    int N, Q; cin >> N >> Q;
    vector<int> A(N);
    for (int i = 0; i < N; i++) cin >> A[i];

    // 賢いやり方不明。
    // ある程度大きなKに対してはMoを適用し、小さなKに対しては空間N√Nで計算しながらやる。
    vector<int> L(Q), R(Q), K(Q);
    for (int i = 0; i < Q; i++) {
        cin >> L[i] >> R[i] >> K[i];
        L[i]--;
        // 0-indexedで開区間
    }

    vector<int> index(Q), priority(Q);
    const int width = (int) sqrt((double) Q) + 10;

    for (int i = 0; i < Q; i++) {
        index[i] = i;
        priority[i] = L[i] / width;
    }

    sort(index.begin(), index.end(), [&](int x, int y) {
                if (priority[x] == priority[y]) {
                    if ((priority[x]) % 2 == 0) return R[x] < R[y];
                    return R[y] < R[x];
                }
                return priority[x] < priority[y];
            });

    const int max_mod = 100000;
    const int sqrt_N = (int) sqrt((double) N);
    vector<int> count_all(max_mod + 1);
    vector<vector<int>> count(sqrt_N);
    for (int i = 0; i < sqrt_N; i++) count[i] = vector<int>(i, 0);


    vector<int> ans(Q);

    int l = 0, r = 0;
    // [l, r)
    for (int i = 0; i < Q; i++) {
        int idx = index[i];

        while (l < L[idx]) {
            count_all[A[l]]--;
            for (int i = 1; i < sqrt_N; i++) count[i][A[l] % i]--;
            l++;
        }
        while (L[idx] < l) {
            l--;
            count_all[A[l]]++;
            for (int i = 1; i < sqrt_N; i++) count[i][A[l] % i]++;
        }

        while (r < R[idx]) {
            count_all[A[r]]++;
            for (int i = 1; i < sqrt_N; i++) count[i][A[r] % i]++;
            r++;
        }
        while (R[idx] < r) {
            r--;
            count_all[A[r]]--;
            for (int i = 1; i < sqrt_N; i++) count[i][A[r] % i]--;
        }

        if (K[idx] < sqrt_N) {
            ans[idx] = count[K[idx]][0];
        }
        else {
            int cur = K[idx];
            while (cur < max_mod) {
                ans[idx] += count_all[cur];
                cur += K[idx];
            }
        }
    }

    for (auto v : ans) cout << v << "\n";
}
0