UVa 10139 Factovisors

这道题同样的代码用 JAVA 会超时，用 C++ 才通过。

#include <iostream>
#include <vector>
#include <cmath>
using namespace std;

// 素数筛的数组长度
long LEN = 1 << 16;

// 用于保存筛出的素数
vector<long> primes;

// 利用筛法求  2^16 以内的素数，保存到 primes 全局变量
void init() {

	bool filter[LEN];
	for ( int i=0; i<LEN; i++ ) filter[i] = true;

	long max = (long)sqrt(LEN)+1;
	for ( int i=2; i<=max; i++ ) {
		if ( filter[i] ) {
			for ( int j=2*i; j<LEN; j+=i )
				filter[j] = false;
		}
	}

	for ( int i=2; i<LEN; i++ )
		if ( filter[i] ) {
			primes.push_back(i);
		}

}

// 计算素数 p ^ k 是否整除 n!
bool divByP2Kth(long p, long k, long n) {

	if ( k==0 ) return true;

	long cnt = 0;
	for ( long i=p; i<=n; i*=p ) {
		cnt += n/i;
		if ( cnt >= k ) return true;
	}

	return false;
}

// 计算 m 是否整除 n!
bool divide(long m, long n) {

	if ( m==0 )
		return false;
	else if ( m== 1 )
		return true;

	for ( unsigned int i=0; i<primes.size(); i++ ) {

		// 求 m 的素因数分解中素数 p 的指数 k
		long p = primes.at(i), k = 0;
		while ( m%p == 0 ) {
			k++; m /= p;
		}

		// p^k 是否整除 n! ?
		if ( ! divByP2Kth(p,k,n) ) return false;

		// 分解完毕
		if ( m==1 ) break;

	}

	// 如果剩余一个大于 2^16  的素因子
	if ( m>1 && ! divByP2Kth(m,1,n) ) return false;

	return true;

}

int main() {

	init();

	long n, M;
	while (cin >> n >> M) {
		bool ok = divide(M,n);
		if ( ok ) {
			cout << M << " divides " << n << "!" << endl;
		} else {
			cout << M << " does not divide " << n << "!" << endl;
		}
	}
	return 0;
}