結果

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

ソースコード

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, frst_ele, c;
    LL 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;
    while (cl <= p) {
        int d = p / cl; // 商（公差）
        int cr = p / d;
        int first = p % cl;
        int last  = p % cr;
        int c = (cl == cr ? 0 : d);
        LL sum = 1LL * (first + last) * (cr - cl + 1) / 2;
        blocks.push_back({ cl, cr, first, c, sum });
        ls.push_back(cl);
        rs.push_back(cr);
        cl = cr + 1;
    }
    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); // 最左块（部分）
                ans += Calc3(lb + 1, rb - 1); // 中间块（完整）
                ans += Calc2(rb, pp[i].second, 1); // 最右块（部分）
            }
        }
        printf("%lld\n", ans);
    }
    
    return 0;
}
0