大区间求区间欧拉函数（洛谷P3601）

题目传送门
欧拉函数：对于n，φ(n)表示1到n内与n互质的数。

题意，定义一个函数f，f(x)表示x以内与x不互质的数的个数。给定一个区间，求出这个区间内f(i) (l<=i<=r)的总和对666623333取模的结果。

数据范围
1<=l<=r<=1e12，r-l<=1e6

求不互质的数的个数，我们只需要求这个数的欧拉函数，然后f(x)=x-φ（x），问题转化为求φ（x）。

因为r-l<=1e6，那么我们可以尝试求每一个数的欧拉函数，但是这个方法是根号级别的，也就是说时间是(r-l)*r^(1/2)，肯定超时了。

因为欧拉函数是个积性函数，对于两个质数n,m，φ（n*m）=φ（n）*φ（m）；所以我们是否对每个合数都找到两个或若干个不同的质数然后求呢？好像不太好做。

还有一种方法，我们先初始化，所有φ（i）=i，然后把1e6以内的素数都筛出来，然后对于区间内的所有数，我们把这个数以内每个素数的倍数都删去，剩下的就是和i互质的数啦。

Code

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+5;
const int inf=0x7fffffff;
const int mod=1e9+7;
const int eps=1e-6;
typedef long long ll;
typedef unsigned long long ull;
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define pii pair<int,int>
#define int long long
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define endl '\n'
bool ok[N];int p[N],cnt;
void f()
{
    for(int i=2;i<=N;i++)
    {
        if(!ok[i])
            p[++cnt]=i;
        for(int j=1;j<=cnt;j++)
        {
            if(i*p[j]>=N)
                break;
            ok[i*p[j]]=1;
            if(i%p[j]==0)
                break;
        }
    }
}
ll a[N],b[N];
void solve(ll l,ll r)
{
    for(int i=1;i<=cnt&&p[i]*p[i]<=r;i++)
    {
        ll pp=p[i];
        for(int x=(pp-l%pp)%pp;x<=r-l;x+=pp)
        {
            a[x]/=pp,a[x]*=pp-1;
            while(b[x]%pp==0)
                b[x]/=pp;
        }
    }
}
signed main()
{
    IOS;
    int l,r;
    cin>>l>>r;
    f();
    for(int i=0;i<=r-l;i++)
        a[i]=i+l,b[i]=i+l;
    solve(l,r);
    ll res=0;
    for(int i=0;i<=r-l;i++)
    {
        if(b[i]!=1)
            a[i]/=b[i],a[i]*=b[i]-1;
        res=(res+i+l-a[i])%666623333;
    }
    cout<<res<<endl;
}