【题意分析】

  本题中,x被称为反质数,当且仅当没有任意一个严格小于x的正整数的约数个数大于x的约数个数。求不超过N的最大反质数。

【解题思路】

  数据范围中最大的N=2*109

  首先可以证明,不超过N的反质数不会拥有9个以上的不同质因数。因为2*3*5*7*11*13*17*19*23*29=6469693230>6*109>N。

  设某数n=∏piki(pi<pi+1),则其约数个数g(n)=∏(ki+1)。(因为每个质数对约数个数的贡献是相互独立的,质数pi的可能选择方案数为(ki+1),所以可以用乘法原理乘起来)。

  显然,对于相同的顺序序列k,选择越小的pi越优,于是最优选择方案就是选择前9个质因数。

  于是暴力枚举的状态数为∏[logpN],则其至多为[log2N]*[log3N]*[log5N]*[log7N]*[log11N]*[log13N]*[log17N]*[log19N]*[log23N]=3779758080。

  显然直接暴力是无法过的,于是需要一些鲁(吉)棒(丽)或玄(松)学(爷)优化。

  所谓鲁棒优化,就是打表。。先把所有的反质数用上面这个爆搜打出来存在表里,然后二分查找即可。

  打表做法的可行性得益于反质数个数的增长极其缓慢,105的数据范围中只有30个反质数,从下图不难看出。

 

  玄学优化呢,有两种方法:

•方法一:考虑对ki的枚举进行优化。一种朴素的想法是同一个素因数的个数过多一定不利于让答案最优,而且越大的质因数个数应当越少,于是可以面向数据调参,限制ki枚举的上限。

•方法二:部分记忆化,f[i][j]表示j的乘积分配给第i个开始的质数最大能达到的约数个数,然后可以对超出记忆化范围的搜索做下界减枝。

  复杂度O(松)。

【参考代码】

  然而当时这题我只用了玄学优化方法一的弱化版,不知为什么就0ms过了?!

  可能有更加紧确的复杂度分析或者bz的数据有毒。。无论是哪一点请读者指出,不胜感激。

 1 #include<cstdio>
 2 #define REP(I,start,end) for(int I=start;I<=end;I++)
 3 const int prime[16]={1,2,3,5,7,11,13,17,19,23,29,31,37,41,43,47};
 4 long long maxsum, bestnum, n;
 5 void getantiprime(long long num, long long k,long long sum,int limit)
 6 {
 7     int i;
 8     long long temp;
 9     if(sum>maxsum)
10     {
11         maxsum=sum;
12         bestnum=num;
13     }
14     if(sum==maxsum&&bestnum>num)
15     bestnum=num;
16     if(k>15)
17         return;
18     temp=num;
19     REP(i,1,limit)
20     {
21         if(temp*prime[k]>n)
22             break;
23         temp*=prime[k];
24         getantiprime(temp,k+1,sum*(i+1),i);
25     }
26 }
27 int main()
28 {
29     scanf("%lld",&n);
30     getantiprime(1,1,1,50);
31     printf("%lld\n",bestnum);
32     return 0;
33 }
View Code