結果
問題 |
No.2940 Sigma Sigma Div Floor Problem
|
ユーザー |
![]() |
提出日時 | 2024-10-18 23:06:56 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 3 ms / 6,000 ms |
コード長 | 1,360 bytes |
コンパイル時間 | 4,284 ms |
コンパイル使用メモリ | 252,728 KB |
最終ジャッジ日時 | 2025-02-24 21:14:25 |
ジャッジサーバーID (参考情報) |
judge1 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 30 |
ソースコード
#include <stdio.h> #include <atcoder/all> #include <bits/stdc++.h> using namespace std; using namespace atcoder; using mint = modint998244353; #define rep(i,n) for (int i = 0; i < (n); ++i) #define Inf32 1000000001 #define Inf64 1000000000000000000LL mint ans= 0; mint i2; mint sum(long long n){ mint ret = 0; ret = n; ret *= n+1; ret *= i2; return ret; } mint sum(long long l,long long r){ if(l==0)return sum(r); return sum(r)-sum(l-1); } mint i6; mint f2(long long n){ mint ret = n; ret *= n+1; ret *= n*2+1; ret *= i6; return ret; } mint sum2(long long l,long long r){ if(l==0)return f2(r); return f2(r)-f2(l-1); } mint Get(long long n,long long L,long long R,long long v){ mint ans = 0; { mint tt = n; tt *= R-L+1; tt *= v; ans += tt; } { mint tt = sum(L-1,R-1); tt *= v; ans -= tt; } { mint tt = sum(v-1); tt *= sum(L,R); ans -= tt; } return ans; } //[[l,r],v] l<=x<=rは, floor(N/x)=v vector<pair<pair<long long,long long>,long long>> get_floor(long long N){ //vector<pair<pair<long long,long long>,long long>> ret; for(long long i=1;i<=N;i++){ long long v = N/i; long long L = i; i = N/v; long long R = i; ans += Get(N,L,R,v); } return {}; } int main(){ i2 = mint(2).inv(); i6 = mint(6).inv(); long long n; cin>>n; get_floor(n); cout<<ans.val()<<endl; return 0; }