算法训练 素因子去重 c语言,【蓝桥杯】算法训练 素因子去重

算法训练 素因子去重

时间限制：1.0s   内存限制：256.0MB

问题描述

给定一个正整数n，求一个正整数p，满足p仅包含n的所有素因子，且每个素因子的次数不大于1

输入格式

一个整数，表示n

输出格式

输出一行，包含一个整数p。

样例输入

1000

样例输出

10

数据规模和约定

n<=10^12

样例解释：n=1000=2^3*5*3，p=2*5=10

题目还是很好理解的，一看到题，我的第一个思路就是先打印素数表，然后直接用整数的质因子分解来做，结果是耗时15ms，内存984.0KB。关于这两个过程可以戳以前的博客：埃氏筛法打印素数表   算术基本定理：整数质因子分解

1 #include

2 using namespace std;

3

4 bool is_prime[10005];

5 int prime[10005];

6 void sieve(int n){

7 int p = 0;

8 for(int i = 2; i <= n; i++) is_prime[i] = true;

9 for(int i = 2; i <= n; i++){

10 if(is_prime[i]){

11 prime[p++] = i;

12 for(int j = 2*i; j<=n; j+=i) is_prime[j] = false;

13 }

14 }

15 }

16

17 int main(int argc, char *argv[]) {

18 sieve(10005);

19 long long n, ans = 1;

20 cin>>n;

21 for(int i = 0;prime[i]*prime[i]<=n; i++){

22 if(n % prime[i] == 0){

23 ans *= prime[i];

24 while(n % prime[i] == 0) n /= prime[i];

25 }

26 }

27 ans *= n;

28 cout<

29

30 return 0;

31 }

后来发现自己把题目理解复杂了，下面给出我的第二个思路，耗时46ms，内存936.0KB。

1 #include

2 using namespace std;

3

4 int main(int argc, char *argv[]) {

5 long long n, ans = 1;

6 cin>>n;

7 for(int k = 2; n != 1; k++){

8 if(n % k == 0){

9 ans = ans * k;

10 while(n % k == 0) n /= k;

11 }

12 }

13 cout<

14

15 return 0;

16 }

这个方法也很直接，从最小的质数2开始循环，判断当前的数k能不能整除n，如果可以就说明k是n的一个质因数，让ans乘上k，再让n除以k直到n不再包含k为止。每次循环k++。有个巧妙的地方在于，因为每次循环结束得到的新的n都已经不再包含上一个质因子k，所以下一次满足整除条件的k一定还是质数(有一点埃氏筛法的感觉)。比如取n=1000，它最小的质因子是2，在k=2的这一轮循环结束后，新的n为125，已经不再包含2，所以k=4时，n已经不再能被4整除，直接进入下一轮循环k=5。这一轮结束后，n=1，说明所有的质因子已经被找到，循环结束。

第三种方法是在网上看到的题解，我在循环那边优化了一下，也很简单直接，有一个素性判断的is_prime函数，然后直接循环判断质因子。下面是代码，耗时0ms，内存936.0KB。

1 #include

2 using namespace std;

3

4 bool is_prime(long long n){ //判断n是否为素数

5 for(long long i = 2; i*i <=n; i++)

6 if(n%i ==0) return false;

7 return true;

8 }

9

10 int main(int argc, char *argv[]) {

11 long long n, m = 1;

12 cin>>n;

13 for(long long i = 2; i*i<=n; i++){

14 if(n%i == 0 && is_prime(i)){

15 m *= i;

16 while(n%i==0) n/=i;

17 }

18 }

19 m *= n;

20 cout<

21

22 return 0;

23 }

最后一种也大同小异，用到了集合set元素，网上看到的题解。耗时46ms，内存932.0KB。(其实几种都差不多)

1 #include

2 #include

3 using namespace std;

4

5 set s;

6 bool isPrime(long long n){

7 for(int i=2; i*i <= n; i++)

8 if(n%i == 0) return false;

9 return true;

10 }

11

12 int main(){

13 long long n, ans=1;

14 cin>>n;

15 for(int i=2; n!=1; i++)

16 if(n%i==0 && isPrime(i)){

17 s.insert(i);

18 while(n%i==0) n/=i;

19 }

20 for(set::iterator it=s.begin(); it!=s.end(); it++)

21 ans *= (*it);

22 cout<

23

24 return 0;

25 }

标签：prime,int,因子,long,蓝桥,算法,ans,include

来源： https://www.cnblogs.com/Aikoin/p/10371719.html