牛客网 14 n的约数（数论）

题意：中文，就是求1~n中约数最多的数

思路：根据唯一分解定理，我们可以知道一个数n = a1^p1*a2^p2...的约数个数为 (1+p1) * (1+p2)....，那么1~n中约数最多的数我们应该怎么求呢？根据公式我们可以知道我们让p1,p2,,,,pn最大就好了，那么我们现在可以枚举一个素数，枚举到他的次方大于n的时候，我们回溯开始乘下一个素数之后我们记录一下当前枚举素数的指数的最大值就好了，这里有一个小优化，就是说你当前枚举的素数的最大指数一定会比你上一个枚举的素数指数小（证明？显然啊（其实我也不会证，敏锐的数学猜想把。。））上代码把:

#include<iostream>
#include<algorithm>
#include<cstdio>
using namespace std;
long long ans;
long long ansnum;//约数个数
long long n;
long long prime[24]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83};
void dfs(long long p,unsigned long long sum,long long cnt,long long up)//p表示我们现在再乘的素数，sum表示的是我们凑出来的数字
//，cnt表示的是我们当前的素数它的个数 ,up表示的是我们枚举素数的上届 
{
	ansnum = max(ansnum,cnt);
	int te = 1;
	for(long long i = 1 ; i <=up ;)
	{	
		if(sum > n / prime[p] ) break; //这里特别gay。。因为n超级大，所以有可能出现超过unsigned long long 的情况，这个就是判断他是否超过unsigned long long 
		sum*=prime[p]; 
		if(sum<0) break;
		if(sum>n) break;
		i++;
		dfs(p+1,sum,cnt*i,i-1);
	}
}
int main()
{
	int t;
	cin>>t;
	while(t--)
	{
		ansnum = 0;
		cin>>n;
		dfs(0,1,1,60);
		cout<<ansnum<<endl;
	}
}