結果
問題 |
No.2127 Mod, Sum, Sum, Mod
|
ユーザー |
![]() |
提出日時 | 2022-11-19 12:38:49 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 43 ms / 2,000 ms |
コード長 | 1,969 bytes |
コンパイル時間 | 4,610 ms |
コンパイル使用メモリ | 252,996 KB |
最終ジャッジ日時 | 2025-02-08 22:43:05 |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 27 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/all> using namespace std; using namespace atcoder; istream &operator>>(istream &is, modint &a) { long long v; is >> v; a = v; return is; } ostream &operator<<(ostream &os, const modint &a) { return os << a.val(); } istream &operator>>(istream &is, modint998244353 &a) { long long v; is >> v; a = v; return is; } ostream &operator<<(ostream &os, const modint998244353 &a) { return os << a.val(); } istream &operator>>(istream &is, modint1000000007 &a) { long long v; is >> v; a = v; return is; } ostream &operator<<(ostream &os, const modint1000000007 &a) { return os << a.val(); } typedef long long ll; typedef vector<vector<int>> Graph; typedef pair<int, int> pii; typedef pair<ll, ll> pll; #define rep(i,n) for (int i = 0;i < (int)(n); i++) #define all(x) x.begin(), x.end() #define rall(x) x.rbegin(), x.rend() #define my_sort(x) sort(x.begin(), x.end()) #define my_max(x) *max_element(all(x)) #define my_min(x) *min_element(all(x)) #define cout_vector(x) rep(i,x.size()) cout << x[i] << ' '; template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } const int INF = (1<<30) - 1; const ll LINF = (1LL<<62) - 1; const int MOD = 998244353; const int MOD2 = 1e9+7; const double PI = acos(-1); vector<int> di = {1,0,-1,0}; vector<int> dj = {0,1,0,-1}; using mint = modint998244353; mint f1(mint x){ return x * (x+1) / 2; } mint f2(mint x){ return x * (x+1) * (2 * x + 1) / 6; } int main(){ cin.tie(0); ios_base::sync_with_stdio(false); int N,M; cin >> N >> M; mint ans = 1; ans *= N; ans *= (N+1); ans *= M; int i = M; while(i){ int v = N/i; int j = N/(v+1); ans -= mint(2 * N + 2) * mint(v) * (f1(i) - f1(j)); ans += mint(v+1) * mint(v) * (f2(i) - f2(j)); i = j; } ans /= 2; cout << ans << endl; }