U138097 小鱼吃大鱼 埃氏筛
题目描述
小P同学在养殖一种非常凶狠的鱼,而且与其他鱼类不同,这种鱼越大越温顺,反而小鱼最凶残。当两条鱼相遇时, 小鱼会不断撕咬大鱼,每一口都咬下与它自身等重的肉(小鱼保持体重不变),直到大鱼的体重小于这条小鱼(若 两条鱼体重相同,一条鱼会将另一条撕咬殆尽)。
现在池塘中有n条鱼,小P想知道哪一对鱼相遇后,被咬的鱼剩余体重最大。
输入格式
单组测试数据。
第一行包含一个整数n,表示鱼的数量。(1 ≤ n ≤ 2e6) 第二行有n个用空格分开的整数ai 表示第i条鱼的体重(1 ≤ ai ≤ 1e6)。
输出格式
输出一个整数代表结果。
输入输出样例
3
3 4 5
2
2
2 2
0
5
2 1 4 3 5
2
说明/提示
数据范围
对于35%的数据,1≤n≤10,1 ≤ ai ≤ 100
对于55%的数据,1≤n≤10000
对于100%的数据,1 ≤ n ≤ 2e6,1 ≤ ai ≤ 1e6
样例解释
当三条鱼的体重分别为3 4 5时,不同对鱼相遇的结果分别是{3,4}=1 {3,5}=2 {4,5}=1,所以只有第一条跟第三条鱼相 遇时,最后大鱼的体重最大,结果为2
题意
给定n个数,求一对数(i, j),使得i>j, 且i%j最大。
题解
考虑枚举每一个数, 筛出它倍数。
对于这幅图, c%x显然大于a,b的结果,容易想到:对于每一个x的倍数, 我们取小于它的第一个数更新答案
很容易想到通过线性筛来筛倍数,但仔细想并不可以。比如12, 使用线性筛每个数只会筛一次,而12作为2,3,4,6的倍数意义是不一样的:假设12左边是10, 作为3的倍数时答案为1,4的倍数时则是2.
于是可以想到通过埃氏筛来做。
考试的时候由于边界细节只有85pts, 要注意,如果使用桶枚举可能会出现:某个数在i*x 和 (i+1)*x 之间, 而(i+1)*x 超过了枚举边间则会漏掉, 具体处理见代码。
优化
如果不是数据太水, 5e5的数据可以把 O(nlog^2n)的埃氏筛瞬间卡炸。
仔细分析后,大多数算法会被卡的情况是不一样的
数据过大:埃氏筛常数小, 复杂度稳定,但在数据过大的情况下,多一个log的致命伤就体现出来了。
数据密集:二分虽然理论上还多一个log,但是在大多数情况下完虐桶+埃氏筛。但是如果没有去重,数据密集就会被卡飞。
数据松散:刚好相反,没有一个重复的数字会让不加优化的桶卡飞(我就是
优化方式: 把上面几种加起来就好了
综合了几位大佬的写法和自己的代码,得到了一种效果显著的优化。
- 用桶储存,这样可以降低一个二分的log, 省掉一个排序,随便还能去重。
- 预处理一个last数组,储存每个数的上一个数是什么,例如
last[c] = b, last[b] = c , 同样的last[3x] = c;
这样我们在查找3x左边的第一个数时, 直接找last[3x]即可。预处理方法见代码。
3. 反向枚举,判断一下现在枚举的这个数是否大于当前答案, 如果这个数都比答案小, 取模结果肯定更小。对答案没有贡献,直接通过。
4. 将i--改为i = last[i] , 朴素的桶在枚举时通过i--, 然后判断这个数是否存在。然而我们既然已经有last数组了,不如好好利用。直接跳到上一个存在的数显然更优。
通过这些优化以及部分细节,可以显著提升效率。 实现其实也很简单
#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;
#define N 1000006
#define rr register
int n, ans=0, las[N];
bool a[N]; //快读
int read(){
int num=0; char c=getchar();
while(c<'0'||c>'9') c=getchar();
while(c>='0' && c<='9') num = num*10+c-'0', c=getchar();
return num;
} int main(){
n=read();
int maxi=0;
for(rr int i=0; i<n; i++){
int l = read();
a[l] = 1;
if(a[l]) maxi = max(maxi, l);
} rr int ls=-1;
//预处理last数组
for(rr int i=0; i<N; i++){
las[i] = ls;
if(a[i]) ls=i;
}
//i = last[i] 优化近100ms
for(rr int i=las[N-1]; i>=1; i=las[i]){
//判断是否对答案无贡献,优化100ms
if(ans >= i) continue;
//枚举倍数
for(rr int j=2; i*j<N; j++){
if(las[j*i]<=(j-1)*i) continue;
ans = max(ans, las[j*i]-(j-1)*i);
}
//处理边界问题 , 优化近200ms, 没有这一句N需要开到2倍才能AC,
ans = max(ans, maxi%i);
}
printf("%d", ans); return 0;
}
U138097 小鱼吃大鱼 埃氏筛的更多相关文章
- CodeForces - 385C Bear and Prime Numbers (埃氏筛的美妙用法)
Recently, the bear started studying data structures and faced the following problem. You are given a ...
- cf1154G 埃氏筛应用
直接用埃氏筛也可以做,但是这题写起来有点恶臭.. 更加简单的写法是直接枚举gcd=k,然后里面再枚举一次i*k,即找到k两个最小的倍数,看起来复杂度很高,但其实也是埃氏筛的复杂度 因为每次枚举gcd, ...
- 「CF779B」「LOJ#10201.」「一本通 6.2 练习 4」Sherlock and His Girlfriend(埃氏筛
题目描述 原题来自:Codeforces Round #400 B. Sherlock 有了一个新女友(这太不像他了!).情人节到了,他想送给女友一些珠宝当做礼物. 他买了 nnn 件珠宝.第 iii ...
- [JXOI 2018] 游戏 解题报告 (组合数+埃氏筛)
interlinkage: https://www.luogu.org/problemnew/show/P4562 description: solution: 注意到$l=1$的时候,$t(p)$就 ...
- 埃氏筛优化(速度堪比欧拉筛) + 洛谷 P3383 线性筛素数 题解
我们一般写的埃氏筛消耗的时间都是欧拉筛的三倍,但是欧拉筛并不好想(对于我这种蒟蒻) 虽然 -- 我 -- 也可以背过模板,但是写个不会的欧拉筛不如写个简单易懂的埃氏筛 于是就有了优化 这个优化还是比较 ...
- 埃氏筛+线段树——cf731F
从2e5-1依次枚举每个数作为主显卡,然后分段求比它大的数的个数,这里的复杂度是调和级数ln2e5,即埃氏筛的复杂度.. #include<bits/stdc++.h> using nam ...
- 数论(8):min_25 筛(扩展埃氏筛)
min_25 筛介绍 我们考虑这样一个问题. \[ans=\sum_{i = 1}^nf(i)\\ \] 其中 \(1 \le n \le 10^{10}\) 其中 \(f(i)\) 是一个奇怪的函数 ...
- Sirni题解(最小生成树,埃氏筛)(继 Liang-梁)
目录 前言 题意 思路 一些建议 前言 本篇是对Liang-梁的Sirni(最小生成树,埃氏筛)的后继博客. 通篇原文:https://blog.csdn.net/qq_37555704/articl ...
- hdu.5212.Code(莫比乌斯反演 && 埃氏筛)
Code Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) Total Submi ...
随机推荐
- Redis 4.0.2分布式锁的Java实现
简介 Redis分布式锁算法有两种,一种是单个Redis实例下的,一种是多个Redis实例的Redlock算法. 官方推荐Redlock算法,但是这个算法需要比较多的Redis实例而且是完全互相独立, ...
- Python-统计序列中元素
问题1: 随机数列[12,5,8,7,8,9,4,8,5,...] 中出现次数最高的3个元素,他们出现的次数 问题2: 对某英文文章的单词,进行词频统计,找出出现次数最搞得10个单词,他们出现的次数是 ...
- sqli-labs第一关 详解
sqli-labs第一关 方法一:手工注入 方法二:sqlmap工具 两种方式,都可以学学,顺便学会用sqlmap,也是不错的.不多说,我们开始吧 方法一: 来到第一关,图上说我们需要一个数字的参数 ...
- c++模板中set(date st):t(st)中的:符号
转载:https://zhidao.baidu.com/question/618119741512344012.html 半角冒号是构造函数里的初bai始化列表 开始du的标识. 如楼上所述: set ...
- 《To B产品经理进阶》
一.沙漏哟:To B产品技术标准化(全网独家) 经济机器是怎样运行的(超级简单模式理解经济运行规律) <俞军产品方法论>(思维模型.交易模型.经济学.心理学) <深度思考六步法> ...
- IDEA中创建父子工程与maven打包Springboot聚合工程报错程序包不存在问题处理
公司新项目需使用java技术栈,便使用IDEA搭建了一个多SpringBoot项目的聚合工程,因为初次使用,遇到了很多问题,maven打包时各种报错,在网上查了好多终于解决了,为巩固记忆,特作此记录. ...
- SpringBoot整合Shiro+MD5+Salt+Redis实现认证和动态权限管理(上)----筑基中期
写在前面 通过前几篇文章的学习,我们从大体上了解了shiro关于认证和授权方面的应用.在接下来的文章当中,我将通过一个demo,带领大家搭建一个SpringBoot整合Shiro的一个项目开发脚手架, ...
- Django ORM 引发的数据库 N+1 性能问题
背景描述 最近在使用 Django 时,发现当调用 api 后,在数据库同一个进程下的事务中,出现了大量的数据库查询语句.调查后发现,是由于 Django ORM 的机制所引起. Django Obj ...
- 洛谷UVA524 素数环 Prime Ring Problem
标签:搜索与回溯 题目: 从1到20这20个数摆成一个环,要求相邻的两个数的和是一个素数. 算法分析: 非常明显,这是一道回溯的题目.从1开始,每个空位有20种可能,只要填进去的数合法:与前面的数不相 ...
- angular http 跨域访问
1.在跟目录中创建配置文件:proxy.config.json ,文件内容如下: { "/api": { "target": "http://192. ...