USACO Training Section1.5 Superprime Rib

最新推荐文章于 2018-09-24 22:51:47 发布

原创最新推荐文章于 2018-09-24 22:51:47 发布 · 168 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

收录于

USACO Training

当前文章被收录于：

USACO Training

56 篇文章 0 人学习

订阅专栏查看详情

当前文章被以下社区和专栏收录：

本文介绍了一个简单的程序，用于计算超级素数。通过优化算法减少运算量，并使用MillerRabin算法进行素数测试。代码实现了一次编译和一次AC。

简单不过的程序，超级素数的计算。注意最高位可选的数与其他位有所不同，其他位的数字可排除偶数，减少运算量。一次编译，一次AC。

顺便查了下作素数测试的Miller Rabin算法，随机算法学过很久很久了，这算是经典的一个。

原题链接：
[url]http://ace.delos.com/usacoprob2?a=bFftpewR5eD&S=sprime[/url]


/*
ID: blackco3
TASK: sprime
LANG: C++
*/
#include <fstream>
#include <math.h>
using namespace std;

int prime[]={2,3,5,7}, n_prime=4 ;
int org[]={1,3,5,7,9}, n_org=5 ;
int n_digit ;
ofstream fout("sprime.out");

inline int is_prime(int num) {
	for( int i=2; i*i<=num; i++)
		if( !(num%i) )
			return 0;
	return 1;
}

void get_sprime(int rank, int cur)
{
	if( rank==n_digit ){
		fout << cur << endl ;
		return ;
	}
	for( int i=0; i<n_org; i++ ){
		int next = cur*10+org[i] ;
		if( is_prime(next) )
			get_sprime(rank+1,next);
	}
}

int main()
{	
	ifstream fin("sprime.in");
	fin >> n_digit ;
	int cur_sum ;
	for( int i=0; i<n_prime; i++){
		cur_sum=prime[i] ;
		get_sprime(1,cur_sum);	
	}

	return 0;
}