poj 2689 Prime Distance

    xiaoxiao2022-06-29  29

    题意 给两个数m,n求m~n间的任意相邻的两个素数的最短距离和最大距离,若不存在两个素数输出

    There are no adjacent primes.

    这个题坑的我WA了十多次,千万注意数据范围:m-n<1000000; 然而m和n均在int范围内所以必须采用映射比如说数m~n之间的任意一个数i都可以用i-m存到数组里,若求的是2000000~3000000,则int数组b;b[0]记录的就是2000000,b[1]记录的就是2000001;

    解题的大致思路:因为两数之差不超过一千万,所以在m~n之间的素因子必定在根号214783648<50000之下,所以一开始只要把小于50000的素数筛选出来即可,再通过这些素数标记m~n之间的合数,剩下的就是这个区间里的素数了。

    代码:

    #include<cmath> #include<cstdio> #include<cstdlib> #include<cstring> #include<iostream> #include<algorithm> using namespace std; long long zd,zx,a[50100],b[1000010],c[50100]; int main(){ long long i,j,k,m,n,end1,end2,sum=0,end3,end4,u; long long p=floor(sqrt(50010)); for(i=2;i<=p;i++){ if(a[i])continue; for(j=i*i;j<=50010;j+=i){ a[j]=1; } } for(i=2;i<=50010;i++)if(!a[i])c[++sum]=i; while(scanf("%lld%lld",&m,&n)!=EOF){ if(m==1)m=2; long long zd=0; long long zx=999999999; long long cnt=-1; memset(b,0,sizeof(b)); for(i=1;i<=sum;i++) { long long l,r; if(m%c[i]==0)l=m/c[i]; else l=m/c[i]+1; r=n/c[i]; for(j=l;j<=r;j++) if(j>1)b[-1*m+c[i]*j]=1; } for(i=m;i<=n;i++) if(b[i-m]==0)b[++cnt]=i; if(cnt<1)printf("There are no adjacent primes.\n"); else { for(i=0;i<cnt;i++){ if(zd<b[i+1]-b[i]){ zd=b[i+1]-b[i]; end1=b[i]; end2=b[i+1]; } if(zx>b[i+1]-b[i]){ zx=b[i+1]-b[i]; end3=b[i]; end4=b[i+1]; } } printf("%lld,%lld are closest, %lld,%lld are most distant.\n",end3,end4,end1,end2); } } return 0; }
    转载请注明原文地址: https://ju.6miu.com/read-1124863.html

    最新回复(0)