結果

問題 No.1548 [Cherry 2nd Tune B] 貴方と私とサイクルとモーメント
ユーザー りあん
提出日時 2021-06-11 22:54:37
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 204 ms / 4,500 ms
コード長 2,617 bytes
コンパイル時間 2,499 ms
コンパイル使用メモリ 214,920 KB
最終ジャッジ日時 2025-01-22 06:43:21
ジャッジサーバーID
(参考情報)
judge5 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 42
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include<bits/stdc++.h>
#include<atcoder/modint>
#include<atcoder/lazysegtree>
using namespace std;
using mint = atcoder::modint998244353;
using S = pair<mint, int>;
S op(S a, S b) {
return { a.first + b.first, a.second + b.second };
}
S e() {
return S(0, 0);
}
S mapping(int a, S b) {
return a == -1 ? b : S(mint(a) * b.second, b.second);
}
int composition(int a, int b) {
return a == -1 ? b : a;
}
int id() {
return -1;
}
vector<atcoder::lazy_segtree<S, op, e, int, mapping, composition, id>> sg(5);
mint prod(int type, const vector<pair<int, int>>& segs) {
mint res(0);
for (auto& s : segs) {
res += sg[type].prod(s.first, s.second).first;
}
return res;
}
mint calc(int type, const vector<pair<int, int>>& segs) {
mint s0 = prod(0, segs);
mint s1 = prod(1, segs);
mint m = s1 / s0;
if (type == 1) {
return (s1 - s0 * m) / s0;
}
else if (type == 2) {
mint s2 = prod(2, segs);
return (s2 - s1 * m * 2 + s0 * m * m) / s0;
}
else if (type == 3) {
mint s2 = prod(2, segs);
mint s3 = prod(3, segs);
return (s3 - s2 * m * 3 + s1 * m * m * 3 - s0 * m * m * m) / s0;
}
else {
mint s2 = prod(2, segs);
mint s3 = prod(3, segs);
mint s4 = prod(4, segs);
return (s4 - s3 * m * 4 + s2 * m * m * 6 - s1 * m * m * m * 4 + s0 * m * m * m * m) / s0;
}
}
int main() {
cin.tie(0);
ios::sync_with_stdio(0);
int n;
cin >> n;
vector<int> a(n);
for (int i = 0; i < n; ++i) {
cin >> a[i];
}
vector<S> a2(n, {1, 1});
for (int i = 0; i < 5; ++i) {
sg[i] = atcoder::lazy_segtree<S, op, e, int, mapping, composition, id>(a2);
for (int j = 0; j < n; ++j) {
a2[j].first *= a[j];
}
}
int q;
cin >> q;
for (int _ = 0; _ < q; ++_) {
int t, u, v, w;
cin >> t >> u >> v >> w;
--u;
--v;
--w;
vector<pair<int, int>> segs;
if (u > v) swap(u, v);
if (u < w && w < v) {
segs.emplace_back(u, v + 1);
}
else {
segs.emplace_back(v, n);
segs.emplace_back(0, u + 1);
}
if (t == 0) {
int b;
cin >> b;
mint mb(1);
for (int i = 0; i < 5; ++i) {
for (auto& s : segs) {
sg[i].apply(s.first, s.second, mb.val());
}
mb *= b;
}
}
else {
cout << calc(t, segs).val() << '\n';
}
}
return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0