【UOJ188】 Sanrd【类min_25筛】

题意:设f(i)f(i)f(i)表示iii的不严格次大质因子(没有为000),求∑i=lrf(i)\sum_{i=l}^rf(i)i=lrf(i)

l≤r≤1011l\leq r\leq10^{11}lr1011

这种和质因数有关的奇奇怪怪的函数的前缀和可以试试魔改min_25筛

S(n,j)=∑i=2n[minp(i)>pj]f(i)S(n,j)=\sum_{i=2}^n[minp(i)>p_j]f(i)S(n,j)=i=2n[minp(i)>pj]f(i)

枚举最小的质因子pkp_kpk以及次数eee

如果pkp_kpk不是次大质因子,直接递归到S(⌊npke⌋,k)S(\lfloor\frac{n}{p_k^e}\rfloor,k)S(pken,k)

否则考虑pkp_kpk的贡献

如果是严格次大,那么枚举最大的质因子

否则就是[e>1][e>1][e>1]

S(n,j)=∑k=j+1pk≤n∑e=1pke≤n(S(⌊npke⌋,k)+∑i=pk+1⌊npke⌋[i∈prime]+[e>1])S(n,j)=\sum_{k=j+1}^{p_k\leq\sqrt n}\sum_{e=1}^{p_k^e\leq n}(S(\lfloor\frac{n}{p_k^e}\rfloor,k)+\sum_{i=p_k+1}^{\lfloor\frac{n}{p_k^e}\rfloor}[i\in prime]+[e>1])S(n,j)=k=j+1pkne=1pken(S(pken,k)+i=pk+1pken[iprime]+[e>1])

注意pk+1p_k+1pk+1可能大于⌊npke⌋\lfloor\frac{n}{p_k^e}\rfloorpken,要特判一下

