結果

問題 No.752 mod数列
コンテスト
ユーザー kakira9618
提出日時 2018-11-09 22:00:14
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
TLE  
実行時間 -
コード長 2,000 bytes
コンパイル時間 1,355 ms
コンパイル使用メモリ 163,956 KB
実行使用メモリ 16,256 KB
最終ジャッジ日時 2024-11-21 05:57:01
合計ジャッジ時間 66,354 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 10 TLE * 21
権限があれば一括ダウンロードができます

ソースコード

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

ll solve(ll P, ll l, ll r) {
    set<ll> D;
    for(ll n = 1; n * n <= P; n++) {
        D.insert(-n);
        D.insert(-(P / n));
    }

    ll ans = (r - l + 1) * P;
    for(auto d : D) {
        int ll, rr;
        d *= -1;

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

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


        ll = max(l, 1LL * ll);
        rr = min(r, 1LL * rr);
        if (ll > rr) continue;

        ans -= d * (ll + rr) * (rr - ll + 1) / 2;
    }
    return ans;
}

void solve() {
    int P, Q;
    cin >> P >> Q;
    for (int q = 0; q < Q; q++) {
        int l, r;
        cin >> l >> r;
        cout << solve(P, l, r) << endl;
    }
}

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