[SDOI2005]反素数ant 洛谷p1463

题目描述

对于任何正整数x，其约数的个数记作g(x)。例如g(1)=1、g(6)=4。

如果某个正整数x满足：g(x)>g(i) 0<i<x，则称x为反质数。例如，整数1，2，4，6等都是反质数。

现在给定一个数N，你能求出不超过N的最大的反质数么？

输入输出格式

输入格式：

一个数N（1<=N<=2,000,000,000）。

输出格式：

不超过N的最大的反质数。

输入输出样例

输入样例#1：

1000

输出样例#1：

840

思路：

这道题实际上是求1~n中因数最多的数中最小的。

因为我们知道，任意一个大于1的整数都可以拆分成多个质数（或质数的幂）相乘

所以我们可以用前几个质数（10几个就够了【标签的打表题应该就是在这里打表了】）构建反素数。

为什么是十几个呢？我们看一下题目的数据范围，一个数N（1<=N<=2,000,000,000）。

【而前9个质数2*3*5*7*11*13*17*19*23*29=6,469,693,230>2,000,000,000】所以实际上只需要前10个质数就可以啦

那么我们每个数S必然可以分解为S=2^x1+3^x2+...+p^xn(p仍旧为质数）

因数的总数cnt=(x1+1)*(x2+1)*(x3+1)*……

我去网上找了个证明：http://blog.csdn.net/baidu_36786519/article/details/53727982

我们需要构建其中因数最多的。因数最多的若有多个则取最小的其中一个并且输出。

但是x1,x2,x3这样的的数怎么枚举呢？

我们不难发现，x1,x2,x3是严格不上升的。为什么呢？

证明：若存在两个数，a^x1+b^x2,和a^x2+b^x1。a<b,x1>x2。两个数的因数个数相同，但是第一个比较小，所以选择第一个。

所以我们这题就大概这么结束了！

#include<iostream>
#include<cstdio>
using namespace std;
int prime[13]={0,2,3,5,7,11,13,17,19,23,29,31,33},used[15];
long long n,best,ans=10000000000,Pow[13][40],TIMES=0;
void MakeTable()
{
	int i,j;
	for(i=1;i<=12;i++){
		Pow[i][1]=prime[i];
		for(j=2;Pow[i][j-1]*prime[i]<=n;j++){
			Pow[i][j]=Pow[i][j-1]*prime[i];
		}
	}
}
void dfs(int now,long long num,long long a)
{
	int i;
	TIMES++;
//	if(TIMES>10) return;
	if(now>11) return;
	if(num>best||(num==best&&a<ans)){
		best=num;
		ans=a;
	}
	for(i=1;i<=used[now-1];i++){
		used[now]=i;
		if(a*Pow[now][i]<=n) dfs(now+1,num*(i+1),a*Pow[now][i]);
		else break;
	}
}
int main()
{
	int i;
	scanf("%lld",&n);
	MakeTable();
	for(i=1;;i++){
		if(!Pow[1][i]) break;
	}
	used[0]=i-1;
	dfs(1,1,1);
	cout<<ans<<endl;
	return 0;
}