poj 1091 (容斥原理)

本文介绍了一种利用质因数分解和快速幂算法解决特定数学问题的方法。通过分解给定整数m的质因数,并使用快速幂算法计算m的n次方,再结合组合数学原理去除重复项,最终求得精确结果。文章详细展示了算法实现过程及代码示例。

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

#include <iostream>
#include <cstdio>
#include <iomanip>
#include <string>
#include <cstdlib>
#include <cstring>
#include <queue>
#include <set>
#include <vector>
#include <map>
#include <algorithm>
#include <cmath>
#include <stack>
#include <stdlib.h>
#include <stdio.h>

#define INF 0x3f3f3f3f
#define LINF 0x3f3f3f3f3f3f3f3f
#define ll long long
#define ull unsigned long long
#define uint unsigned int
#define l(x) x<<1
#define r(x) x<<1|1
#define ms(a,b) memset(a,b,sizeof(a))

using namespace std;

int n, m;
vector<int> x;
void getprime(int m)
{
	//得到m的质因子
	x.clear();
	for (int i = 2; i*i <= m; ++i)
		if (m%i == 0)
		{
			x.push_back(i);
			while (m%i == 0) m /= i;
		}
	if (m != 1) x.push_back(m);
}
int main() {
	while (scanf("%d%d", &n, &m) == 2) {
		getprime(m);//分解m
		ll ans = pow(m, n);//ans为m^n
		for (int i = 1; i < (1 << x.size()); ++i) {//m素因子个数为t。则i<2^t。即m素因子集合的子集个数
			int k = 1, cnt = 0;
			for (int j = 0; j < x.size(); ++j)
				if (i&(1 << j))
					k *= x[j], ++cnt;//k为取出的这个子集的值的积。cnt表示子集中元素的个数
			if (cnt & 1) //等价于cnt%2==1。取出个数为奇数时减,取出个数为偶数时加
				ans -= pow(m / k, n);
			else
				ans += pow(m / k, n);
		}
		printf("%lld\n", ans);
	}
	return 0;
}    
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值