bzoj3309 DZY Loves Math

最新推荐文章于 2019-01-10 16:22:00 发布

sdfzyhx

最新推荐文章于 2019-01-10 16:22:00 发布

阅读量568

点赞数

分类专栏：数学 bzoj

本文链接：https://blog.csdn.net/sdfzyhx/article/details/72854335

版权

bzoj 同时被 2 个专栏收录

251 篇文章 0 订阅

订阅专栏

数学

229 篇文章 0 订阅

订阅专栏

不妨设 $n\ge m$ 。

= = \sum i = 1 n \sum j = 1 m f (gcd (i, j)) \sum d = 1 n f (d) \sum x = 1 ⌊ n d ⌋ μ (x) ⌊ n d x ⌋ ⌊ m d x ⌋ \sum D = 1 n \sum d | D f (d) μ (D d) ⌊ n d x ⌋ ⌊ m d x ⌋

$\begin{align} & \sum_{i=1}^n\sum_{j=1}^mf(\gcd(i,j)) \\ = & \sum_{d=1}^nf(d)\sum_{x=1}^{\lfloor\frac n d\rfloor}\mu(x)\lfloor\frac n{dx}\rfloor\lfloor\frac m{dx}\rfloor \\ = & \sum_{D=1}^n\sum_{d|D}f(d)\mu(\frac D d)\lfloor\frac n{dx}\rfloor\lfloor\frac m{dx}\rfloor \end{align}$
记

g=f∗μ $g=f*\mu$ ，只要我们线性筛出

g $g$ ，就可以

O(n+Tn‾‾√) $O(n+T \sqrt n)$ 解决问题。
设

g(n)=∑d|nf(d)μ(nd) $g(n)=\sum_{d|n}f(d)\mu(\frac n d)$ ，

n=∏ki=1pqii $n=\prod _{i=1}^kp_i^{q_i}$ ，

d=∏ki=1pixi $d=\prod_{i=1}^kp_i{x_i}$ 。显然

xi=qi $x_i=q_i$ 或者

xi=qi−1 $x_i=q_i-1$ ，否则

μ(nd) $\mu(\frac nd)$ 一定为

0 $0$ 。如果

qi $q_i$ 不全相同，把

pi $p_i$ 分成

qi $q_i$ 最大和

qi $q_i$ 不是最大两部分，那么

f(d) $f(d)$ 只取决于前一部分的取值。而对于前一部分的每种选法，我们只关心后一部分选择的

xi=qi−1 $x_i=q_i-1$ 的项的个数的奇偶。显然选奇数个和选偶数个的方案数是相同的，所以

g(n)=0 $g(n)=0$ 。而当

qi $q_i$ 全相同时，我们也可以得到类似的结论，只在所有的

xi $x_i$ 均取到

qi−1 $q_i-1$ 时除外，这时

f(d)=qi−1 $f(d)=q_i-1$ （其他时候都有

f(d)=qi $f(d)=q_i$ ）。结合

μ $\mu$ ，我们得到

g(n)=(−1)k+1 $g(n)=(-1)^{k+1}$ 。于是我们在线性筛的时候维护每个数最小质因数的指数和幂，不难讨论一下求出答案。

#include<cstdio>
#include<algorithm>
using namespace std;
#define LL long long
const int maxn=10000010;
int mnp[maxn],mnq[maxn],g[maxn],sg[maxn],prm[maxn],tot;
void init()
{
    for (int i=2;i<=maxn;i++)
    {
        if (!mnp[i])
        {
            mnp[i]=i;
            mnq[i]=1;
            g[i]=1;
            prm[++tot]=i;
        }
        for (int j=1;j<=tot&&(LL)i*prm[j]<=maxn;j++)
            if (i%prm[j])
            {
                mnp[i*prm[j]]=prm[j];
                mnq[i*prm[j]]=1;
                if (mnq[i]==1) g[i*prm[j]]=-g[i];
            }
            else
            {
                mnp[i*prm[j]]=mnp[i]*prm[j];
                mnq[i*prm[j]]=mnq[i]+1;
                if (i/mnp[i]==1) g[i*prm[j]]=1;
                else if (mnq[i/mnp[i]]==mnq[i*prm[j]]&&mnp[i/mnp[i]]!=mnp[i*prm[j]]) g[i*prm[j]]=-g[i/mnp[i]];
                break;
            }
        sg[i]=sg[i-1]+g[i];
    }
}
LL solve(int n,int m)
{
    LL ret=0;
    for (int i=1,j;i<=n&&i<=m;i=j+1)
    {
        j=min(n/(n/i),m/(m/i));
        ret+=(LL)(sg[j]-sg[i-1])*(n/i)*(m/i);
    }
    return ret;
}
int main()
{
    init();
    int T,n,m;
    scanf("%d",&T);
    while (T--)
    {
        scanf("%d%d",&n,&m);
        printf("%lld\n",solve(n,m));
    }
}