UOJ #188. 【UR #13】Sanrd
2018-08-04 19:32
204 查看
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 ,c ,n,p ;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; }
相关文章推荐
- UOJ187[UR#13]Ernd(斜率优化DP)
- uoj#187. 【UR #13】Ernd
- UOJ#49. 【UR #3】铀仓库
- UOJ 48 [UR #3]核聚变反应强度
- UOJ 49 [UR #3]铀仓库
- UOJ#152. 【UR #10】汉诺塔
- UOJ【UR #12】实验室外的攻防战
- UOJ 74 [UR #6]破解密码
- 【UOJ】#241. 【UR #16】破坏发射台
- UOJ #60 [UR #5] 怎样提高智商
- uoj #152. 【UR #10】汉诺塔
- UOJ 13|【UER #1】跳蚤OS|字典树
- 【UOJ】#242. 【UR #16】破坏蛋糕
- UOJ 49 【UR #3】铀仓库
- [UOJ]#33. 【UR #2】树上GCD
- UOJ #242. 【UR #16】破坏蛋糕
- 【UOJ169】【UR #11】元旦老人与数列
- uoj #118. 【UR #8】赴京赶考 水题
- UOJ 13 [UER #1]跳蚤OS
- UOJ 241. 【UR #16】破坏发射台 [矩阵乘法]