結果

問題 No.2395 区間二次変換一点取得
ユーザー kwm_tkwm_t
提出日時 2023-07-28 22:02:51
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 90 ms / 2,000 ms
コード長 2,437 bytes
コンパイル時間 1,894 ms
コンパイル使用メモリ 196,884 KB
最終ジャッジ日時 2025-02-15 20:25:01
ジャッジサーバーID
(参考情報)
judge3 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

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

#include <bits/stdc++.h>
//#include <atcoder/all>
using namespace std;
//using namespace atcoder;
//using mint = modint1000000007;
//const int mod = 1000000007;
//using mint = modint998244353;
//const int mod = 998244353;
//const int INF = 1e9;
//const long long LINF = 1e18;
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define rep2(i,l,r)for(int i=(l);i<(r);++i)
#define rrep(i, n) for (int i = (n-1); i >= 0; --i)
#define rrep2(i,l,r)for(int i=(r-1);i>=(l);--i)
#define all(x) (x).begin(),(x).end()
#define allR(x) (x).rbegin(),(x).rend()
#define endl "\n"
#define P pair<int,int>
template<typename A, typename B> inline bool chmax(A & a, const B & b) { if (a < b) { a = b; return true; } return false; }
template<typename A, typename B> inline bool chmin(A & a, const B & b) { if (a > b) { a = b; return true; } return false; }
template<class S, S op(S l, S r), S e()>
struct DualSegtree {
private:
int mn = 0;
int mN = 0;
int mlogN = 0;
std::vector<S> A;
void propagate(int i) {
if (i >= mN) return;
A[i * 2] = op(A[i * 2], A[i]);
A[i * 2 + 1] = op(A[i * 2 + 1], A[i]);
A[i] = e();
}
void full_propagate(int i) {
for (int d = mlogN; d > 0; d--) propagate(i >> d);
}
public:
DualSegtree(int n = 0) {
mn = n;
mlogN = 0; mN = 1;
while (mN < n) { mN *= 2; mlogN++; }
A.assign(mN * 2, e());
}
void point_init(int p) {
assert(0 <= p); assert(p < mn);
p += mN;
full_propagate(p);
A[p] = e();
}
S get(int p) {
assert(0 <= p); assert(p < mn);
p += mN;
full_propagate(p);
return A[p];
}
void apply(int l, int r, S x) {
l += mN;
r += mN;
if (l != mN) full_propagate(l - 1);
if (r != mN * 2) full_propagate(r);
while (l < r) {
if (l & 1) { A[l] = op(A[l], x); l++; }
if (r & 1) { r--; A[r] = op(A[r], x); }
l /= 2; r /= 2;
}
}
};
//
using S = long long;
S op(S x, S f) { return { x + f }; }
S e() { return { 0 }; }
int main() {
ios::sync_with_stdio(false);
cin.tie(nullptr);
int n, b, q; cin >> n >> b >> q;
DualSegtree<S, op, e>seg(n);
vector<long long>x(q + 1), y(q + 1), z(q + 1);
x[0] = 1, y[0] = 1, z[0] = 1;
rep(i, q) {
x[i + 1] = (x[i] + 1) % b;;
y[i + 1] = (3LL * y[i] + 2LL * x[i + 1] * z[i]) % b;
z[i + 1] = (3LL * z[i]) % b;
}
rep(i, q) {
int l, m, r; cin >> l >> m >> r;
l--, m--, r--;
seg.apply(l, r + 1, 1);
auto val = seg.get(m);
cout << x[val] << " " << y[val] << " " << z[val] << endl;
}
return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0