社区讨论
丑陋代码,在线求调,qwq
P1835素数密度参与者 3已保存回复 2
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 2 条
- 当前快照
- 1 份
- 快照标识符
- @lo278zrq
- 此快照首次捕获于
- 2023/10/23 09:08 2 年前
- 此快照最后确认于
- 2023/11/03 09:23 2 年前
CPP
#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<string>
#include<cstring>
#include<algorithm>
using namespace std;
const int MAX=50010;
int f[MAX];
bool b[MAX];
bool a[1000001];
void zyz()
{
int i,j,k,t=0;
bool p;
memset(f,0,sizeof(f));
memset(b,false,sizeof(b));
t=1; f[1]=2; i=3; b[2]=true;
while(i<=50000)
{
k=int(sqrt(i));
p=true;
for(j=2; j<=k; j++)
{
if(i%j==0)
{
p=false;
break;
}
}
if(p)
{
t++;
f[t]=i;
b[i]=true;
}
i+=2;
}
}
int main()
{
int l,r,ans,i,j,t,k;
cin>>l>>r;
zyz();
ans=0;
if(r<=50000)
{
for(int i=1; i<=r; i++)
{
if(b[i]) ans++;
}
}
else
{
if(l<=50000)
{
for(int i=1; i<=50000; i++)
{
if(b[i]) ans++;
}
l=50001;
}
k=int(sqrt(r));
i=1;
while(k>f[i]) i++;
k=i;
memset(a,true,sizeof(a));
for(int i=1; i<=k; i++)
{
for(int j=l/f[i]; j<=r/f[i]; j++)
{
if((f[i]*j<=r)&&(f[i]*j>=l))
{
a[f[i]*j-l]=false;
}
}
}
for(int i=0; i<=r-l; i++)
{
if(a[i]) ans++;
}
}
cout<<ans<<endl;
return 0;
}
回复
共 2 条回复,欢迎继续交流。
正在加载回复...