結果

問題 No.752 mod数列
コンテスト
ユーザー 梧桐
提出日時 2026-01-24 17:27:42
言語 C++17
(gcc 15.2.0 + boost 1.89.0)
結果
WA  
実行時間 -
コード長 2,880 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 925 ms
コンパイル使用メモリ 100,016 KB
実行使用メモリ 7,852 KB
最終ジャッジ日時 2026-01-24 17:27:47
合計ジャッジ時間 4,873 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2 WA * 1
other AC * 6 WA * 25
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <iostream>
#include <algorithm>
#include <vector>

using namespace std;

typedef long long LL;
typedef pair<int, int> PII;

const int N = 100010;

struct Block {
    int l, r;
    LL frst_ele, c, sum;
};

int p, q;
PII pp[N];
vector<Block> blocks;
LL ps[N];
vector<int> ls, rs;

LL Calc1(int id, int l, int r) {
    int ele = blocks[id].frst_ele, c = blocks[id].c, cl = blocks[id].l;
    int l_val = ele - c * (l - cl);
    int r_val = ele - c * (r - cl);
    return 1LL * (l_val + r_val) * (r - l + 1) / 2;
}

LL Calc2(int id, int pos, int op) {
    int ele = blocks[id].frst_ele, c = blocks[id].c, cr = blocks[id].r, cl = blocks[id].l;
    if (op == 0) { // 左块
        int l_val = ele - c * (pos - cl);
        int r_val = ele - c * (cr - cl);
        if (c == 0) return 1LL * l_val * (cr - pos + 1);
        return 1LL * (l_val + r_val) * (cr - pos + 1) / 2;
    } else {
        int l_val = ele;
        int r_val = ele - c * (pos - cl);
        if (c == 0) return 1LL * l_val * (pos - cl + 1);
        return 1LL * (l_val + r_val) * (pos - cl + 1) / 2;
    }
}

LL Calc3(int lb, int rb) {
    return ps[rb] - ps[lb - 1];
}

int main() {
    // freopen("mod.in", "r", stdin);
    // freopen("mod.out", "w", stdout);

    scanf("%d%d", &p, &q);
    for (int i = 1; i <= q; ++i) scanf("%d%d", &pp[i].first, &pp[i].second);

    int cl = 1, cr = 1;
    for (int i = 0; cr < p; ++i) {
        if (i != 0) {
            cl = cr + 1;
            cr = p / (p / cl);
        }
        int f_ele = p % cl;
        int l_ele = p % cr;
        if (cl != cr) {
            blocks.push_back({ cl, cr, f_ele, (p % cl) - (p % (cl + 1)), (f_ele + l_ele) * (cr - cl + 1) / 2 });
        } else {
            blocks.push_back({ cl, cr, l_ele, 0, (f_ele + l_ele) * (cr - cl + 1) / 2 });
        }
        ls.push_back(cl), rs.push_back(cr);
    }
    ps[0] = blocks[0].sum;
    for (int i = 1; i < blocks.size(); ++i) {
        ps[i] = ps[i - 1] + blocks[i].sum;
    }

    for (int i = 1; i <= q; ++i) {
        LL ans = 0LL;
        if (pp[i].second > p) {
            ans += 1LL * (pp[i].second - max(pp[i].first, p + 1) + 1) * p;
            pp[i].second = min(pp[i].second, p);
        }
        if (pp[i].first <= p) {
            int lb = upper_bound(ls.begin(), ls.end(), pp[i].first) - ls.begin() - 1;
            int rb = lower_bound(rs.begin(), rs.end(), pp[i].second) - rs.begin();
            if (lb == rb) {
                ans += Calc1(lb, pp[i].first, pp[i].second); // 一个块内
            } else {
                ans += Calc2(lb, pp[i].first, 0); // 最左块（部分）
                if (rb >= lb + 2) {
                    ans += Calc3(lb + 1, rb - 1); // 中间块（完整）
                }
                ans += Calc2(rb, pp[i].second, 1); // 最右块（部分）
            }
        }
        printf("%lld\n", ans);
    }
    
    return 0;
}
0