参见百度百科:http://baike.baidu.com/link?url=6CBrDsXMVeF3zmNwL_ugd-OEeVuVyxRFJVRnt59BJxSPiLVAzc7dEqvlxhl4topCBRE5FjeRoO-JNO3HRLU1Da
其中反素数有两个重要的性质:
<1>性质一:反素数的素因子分解中素因子必须保证是从小到大连续的. PS:因为如果不连续我们总能构造一个比其小的数,而且该数的约数的个数和其相同
<2>性质二:反素数的素因子分解中素因子的指数是单调不增的. PS: a = p1^a1*p2^a2*p3^a3...pk^ak 假设ak-1<ak,我们完全构造得到另一个数b = a / (pk^ak*pk-1^ak-1)*(pk-1^ak*pk^ak-1)这样一来有约数个数的计算公式可得a和b的约数个数是相同的,但是b<a,因此性质二是满足的.
可以依据上述的两条性质对搜索的过程进行相应的剪枝:
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
const int MAXM = 1010; //反素数的个数
const int MAXN = 500010; //求解该范围内的反素数
ll bestNum; //约数最多的数
ll bestSum; //约数最多的约数的个数
ll rprime[MAXM][2]; //2*3*5*7*11*13*17>MAXN
ll prime[14] = {2, 3, 5, 7, 11, 13,