結果

問題 No.752 mod数列
ユーザー kakira9618kakira9618
提出日時 2018-11-09 23:11:44
言語 C++11
(gcc 11.4.0)
結果
WA  
実行時間 -
コード長 2,900 bytes
コンパイル時間 1,986 ms
コンパイル使用メモリ 169,088 KB
実行使用メモリ 13,768 KB
最終ジャッジ日時 2024-11-21 06:44:40
合計ジャッジ時間 44,483 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 5 ms
13,636 KB
testcase_01 AC 5 ms
12,836 KB
testcase_02 AC 4 ms
13,636 KB
testcase_03 AC 3 ms
13,088 KB
testcase_04 AC 2 ms
10,276 KB
testcase_05 AC 5 ms
10,276 KB
testcase_06 AC 3 ms
13,636 KB
testcase_07 AC 5 ms
10,148 KB
testcase_08 AC 5 ms
10,144 KB
testcase_09 AC 4 ms
10,144 KB
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 AC 252 ms
12,708 KB
testcase_16 TLE -
testcase_17 TLE -
testcase_18 AC 1,074 ms
13,768 KB
testcase_19 TLE -
testcase_20 TLE -
testcase_21 TLE -
testcase_22 TLE -
testcase_23 TLE -
testcase_24 TLE -
testcase_25 TLE -
testcase_26 TLE -
testcase_27 TLE -
testcase_28 TLE -
testcase_29 TLE -
testcase_30 WA -
testcase_31 AC 2 ms
6,816 KB
testcase_32 AC 2 ms
6,816 KB
testcase_33 AC 78 ms
11,552 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;
#define pb push_back
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define rep(i,n) for(int i=0;i<(n);i++)
constexpr int MOD = 1000000007;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
constexpr int dx[] = {1, 0, -1, 0, 1, 1, -1, -1};
constexpr int dy[] = {0, -1, 0, 1, 1, -1, -1, 1};

template <typename T> ostream &operator<<(ostream &os, const vector<T> &vec){os << "["; for (const auto &v : vec) {os << v << ","; } os << "]"; return os; }
template <typename T, typename U> ostream &operator<<(ostream &os, const pair<T, U> &p){os << "(" << p.first << ", " << p.second << ")"; return os;}

// \sum_{n = L_i ^ R_i} P - [P / n] * n
// (R_i - L_i * 1) * P - \sum_{n = L_i ^ R_i} [P / n] * n
set<ll> D;
set<ll>::iterator it;

ll solve(ll P, ll l, ll r) {

    ll ans = (r - l + 1) * P;
    for(; it != D.end(); it++) {
        ll l2, r2;
        ll d = *it;
        d *= -1;

        ll ng = P + 1;
        ll ok = 0;
        while(ng - ok > 1) {
            ll c = (ok + ng) / 2;
            if (P / c >= d) {
                ok = c;
            } else {
                ng = c;
            }
        }
        r2 = ok;

        ng = P + 1;
        ok = 0;
        while(ng - ok > 1) {
            ll c = (ok + ng) / 2;
            if (P / c > d) {
                ok = c;
            } else {
                ng = c;
            }
        }
        l2 = ng;

        if (l2 > r) break;
        l2 = max(l, l2);
        r2 = min(r, r2);
        if (l2 > r2) continue;

        ans -= d * (l2 + r2) * (r2 - l2 + 1) / 2;
    }
    if (it != D.begin()) it--;
    return ans;
}

ll solveNaive(ll P, ll l, ll r) {
    ll ret = 0;
    for(ll i = l; i <= r; i++) {
        ret += P % i;
    }
    return ret;
}

void solve() {
    ll P, Q;
    cin >> P >> Q;

    for(ll n = 1; n * n <= P; n++) {
        D.insert(-n);
        D.insert(-(P / n));
    }
    it = D.begin();

    ll d = (ll)sqrt(P);
    vector<ll> B((P + 1 + d - 1) / d);
    for (int i = 0; i < B.size(); i++) {
        ll l = i * d;
        ll r = (i + 1) * d - 1;
        B[i] = solve(P, max(1LL, l), r);
    }

    for (int q = 0; q < Q; q++) {
        ll l, r;
        cin >> l >> r;
        ll ans = 0;
        if (r > P) {
            ans += (r - P) * P;
            r = P;
        }

        if (r - l <= 3 * d) {
            cout << solveNaive(P, l, r) << endl;
            continue;
        }
        
        for(; l % d != 0; l++) {
            ans += P % l;
        }
        for(; r % d != 0; r--) {
            ans += P % r;
        }
        ans += P % r;


        for(int b = l / d; b < r / d; b++) {
            ans += B[b];
        }
        cout << ans << endl;
    }
    
}

int main() {
    std::cin.tie(0);
    std::ios::sync_with_stdio(false);
    cout.setf(ios::fixed);
    cout.precision(16);
    solve();
    return 0;
}
0