3529: [Sdoi2014]数表|莫比乌斯反演|树状数组

最新推荐文章于 2020-09-11 10:58:10 发布

ws_yzy

最新推荐文章于 2020-09-11 10:58:10 发布

阅读量529

点赞数

分类专栏：莫比乌斯反演文章标签：莫比乌斯反演

本文链接：https://blog.csdn.net/ws_yzy/article/details/50916993

版权

莫比乌斯反演专栏收录该内容

6 篇文章 0 订阅

订阅专栏

以下式子均设 $n\leq m$
先设

M (x) = \sum d | x d

$M(x)=\sum_{d|x}d$
在不考虑a的限制下

A n s = \sum i = 1 n \sum j = 1 n M (g c d (i, j))

$Ans=\sum_{i=1}^{n}\sum_{j=1}^nM(gcd(i,j))$

= \sum i = 1 n M (i) * s u m (i)

$=\sum_{i=1}^{n}M(i)*sum(i)$

sum(i) $sum(i)$ 表示

gcd(x,y)=i $gcd(x,y)=i$ 的

(x,y) $(x,y)$ 的个数
显然这个反演一下就可以得到

s u m (i) = \sum x = 1 ⌊ n i ⌋ u (i) * ⌊ n i * x ⌋ * ⌊ m i * x ⌋

$sum(i)=\sum_{x=1}^{\lfloor\frac{n}{i}\rfloor}u(i)*\lfloor\frac{n}{i*x}\rfloor*\lfloor\frac{m}{i*x}\rfloor$

A n s = \sum i = 1 n M (i) * \sum x = 1 ⌊ n i ⌋ u (i) * ⌊ n i * x ⌋ * ⌊ m i * x ⌋

$Ans=\sum_{i=1}^{n}M(i)*\sum_{x=1}^{\lfloor\frac{n}{i}\rfloor}u(i)*\lfloor\frac{n}{i*x}\rfloor*\lfloor\frac{m}{i*x}\rfloor$

= \sum T = 1 n ⌊ n T ⌋ * ⌊ m T ⌋ * (\sum d | T u (d) * M (T d))

$=\sum_{T=1}^{n}\lfloor\frac{n}{T}\rfloor*\lfloor\frac{m}{T}\rfloor*\big(\sum_{d|T}u(d)*M(\frac{T}{d})\big)$
这样只需要线性筛搞出后面的这个

∑d|Tu(d)∗M(Td) $\sum_{d|T}u(d)*M(\frac{T}{d})$ 就可以

n−√ $\sqrt n$ 处理每次询问了
然后再回来考虑a的限制，只有

M(x)≤a $M(x)\leq a$ 才会对答案有贡献，所以可以对询问按

a $a$ 的大小排序，每次都把小于当前的a的

M(x) $M(x)$ 加入到树状数组中，这样离线处理每次询问
然后处理模数时中途直接全部

int $int$ 最后在

and $and$ 一下

int $int$ 的最大值就可以

#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<vector>
#include<cstdio>
#include<queue>
#include<cmath>
#include<set>
#include<map>
#define ll long long
#define N 100005
using namespace std;
int sc()
{
    int i=0,f=1; char c=getchar();
    while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
    return i*f;
}
struct W{int n,m,p,a;}a[N];
int low[N],A[N],prime[N],tr[N];
int M[N],u[N],ans[N],pos[N];
void pre()
{
    int top=0;M[1]=u[1]=1;
    for(int i=2;i<N;i++)
    {
        if(!A[i])
        {
            u[low[i]=prime[++top]=i]=-1;
            M[i]=i+1;
        }
        for(int j=1;i*prime[j]<N;j++)
        {
            A[prime[j]*i]=1;
            if(i%prime[j]==0)
            {
                low[i*prime[j]]=low[i]*prime[j];
                if(low[i]==i)
                    M[i*prime[j]]=M[i]+i*prime[j];
                else M[i*prime[j]]=M[i/low[i]]*M[low[i]*prime[j]];
                break;
            }
            u[prime[j]*i]=-u[i];
            low[prime[j]*i]=prime[j];
            M[prime[j]*i]=M[i]*(prime[j]+1);
        }
    }
}   
bool cmp(W a,W b){return a.a<b.a;}
bool Cmp(int a,int b){return M[a]<M[b];}
void change(int x,int f)
{
    for(;x<N;x+=x&-x)tr[x]+=f;
}
int ask(int x)
{
    int ans=0;
    for(;x;x-=x&-x)ans+=tr[x];
    return ans;
}
int query(int n,int m)
{
    if(n>m)swap(n,m);
    int ans=0;
    for(int i=1,last;i<=n;i=last+1)
    {
        last=min(n/(n/i),m/(m/i));
        ans+=(n/i)*(m/i)*(ask(last)-ask(i-1));
    }
    return ans;
}
int main()
{
    pre();
    int W=sc();
    for(int i=1;i<=W;i++)
        a[i].n=sc(),a[i].m=sc(),a[i].a=sc(),a[i].p=i;
    sort(a+1,a+W+1,cmp);
    for(int i=1;i<N;i++)pos[i]=i;
    sort(pos+1,pos+N,Cmp);
    int now=1;
    for(int i=1;i<=W;i++)
    {
        while(now<N&&M[pos[now]]<=a[i].a)
        {
            int x=pos[now++];
            for(int j=x;j<N;j+=x)
                change(j,M[x]*u[j/x]);
        }
        ans[a[i].p]=query(a[i].m,a[i].n);
    }
    for(int i=1;i<=W;i++)
        printf("%d\n",ans[i]&0x7FFFFFFF);
    return 0;
}