题解 P7586 [COCI2012-2013#1] SNAGA

合集 - COCI(19)

1.题解 P7763 [COCI2016-2017#5] Ronald2023-08-26 2.题解 P6497 [COCI2016-2017#2] Prosječni2023-08-24 3.题解 P7537 [COCI2016-2017#4] Rima2023-08-25 4.题解 P8017 [COCI2013-2014#4] UTRKA2023-08-23 5.题解 P7751 [COCI2013-2014#2] PUTNIK2023-08-22 6.题解 P7752 [COCI2013-2014#2] PALETA2023-08-22

7.题解 P7586 [COCI2012-2013#1] SNAGA2023-08-17

8.题解 P6485 [COCI2010-2011#4] PROSJEK2023-08-12 9.题解 P6447 [COCI2010-2011#1] ŽABE2023-08-11 10.题解P6677 [COCI2019-2020#2] Checker2023-03-27 11.题解P8073 [COCI2009-2010#7] BAKICE2023-03-01 12.题解P8084 [COCI2011-2012#4] BROJ2023-02-26 13.题解P6370 [COCI2006-2007#6] KAMEN2022-12-29 14.题解 P7165 [COCI2020-2021#1] Papričice2024-01-01 15.题解 P9911 [COCI 2023/2024 #2] Kuglice2024-01-24 16.题解 P6356 [COCI2007-2008#3] CUDAK2024-01-29 17.题解 P7309 [COCI2018-2019#2] Kocka2024-01-30 18.题解 P6548 [COCI2010-2011#2] IGRA2024-01-30 19.题解 P6491 [COCI2010-2011#6] ABECEDA2024-01-31

洛谷。

题意

应该好理解，~~没理解的话可以看一下下面暴力的代码。~~

分析

观察这种题，范围 $3 \leq A < B \leq 10^{17}$ ，首先就只有几种算法了（至少我没学过多少种）：倍增，矩阵快速幂，数学，数位 DP，循环，或一些特殊情况。

而这题，他有一个不能整除的性质，看来只有数学来解决了，先确定了一个大的方向。

先打一个暴力，打个小表来分析一下性质：

 inline int js(int x) {
    if(x==2) return 1;
    for(int j=2;j<=x;++j) if(x%j) return 1+js(j);
}

输出前一百个观察一下：

再随机几个我们发现，所有的 $l e n$ 都是小于等于 $4$ 的。

由此，我们再次分析一下我们的操作，可以得到几个性质（令 $f_{i}$ 表示 $i$ 的最小的不能整除的数）：

奇数的 $l e n$ 必然为 $2$ 。
所有小于 $f_{i}$ 的数都是 $i$ 的因数。

看起来有点像废话？所有小于 $f_{i}$ 的数都是 $i$ 的因数，这也就导致，我们的前 $f_{i}$ 个数的 lcm。

输出一下前缀的 lcm：

当我们的 $i$ 大于 $42$ 时，就已经超过了我们的范围了，因此，我们的 $f_{i} \leq 42$ ，这是第一个衍生性质。

再分析一下 $f_{i}$ 在数字上的性质，进而判断 $l e n_{i}$ 的数值。

先给出结论：

$i$ 是一个奇数，那么 $l e n_{i} = 2$ ；
$f_{i}$ 是一个奇数，那么 $l e n_{i} = 3$ ；
$f_{i}$ 是一个偶数，那么 $l e n_{i} = 4$ 。

这三个结论中的前两个都应该是很显然的，着重分析一下第 3 个。

因为 $f_{i}$ 是一个偶数，我们令 $f_{i} = 2 \times k$ ，那么 $k$ 也就已经判断过，是 $i$ 的因数。
那为什么我们的 $2 \times k$ 不是其因数了呢？是因为我们增加了一个 $2$ 的因数。
因此，我们对其造成关键的实际上是 $2^{q}$ ，（ $q$ 为 $2 \times k$ 的 $2$ 的最大幂次）。
而 $2^{q} \leq 2 \times k$ ，由此，我们的 $f_{i}$ 应当就是 $2^{q}$ ，而 $f_{2^{q}} = 3$ 。由此，得出我们的结论 3。

最后轮到了如何求答案了，我们分开处理三种数，可能不太好算，我们令其一部分放在一起考虑。

我们令基础的贡献为 3，那么，我们只要减去奇数的个数，就可以处理出 1，2 的贡献，同时，计算了一部分的 3，我们只需要加上第三种数字的个数即可。

对于这一部分，我们又要分成 4 种，根据 $f_{i}$ 的种类分类，因为 $f_{i} \leq 42$ ， $f_{i} = 2^{q}$ ，所以此时的 $f_{i}$ 只有四种可能：4、8、16、32。

要使 $f_{i}$ 是 $x$ ，我们需要使 $i$ 使 $x$ 前面所有数的 lcm 的因数，同时不被 $x$ 整除，用一个小小容斥即可。

 for(int i=2;i<=5;++i) {
    res+=x/num[i];//num为预处理的前缀lcm
    ll y=num[i]*(1<<i)/__gcd(num[i],(1<<i));
    if(y<0) continue;//防止爆
    res-=x/y;
}

posted @ 2023-08-17 20:52 djh0314 阅读(11) 评论(0) 编辑收藏举报来源

刷新页面返回顶部

登录后才能查看或发表评论，立即登录或者逛逛博客园首页

相关博文：

· 题解 P7752 [COCI2013-2014#2] PALETA

· 题解 P6485 [COCI2010-2011#4] PROSJEK

· P8245 [COCI2013-2014#3] PAROVI & ZLOJ 练习62 D

· CF edu 173 刷题笔记

· 洛谷P10136 暨 USACOJan2024S T3 题解

阅读排行：
· 阿里最新开源QwQ-32B，效果媲美deepseek-r1满血版，部署成本又又又降低了！
· 单线程的Redis速度为什么快？
· SQL Server 2025 AI相关能力初探
· AI编程工具终极对决：字节Trae VS Cursor，谁才是开发者新宠？
· 展开说说关于C#中ORM框架的用法！

蒟蒻一枚

我在时光斑驳深处，聆听到花开的声音。

题解 P7586 [COCI2012-2013#1] SNAGA

题意

分析

公告

一言（ヒトコト）

常用链接

合集

随笔分类

随笔档案

	inline int js(int x) {
	if(x==2) return 1;
	for(int j=2;j<=x;++j) if(x%j) return 1+js(j);
	}

	for(int i=2;i<=5;++i) {
	res+=x/num[i];//num为预处理的前缀lcm
	ll y=num[i]*(1<<i)/__gcd(num[i],(1<<i));
	if(y<0) continue;//防止爆
	res-=x/y;
	}