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

反素数:

如果一个自然数比所有比它小的自然数的约数个数都要多,那么我们就称这个数为一个反素数。例如,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;
}


 

 

 

 

 

posted @ 2013-05-01 22:51  javawebsoa  Views(273)  Comments(0Edit  收藏  举报