博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
poj1284(欧拉函数+原根)
阅读量:6509 次
发布时间:2019-06-24

本文共 1275 字,大约阅读时间需要 4 分钟。

题目链接:https://vjudge.net/problem/POJ-1284

题意:给定奇素数p,求x的个数,x为满足{(xi mod p)|1<=i<=p-1}={1,2,...,p-1}。

思路:题目的实质就是问p有多少原根。

  下面是百度对原根的解释:

    设m是正整数,a是整数,若a模m的阶等于φ(m),则称a为模m的一个原根。(其中φ(m)表示m的欧拉函数)
    假设一个数g是P的原根,那么g^i mod P的结果两两不同,且有 1<g<P, 0<i<P,归根到底就是g^(P-1) = 1 (mod P)当且仅当指数为P-1的时候成立.(这里P是素数).
    简单来说,g^i mod p ≠ g^j mod p (p为素数)//这句话就是满足的条件。
    其中i≠j且i, j介于1至(p-1)之间
    这个题目就是将原根的定义解释了一遍。
  
  有两点重要的原根性质
    1. 模m有原根的充要条件是m= 1,2,4,p,2p,p^n,其中p是奇质数,n是任意正整数。
    2. 当模m有原根时,它有φ(φ(m))个原根。

  某大牛的证明(没看懂...QAQ):

    {xi%p | 1 <= i <= p - 1} = {1,2,...,p-1} 等价于 {xi%(p-1) | 1 <= i <= p - 1} = {0,1,2,...,p-2},即为(p-1)的完全剩余系  

    若x,x2...x(p-1)是(p-1)的完全剩余系,

    根据定理,可以推出若gcd(x, p-1) = 1时, (1,x,...,x(p-2))也是(p-1)的完全剩余系

    因为若xi != xj (mod p-1),那么x*xi != x*xj (mod p-1),与条件m矛盾,所以 xi = xj (mod p-1),

    由此可以确定答案为eu(p-1)。

  知道答案是eu(p-1),代码就很好实现了,筛法打表65525以内的数的欧拉函数即可。

AC代码:

#include
using namespace std;int eu[66000],p;void eular(){ for(int i=2;i<=65536;++i) if(!eu[i]) for(int j=i;j<=65536;j+=i){ if(!eu[j]) eu[j]=j; eu[j]=eu[j]/i*(i-1); }}int main(){ eular(); while(~scanf("%d",&p)){ printf("%d\n",eu[p-1]); } return 0;}

 

转载于:https://www.cnblogs.com/FrankChen831X/p/10827190.html

你可能感兴趣的文章
汇编语言的应用
查看>>
device platform 相应的表
查看>>
php des 加密解密实例
查看>>
【Mac】Mac键盘实现Home, End, Page UP, Page DOWN
查看>>
实战使用Axure设计App,使用WebStorm开发(1) – 用Axure描述需求
查看>>
安德鲁斯----多媒体编程
查看>>
[zz]在linux中出现there are stopped jobs 的解决方法
查看>>
Delphi下实现全屏快速找图找色 一、数据提取
查看>>
查询表字段信息
查看>>
关于机器学习的最佳科普文章:《从机器学习谈起》
查看>>
dxFlowChart运行时调出编辑器
查看>>
NET Framework 3.0 (WinFX) RTM发布
查看>>
图片拼接器
查看>>
C++ TinyXml操作(含源码下载)
查看>>
中断小笔记
查看>>
C#委托、事件、消息(入门级)
查看>>
FreeBinary 格式说明
查看>>
使用Spring Cloud和Docker构建微服务
查看>>
NB-IoT的成功商用不是一蹴而就
查看>>
九州云实战人员为您揭秘成功部署OpenStack几大要点
查看>>