This documentation is automatically generated by online-judge-tools/verification-helper
View the Project on GitHub beet-aizu/library
// verification-helper: PROBLEM https://judge.yosupo.jp/problem/point_set_range_composite #include<bits/stdc++.h> using namespace std; #define call_from_test #include "../../mod/mint.cpp" #include "../../math/linearfunction.cpp" #include "../../segtree/basic/ushi.cpp" #undef call_from_test signed main(){ cin.tie(0); ios::sync_with_stdio(0); int n,q; cin>>n>>q; vector<int> as(n),bs(n); for(int i=0;i<n;i++) cin>>as[i]>>bs[i]; using M = Mint<int, 998244353>; using T = LinearFunction<M>; auto f=[&](T x,T y){return x*y;}; T ti=T::mul_identity(); SegmentTree<T> seg(f,ti); vector<T> vt; for(int i=0;i<n;i++) vt.emplace_back(as[i],bs[i]); seg.build(vt); for(int i=0;i<q;i++){ int t; cin>>t; if(t==0){ int p,c,d; cin>>p>>c>>d; seg.set_val(p,T(M(c),M(d))); } if(t==1){ int l,r,x; cin>>l>>r>>x; cout<<seg.query(l,r)(x)<<"\n"; } } cout<<flush; return 0; }
#line 1 "test/yosupo/point_set_range_composite.test.cpp" // verification-helper: PROBLEM https://judge.yosupo.jp/problem/point_set_range_composite #include<bits/stdc++.h> using namespace std; #define call_from_test #line 1 "mod/mint.cpp" #line 3 "mod/mint.cpp" using namespace std; #endif //BEGIN CUT HERE template<typename T, T MOD = 1000000007> struct Mint{ inline static constexpr T mod = MOD; T v; Mint():v(0){} Mint(signed v):v(v){} Mint(long long t){v=t%MOD;if(v<0) v+=MOD;} Mint pow(long long k){ Mint res(1),tmp(v); while(k){ if(k&1) res*=tmp; tmp*=tmp; k>>=1; } return res; } static Mint add_identity(){return Mint(0);} static Mint mul_identity(){return Mint(1);} Mint inv(){return pow(MOD-2);} Mint& operator+=(Mint a){v+=a.v;if(v>=MOD)v-=MOD;return *this;} Mint& operator-=(Mint a){v+=MOD-a.v;if(v>=MOD)v-=MOD;return *this;} Mint& operator*=(Mint a){v=1LL*v*a.v%MOD;return *this;} Mint& operator/=(Mint a){return (*this)*=a.inv();} Mint operator+(Mint a) const{return Mint(v)+=a;} Mint operator-(Mint a) const{return Mint(v)-=a;} Mint operator*(Mint a) const{return Mint(v)*=a;} Mint operator/(Mint a) const{return Mint(v)/=a;} Mint operator+() const{return *this;} Mint operator-() const{return v?Mint(MOD-v):Mint(v);} bool operator==(const Mint a)const{return v==a.v;} bool operator!=(const Mint a)const{return v!=a.v;} static Mint comb(long long n,int k){ Mint num(1),dom(1); for(int i=0;i<k;i++){ num*=Mint(n-i); dom*=Mint(i+1); } return num/dom; } }; template<typename T, T MOD> ostream& operator<<(ostream &os,Mint<T, MOD> m){os<<m.v;return os;} //END CUT HERE #ifndef call_from_test signed main(){ return 0; } #endif #line 1 "math/linearfunction.cpp" #line 3 "math/linearfunction.cpp" using namespace std; #endif //BEGIN CUT HERE // a * x + b template<typename T> struct LinearFunction{ T a,b; LinearFunction():a(0),b(0){} LinearFunction(T a,T b):a(a),b(b){} using LF = LinearFunction; static LF add_identity(){return LF(T(0),T(0));} static LF mul_identity(){return LF(T(1),T(0));} LF& operator+=(const LF &o){ a+=o.a;b+=o.b; return *this; } LF& operator*=(const LF &o){ a=a*o.a; b=b*o.a+o.b; return *this; } LF operator+(const LF &o)const{return LF(*this)+=o;} LF operator*(const LF &o)const{return LF(*this)*=o;} T operator()(T x) const{return a*x+b;} }; //END CUT HERE #ifndef call_from_test //INSERT ABOVE HERE signed main(){ return 0; } #endif #line 1 "segtree/basic/ushi.cpp" #line 3 "segtree/basic/ushi.cpp" using namespace std; #endif //BEGIN CUT HERE template <typename T> struct SegmentTree{ using F = function<T(T,T)>; int n; F f; T ti; vector<T> dat; SegmentTree(F f,T ti):f(f),ti(ti){} void init(int n_){ n=1; while(n<n_) n<<=1; dat.assign(n<<1,ti); } void build(const vector<T> &v){ int n_=v.size(); init(n_); for(int i=0;i<n_;i++) dat[n+i]=v[i]; for(int i=n-1;i;i--) dat[i]=f(dat[(i<<1)|0],dat[(i<<1)|1]); } void set_val(int k,T x){ dat[k+=n]=x; while(k>>=1) dat[k]=f(dat[(k<<1)|0],dat[(k<<1)|1]); } T query(int a,int b){ if(a>=b) return ti; T vl=ti,vr=ti; for(int l=a+n,r=b+n;l<r;l>>=1,r>>=1) { if(l&1) vl=f(vl,dat[l++]); if(r&1) vr=f(dat[--r],vr); } return f(vl,vr); } template<typename C> int max_right(int s,C &check,T &acc,int k,int l,int r){ if(l+1==r){ acc=f(acc,dat[k]); return check(acc)?-1:k-n; } int m=(l+r)>>1; if(m<=s) return max_right(s,check,acc,(k<<1)|1,m,r); if(s<=l and check(f(acc,dat[k]))){ acc=f(acc,dat[k]); return -1; } int vl=max_right(s,check,acc,(k<<1)|0,l,m); if(~vl) return vl; return max_right(s,check,acc,(k<<1)|1,m,r); } // max t s.t. check(query(s,t)) // O(\log N) template<typename C> int max_right(int s,C &check){ assert(s<n and check(ti) and not check(query(s,n))); T acc=ti; return max_right(s,check,acc,1,0,n); } }; //END CUT HERE #ifndef call_from_test signed main(){ return 0; } #endif #line 10 "test/yosupo/point_set_range_composite.test.cpp" #undef call_from_test signed main(){ cin.tie(0); ios::sync_with_stdio(0); int n,q; cin>>n>>q; vector<int> as(n),bs(n); for(int i=0;i<n;i++) cin>>as[i]>>bs[i]; using M = Mint<int, 998244353>; using T = LinearFunction<M>; auto f=[&](T x,T y){return x*y;}; T ti=T::mul_identity(); SegmentTree<T> seg(f,ti); vector<T> vt; for(int i=0;i<n;i++) vt.emplace_back(as[i],bs[i]); seg.build(vt); for(int i=0;i<q;i++){ int t; cin>>t; if(t==0){ int p,c,d; cin>>p>>c>>d; seg.set_val(p,T(M(c),M(d))); } if(t==1){ int l,r,x; cin>>l>>r>>x; cout<<seg.query(l,r)(x)<<"\n"; } } cout<<flush; return 0; }