P1463 [POI2001] [HAOI2007] 反素数

[POI2001] [HAOI2007] 反素数

题目描述

对于任何正整数 $x$，其约数的个数记作 $g(x)$。例如 $g(1)=1$，$g(6)=4$。

如果某个正整数 $x$ 满足：$\forall 0<i<x$，都有 $g(x)>g(i)$，则称 $x$ 为反质数。例如，整数 $1,2,4,6$ 等都是反质数。

现在给定一个数 $N$，你能求出不超过 $N$ 的最大的反质数么？

输入格式

一个数 $N$。

输出格式

不超过 $N$ 的最大的反质数。

样例 #1

样例输入 #1

1000

样例输出 #1

840

提示

$1\le N\le 2\times 10^9$

分析：

由于质因数分解定理：$$设x=\prod _{i=1}^m{p}_i^{k_i},p_i\in P$$
则根据乘法原理：$$g(x)=\prod _{i=1}^m(k_i+1)$$
$由此可以看出反素数x的p_i是单增的，那么k_i一定是不升的$

证明如下：

证明：

$$即x=\prod _{i=1}^m{p}_i^{k_i},p_i\in P中$$
$存在p_1<p_2<\dots <p_m,k_1\ge k_2\ge \dots \ge k_m$
$设k_1<k_2，x_1=p_1^{k2}\ast p_2^{k1}\dots \ast p_m^{km}=(p_1{p}_2{)}^{k1}\ast p_1^{k2-k1}\ast ...$
$又因为x_2=p_1^{k1}\ast p_2^{k2}\ast ...\ast p_m^{km}=(p_1{p}_2{)}^{k1}\ast p_2^{k2-k1}\ast \dots$
$可知x_1<x_2且约数个数相等，所以x_2必然不是反素数$

code:

#include<bits/stdc++.h>
using namespace std;
#define i64 long long 
#define pii pair<i64,i64> 
#define il inline
#define re register
const i64 MAXN=1e7+5;
int n;
i64 p[]={0,2,3,5,7,11,13,17,19,23,29};//素数表
i64 tot;
pii a[MAXN];

//dfs
il void dfs(i64 k,i64 last,i64 pi,i64 g)//k：第几个质数   last：上一个质数的指数   pi：当前的因数   g：当前约数的个数 
{
	a[++tot]=make_pair(pi,g);//存(i,g(i)) 
	if(k>9) return;
	for(re int j=1;j<=last;j++)//指数不升 
	{
		if(pi*p[k]>n) break;
		pi*=p[k];
		dfs(k+1,j,pi,g*(j+1));
	}
}
//

int main()
{
	cin>>n;
	dfs(1,100,1,1);
	sort(a+1,a+1+tot);
	int ans=a[1].first,maxn=a[1].second;
	for(re int i=2;i<=tot;i++) if(a[i].second>maxn) {ans=a[i].first;maxn=a[i].second;}
	cout<<ans; 
	return 0;
}

ps:

由于数据量只有2e9，而前10个素数的乘积已经大于2e9了，所以只考虑前10个素数即可

end