Primitive Roots(Poj1284)(欧拉函数+原根)

Primitive Roots
Time Limit: 1000MS Memory Limit: 10000K
Total Submissions: 3462 Accepted: 2025

Description

We say that integer x, 0 < x < p, is a primitive root modulo odd prime p if and only if the set { (x i mod p) | 1 <= i <= p-1 } is equal to { 1, ..., p-1 }. For example, the consecutive powers of 3 modulo 7 are 3, 2, 6, 4, 5, 1, and thus 3 is a primitive root modulo 7. 
Write a program which given any odd prime 3 <= p < 65536 outputs the number of primitive roots modulo p. 

Input

Each line of the input contains an odd prime numbers p. Input is terminated by the end-of-file seperator.

Output

For each p, print a single number that gives the number of primitive roots in a single line.

Sample Input

23
31
79

Sample Output

10
8

24

原根求的是 { (x^i mod p) | 1 <= i <= p-1 } is equal to { 1, ..., p-1 }. 

就是一个数x的i次方 (1<=i<=p-1) mod  p等于  a(1<=a<=p-1),即求euler(p),euler(p)

即欧拉函数:求小于或等p且与p互质的个数。

推导原根的过程还没搞懂~ 以后补上推导过程。先练习下欧拉函数使用!

#include<stdio.h>
int euler(int n)
{
	int ans=n;
	for(int i=2;i*i<=n;i++)
	{
		if(n%i==0)
		{
			ans=ans/i*(i-1);
			n/=i;
		}
		while(n%i==0)
		{
			n/=i;
		}
	}
	if(n>1)
	    return ans/n*(n-1);
	else
	    return ans; 
}
int main()
{
	int n;
	while(scanf("%d",&n)!=EOF)
	{
		printf("%d\n",euler(n-1));//或者euler(euler(n))也对。 
	} 
} 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是一个 C++ 函数,它可以用于计算大素数的原根: ``` #include <iostream> #include <vector> #include <algorithm> #include <cmath> using namespace std; // 判断一个数是否为素数 bool is_prime(int n) { if (n <= 1) { return false; } int sqr = sqrt(n); for (int i = 2; i <= sqr; i++) { if (n % i == 0) { return false; } } return true; } // 计算 a 的 b 次方对 p 取模的结果 int power_mod(int a, int b, int p) { int res = 1; while (b) { if (b & 1) { res = (long long)res * a % p; } a = (long long)a * a % p; b >>= 1; } return res; } // 获取大素数 p 的原根 int get_primitive_root(int p) { if (!is_prime(p)) { return -1; } vector<int> factors; int phi = p - 1; int tmp = phi; for (int i = 2; i <= tmp / i; i++) { if (tmp % i == 0) { factors.push_back(i); while (tmp % i == 0) { tmp /= i; } } } if (tmp > 1) { factors.push_back(tmp); } for (int g = 2; g <= p; g++) { bool flag = true; for (int i = 0; i < factors.size(); i++) { if (power_mod(g, phi / factors[i], p) == 1) { flag = false; break; } } if (flag) { return g; } } return -1; } int main() { int p; cout << "请输入一个大素数 p:" << endl; cin >> p; int g = get_primitive_root(p); if (g == -1) { cout << "无法找到 " << p << " 的原根!" << endl; } else { cout << p << " 的一个原根为:" << g << endl; } return 0; } ``` 该函数首先判断输入的 p 是否为素数,如果不是,则输出错误信息并返回 -1。接着,它将 p - 1 分解为质因数的积,并逐个尝试每个可能的原根。对于每个可能的原根 g,它计算 g 的 phi(p) / factors[i] 次方对 p 取模的结果,并检查是否等于 1。如果对于所有质因子 i,g 的 phi(p) / factors[i] 次方对 p 取模的结果都不等于 1,则 g 是 p 的原根。如果无法找到 p 的原根,则输出错误信息并返回 -1。 注意,计算原根的时间复杂度为 O(p log^2 p),因此该函数对于较大的素数可能需要较长的计算时间。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值