反素数[BZOJ1053]

题面描述

传送门

思路

三条引理:

1. 设$m$是1~N中约数个数最多的数中最小的一个。根据m的定义，m显然满足:

   1. $\forall x<m,g(x)<g(m)$
   2. $\forall x>m,g(x)\le g(m)$

   根据反质数定义，1说明了$m$是反质数，2说明了大于$m$的数都不是反质数，

   故$m$即为所求。

2. $1$~$N$中任何数的不同质因子都不会超过$10$个，且所有质因子的指数总和不超过$30$.
   证明：
   因为最小的$11$个质数的乘积$2\ast 3\ast 5\ast 7\ast 13\ast 17\ast 19\ast 23\ast 29\ast 31>2\ast 10^9$
   所以$N\le 2\ast 10^9$不可能多于有$10$个不同的质因子。
   即使仅包含最小质数$2$，$2^{31}>2\ast 10^9$，故指数总和不超过$30$。

3. $\forall x\in [1,N]$,$x$为反质数的充要条件是: x分解质因数可以写作:
   $$2^{c_1}\ast 3^{c_2}\ast 5^{c_3}\ast ......\ast 29^{c_{10}}，c_1\ge c_2\ge c_3\ge ......\ge c_{10}\ge 0$$证明：若$n$存在一大于$29$的质因子$p^k$，由引理1可知，存在整数$m$，约数个数与$n$相同，且小于$n$.
   若$c_i<c_j,i<j$，同理，亦存在整数比$n$小.

综上所述，
我们可以用DFS，尝试确定前十个质数的质数，满足引理条件，总乘积小于N，同时记录约束个数。
每搜索出一个符合条件的数，用引理1更新答案。

AC code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#define ll long long
using namespace std;
const int N=12;
const int inf=2e9;
const int prime[N]={0,2,3,5,7,11,13,17,19,23,29};
int c[N],n;ll ans;int ansc;
int calc(int len)
{
	int k=1;
	for(int i=1;i<=len;i++)k*=(c[i]+1);
	return k;
}
ll get(int now,int len)
{
	ll k=1;
	for(int i=1;i<=len;i++){if(k*prime[now]>n)return 0;k*=prime[now];}
	return k;
}
void dfs(int last,int sum,ll s)
{
	if(last>=10)return ;
	if(sum>ansc)ansc=sum,ans=s;
	if(sum==ansc&&ans>s)ans=s;
	for(int i=c[last];i>=0;i--)
	{
		ll k=get(last+1,i);
		if(k==0||s*k>n)continue;
		c[last+1]=i;
		dfs(last+1,calc(last+1),s*k);
		c[last+1]=0;
	}
}
int main()
{
	c[0]=30;
	scanf("%d",&n);ans=1,ansc=1;
	dfs(0,1,1);
	printf("%lld\n",ans);
	return 0;
}