P4213 杜教筛

首先引入这三个函数：

在引入三个常见的迪利克雷卷积：

下面进入正题，所谓杜教筛，其实就是一个求积性函数前缀和的定理，如果要求的前缀和数量太大，比如前1e9个数的和，就没法通过累加解决，所以就有了我们的杜教筛。杜教筛主要的核心就是一个公式：

这个题的题意是要求莫比乌斯函数u和欧拉函数phi的前缀和，n的大小在2^31以内。
如何利用杜教筛解决这个问题呢？
我们需要做的就是寻找合适的函数帮助我们凑出杜教筛的核心公式，S(n)毋庸置疑就是sigma_u和sigma_phi了，我们只需要找合适的函数来代替g和h即可，函数我也给出了，就是我前面给出的三个函数，那么这个题我采用的是b站大佬视频中采用的：

我在做这个题的时候遇到tle的问题，因为递归迭代的次数太多了，所以就听了实验室xyx巨佬的方法，把前1e7个phi和u用线性筛筛出来，求出前缀和保存在数组里，这样返回1e7以内的数的时候就可以直接返回，大于1e7的就记忆化搜索用哈希map保存，妙哉。
下面附上我的代码：

//#pragma GCC optimize(3,"Ofast","inline")
#include<unordered_map>
#include<unordered_set>
#include<cstdio>
#include<iostream>
#include<cmath>
#include<functional>
#include<cstring>
#include<string>
#include<cstdlib>
#include<queue>
#include<map>
#include<algorithm>
#include<set>
#include<stack>
#include<vector>
#include<sstream>
#include<list>
using namespace std;
typedef long long ll;
const int INF=0x3f3f3f3f;
const int maxn= 2e7+10;
const ll mod=6662333;
const int N=210000;
unordered_map<int,ll>m1;
unordered_map<int,ll>m2;
int prime[maxn],isprime[maxn];
ll phi[maxn];
int num;
ll mob[maxn];
void getphi(int n)
{
    memset(prime,0,sizeof(prime));
    memset(isprime,1,sizeof(isprime));
    memset(mob,0,sizeof(mob));
    phi[1]=1;
    mob[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(isprime[i])
        {
            phi[i]=i-1;
            mob[i]=-1;
            prime[num++]=i;
        }
        for(int j=0;j<num&&i*prime[j]<=n;j++)
        {
            isprime[i*prime[j]]=0;
            if(i%prime[j]==0){
                phi[i*prime[j]]=phi[i]*prime[j];
                mob[i*prime[j]]=0;
                break;
            }
            else
            {
                mob[i*prime[j]]=-mob[i];
                phi[i*prime[j]]=phi[i]*phi[prime[j]];
            }
        }
    }
}
ll s1(ll x)
{
    if(x==1)return 1;
    if(x<=10000000)return mob[x];
    if(m1[x])return m1[x];
    ll sum=0;
    ll l,r;
    for(l=2;l<=x;l=r+1)
    {
        r=x/(x/l);
        sum+=(r-l+1)*s1(x/l);
    }
    return m1[x]=1-sum;
}
ll s2(ll x)
{
    if(x==1)return 1;
    if(x<=10000000)return phi[x];
    if(m2[x])return m2[x];
    ll sum=0;
    ll l,r;
    for(l=2;l<=x;l=r+1)
    {
        r=x/(x/l);
        sum+=(r-l+1)*s2(x/l);
    }
    return m2[x]=(x+1)*x/2-sum;
}
int main()
{
    int t;
    getphi(10000000);
    cin>>t;
    ll x;
    for(int i=2;i<=10000000;i++)
    {
        mob[i]+=mob[i-1];
        phi[i]+=phi[i-1];
    }
    while(t--)
    {
        cin>>x;
        cout<<s2(x)<<" "<<s1(x)<<endl;
    }
    return 0;
}

本人在做的时候，先是筛了1e5，只过了6个数据，因为有种打表的感觉，莫名觉得很爽，所以就增加到1e6，过了9个，还不够，我就露出了阴间的笑容，上1e7，ac！