原题传送门

这道题还算简单

我们要求的期望值:

$$\frac{\sum_{i=l}r\sum_{j=l}rdis[i][j]}{C_{r-l+1}^{2}}$$

当然是上下两部分分别求,下面肥肠容易 ,问题在于如何求上面的

我们珂以把上面的换一个形式(枚举每段路会走几次):

$$\sum_{i=l}^ra[i](r-i+1)(i-l+1)$$

化简一下这个式子:

$$(r-l+1-rl)sum1+(r+l)*sum2-sum3$$

其中\(sum1=\sum_{i=l}^ra[i]\),\(sum2=\sum_{i=l}^ra[i]*i\),\(sum3=\sum_{i=l}^ra[i]*i^2\)

这样就珂以用线段树维护了

#include <bits/stdc++.h>
#define ll long long
#define N 100005
#define getchar nc
using namespace std;
inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline ll read()
{
register ll x=0,f=1;register char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
return x*f;
}
inline void write(register ll x)
{
if(!x)putchar('0');if(x<0)x=-x,putchar('-');
static int sta[25];register int tot=0;
while(x)sta[tot++]=x%10,x/=10;
while(tot)putchar(sta[--tot]+48);
}
struct node{
ll sum[6],tag;
}tr[N<<3];
int n,m;
ll sum1,sum2,sum3;
inline void build(register int x,register int l,register int r)
{
if(l==r)
{
tr[x].sum[4]=1ll*l;
tr[x].sum[5]=1ll*l*l;
return;
}
int mid=l+r>>1;
build(x<<1,l,mid);
build(x<<1|1,mid+1,r);
tr[x].sum[4]=tr[x<<1].sum[4]+tr[x<<1|1].sum[4];
tr[x].sum[5]=tr[x<<1].sum[5]+tr[x<<1|1].sum[5];
}
inline void work(register int x,register int l,register int r,register ll val)
{
tr[x].sum[1]+=1ll*(r-l+1)*val;
tr[x].sum[2]+=tr[x].sum[4]*val;
tr[x].sum[3]+=tr[x].sum[5]*val;
tr[x].tag+=val;
}
inline void pushdown(register int x,register int l,register int r)
{
int mid=l+r>>1;
work(x<<1,l,mid,tr[x].tag);
work(x<<1|1,mid+1,r,tr[x].tag);
tr[x].tag=0;
}
inline void pushup(register int x)
{
tr[x].sum[1]=tr[x<<1].sum[1]+tr[x<<1|1].sum[1];
tr[x].sum[2]=tr[x<<1].sum[2]+tr[x<<1|1].sum[2];
tr[x].sum[3]=tr[x<<1].sum[3]+tr[x<<1|1].sum[3];
}
inline void change(register int x,register int l,register int r,register int L,register int R,register ll val)
{
if(L<=l&&r<=R)
{
work(x,l,r,val);
return;
}
if(tr[x].tag)
pushdown(x,l,r);
int mid=l+r>>1;
if(L<=mid)
change(x<<1,l,mid,L,R,val);
if(R>mid)
change(x<<1|1,mid+1,r,L,R,val);
pushup(x);
}
inline void query(register int x,register int l,register int r,register int L,register int R)
{
if(L<=l&&r<=R)
{
sum1+=tr[x].sum[1];
sum2+=tr[x].sum[2];
sum3+=tr[x].sum[3];
return;
}
if(tr[x].tag)
pushdown(x,l,r);
int mid=l+r>>1;
if(L<=mid)
query(x<<1,l,mid,L,R);
if(R>mid)
query(x<<1|1,mid+1,r,L,R);
}
inline ll gcd(register ll a,register ll b)
{
return !b?a:gcd(b,a%b);
}
int main()
{
n=read(),m=read();
build(1,1,n);
while(m--)
{
char ch=getchar();
while(ch!='C'&&ch!='Q')
ch=getchar();
ll l=read(),r=read()-1;
if(ch=='C')
{
ll v=read();
change(1,1,n,l,r,v);
}
else
{
sum1=sum2=sum3=0;
query(1,1,n,l,r);
ll ansa=(r-l+1-r*l)*sum1+(r+l)*sum2-sum3;
ll ansb=(r-l+2)*(r-l+1)/2;
ll g=gcd(ansa,ansb);
write(ansa/g),putchar('/'),write(ansb/g),puts("");
}
}
return 0;
}

