Problem

Description

\(Y901\) 高速公路是一条重要的交通纽带,政府部门建设初期的投入以及使用期间的养护费用都不低,因此政府在这条高速公路上设立了许多收费站。

\(Y901\) 高速公路是一条由 \(n-1\) 段路以及 \(n\) 个收费站组成的东西向的链,我们按照由西向东的顺序将收费站依次编号为 \(1\sim n\) ,从收费站 \(i\) 行驶到 \(i+1\) (或从 \(i+1\) 行驶到 \(i\) )需要收取 \(V_i\) 的费用。高速路刚建成时所有的路段都是免费的。

政府部门根据实际情况,会不定期地对连续路段的收费标准进行调整,根据政策涨价或降价。

无聊的小 \(A\) 同学总喜欢研究一些稀奇古怪的问题,他开车在这条高速路上行驶时想到了这样一个问题:对于给定的\(l,r(l< r)\) ,在第 \(l\) 个到第 \(r\) 个收费站里等概率随机取出两个不同的收费站 \(a\) 和 \(b\),那么从 \(a\) 行驶到 \(b\) 将期望花费多少费用呢?

Input Format

第一行 \(2\) 个正整数 \(n,m\),表示有 \(n\) 个收费站,\(m\) 次调整或询问

接下来 \(m\) 行,每行将出现以下两种形式中的一种

\(C\ \ l\ \ r\ \ v\) 表示将第 \(l\) 个收费站到第 \(r\) 个收费站之间的所有道路的通行费全部增加 \(v\)

\(Q\ \ l\ \ r\) 表示对于给定的 \(l,r\),要求回答小 \(A\) 的问题

所有 \(C\) 与 \(Q\) 操作中保证 \(1\le l< r\le n\)

Output Format

对于每次询问操作回答一行,输出一个既约分数

若答案为整数 \(ans\) ,输出 \(ans/1\)。

Sample

Input

4 5
C 1 4 2
C 1 2 -1
Q 1 2
Q 2 4
Q 1 4

Output

1/1
8/3
17/6

Explanation

所有 \(C\) 操作中的 \(v\) 的绝对值不超过 \(10000\)

在任何时刻任意道路的费用均为不超过 \(10000\) 的非负整数

所有测试点的详细情况如下表所示

Range

\(Test\) \(n\) \(m\)
1 =10 =10
2 =100 =100
3 =1000 =1000
4 =10000 =10000
5 =50000 =50000
6 =60000 =60000
7 =70000 =70000
8 =80000 =80000
9 =90000 =90000
10 =100000 =100000

Algorithm

线段树

Mentality

很显然,这题和期望概率毫无关系,只需要求出总贡献再去除以 \(C_{r-l+1}^2\) 即可。

那么如何求总贡献呢?

考虑分开计算贡献,我们最后的总贡献肯定是:

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

仔细凝视这个式子肯定是看不出什么的还会显得你很愚蠢 ,动手才是正途。我们先把这两个括号拆开看看:

\[V[i]*(i-l+1)*(r-i+1)=V[i]*(r*i-i^2+i-l*r+l*i-l+r-i+1)\\
=V[i]*(r-l+1-r*l)+V[i]*i*(r+l)-V[i]*i^2
\]

好的,既然我们拆成了这样,那么我们的 \(\sum\) 就成功地被一拆三了!由于我们已经把包含 \(i\) 的项都单独放在了外面,剩下的都是定值,那么我们的总贡献就改为了这样的三个 \(\sum\) 之和:

\[((r-l+1-r*l)*\sum V[i])+((r+l)*\sum V[i]*i)+(\sum V[i]*i^2)
\]

由于 \(r-l+1-r*l\) 和 \(r+l\) 都很普通,我们只需要用线段树维护三个含 \(i\) 的值就好啦!

至于这三个值怎么维护呢?看着满阔怕,但是由于 \(\sum\) 中的每一项都是单独算的,那么不难发现,我们对一个区间都加上 \(v\) 时,这段区间上三个值加上的贡献肯定分别为:\(v*(r-l+1)\) ,\(v*\sum i\),\(v*\sum i^2\) 。那么我们只需要预处理一下 \(i\) 和 \(i^2\) 的前缀和即可!

