結果
問題 | No.2940 Sigma Sigma Div Floor Problem |
ユーザー |
![]() |
提出日時 | 2024-10-18 22:45:49 |
言語 | C (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 2 ms / 6,000 ms |
コード長 | 1,101 bytes |
コンパイル時間 | 2,143 ms |
コンパイル使用メモリ | 30,848 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2024-10-18 22:54:35 |
合計ジャッジ時間 | 8,389 ms |
ジャッジサーバーID (参考情報) |
judge6 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 29 TLE * 1 |
ソースコード
#include<stdio.h> long long int modpow(long long int a, long long int n, long long int p) { long long int res = 1; for (; n > 0; n /= 2, a = a * a % p) if (n % 2 > 0) res = res * a % p; return res; } int main() { long long int n; scanf("%lld", &n); const long long int p = 998244353; long long int ans = 0; long long int i; long long int c, v; for (i = 1; i * i <= n; i++) { v = (n - i + 1) / i; ans += v % p * (v % p + 1) / 2 % p * i % p; ans %= p; ans += ((n - i + 1) - (i * v)) % p * ((v + 1) % p) % p; ans %= p; } long long int t = i; long long int min, mid, max; long long int inv2 = modpow(2, p - 2, p); for (i = t; i > 0; i--) { min = t - 1; max = n + 1; while (max - min > 1) { mid = (max + min) / 2; if (n / mid > i) min = mid; else max = mid; } ans += (min - t + 1) % p * (min + t) % p * inv2 % p * i % p; ans %= p; c = (n - min) - (max * (i - 1)); if (c > 0) { ans += ((c / i + 1) % p * c % p - (c / i) % p * (c / i + 1) % p * inv2 % p * i % p + p) % p * i % p; ans %= p; } } printf("%lld\n", ans); return 0; }