7-178 吸血鬼素数

1994年，美国数学家皮科夫在一篇文章中首次提出了吸血鬼数。如果一个2n（n是自然数）位的自然数等于自己各个数字任意组成的两个n位数的乘积，那么这个自然数就叫吸血鬼数。其中这两个n位数都被称为獠牙或尖牙。比如1260=21×60，所以1260是一个吸血鬼数。注意两个尖牙的个位数不能同时为0，例如虽然126000=210×600，但126000却不是吸血鬼数。在此基础上，2002年，数学家里维拉定义了吸血鬼素数。如果吸血鬼数的两个尖牙都是素数，则这两个尖牙都叫吸血鬼素数。例如，117067=167×701，而167和701都是素数，所以167和701都是3位吸血鬼素数。

输入一个正整数n（小于1000），请判断它是否吸血鬼素数。如果是，打印YES，否则打印NO。

输入格式:

一个小于1000的正整数n

输出格式:

YES或者NO

输入样例:

701

输出样例:

YES

代码长度限制

16 KB

C (gcc)

时间限制

400 ms

内存限制

64 MB

其他编译器

时间限制

400 ms

内存限制

64 MB

代码

首先输入一个正整数n，然后检查n是否为素数。如果n是素数，则遍历所有2n位数中每一对因子x和y（即x * y = n * 10^num_digits(n) + m），其中m是与x和y都为n位数相对应的另一个n位数，判断它们是否符合“吸血鬼数字”的定义：

• x和y包含了n中所有数字（但可能重复出现），且由它们组成的新数字的位数等于2n。
• 对于新数字中的每一种不同的数字d，x、y中包含此数字的次数之和等于新数字中这个数字的出现次数。

如果找到符合条件的一对x和y（或者说找到了吸血鬼数字），则输出“YES”并停止程序；否则输出“NO”。

思路

#include <stdio.h>

//判断某数x是否为素数，如果是返回1，否则返回0
int ss(int x)
{
    int z = 1;
    int i;

    if (x < 2) {
        z = 0; //0和1都不是素数
    }
    else if (x > 2) {
        if (x % 2 == 0) { //偶数一定不是素数
            z = 0;
        }
        for (i = 3; i <= sqrt(x) && z == 1; i += 2) { //从3开始枚举奇数，逐个判断
            if (x % i == 0) {
                z = 0;
                break; //如果发现因子，则直接终止循环
            }
        }
    }

    return z;
}

//计算给定整数x的位数
int weishu(int x)
{
    int z = 0;

    if (x != 0) {
        while (x != 0) { //从低位开始逐位统计
            z++;
            x /= 10;
        }
    }
    else {
        z = 1; //0的位数为1，默认数值为1
    }

    return z;
}

int main(void)
{
    int n, m;
    int nm[10], is[10]; //用于计数数字出现次数的数组
    int h; //n的位数
    int max_h; //枚举范围上限
    int i, j; //循环变量
    int z = 0; //表示找到了吸血鬼数字
    scanf("%d", &n);
    //如果n是素数，则继续检查它是否是吸血鬼数字组成因子
    if (ss(n) == 1) { 
        h = weishu(n); //计算n的位数
        max_h = (int)pow(10, 2 * h); //计算枚举范围上限
        for (i = (int)pow(10, 2 * h - 1); i < max_h; i++) { //i取值范围满足x和y均为n位数，并可以减小至2*h-1位
            if (i % n == 0) {
                m = i / n; //求出y的值
                if (weishu(m) == h && ss(m) == 1) { //如果y的位数也是n，且y是素数，就检查i和m是否符合吸血鬼数字的定义
                    //printf("%d %d %d\n",i,n,m);
                    for (j = 0; j < 10; j++) { //分别统计i、n和m中每个数字出现的次数
                        nm[j] = 0; //初始化数组
                        is[j] = 0;
                    }
                    while (m > 0) { //处理m，统计其中各个数字的出现次数
                        nm[m % 10]++;
                        m /= 10;
                    }
                    m = n;
                    while (m > 0) { //处理n，统计其中各个数字的出现次数
                        nm[m % 10]++;
                        m /= 10;
                    }
                    m = i;
                    while (m > 0) { //处理i，统计其中各个数字的出现次数
                        is[m % 10]++;
                        m /= 10;
                    }
                    z = 1; //假设i和m构成了吸血鬼数字，然后检查是否各个数字都出现了一次
                    for (j = 0; j < 10; j++) {
                        if (nm[j] != is[j]) {
                            z = 0; //如果有一个数字出现次数不同，则认为i和m不是吸血鬼数字
                            break; //终止循环
                        }
                    }
                    if (z == 1) { //找到了吸血鬼数字
                        printf("YES");
                        break; //终止循环
                    }
                }
            }
        }
        if (z == 0) { //未找到吸血鬼数字
            printf("NO");
        }
    }
    else { //如果n不是素数，则直接输出“NO”
        printf("NO"); 
    }

    return 0;
}