[FOI2020]楼房搭建

89 篇文章 0 订阅

题目

题目概要
有一个长度为 n n n 的数组,初始时每个位置都是 0 0 0 。每次操作可以让相邻的两个元素 + 1 +1 +1 + 2 +2 +2(可以是左边 + 1 +1 +1,也可以是右边)。问最少需要操作多少次,使得操作后第 i i i 个元素不小于阀值 h i h_i hi

数据范围与提示
n ⩽ 1 0 6 n\leqslant 10^6 n106 h i ⩽ 1 0 6 h_i\leqslant 10^6 hi106

思路

由于一次操作必然使得总和 + 3 +3 +3,问题变为最小化 “浪费”。为了成全 ( i + 1 ) (i{+}1) (i+1) 而牺牲 i i i 是不明智的:假设 i i i 爆表了,它可以去掉 ⟨ 1 , 2 ⟩ \langle 1,2\rangle 1,2,或者将 ⟨ 2 , 1 ⟩ \langle 2,1\rangle 2,1 换成 ⟨ 1 , 2 ⟩ \langle 1,2\rangle 1,2 。因此从左到右扫描,总是只需要将当前数字填满到 h i h_i hi

普通贪心是,同时让下一个数尽可能小,这样可以避免爆表。但这东西可以隔山打牛:下一个数太小,会导致下下个数过大。

那么只可能是 反悔贪心 嘛。显而易见的反悔方法是,把 ⟨ 2 , 1 ⟩ \langle 2,1\rangle 2,1 换成两个 ⟨ 1 , 2 ⟩ \langle 1,2\rangle 1,2,这可以让下一个数 + 3 +3 +3

但我们还要 对 “反悔操作” 进行反悔。考虑其影响:若不使用反悔操作,则下一个数的 + 3 +3 +3 需要用 3 × ⟨ 1 , 2 ⟩ 3\times\langle 1,2\rangle 3×1,2 ⟨ 1 , 2 ⟩ + ⟨ 2 , 1 ⟩ \langle 1,2\rangle+\langle 2,1\rangle 1,2+2,1 填补,这会使得下下个数 + 3 +3 +3 + 6 +6 +6 。也就是说:每撤销一个 + 3 +3 +3 操作,就可以让再下一个数 + 3 +3 +3 一或两次。

除了反悔操作,遵循下一个数尽量小的原则,优先用 ⟨ 2 , 1 ⟩ \langle 2,1\rangle 2,1 再用 ⟨ 1 , 2 ⟩ \langle 1,2\rangle 1,2 。按上面所说的, ⟨ 2 , 1 ⟩ \langle 2,1\rangle 2,1 可以带来一次反悔,而 ⟨ 1 , 2 ⟩ \langle 1,2\rangle 1,2 不可。

所以就做完了。正确性似乎是显然的——每一步,我们都尽量最小化了 “浪费”,并且给出了反悔的方式。

复杂度 O ( n ) \mathcal O(n) O(n),并且代码极其简单。

代码

代码实现的时候, + 6 +6 +6 看成两个 + 3 +3 +3 的机会就行了。

