結果
| 問題 |
No.1897 Sum of 2nd Max
|
| コンテスト | |
| ユーザー |
miscalc
|
| 提出日時 | 2022-02-14 12:07:48 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 202 ms / 2,000 ms |
| コード長 | 843 bytes |
| コンパイル時間 | 2,196 ms |
| コンパイル使用メモリ | 198,124 KB |
| 最終ジャッジ日時 | 2025-01-27 23:02:47 |
|
ジャッジサーバーID (参考情報) |
judge3 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 30 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
#include <atcoder/modint>
using namespace atcoder;
using mint = modint998244353;
mint sum0(mint a, int n) // sum[i = 0 .. n - 1] a^i
{
if (a == 1)
return n;
return (1 - a.pow(n)) / (1 - a);
}
mint sum1(mint a, int n) // sum[i = 0 .. n - 1] ia^i
{
if (a == 1)
return mint(n) * (n - 1) / 2;
return (sum0(a, n) - 1 - (n - 1) * a.pow(n)) / (1 - a);
}
mint P(mint x, int n, int k)
{
mint y = (x - 1) / x;
mint S = sum0(y, n - 1), T = sum1(y, n - 1);
mint res1 = x.pow(n - 2) * (k - x) * ((n - 1) * S - T); // i < j
mint res2 = x.pow(n - 2) * (k - x + 1) * (S + T); // i > j
return res1 + res2;
}
int main()
{
int N, K;
cin >> N >> K;
mint ans = 0;
for (int x = 1; x <= K; x++)
{
mint tmp = P(x, N, K);
ans += x * tmp;
}
cout << ans.val() << endl;
}
miscalc