CCF刷题计划——因子化简

计算机软件能力认证考试系统

我的蒟蒻暴力法:

很暴力的算法,但是过了。主要就是想着将所有的素数都先dia出来,然后将n对每一个素数都进行一次相爱相杀,同时记录用上的素数和指数。顺便学了个素数筛,之所以把全部素数都提前弄出来,其实就是为了防止多次重复计算素数。其实看到范围限制的时候就可以放弃这个策略了,因为全部的测试数据满足:1<n<=1e10 且 1<k,q≤10,相当于最多只重复个10次而已。

cpp 复制代码
#include <iostream>
#include <cmath>
#include <vector>
using namespace std;
const long long N=1e8;
#define ll long long

ll n;
int k,q;

vector<bool>statue(N,true);	//存放,是否为素数,默认全为素数 
vector<int>primes;

void getVec(vector<pair<int,int>>&vec,ll n)
{
	for(int i=0;i<primes.size();i++)
	{
		int a=primes[i];
		int b=0;
		if(n%a==0)
		{
			while(n%a==0)	//如果是倍数
			{
				n/=a;
				b++;
			}
			vec.push_back({a,b});	//将素数和指数加进去 
		}
	}
}

//欧拉筛,还是现学的
void euler_sieve()	//这里使用欧拉筛将所有的素数罗列 
{
	statue[0]=statue[1]=false;	//默认0,1都不是素数
	
	for(ll i=2;i<=N;i++)	//从2到N遍历一遍
	{
		if(statue[i])	//如果当前是素数
			primes.push_back(i);	//存起来 
		
		for(ll j=0;j<primes.size();j++) //不论是不是素数,都要将该数i和素数集都乘一次
		{
			if(i*primes[j]>N)	break;	//如果超出标记范围,就跳出 
			
			statue[i*primes[j]]=false;	//标记合数 
			
			if(i%primes[j]==0)	break;	//如果标记到最小素倍数(i是primes[j])那就溜啦
			//这一步是为了不重复标记 
		 } 
	}
}

int main()
{
	cin>>q;
	euler_sieve();
	while(q--)
	{
		cin>>n>>k;
		ll div=1;
		vector<pair<int,int>>vec;
		getVec(vec,n);
		for(int i=0;i<vec.size();i++)
		{
			if(vec[i].second<k)	//逆向,除去 
				div*=pow(vec[i].first,vec[i].second);
		}
		cout<<n/div<<endl; 
	}
	return 0;
}

下面是另一篇题解[CCF-CSP------因子化简)

他没有用欧拉筛选出所有的素数,而是在遍历过程中判断是否为素数,相当于而言减少很多不必要的计算。但是借此机会学习我上面提到的欧拉筛也是一桩美事~

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
 
//判断素数
int isPrime(int x) {
	for (int i = 2; i * i <= x; i++) {
		if (x % i == 0) {
			return 0;
		}
	}
	return 1;
}
 
void query(long long n, int k) {
	long long temp = n;
	long long final = 1;
	for(int i=2; i * i <= n; i++) {
		if(isPrime(i)){
			int num = 0;
			while(n%i==0){
				n/=i; 
				num++;
			}
//			printf("%d\n", num);
			if(num>=k){
				while(num--){
					final*=i;
				}
			}
		}
	}
	if(final==temp){
		printf("%lld\n", temp);
	}else if(final==1){
		printf("1\n");
	}else{
		printf("%lld\n", final);
	}
}
 
int main() {
	int q;
	scanf("%d", &q);
	for(int i=0; i<q; i++) {
		long long n;
		int k;
		scanf("%lld%d", &n, &k);
		query(n, k);
	}
	return 0;
}
 
相关推荐
刚学HTML36 分钟前
leetcode 05 回文字符串
算法·leetcode
Yan.love1 小时前
开发场景中Java 集合的最佳选择
java·数据结构·链表
AC使者1 小时前
#B1630. 数字走向4
算法
冠位观测者1 小时前
【Leetcode 每日一题】2545. 根据第 K 场考试的分数排序
数据结构·算法·leetcode
古希腊掌管学习的神2 小时前
[搜广推]王树森推荐系统笔记——曝光过滤 & Bloom Filter
算法·推荐算法
qystca2 小时前
洛谷 P1706 全排列问题 C语言
算法
浊酒南街2 小时前
决策树(理论知识1)
算法·决策树·机器学习
就爱学编程2 小时前
重生之我在异世界学编程之C语言小项目:通讯录
c语言·开发语言·数据结构·算法
学术头条2 小时前
清华、智谱团队:探索 RLHF 的 scaling laws
人工智能·深度学习·算法·机器学习·语言模型·计算语言学
Schwertlilien3 小时前
图像处理-Ch4-频率域处理
算法