然后用min_25的方法预处理出质数个数的前缀和即可

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#include <cmath>
#define MAXN 1000005
using namespace std;
const int N=1e6;
int np[MAXN],pl[MAXN],cnt;
void init()
{np[1]=1;for (int i=2;i<=N;i++){if (!np[i]) pl[++cnt]=i;for (int j=1,x;(x=i*pl[j])<=N;j++){np[x]=1;if (i%pl[j]==0) break;}}
}
typedef long long ll;
ll val[MAXN],n,m;
int key[MAXN],yek[MAXN],tot;
inline int getkey(ll x){return x<=m? key[x]:yek[n/x];}
ll g[MAXN];
ll S(ll n,int j)
{if ((ll)pl[j]*pl[j]>n) return 0;ll sum=0;for (int k=j+1;(ll)pl[k]*pl[k]<=n&&k<=cnt;k++)for (ll e=1,v=pl[k];v<=n;e++,v*=pl[k])sum+=S(n/v,k)+pl[k]*(max(0ll,g[getkey(n/v)]-k)+(e>1));return sum;
}
ll solve(ll N)
{m=sqrt(n=N);tot=0;for (ll l=1,r;l<=n;l=r+1){r=n/(n/l);val[++tot]=n/l;if (val[tot]<=m) key[val[tot]]=tot;else yek[n/val[tot]]=tot;}for (int i=1;i<=tot;i++) g[i]=val[i]-1;for (int j=1;j<=cnt;j++)for (int i=1;(ll)pl[j]*pl[j]<=val[i];i++)g[i]-=g[getkey(val[i]/pl[j])]-j+1;return S(n,0);
}
int main()
{init();ll l,r;cin>>l>>r;cout<<solve(r)-solve(l-1);return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314989.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

小白开学Asp.Net Core 《六》 —— 探究.Net Core 跨平台的奥秘

1、写这篇文章的初衷有好多朋友反馈看不懂我写的开源的一个练手项目&#xff08;GitHub&#xff1a;https://github.com/AjuPrince/Aju.Carefree&#xff09;也有好多人都希望我能写一些简单的入门的文章&#xff0c;记得前几天在群里有人问为什么 .Net Core 能跨平台&#xff…

Codeforces Round #636 (Div. 3) D. Constant Palindrome Sum 思维 + 差分

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先有一个显然的性质就是每组操作最多不会超过两次。 很容易想到一个很暴力的思路&#xff0c;就是枚举x∈[1,2∗k]x \in [1,2*k]x∈[1,2∗k]&#xff0c;让后判断一下每组需要操作几次取…

C#8.0: 在 LINQ 中支持异步的 IAsyncEnumerableT接口

C# 8.0中&#xff0c;提供了一种新的IAsyncEnumerable<T>接口&#xff0c;在对集合进行迭代时&#xff0c;支持异步操作。比如在读取文本中的多行字符串时&#xff0c;如果读取每行字符串的时候使用同步方法&#xff0c;那么会导致线程堵塞。IAsyncEnumerable<T>可…

SCOI2020游记

考试期间发现只有T1看起来可做&#xff0c;脑子发麻硬刚5小时&#xff0c;最后还是只交了40分暴力 考完后用一节数学课口胡出了T2 100pts和T3 80pts&#xff0c;并且发现只有T1不会做 退役了&#xff0c;告辞

Codeforces Round #636 (Div. 3) E. Weights Distributing 思维 + bfs

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; n≤2e5,m≤2e5n\le2e5,m\le2e5n≤2e5,m≤2e5 思路&#xff1a; 怎么感觉每场div3div3div3都有一个巧妙的图论题。 首先如果只有两个点的话&#xff0c;肯定是一次bfsbfsbfs之后取前disdisdis小的www作为边权…

解读 Microsoft.NET.Sdk 的源码,你能定制各种奇怪而富有创意的编译过程

在 csproj 中&#xff0c;Project 中的 Sdk 属性是 MSBuild 15.0 开始支持的&#xff0c;也就是 Visual Studio 2017 开始支持。有了 Sdk 属性的存在&#xff0c;MSBuild 编译过程能够扩展得非常强大&#xff0c;而不止是过去 Import 的一个 props 和 targets 文件。本文将介绍…

【省选2020A卷】作业题【矩阵树】【扩域】【莫比乌斯反演】

传送门 为什么世界上会有这么傻的题啊……我佛了 很显然就是矩阵树强行缝合莫反 设f(n)f(n)f(n)表示所有边权gcd⁡\gcdgcd为nnn的生成树权值和&#xff0c;g(d)g(d)g(d)表示所有边权都是 ddd的倍数的生成树权值和 g(d)∑d∣nf(n)g(d)\sum_{d\mid n}f(n)g(d)d∣n∑​f(n) f(d…

Libre OJ 「BalticOI 2013」非回文数 数位dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 问[l,r][l,r][l,r]内有多少个数是非回文数&#xff0c;即数字中不存在连续几个数为回文数。 l,r≤1e18l,r\le1e18l,r≤1e18 思路&#xff1a; 这么大的范围很明显数位dpdpdp了&#xff0c;容易知道当一个数…

《ASP.NET Core 微服务实战》译者序

最近&#xff0c;我将《ASP.NET Core 微服务实战》一书由英文翻译为中文。这本书是由清华大学出版社引进的&#xff0c;目前还处于最后的排版校对过程中&#xff0c;现将该书的译者序发表于此。以下为译者译全文&#xff1a;“微服务”的概念在 2014 年正式提出之后&#xff0c…

【杭电多校2020】Distinct Sub-palindromes【结论】【构造】

题意&#xff1a;长度为nnn的 本质不同的回文子串个数最少 的小写字母组成的字符串 的个数。 n≤109n\leq 10^9n≤109 回文子串个数至少为不同字母的数量 显然 最少回文子串个数为min⁡(n,3)\min(n,3)min(n,3) n≤2n\leq2n≤2显然 n≥3n\geq3n≥3时 构造abcabcabc...\textt…

.NET 云原生技术使用调查

.NET Core 3.0 即将发布RC&#xff0c;.NET Core 3.0默认更好的支持Docker资源限制&#xff0c;官方团队也在努力让.NET Core成为真正的容器运行时&#xff0c;使其在低内存环境中具有容器感知功能并高效运行。GC堆限制.NET Core减少了CoreCLR默认使用的内存&#xff0c;如G0代…

【杭电多校2020】Minimum Index【Lyndon Word】

题意&#xff1a;给定字符串&#xff0c;求所有前缀的最小后缀。 n≤2107n\leq 2\times10^7n≤2107 最小后缀就是Lyndon分解的最后一段。而Duval本质上是可以重复修改的增量算法&#xff0c;所以是可以做的。 记ansians_iansi​为前缀iii的最小后缀。设维护未确定的循环节的指…

FZU - 2042 The Mad Mathematician 数位dp + 算贡献

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 求sumsumsum。 a,b,c,d,e≤1e18a,b,c,d,e\le1e18a,b,c,d,e≤1e18 思路&#xff1a; 这是一篇无从考究的题解&#xff0c;因为fzu现在进不去。 看到这种题直接考虑数位dpdpdp&#xff0c;对于[A,B],[C,D][…

「PowerBI」分析服务多维数据结构重回关系数据库的一大绝招

在过往Excel催化剂定位的轻量级ExcelPowerBIDesktop的解决方案中&#xff0c;已经做过了近乎完美的PowerBIDesktop数据模型数据导出到Excel工作表的应用&#xff0c;这也是个人版数据应用的重要的刚需功能&#xff0c;由Excel催化剂免费奉献给PowerBI社区作完结。此篇再次带领大…

【杭电多校2020】Fibonacci Sum【斐波拉契通项】【推式子】

题意&#xff1a;设FiF_iFi​为斐波拉契数列&#xff0c;求 ∑i0N(FiC)k\sum_{i0}^N(F_{iC})^ki0∑N​(FiC​)k 模109910^991099 N,C≤1018,k≤105N,C\leq10^{18},k\leq10^5N,C≤1018,k≤105 把斐波拉契暴力拆开 FiC15[(152)iC−(1−52)iC]F_{iC}\frac 1 {\sqrt{5}}[(\frac {…

Codeforces Round #257 (Div. 1) D. Jzzhu and Numbers 高维前缀和 + 容斥

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 完全想不到容斥啊&#xff0c;看了半天也没看懂渍渍渍。 定义f[i]f[i]f[i]表示iii的超集个数&#xff0c;那么选择的方案就是2f[i]−12^{f[i]}-12f[i]−1了&#xff0c;因为不能一个不选所…

花5分钟时间来了解一下高性能网关Kong会有意外收获

前言前几天开源发布了 Kong.Net 项目&#xff0c;收到了大量园友的反馈&#xff0c;开源当天就突破了 100 个star &#xff0c;可喜可贺&#xff0c;但是从侧面也说明&#xff0c;我们 .NetCore 阵营真的非常需要拥抱开源&#xff0c;应该敞开心扉&#xff0c;集众家之长&#…

【杭电多校2020】Total Eclipse【贪心】【并查集】

题意&#xff1a;nnn个点mmm条边的无向图&#xff0c;每个点有一个正点权&#xff0c;每次选择一个连通子图&#xff0c;将里面的权值都减111。求所有点权为000的最小步数。 T≤10,n≤105,m≤2105T\leq 10,n\leq 10^5,m\leq2\times10^5T≤10,n≤105,m≤2105 考虑一个贪心&…

Codeforces Round #725 (Div. 3) G. Gift Set 二分

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 有两种物品分别有x,yx,yx,y个&#xff0c;每次可以从一个拿出aaa个&#xff0c;另一个拿出bbb个分成一组&#xff0c;问最多能分成多少组。 思路&#xff1a; 这个题有一个显然的单调性&#xff0c;所以二…

长沙.NET技术社区·设计到实现

方案最终夭折回到一年前的四月初&#xff0c;也是这样一个气象更新&#xff0c;春意盎然的一年中最好的时光。 在长沙.NET技术圈中&#xff0c;有一个不知名的QQ群&#xff0c;叫做放肆.NET技术交流群&#xff08;群号&#xff1a;44329330&#xff09;&#xff0c;在这个群里面…