【题解】Luogu P2221 [HAOI2012]高速公路的更多相关文章

  1. luogu P2221 [HAOI2012]高速公路题解

    题面 很套路的拆式子然后线段树上维护区间和的题.一般都是把式子拆成区间内几个形如\(\sum i*a_i, \sum i^2 * a_i\)的式子相加减的形式. 考虑一次询问[l,r]的答案怎么算: ...

  2. P2221 [HAOI2012]高速公路(线段树)

    P2221 [HAOI2012]高速公路 显然答案为 $\dfrac{\sum_{i=l}^r\sum_{j=l}^{r}dis[i][j]}{C_{r-l+1}^2}$ 下面倒是挺好算,组合数瞎搞 ...

  3. [Luogu 2221] HAOI2012 高速公路

    [Luogu 2221] HAOI2012 高速公路 比较容易看出的线段树题目. 由于等概率,期望便转化为 子集元素和/子集个数. 每一段l..r中,子集元素和为: \(\sum w_{i}(i-l+ ...

  4. 洛谷 P2221 [HAOI2012]高速公路

    链接: P2221 题意: 有 \(n(1\leq n\leq 10^5)\) 个点,从第 \(i(1\leq i< n)\) 个点向第 \(i+1\) 个点连有边.最初所有边长 \(v_i\) ...

  5. 洛谷P2221 [HAOI2012]高速公路

    线段树 #include<cstdio> #include<cstdlib> #include<algorithm> #include<cstring> ...

  6. P2221 [HAOI2012]高速公路

    思路 考虑每一条边的贡献,然后推式子 \[ \begin{align}&\sum_{i}V_i\times(R-i+1)\times(i-L+1)\\=&\sum_{i}V_i\lef ...

  7. 洛谷P2221 [HAOI2012]高速公路(线段树+概率期望)

    传送门 首先,答案等于$$ans=\sum_{i=l}^r\sum_{j=i}^r\frac{sum(i,j)}{C_{r-l+1}^2}$$ 也就是说所有情况的和除以总的情况数 因为这是一条链,我们 ...

  8. BZOJ2752: [HAOI2012]高速公路(road)

    2752: [HAOI2012]高速公路(road) Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 608  Solved: 199[Submit][ ...

  9. 【线段树】BZOJ2752: [HAOI2012]高速公路(road)

    2752: [HAOI2012]高速公路(road) Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 1621  Solved: 627[Submit] ...

随机推荐

  1. iOS 如何查看APP的jetsamEvent日志

    1.如何在iPhone上查看 设置-通用-分析-分析数据- JetsamEvent+日志 打头的系统日志. 2.如何在Mac 上查看此类分析日志 1.手机链接MAC 2.打开iTunes,点开手机图标 ...

  2. Python 3+selenium+unittest+HTMLTestRunner生成测试报告

    一.下载HTMLTestRunner.py,解压,将它放到 python安装路径的site-packages目录下 https://pan.baidu.com/s/1epWlibxbxWlNoIcxL ...

  3. Linux 部署 xxl-job 注意问题

    问题:Failed to create parent directories for [/data/applogs/xxl-job/xxl-job-admin.log][原因:权限不足] 启动终端: ...

  4. scrapy爬虫框架和selenium的配合使用

    scrapy框架的请求流程 scrapy框架? Scrapy 是基于twisted框架开发而来,twisted是一个流行的事件驱动的python网络框架.因此Scrapy使用了一种非阻塞(又名异步)的 ...

  5. linux驱动调试记录

    linux驱动调试 linux 目录 /proc 下面可以配置驱动的调试信息,比如给proc目录的自己定制的驱动的一文件设置一个变量,然后驱动程序跟了proc的参数值来配置调试级别.类似于内核调试的级 ...

  6. webpack.config.js配置遇到Error: Cannot find module '@babel/core'&&Cannot find module '@babel/plugin-transform-react-jsx' 问题

    下文是网上找到的方法,是因为版本冲突的原因,参照后安装7版本解决 cnpm install -D babel-loader@ babel-core babel-preset-env 一. 问题描述 在 ...

  7. [svc]简单理解什么是rpc调用?跟restapi有何区别?

    什么是rpc调用 restapi调用方式是对数据的crud. 常见的我们写flash写个api,或者借助django drf写个标准的resetapi,一个url可以借助httpget post pu ...

  8. 【Linux】Mac PD set centos static ip

    2,修改Centos的网络设置. (1)进入脚本. vi /etc/sysconfig/network-scripts/ifcfg-eth0 My Mac ip: # 从dhcp改成static BO ...

  9. sping Bean 的生命周期是如何被管理

    1. 实例化一个Bean,也就是我们通常说的new 2. 按照Spring上下文对实例化的Bean进行配置,也就是IOC注入 3. 如果这个Bean实现了BeanNameAware接口,会调用它实现的 ...

  10. JavaScript Object中的函数assign

    Object函数提供了一个叫做assign的函数,用来合并多个对象. Object.assign(...): 你可以传递多个对象给该函数,这些对象中的自有且可枚举的属性,会被拷贝给第一个对象. var ...