結果
問題 | No.752 mod数列 |
ユーザー |
![]() |
提出日時 | 2018-11-10 09:26:37 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 231 ms / 2,000 ms |
コード長 | 916 bytes |
コンパイル時間 | 1,038 ms |
コンパイル使用メモリ | 96,008 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-22 12:09:48 |
合計ジャッジ時間 | 7,040 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 31 |
ソースコード
#include <cstdio> #include <cstring> #include <iostream> #include <string> #include <cmath> #include <bitset> #include <vector> #include <map> #include <set> #include <queue> #include <deque> #include <algorithm> #include <complex> #include <unordered_map> #include <random> using namespace std; typedef long long int ll; typedef pair<int, int> P; ll p, m; ll ss[33000]; ll sl[33000]; ll solve(ll x){ if(x==0) return 0; if(x<=p/(m+1)) return ss[x]; ll ans=ss[p/(m+1)]; ll k=p/x; ans+=sl[k+1]; ans+=(p*x-k*x*(x+1)/2); ans-=(p*(p/(k+1))-k*(p/(k+1))*(p/(k+1)+1)/2); return ans; } int main() { cin>>p; m=(ll)sqrt((double)p); for(ll i=1; i<=p/(m+1); i++){ ss[i]=ss[i-1]+(p%i); } for(ll i=m; i>=1; i--){ ll s=p*(p/i-p/(i+1))-i*(p/i)*(p/i+1)/2+i*(p/(i+1))*(p/(i+1)+1)/2; sl[i]=sl[i+1]+s; } int q; cin>>q; for(int i=0; i<q; i++){ ll l, r; cin>>l>>r; cout<<solve(r)-solve(l-1)<<endl; } return 0; }