#include <cstdio>
#include <iostream>
#include <cstring>
using namespace std;
inline int readint(){
	int a = 0; char c = getchar(), f = 1;
	for(; c<'0'||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}

int main(){
	int n = readint();
	long long ans = 0;
	int v = 0; // already built
	int chance = 0; // how many repent is allowed
	for(int i=1; i<=n; ++i){
		int h = readint(); ans += h;
		h -= v, v = 0; // what's to do
		if(h <= 0){ // finished
			chance = 0; ans += (-h);
			continue; // waste -h
		}
		int x = min(h/3,chance); // how many +3 is applied
		int y = (h-3*x)>>1; // how many <2,1> is applied
		v = (((h-3*x)&1)<<1)+y; // if <1,2> is applied
		chance = (x<<1)+y; // +3 here can be twice +3 there
	}
	ans += v; // build on virtual n+1
	printf("%lld\n",ans/3);
	return 0;
}

后记

题解中提到了 “建隐式图”、“正权反悔边” 啥的,估计更精妙,我只是一个朝三暮四的人罢了。

另外,忘掉具体的操作,这一点也很重要。我在考虑反悔 + 3 +3 +3 时,固执地认为,只能有一次 ⟨ 1 , 2 ⟩ \langle 1,2\rangle 1,2,然后就把自己给绕进去了……反正,用 d p \tt dp dp 的话来说,不是状态定义的量,你管他咧

  • 2
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
由于您没有说明具体的EEG数据,因此我将提供一个示例数据集,并使用matlab中的eeglab和fieldtrip工具箱来实现EEG全通道数据的功率谱密度图像。 步骤如下: 1. 下载示例数据集:EEGlab提供了一些示例数据集,您可以从这里下载一个:https://sccn.ucsd.edu/eeglab/download.php。下载后将其解压缩到matlab当前工作目录中。 2. 启动EEGlab:在matlab命令窗口中输入eeglab,启动EEGlab工具箱。 3. 导入EEG数据:在EEGlab界面中,依选择“File” -> “Import data” -> “Using EEGLAB functions” -> “From ASCII/float file or MATLAB variable” -> “Sampled data” -> 选择示例数据集中的“s1.mat”文件。导入后,您可以在EEGlab的“EEG History”窗口中看到数据导入的过程。 4. 打开fieldtrip工具箱:在matlab命令窗口中输入ft_defaults,启动fieldtrip工具箱。 5. 将EEGlab的EEG数据转换为fieldtrip格式:在matlab命令窗口中输入ft_eeglab2fieldtrip(EEG),将EEGlab中的EEG数据转换为fieldtrip格式。 6. 计算功率谱密度图像:在matlab命令窗口中输入cfg = []; cfg.method = 'mtmfft'; cfg.output = 'pow'; cfg.taper = 'hanning'; cfg.foi = 1:30; cfg.keeptrials = 'yes'; freq = ft_freqanalysis(cfg, data),计算EEG全通道数据的功率谱密度图像。其中,cfg.method指定采用的方法为mtmfft,cfg.output指定输出功率谱密度,cfg.taper指定使用汉宁窗,cfg.foi指定感兴趣的频率范围为1-30Hz,cfg.keeptrials指定保留每个试验的功率谱密度。最后,将计算得到的功率谱密度存储在freq变量中。 7. 绘制功率谱密度图像:在matlab命令窗口中输入cfg = []; cfg.layout = 'eeg1010'; cfg.colorbar = 'yes'; cfg.xlim = [1 30]; cfg.ylim = [0 10]; ft_topoplotER(cfg, freq),绘制EEG全通道数据的功率谱密度图像。其中,cfg.layout指定使用的通道布局为eeg1010,cfg.colorbar指定绘制颜色条,cfg.xlim和cfg.ylim分别指定x轴和y轴的范围。 完整代码如下: % Step 1: download example dataset from EEGlab % Step 2: start EEGlab eeglab % Step 3: import EEG data % Choose File -> Import data -> Using EEGLAB functions -> From ASCII/float file or MATLAB variable -> Sampled data -> select s1.mat % Step 4: start fieldtrip toolbox ft_defaults % Step 5: convert EEGlab data to fieldtrip format data = ft_eeglab2fieldtrip(EEG); % Step 6: calculate power spectral density cfg = []; cfg.method = 'mtmfft'; cfg.output = 'pow'; cfg.taper = 'hanning'; cfg.foi = 1:30; cfg.keeptrials = 'yes'; freq = ft_freqanalysis(cfg, data); % Step 7: plot power spectral density cfg = []; cfg.layout = 'eeg1010'; cfg.colorbar = 'yes'; cfg.xlim = [1 30]; cfg.ylim = [0 10]; ft_topoplotER(cfg, freq);

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值