P3601 签到题

最新推荐文章于 2024-06-25 05:54:01 发布

wamach

最新推荐文章于 2024-06-25 05:54:01 发布

阅读量346

点赞数

分类专栏： luogu

本文链接：https://blog.csdn.net/wamach/article/details/55106896

版权

luogu 专栏收录该内容

26 篇文章 0 订阅

订阅专栏

按畜牧的话讲，写得这么慢根本没资格发博客了。

按畜牧的话讲，这道题是到讲过的傻逼题。

按序幕的话讲，如果是他出的题做不出就该滚粗，而这套题是初三小哥出的。

真的不想就这么萎掉QAQ。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<queue>
#include<cstdlib>
#include<cmath>
using namespace std;
#define rep(i,j,k) for(i=j;i<=k;++i)
#define per(i,j,k) for(i=j;i>=k;--i)
#define ll long long
#define db double
#define ldb long double
#define pli pair<ll,int>
#define mkp make_pair
#define X first
#define Y second
const ll N=1000005,MOD=666623333;
ll pri[N],cnt,vis[N],ans,phi[N],tmp[N];
int main(){
	ll l,r,i,j,y,rg;
	scanf("%lld%lld",&l,&r);rg=r-l;
	rep(i,2,N-1){
		if(!vis[i])pri[++cnt]=i;
		for(j=1;j<=cnt&&(y=i*pri[j])<N;++j){
			vis[y]=1;if(i%pri[j]==0)break;
		}
	}
	rep(i,0,rg)phi[i]=tmp[i]=i+l;
	rep(i,1,cnt)
		for(j=max((l-1)/pri[i]+1,2LL)*pri[i];j<=r;j+=pri[i]){
			phi[j-l]=phi[j-l]/pri[i]*(pri[i]-1);
			while(tmp[j-l]%pri[i]==0)tmp[j-l]/=pri[i];
		}
	rep(i,0,rg)
		if(tmp[i]>1)phi[i]=phi[i]/tmp[i]*(tmp[i]-1);
	rep(i,0,rg)ans=(ans+l+i-phi[i])%MOD;
	printf("%lld\n",ans);
	return 0;
}