bzoj 4916: 神犇和蒟蒻杜教筛

最新推荐文章于 2023-01-12 15:05:38 发布

_beginend

最新推荐文章于 2023-01-12 15:05:38 发布

阅读量591

点赞数

分类专栏：杜教筛

本文链接：https://blog.csdn.net/qq_33229466/article/details/72868110

版权

杜教筛专栏收录该内容

11 篇文章 0 订阅

订阅专栏

题意

求 $A=\sum_{i=1}^N{\mu (i^2)}和B=\sum_{i=1}^N{\varphi (i^2)}$
n<=1e9,答案模1e9+7

分析

md第一问是什么鬼？？？
第二问想到是杜教筛（因为不会洲阁筛），但是没想到用哪个函数和f卷起来。

其实考虑到 $f(d)=\varphi(d)*d$ ，要把常数消掉的话，最好就是把f和 $g(d)=d$ 卷起来，然后就可以得到

\sum i = 1 n \sum j | i f (j) * g (i j) = \sum i = 1 n \sum j | i g (j) f (i j)

$\sum_{i=1}^n\sum_{j|i}f(j)*g(\frac{i}{j})=\sum_{i=1}^n\sum_{j|i}g(j)f(\frac{i}{j})$
容易发现左边等于

\sum i = 1 n i * \sum j | i φ (j) = \sum i = 1 n i 2 = n * ( n + 1 ) * ( 2 n + 1 ) 6

$\sum_{i=1}^ni*\sum_{j|i}\varphi(j)=\sum_{i=1}^ni^2=\frac{n*(n+1)*(2n+1)}{6}$
右边的话，化简一下可以得到

\sum d = 1 n g (d) * \sum i = 1 ⌊ n d ⌋ f (i) = \sum d = 1 n d * S (⌊ n d ⌋)

$\sum_{d=1}^ng(d)*\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}f(i)=\sum_{d=1}^nd*S(\lfloor\frac{n}{d}\rfloor)$
合并一下

S (n) = n * ( n + 1 ) * ( 2 n + 1 ) 6 - \sum d = 2 n d * S (⌊ n d ⌋)

$S(n)=\frac{n*(n+1)*(2n+1)}{6}-\sum_{d=2}^nd*S(\lfloor\frac{n}{d}\rfloor)$
直接上分块大法即可。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<map>
using namespace std;

typedef long long LL;

const int N=10000005;
const int MOD=1000000007;

int n,s[N],phi[N],prime[N/10],tot;
bool not_prime[N];
map<int,int> w;

void get_prime(int n)
{
    phi[1]=1;
    for (int i=2;i<=n;i++)
    {
        if (!not_prime[i]) prime[++tot]=i,phi[i]=i-1;
        for (int j=1;j<=tot&&i*prime[j]<=n;j++)
        {
            not_prime[i*prime[j]]=1;
            if (i%prime[j]==0)
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            phi[i*prime[j]]=phi[i]*(prime[j]-1);
        }
    }
    for (int i=1;i<=n;i++) s[i]=(s[i-1]+(LL)i*phi[i]%MOD)%MOD;
}

int solve(int n)
{
    if (n<=10000000) return s[n];
    if (w[n]) return w[n];
    int ans=0;
    for (int i=2,last;i<=n;i=last+1)
    {
        last=n/(n/i);
        ans=(ans+(LL)((LL)last*(last+1)/2-(LL)i*(i-1)/2)*solve(n/i)%MOD)%MOD;
    }
    ans=((LL)n*(n+1)%MOD*(n*2+1)%MOD*166666668%MOD-ans+MOD)%MOD;
    w[n]=ans;
    return ans;
}

int main()
{
    puts("1");
    get_prime(10000000);
    scanf("%d",&n);
    printf("%d",solve(n));
    return 0;
}