2021-05-29 acwing筛质数

AcWing 1292. 哥德巴赫猜想

1.题意：

2.题解：

3.ac代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+10;
int primes[N];
bool st[N];
int tot;
void init(int n){
    for(int i=2;i<=n;i++){
        if(!st[i]) primes[tot++]=i;
        for(int j=0;primes[j]<=n/i;j++){
            st[primes[j]*i]=1;
            if(i%primes[j]==0) break;
        }
    }
}
int main(){
    int n;
    init(1000002);
    while(cin>>n&&n){
        int flag=1;
        for(int i=0;primes[i]<n;i++){
            if(!st[n-primes[i]]){
                flag=0;
                cout<<n<<" = "<<primes[i]<<" + "<<n-primes[i]<<endl;
                break;
            }
        }
        if(flag){
            cout<<"Goldbach's conjecture is wrong."<<endl;
        }
    }
}

AcWing 1293. 夏洛克和他的女朋友

1.题意：

2.题解：

当一个数是另一个数的质因子时，二者颜色不同，可以想到二分图，将所有数字分为质数和合数两个部分，质数都互质，总有质数是合数的质因子，质数一个颜色，合数一个颜色，所以当2-n+1存在合数时答案就是2，否则答案是1，并且质数被标记为1，合数是2，因为只有质数时显然不能用第二种颜色

3.ac代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+10;
int primes[N];
bool st[N];
int tot;
void sprimes(int n){
    for(int i=2;i<=n;i++){
        if(!st[i]) primes[tot++]=i;
        for(int j=0;primes[j]<=n/i;j++){
            st[primes[j]*i]=1;
            if(i%primes[j]==0) break;
        }
    }
}
int col[N];
int main(){
    sprimes(100005);
    int n;
    cin>>n;
    if(n>=3){
        cout<<2<<endl;
    }else{
        cout<<1<<endl;
    }
    for(int i=1;i<=n;i++){
        if(!st[i+1]){
            cout<<1<<" ";
        }else{
            cout<<2<<" ";
        }
    }
    cout<<endl;
}

AcWing 196. 质数距离

1.题意：

2.题解：

建议反复做
一开始段错误，一个是因为我把差值d付给了maxp和minp，另一个是因为
ll p=primes[i]; for(ll j=max(2*p,(l+p-1)/p*p);j<=u;j+=p){ st[j-l]=1; }
好像这块p和j不是longlong就不行

首先这题l,r很大，而且每次也无法只求该区间的的质数，线性筛这么大范围铁定不行
想到一个性质：任何一个合数n都有质因子,a,n/a <=根号n，l,r<=2³¹,2³¹大概是1e9，根号后是1e4~1e5,而5000>根号2³¹,所以筛出50000内的质数，遍历l，r，有这些质因子的数都是合数，筛掉

3.ac代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+10;
int primes[N];
bool st[N];
int tot;
void get_primes(int n){
    memset(st,0,sizeof st);
    for(int i=2;i<=n;i++){
        if(!st[i]) primes[tot++]=i;
        for(int j=0;primes[j]<=n/i;j++){
            st[primes[j]*i]=1;
            if(i%primes[j]==0) break;
        }
    }
}
int main(){
    int l,u;
    while(cin>>l>>u){
        tot=0;
        get_primes(50000);
        memset(st,0,sizeof st);
        int cnt=0;
        for(int i=0;i<tot;i++){
            ll p=primes[i];
            for(ll j=max(2*p,(l+p-1)/p*p);j<=u;j+=p){
                st[j-l]=1;
            }
        }
        tot=0;
        for(int i=0;i<=u-l;i++){
            if(!st[i]&&i+l>=2) primes[tot++]=i+l;
        }
        if(tot<2) {
            cout<<"There are no adjacent primes."<<endl;
            continue;
        }
        int maxp=0,minp=0;
        for(int i=0;i<tot-1;i++){
            int d=primes[i+1]-primes[i];
            if(d>primes[maxp+1]-primes[maxp]){
                maxp=i;
            }
            if(d<primes[minp+1]-primes[minp]){
                minp=i;
            }
        }
        printf("%d,%d are closest, %d,%d are most distant.\n",
                primes[minp],primes[minp+1],primes[maxp],primes[maxp+1]);
        
    }
}