HAOI2012高速公路bzoj2752 (线段树,数学)
题目大意:
给定一个长度为n的链,一共m次操作
对于每次操作
\(C\ l\ r\ x\)表示将第l个点到第r个点之间的所有道路的权值增加v
\(Q\ l\ r\)在第l个到第r个点里等概率随机取出两个不同的点a和b,那么从a行驶到b将期望花费多少费用呢
QwQ我们可以考虑将期望分为分子和分母两部分
首先考虑分母,分母就是在\(r-l+1\)个点中选两个点的方案数,也就是\({r-l+1}\choose 2\)
而分子就是总权值了
对于一个在\([l,r]\)的点\(i\)来说
它会被计算\((i-l+1)\times (r-i+1)\)次 (就跟分别在左右两个区间枚举端点是一样的)
那么我们求的东西也就变成了
\(\sum_{i=l}^{r} w[i]\times(i-l+1)\times(r-i+1)\)
我们考虑将它展开:
\(\sum_{i=l}^{r} w[i]\times((r+l)\times i-i\times i +(r-l+1-r\times l))\)
再进一步来一波,对于区间\([l,r]\)我们就是求的是:
\((r+l)\times \sum_{i=l}^{r}w[i]\times i + (r-l+1-r\times l)\times \sum_{i=l}^{r}w[i] + \sum_{i=l}^{r}w[i]*i^2\)
而对于更新,我们也不难找出一些规律了咯
所以对于一个区间,我们只需要维护一个\(w[i],w[i]*i,w[i]*i^2,i^2,i\)的sigma值就能处理更新和合并了!
void up(int root)
{
f[root].si=f[2*root].si+f[2*root+1].si;
f[root].sii=f[2*root].sii+f[2*root+1].sii;
f[root].scostii=f[2*root].scostii+f[2*root+1].scostii;
f[root].scosti=f[2*root].scosti+f[2*root+1].scosti;
f[root].scost=f[2*root].scost+f[2*root+1].scost;
}
void pushdown(int root,int l,int r)
{
ll mid = (l+r) >> 1;
if (add[root])
{
add[2*root]+=add[root];
add[2*root+1]+=add[root];
f[2*root].scost+=add[root]*(mid-l+1);
f[2*root+1].scost+=add[root]*(r-mid);
f[2*root].scosti+=add[root]*f[2*root].si;
f[2*root+1].scosti+=add[root]*f[2*root+1].si;
f[2*root].scostii+=add[root]*f[2*root].sii;
f[2*root+1].scostii+=add[root]*f[2*root+1].sii;
add[root]=0;
}
}
emmmmm
query的时候也记得\(\times\)的时候要乘询问的总区间,而不是当然区间!
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<queue>
#include<map>
#include<vector>
#define ll long long
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while (!isdigit(ch)) {if (ch=='-') f=-1;ch=getchar();}
while (isdigit(ch)) {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*f;
}
const int maxn = 1e5+1e2;
struct Node{
ll scost,scosti,scostii,si,sii;
ll ans;
};
Node f[4*maxn];
ll add[4*maxn];
int n,m;
ll count(int xx,int yy)
{
ll x=xx,y=yy;
return (y-x+1)*(y-x+2)/2;
}
void up(int root)
{
f[root].si=f[2*root].si+f[2*root+1].si;
f[root].sii=f[2*root].sii+f[2*root+1].sii;
f[root].scostii=f[2*root].scostii+f[2*root+1].scostii;
f[root].scosti=f[2*root].scosti+f[2*root+1].scosti;
f[root].scost=f[2*root].scost+f[2*root+1].scost;
}
void pushdown(int root,int l,int r)
{
ll mid = (l+r) >> 1;
if (add[root])
{
add[2*root]+=add[root];
add[2*root+1]+=add[root];
f[2*root].scost+=add[root]*(mid-l+1);
f[2*root+1].scost+=add[root]*(r-mid);
f[2*root].scosti+=add[root]*f[2*root].si;
f[2*root+1].scosti+=add[root]*f[2*root+1].si;
f[2*root].scostii+=add[root]*f[2*root].sii;
f[2*root+1].scostii+=add[root]*f[2*root+1].sii;
add[root]=0;
}
}
void build(int root,int l,int r)
{
if (l==r)
{
f[root].si=(long long)1LL*l;
f[root].sii=(long long) 1LL*l*l;
return;
}
int mid = (l+r) >> 1;
build(2*root,l,mid);
build(2*root+1,mid+1,r);
up(root);
}
void update(int root,int l,int r,int x,int y,ll p)
{
if (x<=l && r<=y)
{
ll len=r-l+1;
f[root].scost+=1LL*p*len;
f[root].scosti+=1LL*p*f[root].si;
f[root].scostii+=1LL*p*f[root].sii;
add[root]+=p;
return;
}
pushdown(root,l,r);
int mid = (l+r) >> 1;
if (x<=mid) update(2*root,l,mid,x,y,p);
if (y>mid) update(2*root+1,mid+1,r,x,y,p);
up(root);
}
ll query(int root,int l,int r,int x,int y)
{
if (x<=l && r<=y)
{
ll xx = x,yy=y;
return (long long)1LL*(xx+yy)*f[root].scosti+(long long)1LL*(yy-xx+1-xx*yy)*f[root].scost-f[root].scostii;
}
pushdown(root,l,r);
int mid = (l+r) >> 1;
ll ans=0;
if (x<=mid) ans+=query(2*root,l,mid,x,y);
if (y>mid) ans+=query(2*root+1,mid+1,r,x,y);
return ans;
}
int main()
{
scanf("%d%d",&n,&m);
n--;
build(1,1,n);
for (int i=1;i<=m;i++)
{
char s[10];
scanf("%s",s+1);
if (s[1]=='C')
{
int x,y;
ll z;
scanf("%d%d%lld",&x,&y,&z);
y--;
update(1,1,n,x,y,z);
}
if (s[1]=='Q')
{
int x,y;
scanf("%d%d",&x,&y);
y--;
ll tmp=query(1,1,n,x,y);
ll tmp1=count(x,y);
printf("%lld/%lld\n",tmp/__gcd(tmp1,tmp),tmp1/__gcd(tmp1,tmp));
}
}
return 0;
}
HAOI2012高速公路bzoj2752 (线段树,数学)的更多相关文章
- BZOJ2752: [HAOI2012]高速公路(road)(线段树 期望)
Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 1820 Solved: 736[Submit][Status][Discuss] Descripti ...
- 【bzoj2752】[HAOI2012]高速公路(road) 线段树
题目描述 Y901高速公路是一条重要的交通纽带,政府部门建设初期的投入以及使用期间的养护费用都不低,因此政府在这条高速公路上设立了许多收费站.Y901高速公路是一条由N-1段路以及N个收费站组成的东西 ...
- 2019.01.14 bzoj2752: [HAOI2012]高速公路(线段树)
传送门 线段树菜题. 题意简述:给一条nnn个点的链,链有边权,支持区间修改边权,查询在一段区间内随机选择不同的起点和终点路径的期望总边权和. 思路:考虑每条边的贡献. 考虑对于一段区间[l,r][l ...
- BZOJ 2752: [HAOI2012]高速公路(road)( 线段树 )
对于询问[L, R], 我们直接考虑每个p(L≤p≤R)的贡献,可以得到 然后化简一下得到 这样就可以很方便地用线段树, 维护一个p, p*vp, p*(p+1)*vp就可以了 ----------- ...
- BZOJ 2752: [HAOI2012]高速公路(road) [线段树 期望]
2752: [HAOI2012]高速公路(road) Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 1219 Solved: 446[Submit] ...
- P2221 [HAOI2012]高速公路(线段树)
P2221 [HAOI2012]高速公路 显然答案为 $\dfrac{\sum_{i=l}^r\sum_{j=l}^{r}dis[i][j]}{C_{r-l+1}^2}$ 下面倒是挺好算,组合数瞎搞 ...
- 洛谷P2221 [HAOI2012]高速公路(线段树+概率期望)
传送门 首先,答案等于$$ans=\sum_{i=l}^r\sum_{j=i}^r\frac{sum(i,j)}{C_{r-l+1}^2}$$ 也就是说所有情况的和除以总的情况数 因为这是一条链,我们 ...
- [luoguP2221] [HAOI2012]高速公路(线段树)
传送门 考虑每一段对答案的贡献 用每一段的左端点来表示当前这一段,那么区间就变成了[1,n-1] 如果询问区间[l,r],其中一个点的位置为x,则它对答案的贡献为(x-l)*(r-x)*s[x](s[ ...
- Bzoj 2752 高速公路 (期望,线段树)
Bzoj 2752 高速公路 (期望,线段树) 题目链接 这道题显然求边,因为题目是一条链,所以直接采用把边编上号.看成序列即可 \(1\)与\(2\)号点的边连得是. 编号为\(1\)的点.查询的时 ...
- [bzoj2752]高速公路 题解(线段树)
2752: [HAOI2012]高速公路(road) Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 2102 Solved: 887[Submit] ...
随机推荐
- VMware ESXi 7.0 U2 SLIC & Unlocker USB 网卡驱动集成镜像 202109更新
2021.08.31 更新:集成 "vmkusb-nic-fling"."net-community" 和 "nvme-community" ...
- C++回调机制
一直对回调机制不是很了解,今天索性搜了很多资料顺便整理一下,进步一点点. 1.Callback方式(回调函数) 什么是回调函数? 简而言之,回调函数就是一个通过函数指针调用的函数.如果你把函数的指针( ...
- zap高性能日志
摘要 日志在整个工程实践中的重要性不言而喻,在选择日志组件的时候也有多方面的考量.详细.正确和及时的反馈是必不可少的,但是整个性能表现是否也是必要考虑的点呢?在长期的实践中发现有的日志组件对于计算资源 ...
- NOIP模拟16:「Star Way To Heaven·God Knows·Loost My Music」
T1:Star Way To Heaven 基本思路: 最小生成树. 假如我们将上边界与下边界看作一个点,然后从上边界经过星星向下边界连边,会发现,他会形成一条线将整个矩形分为左右两个部分. ...
- Vue项目中应用TypeScript
一.前言 与如何在React项目中应用TypeScript类似 在VUE项目中应用typescript,我们需要引入一个库vue-property-decorator, 其是基于vue-class-c ...
- 关于vue-cli的安装
(一):*安装 vue-cli 参考: https://cn.vuejs.org/v2/guide/installation.html https://github.com/vuejs/vue-cli ...
- 树莓派的kodi设置遥控器的方法
首先你需要买一个红外接收器,根据卖家的文档,插到树莓派的GPIO串口上, 我的红外接收器是18入口,17出口, 所以我的config.txt文件设置如下 dtoverlay=lirc-rpi,gpio ...
- CodeForce-811C Vladik and Memorable Trip(动态规划)
Vladik and Memorable Trip CodeForces - 811C 有一个长度为 n 的数列,其中第 i 项为 ai. 现在需要你从这个数列中选出一些互不相交的区间,并且保证整个数 ...
- python爬取疫情数据存入MySQL数据库
import requests from bs4 import BeautifulSoup import json import time from pymysql import * def mes( ...
- JAVA语言程序设计课程评价
紧张的又短暂的一个学期结束了,这个学期也许将成为我人生中一个重要的转折点,作为一名半路出家的选手,在初次了解Java语言时我感到非常的迷茫与不知所措.因为之前很多同学都是通过假期时间在家自学,刚转入新 ...