Choose and divide UVA - 10375

本文介绍了一种计算10000以内素数及其对应指数的方法,并提供了一段C++代码实现。该算法首先找出指定范围内的所有素数,接着计算每个素数在特定数学表达式中的出现次数(即指数),最后利用这些素数及其指数来求解最终数值。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

找出10000以内对应的素数,然后计算出式子中每个素数对应的指数值,然后进行计算即可,具体实现见如下代码:

#include<iostream>
#include<vector>
#include<string>
#include<set>
#include<stack>
#include<queue>
#include<map>
#include<algorithm>
#include<cmath>
#include<iomanip>
#include<cstring>
#include<sstream>
#include<cstdio>
#include<deque>
#include<functional>
using namespace std;

vector<int> prime;
int e[10010];

bool is_prime(int n){
	for (int i = 2; i*i <= n; i++){
		if (n%i == 0) return false;
	}
	return true;
}

void getExp(int num,int d){
	for (int i = 0; i < prime.size(); i++){
		while (num%prime[i] == 0){
			num = num / prime[i];
			e[i] += d;
		}
		if (num == 1) break;
	}
}

void CountExp(int num,int d){
	for (int i = 2; i <= num; i++){
		getExp(i,d);
	}
}

int main(){
	for (int i = 2; i <= 10000; i++){
		if (is_prime(i)) prime.push_back(i);
	}
	int p, q, r, s;
	while (cin >> p >> q >> r >> s){
		memset(e, 0, sizeof(e));
		CountExp(p,1);
		CountExp(s, 1);
		CountExp(r - s, 1);
		CountExp(q, -1);
		CountExp(p - q, -1);
		CountExp(r, -1);
		double ans = 1;
		for (int i = 0; i < prime.size(); i++){
			ans *= pow(prime[i], e[i]);
		}
		printf("%.5lf\n",ans);
	}
	return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值