题目链接

题意

灯有\(m\)个亮度等级,\(1,2,...,m\),有两种按钮:

  1. 每次将亮度等级\(+1\),如\(1\rightarrow 2,2\rightarrow 3,...,m-1\rightarrow m,m\rightarrow 1\)
  2. 初始时有一个设定值\(x\),按下该按钮能够从任意亮度等级到达\(x\)

现给定一个序列\(a[1..n]\)代表亮度从\(a_1\rightarrow a_2\rightarrow a_3\rightarrow ... \rightarrow a_n\),问初始值\(x\)设定为多少的时候能够使总共需按按钮的次数最少?

思路

原始思路

从\(1\)到\(m\)枚举\(x\),对于每个\(x\)遍历一遍序列,计算出操作次数。

比较得到一个最小值。时间复杂度\(O(m*n)\)

本质=>优化

思考一下,上面这个思路本质上是什么?

假设我们有一个数组\(c[1..m]\)用来记录操作次数,那么其初始值为\(0\),遍历一遍序列,对每个初始值\(x\)对应的操作次数\(c[x]\)加上\(f(a_i,a_{i+1},x)\).

最后,\(c[\ ]\)数组的最小值即为答案。

再来考虑一下上面的\(f(a_i,a_{i+1},x)\),这又是什么?

记\(s=a_i,t=a_{i+1}\),不妨设\(s\lt t\)则$$f(s,t,x)=

\begin{cases}

t-s, &x\leq s\ ||\ x\gt t\cr

t-x+1, &s\lt x \leq t

\end{cases}$$

// 高兴一下,是个线性函数啊

变化量如图所示:

1			s s+1		  t t+1		   m
|------------|------------|------------|
( +t-s ) ( +s-x+1 )( +t-s )

仔细观察发现可以发现它可以拆成两部分

1									   m
|--------------------------------------|
( +t-s ) 1 s s+1 t t+1 m
|------------|------------|------------|
( +0 ) ( +t-x+1 )( +0 )

  1. 整体加上\(t-s\)
  2. 对\([s+1,t]\)一段加上一个递减的序列。

整体的加直接记录即可,对一段进行的修改怎么办呢?

仔细观察这个递减的序列,发现其是\({0,-1,-2,...,s-t+1}\),于是可以借助前缀和\(O(1)\)完成这个修改。

// 详见后文

于是我们的优化就完成了,成功从\(O(n*m)\)变成\(O(n)\).

对于\(s\gt t\)的情况,只需将\(s\)加上\(m\)考虑即可,就相当于将每个位置拆成两个。

前缀和大法

对区间\([s,t]\)中的元素分别加上\({1,2,...,n}\),要求操作后的序列,方法为:

a[s] += 1, a[t+1] -= (n+1), a[t+2] += n;

求两次前缀和即得操作后的序列。容易验证。

Code

#include <bits/stdc++.h>
#define maxn 200010
using namespace std;
typedef long long LL;
LL c[maxn], a[maxn];
int main() {
int n, m;
scanf("%d%d", &n, &m);
for (int i = 0; i < n; ++i) scanf("%lld", &a[i]);
LL base = 0;
for (int i = 1; i < n; ++i) {
LL s = a[i-1], t = a[i];
if (t < s) t += m;
base += t - s;
c[s+2] += 1, c[t+1] -= t-s, c[t+2] += t-s-1;
}
for (int i = 1; i <= (m<<1); ++i) c[i] += c[i-1];
for (int i = 1; i <= (m<<1); ++i) c[i] += c[i-1];
LL maxx = 0;
for (int i = 1; i <= m; ++i) maxx = max(maxx, c[i] + c[i+m]);
printf("%lld\n", base - maxx);
return 0;
}

