注意到超级n位(n>1)素数的个位数必然是3或7,则递推的初始条件:a[1] = 3,a[2] = 7, g = 2.
#include<iostream>
#include<stdio.h>
#include<math.h>
using namespace std;
int isPrime(double x); //判断是否为素数
int main(){
cout<<"请确认n值:";
int g; //两位超级素数的个数
int n; //n位超级素数
int s; //n位超级素数总和
int t;
double d,f,a[20000],b[20000],e[20];
cin>>n;
g = 2;
s = 0;
a[1] = 3;a[2] = 7;e[1] = 1;
for(int k=2;k<=n;k++){ //从个位数开始
e[k] = e[k-1]*10;
t = 0;
for(int j=1;j<=9;j++){ //如 个位的大一位 百位,添加一个百位数
for(int i=1;i<=g;i++){
f = j*e[k] + a[i];
if(isPrime(f) == 1){
t++;b[t]=f;
if(k==n) s++;d = f; //统计并记录最大超级素数
}
}
}
g = t;
for(int i=1;i<=g;i++) a[i] = b[i]; //g个k位b[i]赋值给a[i]
}
cout<<"共"<<s<<"个"<<n<<"超级素数"<<endl;
cout<<"最大数为"<<d<<endl;
return 0;
}
//判断是否为素数
int isPrime(double x){
for(int i=2;i<=sqrt(x);i++){
if(fmod(x,i) == 0){
return 0;
}
}
return 1;
}