結果

問題 No.752 mod数列
ユーザー SHIJOU
提出日時 2020-08-20 17:49:32
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 124 ms / 2,000 ms
コード長 1,654 bytes
コンパイル時間 2,323 ms
コンパイル使用メモリ 201,220 KB
最終ジャッジ日時 2025-01-13 04:27:11
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff # 

//#define _GLIBCXX_DEBUG
#include <bits/stdc++.h>
#define rep(i, n) for(int i=0; i<n; ++i)
#define all(v) v.begin(), v.end()
#define rall(v) v.rbegin(), v.rend()
using namespace std;
using ll = int64_t;
using ld = long double;
using P = pair<int, int>;
using vs = vector<string>;
using vi = vector<int>;
using vvi = vector<vi>;
template<class T> using PQ = priority_queue<T>;
template<class T> using PQG = priority_queue<T, vector<T>, greater<T> >;
const int INF = 0xccccccc;
const ll LINF = 922337203685477580LL;
template<typename T1, typename T2>
inline bool chmax(T1 &a, T2 b) {return a < b && (a = b, true);}
template<typename T1, typename T2>
inline bool chmin(T1 &a, T2 b) {return a > b && (a = b, true);}
template<typename T1, typename T2>
istream &operator>>(istream &is, pair<T1, T2> &p) { return is >> p.first >> p.second;}
template<typename T1, typename T2>
ostream &operator<<(ostream &os, const pair<T1, T2> &p) { return os << p.first << ' ' << p.second;}



//head

int p, q, m, u;
map<int, ll> mp;

int main() {
	ios::sync_with_stdio(false);
	cin.tie(0);
	cin >> p >> q;
	m = sqrt(p);
	u = p/m;
	for(int i = 1; i <= u; i++) {
		mp[i] = mp[i-1] + p%i;
	}
	for(int j = m-1; j > 0; j--) {
		auto x = *mp.rbegin();
		int y = p/j;
		mp[y] = x.second + (ll)(y-x.first)*(p%y + p%(x.first+1))/2;
	}

	while(q--) {
		int lr[2];
		cin >> lr[0] >> lr[1];
		lr[0]--;
		ll box[2];
		rep(u, 2) {
			map<int, ll>::iterator lb = mp.lower_bound(lr[u]);
			if(lb->first == lr[u]) box[u] = lb->second;
			else {
				lb = prev(lb);
				box[u] = lb->second + (ll)(lr[u]-lb->first)*(p%lr[u] + p%(lb->first+1))/2;
			}
		}
		cout << box[1] - box[0] << '\n';
	}
}
0