（第十届蓝桥杯国赛）D：求值（数学/暴力）

题意：学习了约数后，小明对于约数很好奇，他发现，给定一个正整数 t，总是可以找到含有 t 个约数的整数。小明对于含有 t 个约数的最小数非常感兴趣，并把它定义为 S t 。
例如 S1 = 1, S2 = 2, S3 = 4, S4 = 6，· · · 。
现在小明想知道，当 t = 100 时，S t 是多少？即 S100 是多少？

答案：45360

首先说明一点，这道题目可以用暴力来做，等会我会附上暴力代码，这道题的数据范围较小，所以暴力能很快求出来，但是我下面会讲述一种通用的数学求解方法。

比如要求Sx,也就是求约束个数为x的最小数，我们先来看一下对于任意一个数W，他的约数个数是多少，假设W的质因子分解式是W=a1^p1+a2^p2+……+an^pn,那么W的约数个数就是(p1+1)*(p2+1)*……*(pn+1),也就是说W的约数个数仅与W的质因子分解式中的质数的幂次有关，而与具体的质数无关，那么根据贪心性质我们显然可以知道如果pi越大，那么我们就应该让他对应的质数越小，这样才会使得最终的值最小，也就是先对pn……pn进行从大到小进行排序，然后依次令2，3，5……这些质数对应于相应的幂次即可求得最优解。就拿该题举个例子：100=(1+1)*(1+1)*(4+1)*(4+1),四个幂次分别为1，1，4，4，从大到小排序为4，4，1，1，对应于2，3，5，7得2^4*3^4*5^1*7^1=45360,这就是利用贪心求解的思路，这样的复杂度要比暴力小很多，而且能求解的范围也会比暴力大很多。

代码（输入x求Sx）：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<vector>
using namespace std;
const int N=10003;
int prime[N],cnt;
bool vis[N];
vector<int>p;
void init()//打素数表 
{
	for(int i=2;i<N;i++)
	{
		if(!vis[i]) prime[++cnt]=i;
		for(int j=1;j<=cnt&&i*prime[j]<N;j++)
		{
			vis[i*prime[j]]=true;
			if(i%prime[j]) break;
		}
	}
}
void f(int x)//求x质因子分解的幂次，存入p中
{
	for(int i=1;prime[i]<=sqrt(x);i++)
	{
		while(x%prime[i]==0)
		{
			p.push_back(prime[i]-1);
			x/=prime[i];
		}
	}
	if(x!=1)  p.push_back(x-1);//考虑x最后剩余一个质数的情况 
}
long long qmi(long long a,long long b)
{
	long long ans=1;
	while(b)
	{
		if(b&1) ans*=a;
		a*=a;
		b>>=1;
	}
	return ans;
}
long long dfs()
{
	long long ans=1;
	sort(p.begin(),p.end());
	for(int i=p.size()-1,j=1;i>=0;i--,j++)
		ans*=qmi(prime[j],p[i]);
	return ans;
}
int main()
{
	init();
	int x;
	cin>>x;
	f(x);//求x质因子分解的幂次，存入p中 
	cout<<dfs();
	return 0;
}

暴力方法：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<vector>
using namespace std;
int f(int x)
{
	int cnt=0;
	for(int i=1;i*i<=x;i++)
	{
		if(x%i==0)
		{
			if(i*i!=x) cnt+=2;
			else cnt++;
		}
	}
	return cnt;
}
int main()
{
	for(int i=100;;i++)
		if(f(i)==100)
		{
			printf("%d",i);
			break;
		}
	return 0;
}