Andrew Stankevich Contest 10 D More Divisors 反素数

最新推荐文章于 2019-07-21 19:34:58 发布

R灬O灬J

最新推荐文章于 2019-07-21 19:34:58 发布

阅读量393

点赞数

分类专栏：素数文章标签： ACM C++ 数学

本文链接：https://blog.csdn.net/R_O_J/article/details/52524444

版权

素数专栏收录该内容

1 篇文章 0 订阅

订阅专栏

今天学了个新东西，反素数，感觉这玩意主要是套模板。

模板：http://www.cnblogs.com/tiankonguse/archive/2012/07/29/2613877.html

以下套用大牛的讲解。

反素数的定义：对于任何正整数，其约数个数记为，例如，如果某个正整数满足：对任意的正整

数，都有，那么称为反素数。

从反素数的定义中可以看出两个性质：

（1）一个反素数的所有质因子必然是从2开始的连续若干个质数，因为反素数是保证约数个数为的这个数尽量小

（2）同样的道理，如果，那么必有

代码：

#include<bits/stdc++.h>
#define MEM(a,x) memset(a,x,sizeof(a));
#define MEMINF(a) memset(a,0x3f,sizeof(a));
using namespace std;
typedef long long LI;
const int MAXN=205;
const int INF=0x3f3f3f3f;
const int MOD=1000000007;
LI n;
LI bestSum,bestNum;
LI prim[15]={2,3,5,7,11,13,17,19,23,29,31,37,43,47};
void getNum(LI num,LI k,LI sum, LI limit) {
  if (num>n) return;
  if (sum>bestSum) {
    bestSum=sum;
    bestNum=num;
  } else if (sum==bestSum&&num<bestNum) {
    bestNum=num;
  }
  if (k>=13) return ;
  for (LI i=1,p=num; i<=limit; ++i) {
    if(p*prim[k]>n)break;
    p*=prim[k];
    getNum(p,k+1,sum*(i+1),i);
  }
}

LI log2(LI n) {
  LI i=0;
  LI p=1;
  while (p<n) {
    p*=2;
    i++;
  }
  return i;
}


int main() {
  freopen("divisors.in","r",stdin);
  freopen("divisors.out","w",stdout);
  cin>>n;
  bestNum=n;
  bestSum=0;
  getNum(1,0,1,log2(n));
  cout<<bestNum<<endl;
}