bzoj3309 DZY Loves Math

最新推荐文章于 2019-02-09 11:14:48 发布

Euryale_

最新推荐文章于 2019-02-09 11:14:48 发布

阅读量369

点赞数

分类专栏：数论文章标签：线性筛莫比乌斯反演

本文链接：https://blog.csdn.net/u010336344/article/details/67633650

版权

数论专栏收录该内容

25 篇文章 0 订阅

订阅专栏

题目链接：bzoj3309
题目大意：
对于正整数n，定义f(n)为n所含质因子的最大幂指数。给定正整数a,b，求 $\sum_{i=1}^{a}\sum_{j=1}^{b}f(gcd(i,j))$ 。T<=10000，1<=a,b<=10^7

题解：
莫比乌斯反演+线性筛
sb的我一开始柿子都画错了
要求

a n s = \sum i = 1 a \sum j = 1 b f (g c d (i, j))

$ans=\sum_{i=1}^{a}\sum_{j=1}^{b}f(gcd(i,j))$
即

\sum i = 1 a \sum j = 1 b \sum d | i, d | j, g c d (i, j) = d f (d)

$\sum_{i=1}^{a}\sum_{j=1}^{b}\sum_{d|i,d|j,gcd(i,j)=d}f(d)$
（蒟蒻就是漏掉了

gcd(i,j)=d $gcd(i,j)=d$ QWQ）
即

\sum d = 1 m i n (a, b) \sum i = 1 ⌊ a d ⌋ \sum j = 1 ⌊ b d ⌋ f (d) [g c d (i, j) = 1]

$\sum_{d=1}^{min(a,b)}\sum_{i=1}^{\lfloor\frac{a}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{b}{d}\rfloor}f(d)[gcd(i,j)=1]$
反演一下：

\sum d = 1 m i n (a, b) f (d) \sum i = 1 ⌊ a d ⌋ \sum j = 1 ⌊ b d ⌋ \sum t | i, t | j μ (t)

$\sum_{d=1}^{min(a,b)}f(d)\sum_{i=1}^{\lfloor\frac{a}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{b}{d}\rfloor}\sum_{t|i,t|j}\mu(t)$
即

\sum d = 1 m i n (a, b) f (d) \sum t = 1 m i n (⌊ a d ⌋, ⌊ b d ⌋) μ (t) \times ⌊ a d t ⌋ \times ⌊ b d t ⌋

$\sum_{d=1}^{min(a,b)}f(d)\sum_{t=1}^{min(\lfloor\frac{a}{d}\rfloor,\lfloor\frac{b}{d}\rfloor)}\mu(t)\times \lfloor\frac{a}{dt}\rfloor\times \lfloor\frac{b}{dt}\rfloor$
修改求和指标，枚举T=dt
则式子变成

\sum T = 1 m i n (a, b) ⌊ n T ⌋ \times ⌊ m T ⌋ \sum d | T f (d) \times μ (T d)

$\sum_{T=1}^{min(a,b)}\lfloor\frac{n}{T}\rfloor\times \lfloor\frac{m}{T}\rfloor\sum_{d|T}f(d)\times \mu(\frac{T}{d})$
嗯..然后我就不会了。
一开始打错的时候以为只用预处理出

f(i) $f(i)$ ，而现在要搞

∑d|Tf(d)×μ(Td) $\sum_{d|T}f(d)\times \mu(\frac{T}{d})$ 啊><。
接下来的看这个吧。我觉得写得最详细了。嗯还是新出炉的。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;
#define N 10000100

bool ispri[N];LL f[N];int ct[N],g[N],pri[N/4],cnt;
int mymin(int x,int y){return (x<y)?x:y;}
void pre(int lim)
{
    cnt=0;f[1]=0;ct[1]=0;
    for (int i=2;i<=lim;i++)
    {
        if (!ispri[i]) {pri[++cnt]=i;f[i]=1;ct[i]=1;g[i]=i;}
        for (int j=1;j<=cnt && i*pri[j]<=lim;j++)
        {
            int k=i*pri[j];
            ispri[k]=true;
            if (i%pri[j]==0)
            {
                ct[k]=ct[i]+1;
                g[k]=g[i]*pri[j];
                int tmp=i/g[i];
                if (tmp==1) f[k]=1;
                else f[k]=(ct[tmp]==ct[k])?-f[tmp]:0;
                break;
            }
            ct[k]=1;g[k]=pri[j];
            f[k]=(ct[i]==1)?-f[i]:0;
        }
    }
    for (int i=2;i<=lim;i++) f[i]+=f[i-1];
}
int main()
{
    //freopen("a.in","r",stdin);
    //freopen("a.out","w",stdout);
    int T,lim,n,m,i,r;LL ans;
    scanf("%d",&T);
    lim=10000000;pre(lim);
    while (T--)
    {
        scanf("%d%d",&n,&m);
        lim=mymin(n,m);ans=0;
        for (i=1;i<=lim;i=r+1)
        {
            r=mymin(n/(n/i),m/(m/i));
            ans+=(f[r]-f[i-1])*(LL)(n/i)*(LL)(m/i);
        }
        printf("%lld\n",ans);
    }
    return 0;
}