BZOJ 1053: [HAOI2007]反素数ant

1053: [HAOI2007]反素数ant

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 3697  Solved: 2162
[Submit][Status][Discuss]

Description

  对于任何正整数x,其约数的个数记作g(x)。例如g(1)=1、g(6)=4。如果某个正整数x满足:g(x)>g(i) 0<i<x,则称x为反质数。例如,整数1,2,4,6等都是反质数。现在给定一个数N,你能求出不超过N的最大的反质数么?

Input

  一个数N(1<=N<=2,000,000,000)。

Output

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

Sample Input

1000

Sample Output

840

题解

根据唯一分解,x=p1^k1*p2^k2*…*pn^kn,所以g(x)=(k1+1)*(k2+1)*…*(kn+1)。

根据题目要找约数个数大于前面所有数的约数个数的最大的数,那么如果k1>k2,可以交换k1,k2找到一个约数个数相等但是小于它的数,不满足题目条件。所以必须满足k1<=k2<=…<=kn。

发现质数从2乘到29已经达到了2000000000,所以枚举2到29的质数就可以了。

代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<iostream>
#define LL long long
using namespace std;
int n,ans,mx;
int prime[15]={0,2,3,5,7,11,13,17,19,23,29};
void dfs(int num,int cnt,int now,int last){
	if(prime[now]==0){
		if(cnt>mx){
			ans=num;
			mx=cnt;
		}
		if(cnt==mx&&num<ans)ans=num;
		return;
	}
	LL temp=1;
	for(int i=0;i<=last&&num*temp<=n;i++){
		dfs(num*temp,cnt*(i+1),now+1,i);
		temp*=prime[now];
	}
}
int main(){
	scanf("%d",&n);
	dfs(1,1,1,32);
	printf("%d\n",ans);
	return 0;
}

转载于:https://www.cnblogs.com/chezhongyang/p/7729635.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值