PAT 2019-3 7-1 Sexy Primes

Description:

Sexy primes are pairs of primes of the form (p, p+6), so-named since "sex" is the Latin word for "six". (Quoted from http://mathworld.wolfram.com/SexyPrimes.html)

Now given an integer, you are supposed to tell if it is a sexy prime.

Input Specification:

Each input file contains one test case. Each case gives a positive integer N (≤).

Output Specification:

For each case, print in a line Yes if N is a sexy prime, then print in the next line the other sexy prime paired with N (if the answer is not unique, output the smaller number). Or if N is not a sexy prime, print No instead, then print in the next line the smallest sexy prime which is larger than N.

Sample Input 1:

47

Sample Output 1:

Yes
41

Sample Input 2:

21

Sample Output 2:

No
23

Keys:

• 素數

Attention:

• 這幾次很愛考素數，奶一口歐拉算法，9月份考試能不能中

Code:

/*
SP數，P和P+6都是素數
給定一個正整數，判斷N是否為SP數，打印與他配對的另一個素數（如果不唯一，輸出較小的一個
如果N不是SP數，打印大於N的最小SP數

1，判斷i-6是否大於n
2，n-6可能小於0,
*/
#include<cstdio>

bool IsPrime(long long n)
{
    if(n <= 1)
        return false;
    for(long long i=2; i*i<=n; i++)
        if(n%i == 0)
            return false;
    return true;
}

int main()
{
#ifdef ONLINE_JUDGE
#else
    freopen("Test.txt", "r", stdin);
#endif // ONLINE_JUDEG

    int n;
    scanf("%d",&n);
    if(IsPrime(n-6) && IsPrime(n))
        printf("Yes\n%d", n-6);
    else if(IsPrime(n) && IsPrime(n+6))
        printf("Yes\n%d", n+6);
    else
    {
        for(int i=n+1; i<1e9; i++)
            if(IsPrime(i-6) && IsPrime(i))
            {
                printf("No\n%d", i-6>n?i-6:i);
                break;
            }
    }

    return 0;
}