[SDOI2014]数表

2022-02-08 11:31:48 阅读：144 来源： 互联网

标签：lfloor frac limits sum rfloor SDOI2014 数表 sigma

题意：数表
有一个表，其中\(val(i,j)=\sigma(gcd(i,j))\)给出\(n,m,a\)求所有\(val\)不超过\(a\)的和。
思路：
假如没有\(a\)的限制：
\(\sum\limits_{d=1}^n\sigma(d)\sum\limits_{i=1}^n\sum\limits_{j=1}^m[gcd(i,j)==d]\)
\(\sum\limits_{d=1}^n\sigma(d)\sum\limits_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum\limits_{j=1}^{\lfloor\frac{m}{d}\rfloor}[gcd(i,j)==d]\)
\(\sum\limits_{d=1}^n\sigma(d)\sum\limits_{k=1}^{\lfloor\frac{n}{d}\rfloor}\lfloor\frac{n}{k*d}\rfloor\lfloor\frac{m}{k*d}\rfloor\mu(k)\)
令\(T=k*d\)
\(\sum\limits_{T=1}^n\lfloor\frac{n}{T}\rfloor\lfloor\frac{m}{T}\rfloor\sum\limits_{d|T}\sigma(d)\mu(\frac{T}{d})\)

为什么不把\(val\)带进去枚举\(gcd(i,j)\)的因数，因为我们还是要保留完整的\(\sigma\)，才能处理后面的\(a\)

令\(g(T)=\sum\limits_{d|T}\sigma(d)\mu(\frac{T}{d})\)
发现这里跟询问的\(n,m\)无关，于是考虑离线按\(a\)排序动态加入\(sigma(d)\),更新\(g(T)\)枚举倍数\(T\)即可。加入复杂度也是\(O(nlogn)\)不过考虑到数论分块要求前缀和+修改我们便用树状数组维护\(g()\)，复杂度\(O(nlog^2n)\)。

code

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+5;
const ll mod=(1ll)<<31;
ll as[N],sg[N],sm[N],c[N],g[N];
int ptot,p[N],miu[N];
bool is_p[N];
struct node {int n,m,a,id;}Q[N];
struct data {int p;ll val;}P[N];
bool cmp2(data u,data v) {return u.val<v.val;}
bool cmp(node u,node v) {return u.a<v.a;}
int lowbit(int u) {return u&(-u);}
void Update(int u,ll w) {for(;u<N;u+=lowbit(u))c[u]+=w;}
ll Sum(int u) {ll res=0;for(;u;u-=lowbit(u))res=(res+c[u])%mod;return res;}
void _xxs() {
    is_p[1]=miu[1]=sg[1]=sm[1]=1;
    for(int i=2;i<N;i++) {
        if(!is_p[i])p[++ptot]=i,sg[i]=1,sm[i]=i+1,miu[i]=-1;
        for(int j=1,x;j<=ptot&&(x=i*p[j])<N;j++) {
            is_p[x]=1;
            if(i%p[j]==0) {miu[x]=0;sg[x]=sg[i];sm[x]=sm[i]*p[j]+1;break;}
            miu[x]=miu[i]*-1;sg[x]=sg[i]*sm[i];sm[x]=p[j]+1;
        }
    }
    for(int i=2;i<N;i++) sg[i]=sg[i]*sm[i]%mod;
//  for(int i=1;i<=10;i++) printf("%lld ",sg[i]);puts("");
}
void init() {
    for(int i=1;i<N;i++)P[i]=(data){i,sg[i]};
    sort(P+1,P+N,cmp2);
}
int main() {
    _xxs();init();
    int q;scanf("%d",&q);
    for(int i=1;i<=q;i++)scanf("%d%d%d",&Q[i].n,&Q[i].m,&Q[i].a),Q[i].id=i;
    sort(Q+1,Q+1+q,cmp);
    int j=0;
    for(int i=1;i<=q;i++) {
        int n=Q[i].n,m=Q[i].m,a=Q[i].a;
        if(n>m)swap(n,m);
        while(j<N&&P[j+1].val<=a) {
            j++;ll phi=P[j].val,x=P[j].p;
            for(int k=x;k<N;k+=x) Update(k,phi*miu[k/x]%mod);
        }
        ll ans=0;
        for(int l=1,r;l<=n;l=r+1) {
            r=min(n/(n/l),m/(m/l));
            ans=(ans+1ll*(n/l)*(m/l)%mod*(Sum(r)-Sum(l-1)))%mod;
        }
        as[Q[i].id]=(ans+mod)%mod;
    }
    for(int i=1;i<=q;i++)printf("%lld\n",as[i]);
    return 0;
}

标签：lfloor,frac,limits,sum,rfloor,SDOI2014,数表,sigma
来源： https://www.cnblogs.com/bestime/p/15870600.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

[SDOI2014]数表