結果

問題 No.1097 Remainder Operation
ユーザー outlineoutline
提出日時 2020-12-14 10:11:21
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 88 ms / 2,000 ms
コード長 2,580 bytes
コンパイル時間 1,008 ms
コンパイル使用メモリ 87,632 KB
実行使用メモリ 25,376 KB
最終ジャッジ日時 2023-10-20 04:54:07
合計ジャッジ時間 3,930 ms
ジャッジサーバーID
(参考情報)
judge13 / judge11
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,348 KB
testcase_01 AC 2 ms
4,348 KB
testcase_02 AC 2 ms
4,348 KB
testcase_03 AC 2 ms
4,348 KB
testcase_04 AC 2 ms
4,348 KB
testcase_05 AC 2 ms
4,348 KB
testcase_06 AC 2 ms
4,348 KB
testcase_07 AC 9 ms
5,016 KB
testcase_08 AC 8 ms
5,016 KB
testcase_09 AC 9 ms
5,016 KB
testcase_10 AC 8 ms
5,016 KB
testcase_11 AC 8 ms
5,016 KB
testcase_12 AC 85 ms
20,492 KB
testcase_13 AC 87 ms
20,492 KB
testcase_14 AC 87 ms
20,488 KB
testcase_15 AC 88 ms
20,500 KB
testcase_16 AC 88 ms
20,496 KB
testcase_17 AC 65 ms
21,628 KB
testcase_18 AC 57 ms
20,300 KB
testcase_19 AC 65 ms
25,376 KB
testcase_20 AC 65 ms
25,376 KB
testcase_21 AC 61 ms
23,080 KB
testcase_22 AC 61 ms
23,080 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;

using ll = long long;

int main(){
    ios::sync_with_stdio(false);
    cin.tie(nullptr);

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

    // 強連結成分分解で周期性を管理
    vector<vector<int>> graph(N), rev(N);
    for(int i = 0; i < N; ++i){
        int j = (i + A[i]) % N;
        graph[i].emplace_back(j);
        rev[j].emplace_back(i);
    }

    vector<int> vs, used(N);
    auto dfs = [&](auto&& self, int from) -> void {
        used[from] = true;
        for(int to : graph[from]){
            if(!used[to])   self(self, to);
        }
        vs.emplace_back(from);
    };
    for(int v = 0; v < N; ++v){
        if(!used[v])    dfs(dfs, v);
    }
    reverse(vs.begin(), vs.end());

    int scc_num = 0;
    vector<int> scc_id(N, -1);
    auto rdfs = [&](auto&& self, int from) -> void {
        scc_id[from] = scc_num;
        for(int to : rev[from]){
            if(scc_id[to] == -1)    self(self, to);
        }
    };
    for(int v : vs){
        if(scc_id[v] == -1){
            rdfs(rdfs, v);
            scc_num += 1;
        }
    }

    vector<vector<int>> scc_group(scc_num);
    for(int v = 0; v < N; ++v){
        scc_group[scc_id[v]].emplace_back(v);
    }

    // ループに入るまでの累積和をとる
    // 頂点数 1 のループに注意!
    int u = 0;
    vector<ll> path_sum(1, 0);
    while(u != graph[u][0] && scc_group[scc_id[u]].size() == 1){
        path_sum.emplace_back(path_sum.back() + A[u]);
        u = graph[u][0];
    }

    // ループの 1 周期分の累積和をとる
    vector<ll> loop_sum(1, 0);
    loop_sum.emplace_back(A[u]);
    int cur = graph[u][0];
    while(cur != u){
        loop_sum.emplace_back(loop_sum.back() + A[cur]);
        cur = graph[cur][0];
    }

    cin >> Q;
    for(int i = 0; i < Q; ++i){
        ll K;
        cin >> K;
        // ループに到達しない場合
        if(K < (ll)path_sum.size()){
            cout << path_sum[K] << '\n';
            continue;
        }

        // ループに到達させる
        int v = u;
        ll ans = path_sum.back();
        K -= (ll)path_sum.size() - 1;

        // 何回ループできるかを計算して和を ans に足し込む
        int sz = scc_group[scc_id[v]].size();
        ll loop_cnt = K / sz;
        ans += loop_sum.back() * loop_cnt;

        // あまり分を ans に足し込む
        ans += loop_sum[K % sz];

        cout << ans << '\n';
    }

    return 0;
}
0