【板子】欧拉函数

最新推荐文章于 2021-09-07 21:58:46 发布

原创最新推荐文章于 2021-09-07 21:58:46 发布 · 225 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

板子专栏收录该内容

5 篇文章

订阅专栏

本文深入探讨了欧拉函数的定义与计算方法，通过一个典型算法问题——求解小于等于n的正整数中与n互质的整数数量，详细介绍了如何使用欧拉函数进行高效计算。文章提供了完整的C++代码实现，并解释了欧拉函数背后的数学原理。

例题

A - Relatives
Given n, a positive integer, how many positive integers less than n are relatively prime to n? Two integers a and b are relatively prime if there are no integers x > 1, y > 0, z > 0 such that a = xy and b = xz.
Input
There are several test cases. For each test case, standard input contains a line with n ≤ 1,000,000,000. A line containing ‘0’ follows the last case.
Output
For each test case there should be single line of output answering the question posed above.

Sample Input
7
12
0
Sample Output
6
4
题目大意：输入n,求出1到n(包括n)中与n互质的整数的个数，注意当n=1时，结果为0。

解题方法：欧拉函数

欧拉函数的定义：
欧拉函数是计算少于或等于n的数中与n互质的数的数目,计算公式为
euler(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…(1-1/pn)
其中p1,p2……pn为x的所有素因数，x是不为0的整数。

#include <iostream>
#include <algorithm>
#include <string>
#include <cstring>
#include <map>
#include <set>
#include <deque>
#include <queue>
#include <vector>
#include <cstdio>
#include <cstdlib>
#include <iomanip>
#include <math.h>
#define ull unsigned long long
#define pi acos(-1)
using namespace std;
typedef long long ll;
const ll inf = 0x3f3f3f3f;
const ll maxn = 1e7+5;
const ll mod = 1e9+7;
ll euler(ll n){//欧拉函数是少于或等于n的数中与n互质的数的数目,计算公式为euler(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…(1-1/pn)其中p1,p2……pn为x的所有素因数，x是不为0的整数
	ll ans=n;
	for(int i=2;i*i<n;i++){
		if(n%i==0) ans=ans/i*(i-1);
		while(n%i==0) n/=i;
	}
	if(n>1) ans=ans/n*(n-1);
	return ans;
}
int main(){
	ll n;
	while(cin>>n&&n){
		cout<<euler(n)<<'\n';
	} 
	return 0;
}