最后和 \(C_{r-l+1}^2\) 同除 \(gcd\) 即可。

啥?代码有哪些坑点?没啥坑点,自己体会 = = 。

Code

#include <cstdio>
#include <iostream>
using namespace std;
#define ls (o << 1)
#define rs ((o << 1) + 1)
#define mid ((l + r) >> 1)
int n, m;
long long x, L, R, ans[4], a[100001], q1[100001], q2[100001], adv[400001],
sum[4][400001];
char opt;
void pushup(int o) {
for (int i = 1; i <= 3; i++) sum[i][o] = sum[i][ls] + sum[i][rs];
}
void pushdown(int o, int l, int r) {
sum[1][ls] += 1ll * adv[o] * (mid - l + 1);
sum[2][ls] += adv[o] * (q1[mid] - q1[l - 1]);
sum[3][ls] += adv[o] * (q2[mid] - q2[l - 1]);
sum[1][rs] += 1ll * adv[o] * (r - mid);
sum[2][rs] += adv[o] * (q1[r] - q1[mid]);
sum[3][rs] += adv[o] * (q2[r] - q2[mid]);
adv[ls] += adv[o], adv[rs] += adv[o];
adv[o] = 0;
}
void add(int o, int l, int r) {
if (l >= L && r <= R) {
adv[o] += x;
sum[1][o] += 1ll * x * (r - l + 1);
sum[2][o] += x * (q1[r] - q1[l - 1]);
sum[3][o] += x * (q2[r] - q2[l - 1]);
return;
}
pushdown(o, l, r);
if (mid >= L) add(ls, l, mid);
if (mid < R) add(rs, mid + 1, r);
pushup(o);
}
void query(int o, int l, int r) {
if (l >= L && r <= R) {
for (int i = 1; i <= 3; i++) ans[i] += sum[i][o];
return;
}
pushdown(o, l, r);
if (mid >= L) query(ls, l, mid);
if (mid < R) query(rs, mid + 1, r);
pushup(o);
}
long long gcd(long long a, long long b) { return !b ? a : gcd(b, a % b); }
int main() {
freopen("2221.in", "r", stdin);
freopen("2221.out", "w", stdout);
cin >> n >> m;
for (int i = 1; i <= n; i++) {
scanf("%lld", &a[i]);
q1[i] = q1[i - 1] + 1ll * i;
q2[i] = q2[i - 1] + 1ll * i * i;
}
while (m--) {
scanf("%s%lld%lld", &opt, &L, &R);
R--;
if (opt == 'C') {
scanf("%lld", &x);
add(1, 1, n);
} else {
ans[1] = ans[2] = ans[3] = 0;
query(1, 1, n);
long long Ans = (R - L + 1 - R * L) * ans[1] + (L + R) * ans[2] - ans[3],
Tot = (R - L + 2) * (R - L + 1) / 2;
long long Gcd = gcd(Ans, Tot);
printf("%lld/%lld\n", Ans / Gcd, Tot / Gcd);
}
}
}