AtCoder Regular Contest 077 E - guruguru 线性函数 前缀和的更多相关文章

  1. AtCoder Regular Contest 077 E - guruguru

    https://arc077.contest.atcoder.jp/tasks/arc077_c 有m个点围成一个圈,按顺时针编号为1到m,一开始可以固定一个位置x,每次操作可以往顺时针方向走一步或直 ...

  2. AtCoder Regular Contest 077 被虐记&题解

    直到\(7:58\)才知道今天\(8:00\)有\(AtCoder\)的菜鸡来写题解啦. C - pushpush 题目: 给定一个长为\(n\)的序列,第\(i\)次操作做如下的事 : 将\(a_i ...

  3. AtCoder Regular Contest 077 D - 11

    题目链接:http://arc077.contest.atcoder.jp/tasks/arc077_b Time limit : 2sec / Memory limit : 256MB Score ...

  4. AtCoder Regular Contest 077 C - pushpush

    题目链接:http://arc077.contest.atcoder.jp/tasks/arc077_a Time limit : 2sec / Memory limit : 256MB Score ...

  5. AtCoder Regular Contest 077

    跟身在国外的Marathon-fan一起打的比赛,虽然最后没出F但还是涨分了. C - pushpush 题意:n次操作,每次往一个序列后面塞数,然后把整个序列翻转. #include<cstd ...

  6. 【arc077f】AtCoder Regular Contest 077 F - SS

    题意 给你一个形如"SS"的串S,以及一个函数\(f(x)\),\(x\)是一个形如"SS"的字符串,\(f(x)\)也是一个形如"SS"的 ...

  7. AtCoder Regular Contest 096

    AtCoder Regular Contest 096 C - Many Medians 题意: 有A,B两种匹萨和三种购买方案,买一个A,买一个B,买半个A和半个B,花费分别为a,b,c. 求买X个 ...

  8. AtCoder Regular Contest 098

    AtCoder Regular Contest 098 C - Attention 题意 给定一个只包含"E","W"字符串,可以花一的花费使他们互相转换.选定 ...

  9. AtCoder Regular Contest 061

    AtCoder Regular Contest 061 C.Many Formulas 题意 给长度不超过\(10\)且由\(0\)到\(9\)数字组成的串S. 可以在两数字间放\(+\)号. 求所有 ...

随机推荐

  1. 整合mybatis和spring时 Error creating bean with name 'sqlSessionFactory' defined in class path resource

    今天在整合mybatis和spring的时候出的错 报错如下 Exception in thread "main" org.springframework.beans.factor ...

  2. MongDB之各种删除操作

    接口IMongDaoDelete: package com.net.test.mongdb.dao; public interface IMongDaoDelete { public void del ...

  3. Flask初学者:配置文件

    如果设置项比较少的话可以使用“app.config['param_name']=value”的形式直接使用,如果需要设置的参数比较多的话,可以单独新建一个配置文件用来存放配置信息,配置文件中的参数需大 ...

  4. CentOS6.5生产环境系统安装

    CentOS 6.5系统安装 1-1 将预先准备的CentOS 6.5安装光盘插入光驱中,开机/重启系统时,系统会进行自检,自检完毕就会出现安装系统时的引导界面,如图1-1所示.1-2 使用键盘方向键 ...

  5. HDU - 5017 Ellipsoid(模拟退火)

    题意 给一个三维椭球面,求球面上距离原点最近的点.输出这个距离. 题解 模拟退火. 把\(z = f(x, y)\)函数写出来,这样通过随机抖动\(x\)和\(y\)坐标就能求出\(z\). 代码 / ...

  6. (PowerDesigner&Sqlite)PD中设计完表后,将其导入数据库中

    本人连接过SQLServer跟SQLite Ⅰ.SQLServer,百度,转一下:http://jingyan.baidu.com/article/7f766daf465e9c4101e1d0d5.h ...

  7. cf965c Greedy Arkady

    呸,大傻逼题,我更傻逼ref #include <iostream> using namespace std; typedef long long ll; ll n, k, m, d, a ...

  8. 矩阵儿快速幂 - POJ 3233 矩阵力量系列

    不要管上面的标题的bug 那是幂的意思,不是力量... POJ 3233 Matrix Power Series 描述 Given a n × n matrix A and a positive in ...

  9. 解决pymysql不能实时查询最新的数据

    #在网上查询到的原因为: InnoDB 的默认隔离级别.它可以防止任何被查询的行被其他事务更改,从而阻止不可重复的读取,而不是 幻读取.它使用中度严格的锁定策略,以便事务内的所有查询都会查看同一快照中 ...

  10. Python+Selenium练习篇之6-利用class name定位元素

    有时候,我们在用firepath(不会的请点这里)查看元素的XPath信息,发现没有可以用来定位的id信息,这个时候我们就需要考虑用其他的可用的来定位元素.本文介绍如何通过元素节点中class nam ...