P1463 [POI2001][HAOI2007] 反素数——解题报告

一、题目链接

P1463 [POI2001][HAOI2007] 反素数

二、题目大意

设函数$f(x)$表示$x$的约数个数。对于一个正整数$x$，若满足任意$0<i<x$，都有$f(x)$大于$f(i)$。那么称$x$为反质数。现给定数$N$，找到$1\sim N$中最大的反质数。

数据范围：$1\le N\le 2\times 10^9$。

三、题目分析

首先，$1\sim N$的最大反质数一定是$1\sim N$中约数个数最多的数中最小的一个。下面给出证明：

• 设这个数为$x$，因为$f(x)$最大是函数$g$的最大值，且$x$是最小的那一个，所以小于$x$的数的$f(i)<f(x)$成立，所以$x$为反质数。
• 因为$f(x)$最大是函数g的最大值，且$x$是最小的那一个，所以大于$x$的数的$f(i)\le f(x)$成立，所以大于$x$的所有数一定不是反质数。

$x$为反质数的必要条件是：$x$分解质因数后可以写成$2^{c_1}\times 3^{c_2}\times 5^{c_3}\times 7^{c_4}\times 11^{c_5}\times 13^{c_6}\times 17^{c_7}\times 19^{c_8}\times 23^{c_9}\times 29^{c_{10}}$。并且$c_1\ge c_2\ge c_3\ge c_4\ge c_5\ge c_6\ge c_7\ge c_8\ge c_9\ge c_{10}$。下面给出证明：

• 反证法：若存在质因数分解中有一项$p^k(p>29)$，则必定有不超过$29$的质因子$p^{\prime }$不能整除$x$。显然$x/p^k\times p^{\prime k}$的约数个数和$x$相同，但$x/p^k\times p^{\prime k}<x$所以根据$1$，发现与定义矛盾，所以质因数一定小于$29$。
• 反证法：若$x$的质因数不是连续的，同上用$p^{\prime k^{\prime }}$代替最大的质因数$p$的项$p^k$得到更小的值，与定义矛盾。
• 反证法：若$x$的质因数的指数不单调递减，则一定有相邻两项$p_1^a\times p_2^b$满足$p_1<p_2$，且$a>b$。而同上替换为$p_1^b\times p_2^a$，一定更小，且约数数目相同。与定义矛盾。

综上所述，我们只需要$DFS$搜索答案即可，因为实际上的搜索范围很小，所以并不会超时。

四、正解程序

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <algorithm>

using namespace std;
typedef long long ll;
const ll maxn=100010; 
ll prime[maxn/10],count1=0,n;
bool vis[maxn];
void oula()
{
	memset(vis,false,sizeof(vis));
	for(ll i=2;i<=maxn-10;i++)
	{
		if(!vis[i])
			prime[++count1]=i;
		for(ll j=1;j<=count1 && prime[j]*i<=maxn-10;j++)
		{
			vis[prime[j]*i]=true;
			if(i%prime[j]==0)
				break;
		}
	}
}
ll ans=0x7fffffff,Max=-1;
void dfs(ll k,ll val,ll last,ll times)
{
    if(k>16)
        return;
    if(Max<times)
    {
        Max=times;
        ans=val;
    }
    else if(Max==times)
        ans=min(val,ans);
    for(ll i=1;i<=last;i++)
    {
        ll temp=val*prime[k];
        if(temp>n)
            return;
        dfs(k+1,temp,i,times*(i+1));
        val*=prime[k];
    }
}
int main()
{
    oula();
    scanf("%lld",&n);
    dfs(1,1,100,1);
    printf("%lld\n",ans);
    
    return 0;
}