Programming Challenges 习题7.6.4

PC/UVa：110704/10139

Factovisors

判断m是否可以整除n!。

这题很简单，从n开始和m进行约分，直到m变成1。还有就是特殊情况特殊处理，m为0的测试用例在uDebug上显示可以整除任何值，但是UVa上就是WA，因为这个错了4次；m为1的也要特殊处理，原因在下面。

这道题反映出了除法是有多么慢！如果不判断质数直接进行辗转相除，uDebug上第二个都是较大数的例子就会算十几分钟。如果去掉代码中注释的那两行，并把循环之前相同的检查去掉，以及约分之后相同的检查也去掉，程序语义上没有变化，但是提交后就是超时。通过30多次坚持不懈的提交，终于知道了原因。在m没有变化的情况下，每次都判断m是不是质数就导致大量无意义的除法运算，所以就超时！！！

因为代码只能这么写，所以m为1也要特殊处理，因为1既不是质数也不是合数。

如果没有超时的限制，同时不单独判断m为1的情况，直接将循环下限变为1，此时当n为0时，会直接返回false，所以接着把循环下限改成0，但此时会判断1是不是质数，根据代码会返回是质数，所以也不对，如果再把1的判断结果返回false，那么辗转相除法会返回1，此时不会进行约分，所以依然返回false。

所以代码只能这么写。

#include <iostream>

using namespace std;

int GCD(int n, int m)
{
	int r = 0;
	if (n < m) n ^= m ^= n ^= m;
	while (m != 0){
		r = n % m;
		n = m;
		m = r;
	}
	return n;
}

bool Prime(int m)
{
	if (m == 2) return true;
	if (m % 2 == 0) return false;
	for (int i = 3; i * i <= m; i += 2)
	{
		if (m % i == 0) return false;
	}
	return true;
}

bool divided(int n, int m)
{
	if (m == 0) return false;
	if (m == 1) return true;
	if (m <= n) return true;
	else if (Prime(m)) return false;
	else{
		int gcd;
		for (int tmp = n; tmp > 1; tmp--)
		{
			//if (m <= tmp) return true;
			//else if (Prime(m)) return false;
			gcd = GCD(m, tmp);
			if (gcd > 1){
				m /= gcd;
				if (m < tmp) return true;
				else if (Prime(m)) return false;
			}
		}
		return false;
	}
}

int main()
{
	int n = 0, m = 0;
	while (cin >> n >> m){
		bool b = divided(n, m);
		if (b) cout << m << " divides " << n << '!' << endl;
		else cout << m << " does not divide " << n << '!' << endl;
	}
	return 0;
}
/*
6 9
6 27
20 10000
20 100000
1000 1009
*/