UOJ #188. 【UR #13】Sanrd

Description

给定 $\sum_{i=l}^r f[i]$
$f[i]=$ 把 $i$ 的每一个质因子都从小到大排列成一个序列($p_i^{c_i}$要出现 $c_i$ 次)后 , 第二大的质因子.
题面

Solution

符合 $Min25$ 筛的处理顺序.
递归处理每个质因子作为次大值时的贡献,和不作为次大值时贡献的方案数 , 预处理一下区间质数个数就行了.

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+10;
ll w[N],c[N],n,p[N];int m=0,sqr,cnt=0;
inline int id(ll x){return x<=sqr?x:m-(n/x)+1;}
inline ll S(ll n,int k){if(n<=p[k])return 0;ll ret=p[k]*(c[id(n)]-c[p[k]]);for(int i=k+1;i<=cnt && n/p[i]>=p[i];i++){for(ll j=n;(j/=p[i])>=p[i];)ret+=S(j,i)+p[i];}return ret;
}
inline ll solve(ll now){m=0,cnt=0,n=now,sqr=sqrt(n);for(ll i=1,j;i<=n;i=j+1)w[++m]=j=n/(n/i),c[m]=j-1;for(int i=2;i<=sqr;i++){if(c[i]==c[i-1])continue;p[++cnt]=i;for(int j=m;w[j]/i>=i;j--)c[j]-=c[id(w[j]/i)]-c[i-1];}return S(n,0);
}
int main(){freopen("pp.in","r",stdin);freopen("pp.out","w",stdout);ll l,r;cin>>l>>r;cout<<solve(r)-solve(l-1);return 0;
}

转载于:https://www.cnblogs.com/Yuzao/p/9270403.html

UOJ #188. 【UR #13】Sanrd相关推荐

uoj#188. 【UR #13】Sanrd（Min_25筛）
题面传送门题解这是一道语文题不难看出,题目所求即为$l$到$r$中每个数的次大质因子我们考虑$Min\_25$筛的过程,设 \[S(n,j)=\sum_{i=1}^nsec_p( ...
uoj #118. 【UR #8】赴京赶考水题
#118. [UR #8]赴京赶考 Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://uoj.ac/problem/118 Description ...
UOJ#31 【UR #2】猪猪侠再战括号序列
传送门http://uoj.ac/problem/31 大家好我是来自百度贴吧的_叫我猪猪侠,英文名叫_CallMeGGBond. 我不曾上过大学,但这不影响我对离散数学.复杂性分析等领域的兴趣:尤其 ...
UOJ #31 【UR #2】猪猪侠再战括号序列
大家好我是来自百度贴吧的_叫我猪猪侠,英文名叫_CallMeGGBond. 我不曾上过大学,但这不影响我对离散数学.复杂性分析等领域的兴趣:尤其是括号序列理论,一度令我沉浸其中,无法自拔.至于OI算法 ...
UOJ#370. 【UR #17】滑稽树上滑稽果动态规划
原文链接www.cnblogs.com/zhouzhendong/p/UOJ370.html 题解首先易知答案肯定是一条链,因为挂在链的最下面肯定比挂在其他节点上赚. 问题被转化成了从一个集合中不断 ...
[UOJ #167]【UR #11】元旦老人与汉诺塔
题目大意:给你一个有$n$个盘子的汉诺塔状态$S$,问有多少种不同的操作方法,使得可以在$m$步以内到达状态$T$.$n,m\leqslant100$ 题解:首先可以知道的是,一个状态最多可以转移到其 ...
uoj#119. 【UR #8】决战圆锥曲线（线段树+复杂度分析）
题解传送门题解然而要我来说我感觉只是个爆搜啊-- //minamoto #include<bits/stdc++.h> #define R register #define ll l ...
UOJ#192. 【UR #14】最强跳蚤
题目大意:给定一颗树,每条边有一个权值,求有多少有序点对使得这两点间路径权值乘起来是完全平方数我们可以给每个素数随机一个权值,然后把每条边的权值分解质因数,把每个素数换成对应的权值然后异或起来,这样 ...
UOJ#31. 【UR #2】猪猪侠再战括号序列 splay
显然可以直接凑左面全是左括号,右面全是右括号的情况. 然后区间翻转就用 splay 模拟好了. splay 的时候一定注意一点: 如果没有调用 find(x),就一定要手动把 x 及其祖先的节点 pu ...

UOJ #188. 【UR #13】Sanrd

Description

Solution

UOJ #188. 【UR #13】Sanrd相关推荐

最新文章

热门文章