結果

問題 No.1838 Modulo Straight
コンテスト
ユーザー hitonanode
提出日時 2021-12-19 17:13:38
言語 C++17(gcc12)
(gcc 12.4.0 + boost 1.90.0)
コンパイル:
g++-12 -O2 -lm -std=c++17 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 219 ms / 2,000 ms
コード長 1,568 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 5,417 ms
コンパイル使用メモリ 267,224 KB
実行使用メモリ 45,828 KB
最終ジャッジ日時 2026-06-24 20:16:21
合計ジャッジ時間 12,865 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 38
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

// Validator
#include <algorithm>
#include <cassert>
#include <iostream>
#include <vector>
using namespace std;
#include "testlib.h"

#include <atcoder/fenwicktree>

int main(int argc, char **argv) {

    registerValidation(argc, argv);
    int M = inf.readInt(2, 10000000);
    inf.readSpace();
    int K = inf.readInt(1, 400000 / M);
    inf.readEoln();

    const int N = M * K;
    vector<int> A(N);
    for (int i = 0; i < N; ++i) {
        A[i] = inf.readInt(0, M - 1);
        if (i + 1 < N) {
            inf.readSpace();
        } else {
            inf.readEoln();
        }
    }
    inf.readEof();

    vector<int> ord(N);
    vector<int> freq(M);
    vector<vector<int>> v2i(M);
    for (int i = 0; i < N; ++i) {
        cin >> A[i];
        ord[i] = A[i] + freq[A[i]] * M;
        v2i[A[i]].push_back(i);
        freq[A[i]]++;
    }
    assert(freq == vector<int>(M, K));

    long long current = 0;
    atcoder::fenwick_tree<int> bit(N);
    for (int i = 0; i < N; ++i) {
        current += bit.sum(ord[i], N);
        bit.add(ord[i], 1);
    }
    long long ret = current;
    vector<vector<int>> prvs(K);
    for (int i = 0; i < N; ++i) prvs.at(ord[i] / M).push_back(i);
    for (auto &v : prvs) sort(v.begin(), v.end());

    for (auto &is : v2i) {
        for (int t = 0; t < K; ++t) {
            int i = is.at(t);
            int nlo = lower_bound(prvs[t].begin(), prvs[t].end(), i) - prvs[t].begin();
            int nhi = M - 1 - nlo;
            current += nhi - nlo;
        }
        ret = min(ret, current);
    }
    cout << ret << '\n';
}
0