xmu 1307一把心酸一把泪啊

最新推荐文章于 2024-07-02 08:45:18 发布

vsooda

最新推荐文章于 2024-07-02 08:45:18 发布

阅读量769

点赞数

CC 4.0 BY-SA版权

分类专栏： xmu

本文链接：https://blog.youkuaiyun.com/vsooda/article/details/7316001

xmu 专栏收录该内容

8 篇文章

订阅专栏

本文详细介绍了使用C++实现的筛法求质数算法，并结合序列求和的应用实例。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<stdio.h>
#include<string.h>
#include<math.h>
int a[1000010];
#define MAXN 1000010
unsigned int plist[60000],pcount;
unsigned int isprime[(MAXN>>5)+1];
#define setbitzero(a) (isprime[(a)>>5]&=(~(1<<((a)&31))))
#define setbitone(a) (isprime[(a)>>5]|=(1<<((a)&31)))
#define ISPRIME(a) (isprime[(a)>>5]&(1<<((a)&31)))
void initprime(){
    int i,j,m;
    int t=(MAXN>>5)+1;
    for(i=0;i<t;++i)isprime[i]=2863311530;
    plist[0]=2;setbitone(2);setbitzero(1);
    m=(int)sqrt(MAXN);
    for(pcount=1,i=3;i<=m;i+=2)
        if(ISPRIME(i))
            for(plist[pcount++]=i,j=i<<1;j<=MAXN;j+=i)
                setbitzero(j);
    if(!(i&1))++i;
    for(;i<=MAXN;i+=2)
		if(ISPRIME(i))
		{
			plist[pcount++]=i;
			a[i]=i;
		}
}


int main()
{
	int t,i,j,temp;
	scanf("%d",&t);
	memset(a,0,sizeof(a));
	initprime();
	for(i=2;i<=1000010;i++)
	{
		if(a[i]==0)
		{
			temp=i;
			for(j=2;j<=temp;j++)
			{
				if(temp%j==0)
				{
					a[temp]=a[temp/j]+j;
					temp=temp/j;
					break;
				}
			}
		}	
	}
	while(t--)
	{
		int l,r;
		long long sum=0;
		scanf("%d%d",&l,&r);
		for(i=l;i<=r;i++)
			sum+=a[i];
		printf("%lld\n",sum);
	}
	return 0;

}