您的位置:首页 > 其它

反素数 模板 求因子的个数为n的最小的数是什么

2013-05-01 22:51 274 查看
反素数:

如果一个自然数比所有比它小的自然数的约数个数都要多,那么我们就称这个数为一个反素数。例如,1、2、4、6、12和24都是反素数。

不大于n的最大的反素数模板:

#include <math.h>
#include <iostream>
using namespace std;
#define  ll __int64
const ll prime[16] = {1, 2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47};
ll maxsum, bestnum, n;
void dfs(ll num, ll k, ll sum, ll limit)
//num:当前枚举到的数,k:枚举到的第k大的质因子;sum:该数的约数个数;limit:质因子个数上限;
{
int i;
ll tmp;
if (sum > maxsum)
{
maxsum = sum;
bestnum = num;//如果约数个数更多,将最优解更新为当前数;
}
if (sum == maxsum && bestnum > num)
{//如果约数个数相同,将最优解更新为较小的数;
bestnum = num;
}
if (k > 15) return ;
tmp = num;
for (i = 1; i <= limit; ++i)
{//开始枚举每个质因子的个数
if (tmp * prime[k] > n) {
break;
}
tmp *= prime[k];//累乘到当前数
dfs(tmp, k + 1, sum * (i + 1), i);//继续下一步搜索
}
}
int main()
{
ll i, j;
int n;
cin>>n;
dfs(1, 1, 1, 50);
printf("%lld\n", bestnum);
return 0;
}


下面的是打表法 可以求1-1000 甚至更多的反素数 (只要改下数组大小就好)

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
typedef __int64 lld;
lld p[1010];//p[i] 表示为因子个数为i的最小整数是什么
lld prime[30]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53};
int  maxn;//最大的可能的值
void getartprime(lld cur,int cnt,int limit,int k)
{
//cur:当前枚举到的数;
//cnt:该数的因数个数;
//limit:因数个数的上限;2^t1*3^t2*5^t3……t1>=t2>=t3…… 一般取50就够了
//第k大的素数
//if(cur>((lld)1<<60) || cnt>150) return ;
if(cur>maxn) return ;//如果当前的数大于我们要求的 最大的数 maxn 就寻找完毕了
if(p[cnt]!=0 && p[cnt]>cur)//当前的因数个数已经记录过且当时记录的数比当前枚举到的数要大,则替换此因数个数下的枚举到的数
p[cnt]=cur;
if(p[cnt]==0)//此因数个数的数还没有出现过,则记录
p[cnt]=cur;
lld temp=cur;
for(int i=1;i<=limit;i++)//枚举数
{
temp=temp*prime[k];
if(temp>maxn) return;
getartprime(temp,cnt*(i+1),i,k+1);

}
}
int main()
{
getartprime(1,1,50,0);
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: