[Cogs]-693-Antiprime数

693. Antiprime数 – COGS.

//部分引用自byvoid.com，代码原创

如果一个数X是反质数，则它的约数的个数大于所有Y(X>Y)的约数的个数。也就是说反质数是最小的具有相同个约数的数。寻找不大于N的最

大反质数问题可以转化成，寻找不大于N的约数个数最多的最小正整数。

求一个数的约数个数可以用乘法原理，例如75=3^1×5^2，则75有(1+1)(2+1)=6个约数。对于这道题我们可以采取搜索的方法。按照质因数从小到大依次枚举指数，找出最多的约数个数情况下最小的正整数。

设A = ∏p[i]^a[i]    (1≤i≤k)    (字体显示有问题，∏表示求从…到…的积)

A的约数个数 = ∏(1 + a[i])      ,(1≤i≤k)

A的所有约数和 =  ∑x   (gcd(A,x)==x)  =  ∏∑p[j]^i  ,(0≤i≤a[j]，1≤j≤k)

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <vector>
#include <queue>
#include <algorithm>
  
#define inf 0xefffffff
  
using namespace std;
  
const int p[10] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29};
long long h, mxdivi, ans;
  
void dfs(int pi, long long divi, long long num)
//第pi个质因子，有divi个约数，数字为num
{
 if (divi > mxdivi)
  mxdivi = divi, ans = num;
 if (divi == mxdivi && num < ans)
  ans = num;
 for (int i = 0; num <= h && pi < 9; num *= p[pi], i++)
  dfs(pi + 1, divi * (i + 1), num);
}
  
int main()
{
 freopen("antip.in", "r", stdin);
 freopen("antip.out", "w", stdout);
  
 scanf("%lld", &h);
 dfs(0, 1, 1);
 printf("%lld", ans);
 return 0;
}