杜教筛学习笔记

QwQ一个不会莫比乌斯反演的蒟蒻来写杜教筛的博客了

这里写图片描述

这里写图片描述
这个是杜教筛的一般形式

中间那个先枚举几倍,实际上相当于把令 i=kd i = k ∗ d 然后进行k和d枚举

这么空说怎么好理解

我们来引入两道例题吧

51nod 51 n o d 莫比乌斯函数之和

i=1nμ(i) ∑ i = 1 n μ ( i )

一看这个题,貌似没什么头绪呀。

我们可以现推一下

因为 μ1=e μ ∗ 1 = e (或者写成 d|nμ(d)=[n==1] ∑ d | n μ ( d ) = [ n == 1 ]

所以

μ(n)=e(n)d|n,d!=nμ(d) μ ( n ) = e ( n ) − ∑ d | n , d ! = n μ ( d )

那么

ans=i=1n(e(i)d|i,d!=iμ(i)) a n s = ∑ i = 1 n ( e ( i ) − ∑ d | i , d ! = i μ ( i ) )

因为 ni=1e(i)=1 ∑ i = 1 n e ( i ) = 1

所以

ans=1i=1nd|i,d!=iμ(i) a n s = 1 − ∑ i = 1 n ∑ d | i , d ! = i μ ( i )

我们令 i=kd i = k ∗ d ,然后分别枚举k和d

ans=1k=2nd=1nkμ(i) a n s = 1 − ∑ k = 2 n ∑ d = 1 ⌊ n k ⌋ μ ( i )

这里k从2开始枚举的原因是因为 d!=i d ! = i

到这里我们能发现对于 nk n k 可以整除分块 且 nkd=1μ(i) ∑ d = 1 ⌊ n k ⌋ μ ( i ) 这个式子实际上是一个子问题,可以通过递归求值,只需要记忆化一下,就可以解决了

那么到这里,我们杜教筛的大致思路也就出来了

1.将一些小数的ans值筛出来,然后记忆化
2.对于一个数 x x ,我们可以进行分块,然后递归求解

直接上代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<unordered_map> 
using namespace std;

inline long long read()
{
  long long x=0,f=1;char ch=getchar();
  while (!isdigit(ch)) {if (ch=='-') f=-1;ch=getchar();}
  while (isdigit(ch)) {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
  return x*f;
}

const int maxn = 4700000;

int prime[maxn],check[maxn];
int mu[maxn];
int tot;
unordered_map<long long,int> mp;

void init(int n)
{
    mu[1]=1;
    check[1]=1;
    for (int i=2;i<=n;i++)
    {
        if (!check[i])
        {
            prime[++tot]=i;
            mu[i]=-1;
        }
        for (int j=1;j<=tot;j++)
        {
            if (i*prime[j]>n) break;
            check[i*prime[j]]=1;
            if (i%prime[j]==0) break;
            mu[i*prime[j]]=-mu[i];
        }
    }
    for (int i=2;i<=n;i++) mu[i]+=mu[i-1]; 
}

long long l,r;

long long dfs(long long x)
{
    if (x<=maxn) return mu[x];
    if (mp[x]) return mp[x];
    long long ans=1;
    for (long long i=2,j=0;i<=x;i=j+1)
    {
        j=x/(x/i);
        ans=ans-dfs(x/i)*(j-i+1);
    }
    mp[x]=ans;
    return ans;
}

int main()
{
  init(maxn);
  l=read(),r=read();
  cout<<dfs(r)-dfs(l-1); 
  return 0;
}

第二个问题是51nod 欧拉函数之和

其实这两个问题是差不多的

针对这个问题,我们要求的是

i=1nϕ(i) ∑ i = 1 n ϕ ( i )

还是同样
因为 id=1ϕ i d = 1 ∗ ϕ (或者写成 d|nϕ(d)=n ∑ d | n ϕ ( d ) = n

那么

ϕ(n)=idd|n,d!=nϕ(d) ϕ ( n ) = i d − ∑ d | n , d ! = n ϕ ( d )

那我们要求的

ans=i=1n(idd|n,d!=nϕ(d)) a n s = ∑ i = 1 n ( i d − ∑ d | n , d ! = n ϕ ( d ) )

ans=(n+1)n2i=1nd|i,d!=iϕ(d) a n s = ( n + 1 ) n 2 − ∑ i = 1 n ∑ d | i , d ! = i ϕ ( d )

i=kd i = k ∗ d
ans=(n+1)n2k=2nd=1nkϕ(d) a n s = ( n + 1 ) n 2 − ∑ k = 2 n ∑ d = 1 ⌊ n k ⌋ ϕ ( d )

ans=(n+1)n2k=2nans(nk) a n s = ( n + 1 ) n 2 − ∑ k = 2 n a n s ( ⌊ n k ⌋ )

然后就可以和上一道题一样的思路,直接做就好

上代码(注意取膜的时候的一些注意事项)

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#define ll long long
#include<unordered_map>
using namespace std;

inline ll read()
{
  ll x=0,f=1;char ch=getchar();
  while (!isdigit(ch)) {if (ch=='-') f=-1;ch=getchar();}
  while (isdigit(ch)) {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
  return x*f;
}

const int maxn = 6e6+1e2;
const ll mod = 1e9+7;

ll phi[maxn];
ll prime[maxn];
int check[maxn];
ll l,r;
int tot;

ll qsm(ll i,ll j)
{
    ll ans=1;
    while (j)
    {
        if (j&1) ans=ans*i%mod;
        i=i*i%mod;
        j>>=1;
    }
    return ans;
}

ll inv = qsm(2,mod-2);

unordered_map<long long,long long> mp,mp1;

void init(ll n)
{
    phi[1]=1;
    check[1]=1;
    for (int i=2;i<=n;i++)
    {
        if (!check[i])
        {
            prime[++tot]=i;
            phi[i]=i-1;
        }
        for (int j=1;j<=tot;j++)
        {
            if (i*prime[j]>n) break;
            check[i*prime[j]]=1;
            if (i%prime[j]==0)
            {
                phi[i*prime[j]]=phi[i]*prime[j]%mod;
                break;
            }
            else
            {
                phi[i*prime[j]]=phi[i]*(prime[j]-1)%mod;
            }
        }
    }
    for (int i=2;i<=n;i++) phi[i]=(phi[i]+phi[i-1])%mod;
}

ll dfs(ll x)
{
    if (x<=maxn) return phi[x]%mod;
    if (mp1[x]) return mp[x]%mod;
    long long ans=x%mod*(x%mod+1)%mod*inv%mod;
    for (ll i=2,j=0;i<=x;i=j+1)
    {
        j=x/(x/i);
        ans=(ans-(j-i+1)%mod*dfs(x/i)%mod+mod)%mod;
    }
    ans=(ans%mod+mod)%mod;
    mp[x]=ans;
    mp1[x]=1;
    return ans;
}

int main()
{
  init(maxn-10);
  l=read();
  cout<<dfs(l)%mod;
  return 0;
}
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值