CSP-因子化简 满分题解

题目回顾

因子化简https://sim.csp.thusaac.com/contest/32/problem/1

C++满分题解

思路概述

         任意一个正整数n，有若干素因数，n可以表示为:

可采用数据结构map<int, int>来存储以上表示，key代表底数，value代表次数。

核心步骤：判断素数

        流程框图中的核心步骤是判断素数，若不是素数则返回一个素因数。这一步是时间复杂度的关键所在。这里可以使用一个结论来提速：

如果一个正整数n(n>3)是素数，则 n mol 6 余 1 或 n mol 6 余 5

        简要证明如下：

        若 n mol 6 余 2 or 4 or 0，n一定是偶数，若n mol 6 余 3，则n是3的倍数，有素因数3。故若一个大于3的正整数是素数，他与6相除所得余数一定是1 or 5；

        由此可设计核心步骤算法：判断一个数n是否为素数，若是返回0，否则返回一个素因数。

        在otherwise情况下，另i=2，累加至sqrt(n)，若存在能够整除n的i，则n不是素数，并且i是n的素因数。否则，n是素数。 这里有一个值得讨论的点：为什么i一定是n的素因数？

        可以假设在 i 的累加过程中第一个能够整除 n 的数如果是合数，那么这个i的因数 p 一定能够整除 n，并且这个因数 p<i ，所以第一个能整除 n 的 i 一定是素数。

完整代码

#include<iostream>
#include<stdio.h>
#include<map>
#include<fstream>
#include<string>
#include<math.h>
using namespace std;

int isSushu(long long num) {
	if (num == 2 || num == 3)
		return 0;
	else if (num % 6 == 2 || num % 6 == 4) {
		return 2;
	}
	else if (num % 6 == 3 || num % 6 == 0) {
		return 3;
	}
	else {
		long long end = sqrt(num);
		for (int i = 2; i <= end; i++) {
			if (num % i == 0)
				return i;
		}
		return 0;
	}
}

long long simplify(long long n, int k) {
	map<int, int> dishu_cishu;
	long long num = n;
	int dishu = isSushu(num);
	while (dishu != 0) {
		dishu_cishu[dishu] += 1;
		num = num / dishu;
		dishu = isSushu(num);
	}
	dishu_cishu[num] += 1;
	map<int, int>::iterator it;
	long long result = 1;
	for (it = dishu_cishu.begin(); it != dishu_cishu.end(); it++) {
		if (it->second >= k) {
			result *= pow(it->first, it->second);
		}
	}
	return result;
}

int main() {
	int q, k;
	long long n;
	/*ifstream infile;
	infile.open("D://C,C++程序设计/MyProject/MyProject/input.txt", ios::in);
	if (!infile.is_open())
	{
		cout << "读取文件失败" << endl;
		return 0;
	}*/
	//infile >> q;
	cin >> q;
	while (q != 0) {
		q--;
		//infile >> n;
		//infile >> k;
		cin >> n >> k;
		cout << simplify(n, k) << endl;
	}
}