题目链接:https://www.acwing.com/problem/content/description/200/
本作品采用知识共享署名-相同方式共享 4.0 国际许可协议进行许可。
对于任何正整数x,其约数的个数记作g(x),例如g(1)=1、g(6)=4。
如果某个正整数x满足:对于任意的小于x的正整数 i,都有g(x)>g(i) ,则称x为反素数。
例如,整数1,2,4,6等都是反素数。
现在给定一个数N,请求出不超过N的最大的反素数。
输入格式
一个正整数N。
输出格式
一个整数,表示不超过N的最大反素数。
数据范围
1≤N≤2∗109
输入样例:
1000
输出样例:
840
分析:
要不是有书,给我三十年,我也做不出这道题。 好菜呀。 关键有书我还能wa。
讲一下自己的理解吧。
明确几点概念:
正整数N被唯一分解为N=p1^c1* p2^2 * ~~pn ^cn
我们可得:
此正整数约数的个数为:
(c1 + 1) * (c2 +1) ~ ~(cn + 1);
第一反应求出1到N里面所有数的约数个数,并从左到右进行扫描找到最大的约数个数,如果最大数有多个,那应该是数小的那一个。到此方法必超时。
但可得引理一:
1到N中最大的反质数,就是1到N中约数个数最多的数中的最小一个数。
引理二:
看N的范围为1到2e9.那么1到N中任何数的不同质因子个数一定不会超过十个
eg:
2 * 3 * 5 * 7 * 11 * 13* 17 * 19 * 23 * 29 * 31 > 2e9.
那么同理,我们可的最高次幂不会超过31
eg:
2^31 > 2e9
引理三:
我们要求的那个数势必是下面形式:
2^c1 * 3^c2 * ~~* 31^cn。并且c1>= 2 >= cn
证明:不是非上升幂。那么我们通过交换质因子的方法,可以得到一个比他小的数。
最后在跑2e9的时候超时了。原因归咎是sum爆了longlong,成为了负数,所以特别加了一个判定是否为负数的情况。
#include"string.h"
#include"math.h"
#include"algorithm"
using namespace std;
typedef long long ll;
ll N,m,max_cnt;
ll prime[12] = {0,2,3,5,7,11,13,17,19,23,29,31};
void dfs(int id,int x,ll sum,ll cnt)
{
if(sum > N || sum < 0) return ;
//printf("sum = %lld m = %lld\n",sum,m);
if(x == 0 || id > 11)
{
if(cnt > max_cnt)
{
m = sum; max_cnt = cnt;
}
else if(cnt == max_cnt)
m = min(m,sum);
return ;
}
dfs(id,x - 1,sum,cnt);
dfs(id + 1,x,sum * (ll)pow(prime[id],x),cnt * (x + 1));
}
int main()
{
scanf("%lld",&N);
m = N + 1;
max_cnt = 1;
dfs(1,30,1LL,1LL);
printf("%lld\n",m);
}