#include #include #include using namespace std; using int64 = long long; struct LazySegmentTree { private: int64 n; vector node, lazy; public: LazySegmentTree(vector v) { int64 sz = (int64)v.size(); n = 1; while(n < sz) n *= 2; node.resize(2*n-1); lazy.resize(2*n-1, 0); for(int64 i=0; i=0; i--) node[i] = node[i*2+1] + node[i*2+2]; } // k 番目のノードについて遅延評価を行う void eval(int64 k, int64 l, int64 r) { // 遅延配列が空でない場合、自ノード及び子ノードへの // 値の伝播が起こる if(lazy[k] != 0) { node[k] += lazy[k]; // 最下段かどうかのチェックをしよう // 子ノードは親ノードの 1/2 の範囲であるため、 // 伝播させるときは半分にする if(r - l > 1) { lazy[2*k+1] += lazy[k] / 2; lazy[2*k+2] += lazy[k] / 2; } // 伝播が終わったので、自ノードの遅延配列を空にする lazy[k] = 0; } } void add(int64 a, int64 b, int64 x, int64 k=0, int64 l=0, int64 r=-1) { if(r < 0) r = n; // k 番目のノードに対して遅延評価を行う eval(k, l, r); // 範囲外なら何もしない if(b <= l || r <= a) return; // 完全に被覆しているならば、遅延配列に値を入れた後に評価 if(a <= l && r <= b) { lazy[k] += (r - l) * x; eval(k, l, r); } // そうでないならば、子ノードの値を再帰的に計算して、 // 計算済みの値をもらってくる else { add(a, b, x, 2*k+1, l, (l+r)/2); add(a, b, x, 2*k+2, (l+r)/2, r); node[k] = node[2*k+1] + node[2*k+2]; } } int64 getsum(int64 a, int64 b, int64 k=0, int64 l=0, int64 r=-1) { if(r < 0) r = n; if(b <= l || r <= a) return 0; // 関数が呼び出されたら評価! eval(k, l, r); if(a <= l && r <= b) return node[k]; int64 vl = getsum(a, b, 2*k+1, l, (l+r)/2); int64 vr = getsum(a, b, 2*k+2, (l+r)/2, r); return vl + vr; } }; int main() { int64 n,a,b; cin >> n >> a >> b; vector x(n,0LL); for (int64 i=0; i> x[i]; LazySegmentTree rseg(vector(n,0LL)); LazySegmentTree lseg(vector(n,0LL)); cerr<<"calc1"<=0; i--) { auto l = lower_bound(x.begin(), x.end(), x[i]-b); auto r = upper_bound(x.begin(), x.end(), x[i]-a); int64 lindex = l - x.begin(); int64 rindex = r - x.begin(); lseg.add(lindex, rindex, rindex-lindex + lseg.getsum(i,i+1)); } cerr<<"ans"<