Usaco Training Section 1.6 Superprime Rib

让你求n位数中的超级素数(sp)。(如7331,7331、733、73、7都是素数，7331就是sp)

这题我写了三个版本。第一个版本是线性筛，求出所有素数，判断前缀是否为素数。但打完自测时发现会超时，换了第二个版本。第二个版本就是暴力。暴力枚举所有n位数，再判断。打了一半突然想到第三种做法。因为sp有一个性质：一个sp的前缀肯定也是sp。于是我们可以用构造的方法。一位的sp只有2,3,5,7。然后可以构造两位的sp，就在一位的sp之后加上奇数(1,3,5,7,9),再判断是否为素数。然后三位，四位，……n位。由于我前两个版本没有提交，这题一遍过。USACO又发了个有趣的congratulation。

AC代码：

#include<bits/stdc++.h>
using namespace std;

inline int read(){
	int x=0;char c=getchar();
	while(c<'0'||c>'9') c=getchar();
	while(c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x;
}

int a[9][17],b[9];//17:用程序算出n位数中最多只有16个sp 
int main()
{
	freopen("sprime.in","r",stdin);
	freopen("sprime.out","w",stdout);
	int n=read();
	if(n==1){
		printf("2\n3\n5\n7\n");
		return 0;
	}
	int l=1,r=0;
	for(int i=1;i<n;++i) l=l*10,r=r*10+9;
	a[1][1]=2,a[1][2]=3,a[1][3]=5,a[1][4]=7;
	b[1]=4;
	for(int i=2;i<=n;++i){
		for(int j=1;j<=b[i-1];++j){
			for(int k=1;k<=9;k+=2){
				int x=a[i-1][j]*10+k;
				bool f=1;
				for(int p=2;p*p<=x;++p)
					if(x%p==0){f=0;break;}
				if(f) a[i][++b[i]]=x;
			}
		}
	}
	for(int i=1;i<=b[n];++i) printf("%d\n",a[n][i]);
}