hdu2964-Prime Bases

http://acm.hdu.edu.cn/showproblem.php?pid=2964

题意，给你一个整数n,现在要你分解成 n = k1 * ( 2 * 3 * ....*x1 ) + k2 * ( 2 * 3 * .... *x2 ) + ........;其中后面均为素数，且是由最小的2递增相乘；

分析：首先打印素数表；然后使用数组a【】来储存到从第一个素数2到第几个素数乘积。找到第i个小于n，第i +1个大于n的数值 i ；

然后分解n.使用数组b【】来储存a【i】的个数；同理执行多次，知道将n分解完毕；

如下的这种输出格式比较好，大家可以参考下；

参照做的；

 

#include<iostream>

#include<cstdio>

#include<cstring>

#include<cmath>

#include<algorithm>

#include<bitset>

#include<iomanip>



using namespace std;

const int prime[]={2,3,5,7,11,13,17,19,23,29,31,37,39,41,43,47,51,53,57,59};  

__int64 a[21];  

int b[21],rem; 



void init( )

{

	a[ 0 ] = 1 ;

	for( int i =1 ; i < 21 ; ++i )

		a[ i ] = a[ i - 1 ] * prime[ i - 1 ] ;

}



void deal( int n )

{

	int i ;

	for( int i = 0 ; i < 21 ; ++i )

	{

		if( a[ i ] <= n && a[ i + 1 ] > n )

		{

			rem = i ;

			break ;

		}

	}

	memset( b , 0 ,sizeof( b ) ) ;

	for( int i = rem ; i >=0 ; --i )

	{

		b[ i ] = n / a[ i ] ;

		n %= a[ i ] ;

	}

}



void output( int n )

{

	cout << n << " = " ;

	for( int i = 0 ; i <= rem ; ++i )

	{

		

		if( b[ i ] )

		{

			cout << b[ i ] ;

			for( int j = 0 ; j < i ; ++j )

				cout << "*" << prime[ j ] ;

			if( i != rem )

				cout << " + " ;

		}

	}

	cout << endl ;

}



int main( )

{

	init( ) ;

	int n ;

	while( cin >> n , n )

	{

		deal( n ) ;

		output( n );

	}

	return 0 ;

}