首先便于计算,我们可以将每一条边缩成一个点,于是右端点\(r\)需要减\(1\).
从\(a\)行驶到\(b\)将期望花费
\[ \begin{align} ans=\dfrac{Ans}{C^{2}_{r-l+2}} \end{align} \]
我们可以对于每一条边统计它被算了多少次.
\[ \begin{align} \therefore Ans&=\sum^r_{i=l}v[i]*(i-l+1)(r-i+1)\\ &=\sum^r_{i=l}v[i]*[-i^2+(l+r)*i+(1+r)(1-l)]\\ &=-\sum^r_{i=l}v[i]*i^2+(l+r)\sum^r_{i=l}v[i]*i-(r+1)(l-1)\sum^r_{i=l}v[i] \end{align} \]
于是我们只需要在线段树中维护三个值.
\[ \begin{align} tr[0][x]&=\sum^{x.r}_{i=x.l}a[i]\\ tr[1][x]&=\sum^{x.r}_{i=x.l}a[i]*i\\ tr[2][x]&=\sum^{x.r}_{i=x.l}a[i]*i^2\\ \end{align} \]
合并非常简单.
\[ \begin{align} tr[0][x]&=tr[0][x<<1]+tr[0][x<<1|1]\\ tr[1][x]&=tr[1][x<<1]+tr[1][x<<1|1]\\ tr[2][x]&=tr[2][x<<1]+tr[2][x<<1|1]\\ \end{align} \]
但是如何修改呢?
对于\(tr[0]\):
\[ \begin{align} \Delta tr[0][x]&=\sum^{x.r}_{i=x.l}\Delta v[i]\\ &=k\sum^{x.r}_{i=x.l}1\\ &=k*(x.r-x.l+1) \end{align} \]
对于\(tr[1]\):
\[ \begin{align} \Delta tr[1][x]&=\sum^{x.r}_{i=x.l}\Delta v[i]*i\\ &=k\sum^{x.r}_{i=x.l}i\\ &=k*(x.r+x.l)*(x.r-x.l+1)/2 \end{align} \]
对于\(tr[2]\):
\[ \begin{align} \Delta tr[2][x]&=\sum^{x.r}_{i=x.l}\Delta v[i]*i^2\\ &=k\sum^{x.r}_{i=x.l}i^2\\ &=k(\sum^{x.r}_{i=1}i^2-\sum^{x.l-1}_{i=1}i^2) \end{align} \]
相信大家都学过平方相加的公式:
\[ \begin{align} 1^2+2^2+......+n^2=\frac{n(n+1)(2n+1)}{6} \end{align} \]
没学过的可以康康
\[ \begin{align} \therefore\Delta tr[2][x]&=\frac{x.r(x.r+1)(2x.r+1)-(x.l-1)x.l(2x.l-1)}{6} \end{align} \]
最后:
一定要开longlong!
#pragma GCC optimize(3)
#include<bits/stdc++.h>
#define il inline
#define rg register
#define gi read<int>
using namespace std;
typedef long long ll;
const int O = 1e5 + 10;
template<class TT>
il TT read() {
TT o = 0,fl = 1; char ch = getchar();
while (!isdigit(ch) && ch != '-') ch = getchar();
if (ch == '-') fl = -1, ch = getchar();
while (isdigit(ch)) o = o * 10 + ch - '0', ch = getchar();
return fl * o;
}
char ch[2];
int n, m;
ll xx, yy, zz;
class SegmentTree {
private:
ll tr[6][O << 2], lazy[O << 2];
il void pushup(int x) {
tr[0][x] = tr[0][x << 1] + tr[0][x << 1 | 1];
tr[1][x] = tr[1][x << 1] + tr[1][x << 1 | 1];
tr[2][x] = tr[2][x << 1] + tr[2][x << 1 | 1];
}
il void work(int delta, int x, int len) {
tr[0][x] += 1ll * delta * len;
tr[1][x] += 1ll * delta * tr[4][x];
tr[2][x] += 1ll * delta * tr[5][x];
lazy[x] += delta;
}
il void pushdown(int x, int l, int r) {
if(!lazy[x]) return ;
int mid = l + r >> 1;
work(lazy[x], x << 1, mid - l + 1);
work(lazy[x], x << 1 | 1, r - mid);
lazy[x] = 0;
}
public:
il void Build(int x, int l, int r) {
if (l == r) {
tr[4][x] = l;
tr[5][x] = 1ll * l * l;
return ;
}
int mid = l + r >> 1;
Build(x << 1, l, mid);
Build(x << 1 | 1, mid + 1, r);
tr[4][x] = tr[4][x << 1] + tr[4][x << 1 | 1];
tr[5][x] = tr[5][x << 1] + tr[5][x << 1 | 1];
}
il void Modify(int x, int l, int r, int L, int R, int delta) {
if (l > R || L > r) return ;
if (L <= l && r <= R) return work(delta, x, r - l + 1);
pushdown(x, l, r);
int mid = l + r >> 1;
Modify(x << 1, l, mid, L, R, delta);
Modify(x << 1 | 1, mid + 1, r, L, R, delta);
pushup(x);
}
il void Query(int x, int l, int r, int L, int R) {
if (l > R || L > r) return ;
if (L <= l && r <= R) {
xx += tr[0][x];
yy += tr[1][x];
zz += tr[2][x];
return ;
}
pushdown(x, l, r);
int mid = l + r >> 1;
Query(x << 1, l, mid, L, R);
Query(x << 1 | 1, mid + 1, r, L, R);
}
}st;
int main() {
n = gi(), m = gi();
st.Build(1, 1, n);
while (m--) {
scanf("%s", ch);
int l = gi(), r = gi() - 1;
if (ch[0] == 'C') st.Modify(1, 1, n, l, r, gi());
else {
xx = yy = zz = 0;
st.Query(1, 1, n, l, r);
ll X = 1ll * (l + r) * yy + 1ll * (r + 1) * (1 - l) * xx - zz;
ll Y = 1ll * (r - l + 2) * (r - l + 1) >> 1, g = __gcd(X, Y);
X /= g, Y /= g;
printf("%lld/%lld\n", X, Y);
}
}
return 0;
}