牛客小白月赛23 B-阶乘

牛客小白月赛23 B-阶乘

题目链接（Click）
时间限制：C/C++ 1秒，其他语言2秒
空间限制：C/C++ 262144K，其他语言524288K
64bit IO Format: %lld

题目描述

给定一个正整数 p
求一个最小的正整数 n，使得 n! 是 p 的倍数

输入描述:

第一行输入一个正整数 T 表示测试数据组数
接下来 T 行，每行一个正整数 p

输出描述:

输出 T 行，对于每组测试数据输出满足条件的最小的 n

样例输入

4
1
2
4
8

样例输出

1
2
4
4

数据范围

T≤10³ , p≤10⁹

分析

这题目短得可怕，信息量确实巨大，让我想起了当年被数学题支配的恐惧。
一看这数据范围 T≤10³ p≤10⁹ ，再一看这阶乘 n! ，当场去世
不仅时间上捉襟见肘，空间上更是杯水车薪
如果你有那么一秒想过朴素算法，编译器可能都脊背发凉
不过仔细一想，也对，这不就是ACM的尿性嘛…

真·分析

这该如何是好，再来看看题目：求一个最小的正整数 n，使得 n! 是 p 的倍数
n! 是 p 的倍数。倍数，如何定义倍数呢，我们经常用 a%b==0 来判断a是b的倍数
但眼下，我们显然不可能求出n!
那么换个思路，什么才叫做倍数，a是b的整数倍，也就是说a=b*c (c为正整数)
即，a的因子中，刨去构成b的部分，剩下的部分是一个正整数，即可

先对p进行质因数分解，然后将这些因子与n!的因子进行匹配即可
两种情况
第一种是：恰好匹配，比如2匹配2
第二种是：倍数匹配，比如2*2 ，第一个2匹配n!中的2，第二个2匹配n!中的4，(4%2==0)

口说无凭，举个栗子
6!=1×2×3×4×5×6
180=2×2×3×3×5
2->6!的2
2->6!的4
3->6!的3
3->6!的6
5->6!的5

注意：n!中的因子必须要是匹配数的倍数，这样整体才是倍数关系
比如:32=2×2×2×2×2
2->n!的2
2->n!的4
2->n!的4
2->n!的6
2->n!的8 (注意，这里不能匹配6，因为匹配后6%4!=0)

具体实现上要用map保存已匹配数，如m[6]=2，代表n!中的6已经被匹配了2

反正我也讲不清楚了，先上代码…
"Talk is Cheap. Show me the Code."

#include<stdio.h>
#include<stdlib.h>
#include<algorithm>
#include<string.h>
#include<string>
#include<vector>
#include<math.h>
#include<queue>
#include<map>
using namespace std;
int main(void)
{
	int T;
	scanf("%d", &T);
	while (T--)
	{
		int p;
		scanf("%d", &p);
		vector<int>fact;
		int lim = sqrt(p);
		for (int i = 2; i <= lim; i++) {//质因数分解 复杂度O(logn)
			if (p == 1)break;
			if (p % i == 0) {
				do {
					fact.push_back(i);
					p /= i;
				} while (p % i == 0);
			}
		}
		fact.push_back(p);//假如因子含有质数，少不了这一步
		map<int,int>m;
		for (auto n : fact)
			if (m.find(n) == m.end())
				m[n] = n;
			else {
				int i = 2;
				while (true) {//寻找合适匹配
					if (m.find(n * i) == m.end()) {
						m[n * i] = n;
						break;
					}
					else if (m[n * i] <= i && ((n*i)%(m[n*i]*n)==0)) {//必须是倍数关系
						m[n * i] *= n;
						break;
					}
					i++;
				}
			}
		int Max = 0;
		for (auto n : m)
			Max = max(Max, n.first);
		printf("%d\n", Max);
	}
	return 0;
}