大家好，又见面了，我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。
Jetbrains全家桶1年46，售后保障稳定
今天学了数论。。。求原根真的好暴力
设模数为p,
我们把 p − 1 p-1 p−1分解质因数， p 1 , p 2 … p m p_1,p_2 \ldots p_m p1,p2…pm对于每一个 2 ≤ i ≤ p − 1 2 \leq i \leq p-1 2≤i≤p−1 ,判断 i p − 1 p j % p i^{p-1 \over p_j} \%p ipjp−1%p是否为1，如果是，那么这个数就不是原根，否则就是
AC Code
#include<cstdio>
#include<iostream>
#include<cmath>
using namespace std;
int zys[50];
int p;
int o;
int n;
int ksm(int a,int c,int b)
{
long long re=1;
long long t=a;
while(c)
{
if(c&1)re=re%b*t%b;
t=t%b*t%b;
c>>=1;
}
return re;
}
int pd(int x)
{
for(int i=1;i<=o;i++)
{
if(ksm(x,(p-1)/zys[i],p)%p==1)return 0;
}
return 1;
}
int main()
{
scanf("%d",&p);
int pp=p;
p--;
for(int i=2;i<=sqrt(p);i++)
{
if(p==1)break;
if(p%i==0)
{
o++;
zys[o]=i;
while(p%i==0)
{
p/=i;
}
}
}
if(p!=1)
{
o++;
zys[o]=p;
}
p=pp;
for(int i=2;i<p;i++)
{
if(pd(i))
{
printf("%d",i);
break;
}
}
return 0;
}
版权声明：本文内容由互联网用户自发贡献，该文观点仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请联系我们举报，一经查实，本站将立刻删除。
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/215435.html原文链接：https://javaforall.net
