码迷,mamicode.com
首页 > 其他好文 > 详细

[BZOJ 1053] [HAOI 2007] 反素数ant

时间:2014-11-25 23:05:37      阅读:252      评论:0      收藏:0      [点我收藏+]

标签:blog   http   io   ar   os   使用   sp   for   on   

题目链接:BZOJ 1053

 

想一想就会发现,题目让求的 1 到 n 中最大的反素数,其实就是 1 到 n 中因数个数最多的数。(当有多于一个的数的因数个数都为最大值时,取最小的一个)

考虑:对于一个整数 n ,如何求 n 的因数的个数?

  将 n 分解质因数,n = p1^a1 * p2^a2 * p2^a3 * ...... * px^ax 。(其中 pi 为质因数, ai 为质因数的指数)

  那么 n 的因数的个数为 :Π (ai+1)     (使用组合数学的知识很容易看出)

那么我们想要找到 1 到 n 的数中因数最多的一个,我们可以使用质数相乘,直接搜索。

注意最后的答案 Ans = p1^a1 * p2^a2 * ...... * px^ax ,满足当 p1 < p2 < p3 ... < px 时,a1 <= a2 <= a3 ... <= ax 。

因为对于某两个素数 p1, p2 (p1 < p2) ,若 a2 > a1,那么交换 a1, a2 之后,因数的个数并没有改变,相乘得到的数字却减小了,因此更优。

另外重要的一点是,我们只会用到前 10 个素数,因为前 10 个素数相乘就已经超过 n 的最大范围,使用更大的素数是没有意义的,不会更优。

 

代码如下:

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cmath>

using namespace std;

const int INF = 0x3fffffff;
const int Prime[15] = {0, 2, 3, 5, 7, 11, 13, 17, 19, 23, 29};

typedef long long LL;

LL n, Ans, Cnt;

void DFS(LL Num, LL MaxP, int Now, int NowCnt) {
	if (NowCnt > Cnt || (NowCnt == Cnt && Num < Ans)) {
		Ans = Num;
		Cnt = NowCnt;
	}
	if (Now > 10) return;
	LL NowNum = Num;
	for (int i = 1; i <= MaxP; i++) {
		NowNum *= (LL)Prime[Now];
		if (NowNum > n) return;
		DFS(NowNum, i, Now + 1, NowCnt * (i + 1));
	}
}

int main() 
{
	scanf("%lld", &n);
	Ans = 1; Cnt = 1;
	DFS(1, INF, 1, 1);
	printf("%lld\n", Ans);
	return 0;
}

  

[BZOJ 1053] [HAOI 2007] 反素数ant

标签:blog   http   io   ar   os   使用   sp   for   on   

原文地址:http://www.cnblogs.com/F-Magician/p/4121948.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!