[51nod][积性函数][杜教筛]最大公约数之和 V3

最新推荐文章于 2022-01-17 20:21:30 发布

Gzb1128

最新推荐文章于 2022-01-17 20:21:30 发布

阅读量280

点赞数

CC 4.0 BY-SA版权

分类专栏：杜教筛积性函数

本文链接：https://blog.youkuaiyun.com/qq_36993218/article/details/79651495

积性函数同时被 2 个专栏收录

7 篇文章

订阅专栏

杜教筛

6 篇文章

订阅专栏

本文介绍了一道关于杜教筛算法的题目，并分享了作者对该题的解题思路及实现过程。通过反演和筛选phi函数，实现了高效的求解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1237

sol:

本来已经不想做裸题了，但是之前的题总是要调一下才能过，就很不优美
这题终于一遍过了，感觉差不多把调试的坑点给找完了。

杜教筛的精髓似乎也掌握到一点了。
简单反演

$\sum_d^nd\sum_i^{\frac nd}\sum_j^{\frac nd}(i,j)=1$
$\sum_d^nd(2(\sum_j^{\frac nd}\phi(i))-1)$
筛一下phi。这里的phi我已经不是用tangls的那种推导方式了，而是使用了套路的推导方式，感觉确实已经很通用了。

#include<iostream>
#include<string>
#include<cstring>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cmath>
using namespace std;
typedef long long ll;

ll n;

inline int read()
{
    char c;
    bool pd=0;
    while((c=getchar())>'9'||c<'0')
    if(c=='-') pd=1;
    int res=c-'0';
    while((c=getchar())>='0'&&c<='9')
    res=(res<<3)+(res<<1)+c-'0';
    return pd?-res:res;
}
const int N=1e6;
const int pyz=1e9+7;
int inv2;
int phi[N+7],s[N+7];
inline int ksm(int s,int t)
{
    int res=1;
    while(t)
    {
        if(t&1) res=(ll)res*s%pyz;
        s=(ll)s*s%pyz;
        t>>=1;
    }
    return res;
}
inline get_f(ll x)
{
    if(x<=N) return phi[x];
    if(s[n/x]) return s[n/x];
    s[n/x]=(ll)x%pyz*(x%pyz+1)%pyz*inv2%pyz;
    ll a,b;
    for(ll i=2;i<=x;++i)
    {
        a=x/i;
        b=x/a;
        s[n/x]=(s[n/x]-(ll)get_f(a)*((b-i+1)%pyz)%pyz+pyz)%pyz;
        i=b;
    }
    return s[n/x];
}
int prime[N];
bool is[N+7]; 
int main()
{
//  freopen("1237.in","r",stdin);
//  freopen(".out","w",stdout);
    cin>>n;
    inv2=ksm(2,pyz-2);
    ll a,b,ans=0;
    phi[1]=1;
    for(int i=2;i<=N;++i)
    {
        if(!is[i]) prime[++prime[0]]=i,phi[i]=i-1;
        for(int j=1;j<=prime[0]&&i*prime[j]<=N;++j)
        {
            is[i*prime[j]]=1;
            if(!(i%prime[j]))
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            phi[i*prime[j]]=phi[i]*(prime[j]-1);
        }
        phi[i]=(phi[i-1]+phi[i])%pyz;
    }
    for(ll i=1;i<=n;++i)
    {
        a=n/i;
        b=n/a;
        ans=(ans+(i+b)%pyz*((b-i+1)%pyz)%pyz*inv2%pyz*(2*get_f(a)%pyz-1)%pyz)%pyz;
        i=b;
    }
    ans=(ans+pyz)%pyz;
    cout<<ans;
}