日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪(fǎng)問(wèn) 生活随笔!

生活随笔

當(dāng)前位置: 首頁(yè) > 编程资源 > 编程问答 >内容正文

编程问答

UOJ #188. 【UR #13】Sanrd

發(fā)布時(shí)間:2023/12/1 编程问答 41 豆豆
生活随笔 收集整理的這篇文章主要介紹了 UOJ #188. 【UR #13】Sanrd 小編覺(jué)得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.

Description

給定 \(\sum_{i=l}^r f[i]\)
\(f[i]=\)\(i\) 的每一個(gè)質(zhì)因子都從小到大排列成一個(gè)序列(\(p_i^{c_i}\)要出現(xiàn) \(c_i\) 次)后 , 第二大的質(zhì)因子.
題面

Solution

符合 \(Min25\) 篩的處理順序.
遞歸處理每個(gè)質(zhì)因子作為次大值時(shí)的貢獻(xiàn),和不作為次大值時(shí)貢獻(xiàn)的方案數(shù) , 預(yù)處理一下區(qū)間質(zhì)數(shù)個(gè)數(shù)就行了.

#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; }

轉(zhuǎn)載于:https://www.cnblogs.com/Yuzao/p/9270403.html

總結(jié)

以上是生活随笔為你收集整理的UOJ #188. 【UR #13】Sanrd的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò),歡迎將生活随笔推薦給好友。