USACO-Section1.5 Superprime Rib[其他][暴力枚举]

2017-5-30

题目

寻找超级素数，比如7331，本身是素数，733也是，73和7都是，那么它就符合条件。比如四位数的超级素数有：
2333
2339
2393
2399
2939
3119
3137
3733
3739
3793
3797
5939
7193
7331
7333
7393
输入n表示n位数，输出所有可能。

题解

有了上一题的经验，我根本就不想搞什么素数表，也不想什么枚举判断。于是想出了一种逐层搜寻的方法。
一位数的素数有2，3，5，7，这就使得超级素数的首位必须是这四个，后几位的数字都有机会作为需要判断的数字的末尾，故只能是1，3，7，9。所以将第一位的素数乘以10后加上四者之一，仍为素数的放入数组内保存，两位数找完后再将所有结果乘以10后加上四者之一，仍为素数的放入三位数对应的数组保存，以此类推得到所求（最多8位）
尝试情况共5*4^7种。

代码

/*
ID: xcwhkh1
LANG: C 
TASK: sprime
*/
#include <stdio.h>
int pri(int x)
{
    for(int i=3;i*i<=x;i+=2)
        if(x%i==0)
            return 0;
    return 1;
}
int a[8][100000]={0};//max<5*4^7
int main(){
    FILE *fin=fopen("sprime.in","r");
    FILE *fout=fopen("sprime.out","w");
    int n,num;
    int prime[4]={1,3,7,9};//尾部添的数字 
    fscanf(fin,"%d",&n);
    a[0][0]=2;
    a[0][1]=3;
    a[0][2]=5;
    a[0][3]=7;
    for(int i=1;i<=n;i++)//计算到n位数为止 
    {
        num=0;
        for(int pre=0;a[i-1][pre]!=0;pre++)
            for(int j=0;j<4;j++)
                if(pri(a[i-1][pre]*10+prime[j]))//添数字
                    a[i][num++]=a[i-1][pre]*10+prime[j];
    }
    for(int i=0;a[n-1][i]!=0;i++)
        fprintf(fout,"%d\n",a[n-1][i]);
    return 0;
}