SDOI2015约数个数和

最新推荐文章于 2020-04-22 15:48:02 发布

转载最新推荐文章于 2020-04-22 15:48:02 发布 · 223 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/LiGuanlin1124/p/10045966.html

文章标签：

#数据结构与算法

本文深入探讨了数论中反演与整除分块算法的应用，通过具体实例展示了如何利用这些算法解决复杂问题。文章首先介绍了算法的基本原理，随后通过代码实现详细说明了算法的具体应用过程。

题目描述

题解：

有一个式子：

证明先不说了。

然后倒一波反演：

然后整除分块就好了。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define N 50050
#define ll long long
int t,n,m;
int pri[N],cnt,mu[N];
ll sum[N],f[N];
bool vis[N];
void get_mu()
{
    mu[1]=sum[1]=1;
    for(int i=2;i<=50000;i++)
    {
        if(!vis[i])
        {
            pri[++cnt]=i;
            mu[i]=-1;
        }
        for(int j=1;j<=cnt&&i*pri[j]<=50000;j++)
        {
            vis[i*pri[j]]=1;
            if(i%pri[j])mu[i*pri[j]]=-mu[i];
            else break;
        }
        sum[i]=sum[i-1]+mu[i];
    }
}
void get_f()
{
    for(int x=1;x<=50000;x++)
    for(int i=1,nxt;i<=x;i=nxt+1)
    {
        nxt=x/(x/i);
        f[x]+=1ll*(nxt-i+1)*(x/i);
    }
}
int main()
{
    scanf("%d",&t);
    get_mu();
    get_f();
    while(t--)
    {
        scanf("%d%d",&n,&m);
        if(n>m)swap(n,m);
        ll ans = 0;
        for(int i=1,nxt;i<=n;i=nxt+1)
        {
            nxt = min(n/(n/i),m/(m/i));
            ans+=1ll*(sum[nxt]-sum[i-1])*f[n/i]*f[m/i];
        }
        printf("%lld\n",ans);
    }
    return 0;
}