【洛谷】P3601 签到题

最新推荐文章于 2024-07-09 01:32:36 发布

EagleEyeKestrel

最新推荐文章于 2024-07-09 01:32:36 发布

阅读量1k

点赞数 2

CC 4.0 BY-SA版权

分类专栏：洛谷数学

本文链接：https://blog.csdn.net/u014156276/article/details/108395506

洛谷同时被 2 个专栏收录

17 篇文章

订阅专栏

数学

13 篇文章

订阅专栏

洛谷P3601 签到题

题目描述
我们定义一个函数：qiandao(x)为小于等于x的数中与x不互质的数的个数。

这题作为签到题，给出l和r，要求求 $∑i=lrqiandao(i)mod666623333\sum_{i=l}^r qiandao(i) mod 666623333$ .

输入格式
一行两个整数，l、r。

输出格式
一行一个整数表示答案。

输入输出样例
输入 #1

233 2333

输出 #1

说明/提示

对于100%的数据， $1≤l≤r≤1012,r−l≤1061\leq l\leq r\leq 10^{12},r-l\leq 10^6$ .

这题一看就是要求l到r范围内所有x，与x不互质且小于等于x的数个数。一看数据范围，l和r是 $10^{12}$ 范围，暴力不太行，但r-l的范围却是可观的。考虑用数组phi[i]存储与l+i互质且小于等于i的元素个数。最关键的是欧拉函数的公式别忘记了：
$φ(x)=x∏i=1n(1−1pi)\varphi(x)=x\prod_{i=1}^n(1-\frac{1}{p_i})$
其中 $,pnp_1,p_2,\cdots,p_n$ 为x所有质因数。那么接下来就需要求l到r范围内数的所有质因数，一个个遍历不太好，比较好的办法是先预处理出来 $10^6$ 范围内的素数（基础的线性筛），然后遍历这些素数，将i的phi初始化为i，如果i含有这个素因子就减小。代码如下：

now标记素数数组primes当前的位置，isprime用于预处理素数。
phi[i]表示i+l的欧拉函数，remain[i]表示i+l除掉之前的素因子之后还剩下什么。由于处理 $10^6$ 以上的素因子开销太大，而这么大的素因子的个数至多一个，所以如果最后remain[i]>1，就说明这是一个大于 $10^6$ 的素因子。

#include<bits/stdc++.h>
using namespace std;
long l,r,now=1,res=0,mod=666623333;
long isprime[1000010],primes[1000010];
long phi[1000010],remain[1000010];
void build(){
    isprime[1]=0;
    for(int i=2;i<=1e6;i++){
        if(isprime[i]){
            primes[now++]=i;
        }
        for(int j=1;j<now&&primes[j]*i<=1e6;j++){
            isprime[primes[j]*i]=0;
            if(i%primes[j]==0) break;
        }
    }
}
int main(){
    cin>>l>>r;
    memset(isprime,1,sizeof(isprime));
    build();
    for(long i=l;i<=r;i++){
        phi[i-l]=i;
        remain[i-l]=i;
    }
    for(long i=1;i<now&&primes[i]*primes[i]<=r;i++){
        long p=primes[i];
        for(long j=((l-1)/p+1)*p;j<=r;j+=p){
            phi[j-l]=phi[j-l]/p*(p-1);
            while(remain[j-l]%p==0) remain[j-l]/=p;
        }
    }
    for(long i=l;i<=r;i++){
        if(remain[i-l]>1){
            phi[i-l]=phi[i-l]/remain[i-l]*(remain[i-l]-1);
        }
        res=(res+i-phi[i-l])%mod;
    }
    cout<<res;
    return 0;
}