ZOJ 2562 HDU 4228 反素数

反素数：

对于不论什么正整数x,起约数的个数记做g(x).比如g(1)=1,g(6)=4.

假设某个正整数x满足:对于随意i(0<i<x),都有g(i)<g(x),则称x为反素数.

ZOJ 2562 反素数

由于写了POJ 2886的线段树。然后里面有反素数，曾经没遇到过，所以先搞这两题普及一下知识再说。

#include<iostream>

#include<cstdio>

#include<cstring>

#include<algorithm>

#include<map>

#include<queue>

#include<set>

#include<cmath>

#include<bitset>

#define mem(a,b) memset(a,b,sizeof(a))

#define lson i<<1,l,mid

#define rson i<<1|1,mid+1,r

#define INF 510010

#define maxn 400010

using namespace std;

typedef long long ll;

typedef unsigned long long ull;

ll prime[20]={2,3,5,7,11,13,17,19,23,29,31,37,39,41,43,47,53};

ll n;

ll bestcurr;//bestcurr 同样最大因数个数中值最小的数

ll largecnt;//largecnt：n范围内最大的因数个数

void getarcprime(ll curr,int cnt,int limit,int k)

{

    if(curr>n) return ;

    if(largecnt<cnt)//此时枚举到的因数个数比之前记录的最大的因数个数要大，就替换最大因数个数

    {

        largecnt=cnt;

        bestcurr=curr;

    }

    if(largecnt==cnt && bestcurr>curr)//替换最优值

        bestcurr=curr;

    ll temp=curr;

    for(int i=1;i<=limit;i++)

    {

        temp=temp*prime[k];

        if(temp>n) return;

        getarcprime(temp,cnt*(i+1),i,k+1);

    }

}

int main()

{

    while(scanf("%lld",&n)!=EOF)

    {

        bestcurr=0;

        largecnt=0;

        getarcprime(1,1,50,0);

        printf("%lld\n",bestcurr);

    }

    return 0;

}

HDU 4228

这题就是上题的延伸吧，就是求出每一个然后打表。

#include<iostream>

#include<cstdio>

#include<cstring>

#include<algorithm>

#include<map>

#include<queue>

#include<set>

#include<cmath>

#include<bitset>

#define mem(a,b) memset(a,b,sizeof(a))

#define lson i<<1,l,mid

#define rson i<<1|1,mid+1,r

#define INF 510010

#define maxn 400010

using namespace std;

typedef long long ll;

typedef unsigned long long ull;

ll p[1010];

ll prime[30]= {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53};

void getartprime(ll cur,int cnt,int limit,int k)

{

    //cur：当前枚举到的数；

    //cnt：该数的因数个数；

    //limit：因数个数的上限；2^t1*3^t2*5^t3……t1>=t2>=t3……

    //第k大的素数

    if(cur>(1LL<<60) || cnt>150) return ;

    if(p[cnt]!=0 && p[cnt]>cur)//当前的因数个数已经记录过且当时记录的数比当前枚举到的数要大，则替换此因数个数下的枚举到的数

        p[cnt]=cur;

    if(p[cnt]==0)//此因数个数的数还没有出现过，则记录

        p[cnt]=cur;

    ll temp=cur;

    for(int i=1; i<=limit; i++) //枚举数

    {

        temp=temp*prime[k];

        if(temp>(1LL<<60)) return;

        getartprime(temp,cnt*(i+1),i,k+1);

    }

}

int main()

{

    int n;

    getartprime(1,1,75,0);

    for(int i=1; i<=75; i++)

    {

        if(p[i*2-1]!=0 && p[i*2]!=0)

            p[i]=min(p[i*2-1],p[i*2]);

        else if(p[i*2]!=0) p[i]=p[i*2];

        else p[i]=p[i*2-1];

    }

    while(scanf("%d",&n),n)

        printf("%I64d\n",p[n]);

    return 0;

}