【HAOI 2012】高速公路的更多相关文章

  1. HAOI 2012 高速公路

    https://www.luogu.org/problem/show?pid=2221 题目描述 Y901高速公路是一条重要的交通纽带,政府部门建设初期的投入以及使用期间的养护费用都不低,因此政府在这 ...

  2. 大暴力——[HAOI]2012音量调节

    题目:[HAOI]2012音量调节 描述: 问题描述 一个吉他手准备参加一场演出.他不喜欢在演出时始终使用同一个音量,所以他决定每一首歌之前他都要改变一次音量.在演出开始之前,他已经做好了一个列表,里 ...

  3. BZOJ 2750 HAOI 2012 Road 高速公路 最短路

    题意: 给出一个有向图,求每条边有多少次作为最短路上的边(任意的起始点). 范围:n <= 1500, m <= 5005 分析: 一个比较容易想到的思路:以每个点作为起点,做一次SPFA ...

  4. [HAOI 2012]音量调节

    Description 一个吉他手准备参加一场演出.他不喜欢在演出时始终使用同一个音量,所以他决定每一首歌之前他都要改变一次音量.在演出开始之前,他已经做好了一个列表,里面写着在每首歌开始之前他想要改 ...

  5. 解题:HAOI 2012 道路

    题面 这题不开O2怎么过=.= 可能这种有关最短路的计数题做多了就有些感觉了...... 以每个点为基准跑出一张最短路图,然后对每个边$(u,v)$统计两个东西.一个$pre[u]$表示到达$u$这个 ...

  6. BZOJ 2749 HAOI 2012 外星人 数论 欧拉函数

    题意: 给出一个数,给出的形式是其分解质因数后,对应的质因数pi及其次数qi,问对这个数不停求phi,直至这个数变成1,需要多少次.(多组数据) 范围:pi <= 1e5,qi <= 1e ...

  7. [HAOI 2012] Road

    [题目链接] https://www.lydsy.com/JudgeOnline/problem.php?id=2750 [算法] 考虑计算每个点对每条边的贡献 对于每个点首先运行SPFA或Dijks ...

  8. [HAOI 2012] 容易题

    [题目链接] https://www.lydsy.com/JudgeOnline/problem.php?id=2751 [算法] 考虑k = 0的情况 , 根据乘法原理 : Ans = (n * ( ...

  9. [HAOI 2012] 外星人

    [题目链接] https://www.lydsy.com/JudgeOnline/problem.php?id=2749 [算法] 首先 , 每次对一个数x进行操作 , 只会使该数减少一个2的因子 那 ...

随机推荐

  1. CString数组和CStringArray

    CStringArray是编译器定义的类型!可以进行一些(如:访问.增.删.改)等操作. 集中单个字符串的操作使用Cstring,集中一批字符串的管理使用CstringArray. 一个是动态,CSt ...

  2. filter滤镜效果(css3属性)

    <!DOCTYPE html> <html> <head> <style> img { width: 33%; height: auto; float: ...

  3. swiper嵌套小demo(移动端触摸滑动插件)

    swiper(移动端触摸滑动插件) tip:自己敲得Swiper 的小demo,可以复制粘贴看看效果哦. swiper的js包css包下链接地址 :  https://github.com/Clear ...

  4. 目前常用的加密算法有DES(Data Encryption Standard)和IDEA(International Data Encryption Algorithm)国际数据加密算法等,请用工厂方法实现加密算法系统。提交该系统的代码,该系统务必是一个可以能够直接使用的系统,查阅资料完成相应加密算法的实现;

    1.加密算法的类图结构 2.源代码 2.1代码运行截图 2.2代码的目录结构 2.3具体代码 MethodFactory.java package jiami; public interface Me ...

  5. 微信小程序开发笔记02

    今天学习了微信小程序开发用到的语言,wxml与wxss语言基本语法与html和css基本语法相似,学习起来相对简单.在小程序主要的语言是js(javascript,跟准确的说是jqery) ,由于这种 ...

  6. WdatePicker日历添加事件,在任意月改变时处理日期事件

    原由 在做系统时根据要求有时候需要屏蔽掉某些特殊的日期,像周日或者法定假日,以及一些调班的日期:使用WdatePicker可以屏蔽掉周日和大多数法定假日,但像清明或者调班的日期则不好处理. 想法 1: ...

  7. layer[申明:来源于网络]

    layer[申明:来源于网络] 地址:http://layer.layui.com/

  8. [No0000193]Chrome浏览器控制台(console)花式调试

    对前端开发者来说,Chrome Dev Tools(开发者工具,以下简称CDT)是一个不可或缺的开发调试工具,但是你可能只用过console.log(),却不知道console还有很多功能强大的调试方 ...

  9. Improved GAN

    https://www.bilibili.com/video/av9770302/?p=16 从之前讲的basic gan延伸到unified framework,到WGAN 再到通过WGAN进行Ge ...

  10. 记录常用的adb命令

    1.启动adb服务 adb start-server 2.关闭服务 adb kill-server 3.进入shell环境 adb shell 4.安装应用 adb install -r xxx.ap ...