牛客小白月赛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≤103 , p≤109
分析
这题目短得可怕,信息量确实巨大,让我想起了当年被数学题支配的恐惧。
一看这数据范围 T≤103 p≤109 ,再一看这阶乘 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;
}